Hint:
To understand how to work with this type of integral, first consider an integral of Brownian motion:
$$I = \int_{0}^{T} B_t dt$$
The integral makes sense because Brownian motion has almost-surely continuous sample paths. Consider the approximation as a Riemann sum over a partition of $[0,T]$:
$$ S_n = \sum_{k=1}^{n} B_{t_k} (t_k - t_{k-1})$$
Now you can think of $S_n$ as a random variable with a multi-variate normal distribution. Try to establish properties of this random variable such as the mean, variance, and higher moments. You will need to use the fact that the Brownian motion at different times is correlated, i.e.,
$$E(B_{t_1}B_{t_2}) = min(t_1,t_2)$$
Then consider taking limits as $n \rightarrow \infty$ ($\max(t_k-t_{k-1}) \rightarrow 0$).
You will find for example that
$$E(I) = 0, E(I^2) = T^3/3, ...$$
Then you can move on to a more general case where the integrand is a function of the Brownian motion.