I've found in different references that $\int_0^t g(s)dB_s$ is a Wiener integral, where $g(s)$ is a bounded, continuous function and $dB_s$ is the standard Brownian Motion. Also, that $\int_0^t g(s)dB_s \sim \mathcal{N}(0, \int_0^t g^2(s) ds)$. Why is this?
My attempt so far, is to approximate it as a series of normal random variables. Since the increments $B_{t_{i+1}}-B_{t_i}$ have a distribution $\mathcal{N}(0, t_{i+1}-t_i)$ and they are pairwise independent. Then, thinking it as a Riemann integral $\sum_{i=0}^{\infty} g_i(s) (B_{t_{i+1}}-B_{t_i}) \sim \mathcal{N}(0, \sum_{i=0}^{\infty} g_i^2(s)(t_{i+1}-t_i))$ for some $s\in [t_i, t_{i+1}]$. My question is, is it posible to take the limit $\lim_{t_{i+1} \to t_i} \sum_{i=0}^{\infty} g_i^2(s) (t_{i+1}-t_i) = \int_0^t g^2(s)ds$ in this case? Why?
We consider the general case that $g: [0,t] \longrightarrow \mathbb{R}$ is a square integrable function.
Choose a sequence of simple/elementary functions $(g_n)_{n\in\mathbb{N}}$ that converge in $L^2([0,t])$ to $g$. Then, by Ito's Isometry one has $$ \mathbb{E}\left[ \left(\int_0^t g(s) d B_s - \int_0^t g_n(s) d B_s \right)^2 \right] = \mathbb{E}\left[ \left(\int_0^t (g(s)-g_n(s)) d B_s\right)^2\right] = \|g_n-g\|_{L^2}^2 \to 0, $$ as $n\to\infty$. This implies that $\left(\int_0^t g_n(s) d B_s\right)_{n\in\mathbb{N}}$ converges in $\int_0^t g(s) d B_s$ in distribution.
Next, calculate the characteristic function of $\int_0^t g_n(s) d B_s$. As you have seen, $\int_0^t g_n(s) d B_s$ is simply a sum of independent Normal-distributed variables and you should see that this is the characteristic function of a $\mathcal{N}(0,\int_0^t g_n^2(s) d s)$-distributed random variable. Use $L^2$-convergence of $(g_n)_{n\in\mathbb{N}}$ to $g$ to show that the characteristic function of $\int_0^t g_n(s) d B_s$ converges point-wise to the characteristic function of a $\mathcal{N}(0,\int_0^t g^2(s)ds)$-distributed random variable. This implies $$ \operatorname{law}(\int_0^t g_i(s) d B_s) \Longrightarrow \mathcal{N}(0,\int_0^t g^2(s) d s) $$ as $n\to\infty$. Since the limit in distribution is unique one has $$ \operatorname{law}(\int_0^t g(s) d B_s) = \mathcal{N}(0,\int_0^t g^2(s) d s). $$