This paper is concerned with a stability problem on perturbations near a physically important steady state solution of the 3D MHD system. We obtain three major results. The first assesses the existence of global solutions with small initial data. Second, we derive the temporal decay estimate of the solution in the L2-norm, where to prove the result, we need to overcome the difficulty caused by the presence of linear terms from perturbation. Finally, the decay rate in L2 space for higher order derivatives of the solution is established.