A more general version of the result Beenakker cites that allows arbitrary $B$ (not just symmetric):$\newcommand\tr{\operatorname{tr}}\newcommand\E{\mathbb{E}}$
 
Take $x\sim\mathcal{N}(0,I)$ and recall that $\tr(C) = \E[x^*Cx]$. If $B = P \Lambda P^{-1}$ is diagonalizable let $\Lambda^{1/2}$ be a (possibly complex) square root of $\Lambda$. Then using the cyclic property of the trace
$$
\begin{align}
\tr(AB) &= \tr(\Lambda^{1/2} P^{-1}AP\Lambda^{1/2})\\
&= \E[(\Lambda^{1/2}x)^* P^{-1}A P (\Lambda^{1/2}x)]\\
&= \sum_{i\neq j} \lambda_i^{1/2}\lambda_j^{1/2}\E[x_ix_j](P^{-1}A P)_{ij} + \sum_{i} \lambda_i\E[x_i^2](P^{-1}A P)_{ii}\\
&= \sum_{i} \lambda_i(P^{-1}A P)_{ii}.
\end{align}
$$
From here, since $(P^{-1} A P)_{ii}\geq 0$ by the PSD condition, it follows that
$$\tr(AB) \leq \lambda_{\max}\sum_{i} (P^{-1}A P)_{ii} = \lambda_{\max}\tr(P^{-1}AP) = \lambda_{\max}\tr(A)$$
and
$$\tr(AB) \geq \lambda_{\min}\sum_{i} (P^{-1}A P)_{ii} = \lambda_{\min}\tr(P^{-1}AP) = \lambda_{\min}\tr(A)$$
as desired. Density will tell us this also holds for non-diagonalizable $B$.