Key idea. Subtract from X1 the multiple of the sum that is uncorrelated with the sum. For jointly normal variables, uncorrelated means independent.
Let S=X1+X2. The pair (X1,S) is a linear image of the independent normals (X1,X2), so it is jointly normal. We have
Var(S)=9+16=25,Cov(X1,S)=Var(X1)=9.
Put R=X1−259S. Then
Cov(R,S)=9−259⋅25=0.
Since (R,S) is jointly normal and uncorrelated, R is independent of S. Also E[R]=0. Therefore
E[X1∣S]=259S+E[R]=259S.
So each summand takes a share of the sum proportional to its variance.
At S=5 this gives 259⋅5=59.
The conditional variance is Var(R)=9−2581=25144. Hence
X1∣S=5 ∼ N(59,25144).
E[X1∣X1+X2=5]=59