Your line $$E(aXY+aZ+bEY+bEZ)-[(aEX+b)(EY+EZ)]$$ has a mistake: it should be $$E(aXY+a\color{red}{XZ}+bEY+bEZ)-[(aEX+b)(EY+EZ)]\tag1$$ After you fix this, you should be able to collect terms properly: $$ \begin{align} (1)&=aE(XY)+aEXZ+bEY+bEZ-[aEXEY+aEXEZ+bEY+bEZ]\\\ &=aE(XY)-aEXEY +aEXZ-aEXEZ \end{align} $$ since the terms $bEY$ and $bEZ$ drop out.