Let $A$, $B$, $\epsilon_a$ and $\epsilon_b$ be independent normal random variables:
$A \sim N(\mu_A, \sigma^{2}_A)$
$B \sim N(\mu_B, \sigma^{2}_B)$
$\epsilon_a \sim N(0, \sigma^{2}_{\epsilon_A})$
$\epsilon_b \sim N(0, \sigma^{2}_{\epsilon_B})$
Let Z be defined as $Z \equiv A + B + \epsilon_a + \epsilon_b$.
I am interested in computing the expectation: $E(AB|Z)$.
From questions like this I realize that I (most likely) cannot simply write $E(AB|Z) = E(A|Z)E(B|Z)$, which would be easy to solve using the bivariate normal distribution.
How can I go about computing $E(AB|Z)$ in terms of means and variances?