Hi Apeirohedron,
your expectation is $\sigma_x^2 + b^2$, since $E[(b+X)^2|Y] = E[X^2|Y] + 2bE[X|Y] + b^2 = \sigma_x^2 + 2bE[X] + b^2 = \sigma_x^2 + b^2$. $E[X|Y] = E[X]$ holds since $X$ and $Y$ are independent.
In the general case: You know that $X$ and $Z$ are Gaussian. In particular, you know that $X\sim \mathcal N(0,\sigma_x^2)$ and $Y|X \sim \mathcal N(x,\sigma_z^2)$. Therefore, playing with conditional Gaussians, you get $X|Y \sim \mathcal N\left(\frac{y}{\sigma_z^2/\sigma_x^2 + 1},(\sigma_x^{-2}+\sigma_z^{-2})^{-1}\right)$. You can find the formulae for that in
Bishop CM. Pattern Recognition and Machine Learning (Information Science and Statistics). Springer; 2007.
on page 93. Your conditional expectation is, therefore, $E[X|Y] = \frac{y}{\sigma_z^2/\sigma_x^2 + 1}$.