The combined effects of radiation and mass transfer on a steady MHD two-dimensional Marangoni convection flow over a flat surface in presence of Joule heating and viscous dissipation under influence of suction and injection is studied numerically. The general governing partial differential equations are transformed into a set of nonlinear ordinary differential equations by using unique similarity transformation. Numerical solutions of the similarity equations are obtained using the Runge-Kutta method along with shooting technique. The effects of governing parameters on velocity, temperature, and concentration as well as interface velocity, the surface temperature gradient, and the surface concentration gradient were presented in graphical and tabular forms. Comparisons with previously published work are performed and the results are found to be in excellent agreement.