A Finite Element Interior Penalty Method for MHD in heterogeneous domains