Question: Assume $$X_t=e^{-at}\int_0^te^{au}dB_u,$$ where $B$ is a standard Brownian motion. Show that for any $t,s\in \mathbb{R}_+$, we have $$\mathbb{E}\big[(X_s^2-\mathbb{E}[X_s^2])(X_t^2-\mathbb{E}[X_t^2])\big]=2\mathbb{E}[X_sX_t]^2.$$ The proof shouldn't be calculating explicit out the Left hand side and Right hand side and find them equal.
My trail:
First I proved by Isometry that $X_t$ is Gaussian, i.e. $\sim\mathcal{N}(0,\frac{1}{2a}(1-e^{-2at}))$ by $$\mathbb{E}[X_t^2]=\mathbb{E}[(e^{-at}\int_0^te^{au}dB_u)^2]=e^{-2at}\int_0^te^{2au}du=\frac{1}{2a}(1-e^{-2at}).$$
Then $$\mathbb{E}\big[(X_s^2-\mathbb{E}[X_s^2])(X_t^2-\mathbb{E}[X_t^2])\big]=\mathbb{E}[X_s^2X_t^2]-\mathbb{E}[X_s^2]\mathbb{E}[X_t^2].$$ Then I don't know how to proceed.
Hints: