$\newcommand{\E}{\mathbb{E}}$
Unfortunately the upper bound $n^{3/2}$ you obtained cannot be improved. In fact 
$$\E \sup_{Q} \langle W, Q Q^T \rangle = \Theta(n^{3/2}).$$

First of all, note that it is enough to bound this quantity for $Z$ a matrix with independent gaussians (as opposed to a symmetric matrix), since $\sqrt{2} W = Z + Z^T + D$ where $Z$ is a matrix with independent gaussian entries, and $D$ is a diagonal matrix with entries $D_{ii} = (\sqrt{2} - 2)Z_{ii}$. It is fairly easy to bound $\mathbb{E} \sup_Q \langle D, Q Q^T \rangle \leq O(\sqrt{n})$, and since for any $Q$ we have $\langle Z^T, Q Q^T\rangle = \langle Z, Q Q^T \rangle$, we get

$$
\E \sup_Q \langle W, Q Q^T \rangle \approx \E \sup_Q \langle Z, Q Q^T \rangle
$$

Moreover note that by restricting our choice to $\theta_i \in \pm \pi/2$ we can get $Q Q^T$ to be any rank one matrix of form $q q^T$ for $q \in \{\pm 1\}$. On the other hand, for an upper bound we can use a simple triangle inequality to get
$$
\sup_{q \in \{\pm 1\}^n} \langle Z, q q^T \rangle \leq \sup_{Q} \langle Z, Q Q^T \rangle \leq 2 \sup_{q \in [\pm 1]^n} \langle Z, q q^T \rangle,
$$
i.e. it is enough to understand the maximum inner product of a gaussian matrix, with a rank one matrix of form $q q^T$ for $q \in [\pm 1]^n$. Clearly for any realization $Z$ the supremum is obtained at one of the vertices, so in fact we get
$$
\E \sup_{Q} \langle Z, Q Q^T \rangle \approx \E \sup_{q \in \{\pm 1\}^n} \langle Z, q q^T \rangle.
$$

Let us consider now a gaussian process $X_q := \langle Z, q q^T \rangle$, for $q \in \{\pm 1\}^n$. By the Talagrands majoring measure theorem (see for example section 3.4 in [chaining](https://dash.harvard.edu/bitstream/handle/1/34872844/57723140.pdf?sequence=1&isAllowed=yand) for relevant definitions) we have
$$
\E \sup_{q \in \{\pm 1\}^n} X_q \approx \gamma_2(\{\pm 1\}^n, d_X), \tag{1}
$$
where $d_X : \{\pm 1\}^n \times \{ \pm 1\}^n \to \mathbb{R}_{\geq 0}$ is a pseudo metric given by $d_X(u, v) := \sqrt{\E (X_u - X_v)^2}$. 

Note that for $u, v \in \{\pm 1\}^n$, we have 
$$
\mathbb{E} (X_u - X_v)^2 = \|u u^T - v v^T\|_F^2,
$$
where $\|A\|_F^2 := \sum_{ij} A_ij^2$. When $\|u\|_2 = \|v\|_2 = \theta$ an elementary calculation leads to 
$$
\|u u^T - v v^T\|_F = \sqrt{2} \theta \|u - v\|_2,
$$
hence 
$$\gamma_2(\{\pm 1\}^n, d_X) = \gamma_2(\{\pm 1\}^n, 2 \sqrt{n} d_2) = 2 \sqrt{n} \gamma_2(\{\pm 1\}^n, d_2),\tag{2}$$ where $d_2(u,v) = \|u - v\|_2$. 

Finally, applying the majorizing measures theorem again, we know that $\gamma_2(\{\pm 1\}^n, d_2) \approx \E \sup_{q \in \{\pm 1\}} \langle q, G \rangle$, where $G$ is a gaussian vector in $\mathbb{R}^n$.

This gives 
$$\gamma_2(\{\pm 1\}^n, d_2) \approx \E \|G\|_1 \approx n.\tag{3}$$

Combining (1), (2) and (3) gives
$$
\E \sup_{q \in \{\pm 1\}^n} X_q \approx n^{3/2}.
$$