$$
d Y \left(t\right) := d \left[\int_0^t{a \left(s\right)\mathrm{d}W_s}\right]
= a \left(t\right) dW_t
$$
Note that since $Y$ is a driftless process, it is a local martingale, and because $a$ is bounded, a true martingale. Its quadratic variation is given by
$$
\langle Y \left(\cdot\right)\rangle_t = \int_0^t{a^2 \left(s\right)\mathrm{d}s}
$$
by definition of the stochastic integral with respect to the Wiener process.
Using Itō's lemma,
$$
d \left[\left[Y \left(t\right)\right]^2\right] = 2 Y \left(t\right) d Y\left(t\right) + d \langle Y \left(\cdot\right)\rangle_t
$$
Subtracting the differential of the time integral, i.e. $a^2 \left(t\right) \, dt$, removes the drift term due to $d \langle Y \left(\cdot\right)\rangle_t$ and you are done.