A magnetic field, imposed on turbulent flow of an electrically conductive fluid, is known to cause preferential damping of the velocity and its fluctuations in the direction of Lorentz force, thus leading to an increase in stress anisotropy. Based on direct numerical simulations ͑DNS͒, we have developed a model of magnetohydrodynamic ͑MHD͒ interactions within the framework of the second-moment turbulence closure. The MHD effects are accounted for in the transport equations for the turbulent stress tensor and energy dissipation rate-both incorporating also viscous and wall-vicinity nonviscous modifications. The validation of the model in plane channel flows with different orientation of the imposed magnetic field against the available DNS (Reϭ4600,Haϭ6), large eddy simulation (Reϭ2.9ϫ10 4 ,Haϭ52.5,125) and experimental data (Reϭ5.05ϫ10 4 and Reϭ9ϫ10 4 , 0рHaр400), show good agreement for all considered situations.