The effects of energetic particles on MHD instabilities is a key issue in the basic physics that will be important in burning plasma experiments such as ITER. Recently, experiments (JET, JT-60U and DIII-D) have shown that the kinetic effects of energetic particles can play a crucial role in the stability of the m/n = 2/1 tearing mode, where the fraction of energetic particle βfrac ≡ βh/β is high (β = 2μ0P/B2, P is pressure, B is magnetic field and βh is the energetic particle β). For the study of the kinetic effects of energetic particles on resistive MHD stability, using model equilibria based on DIII-D experimental reconstructions, the non-ideal MHD linear stability of cases unstable to the 2/1 mode is investigated including a δf PIC model for the energetic particles coupled to the nonlinear 3D resistive MHD code NIMROD (Kim et al 2008 Phys. Plasmas15 072507). The growth of unstable modes is calculated at a series of β, βfrac and S = τR/τA (the ratio of the resistive time to Alfvén time), spanning from the resistive to the ideal unstable regime of the mode, up to and above βfrac = 0.25, and well into the asymptotic regime of the resistivity. It is observed that energetic particles have significant damping and stabilizing effects at experimentally relevant β, βfrac and S, and less damping and stabilizing effects at the ideal unstable regime. It is also observed that energetic particles excite a real frequency of the 2/1 mode. The growth rates significantly reduce at experimentally relevant β, βfrac and S due to the mode interacting with the trapped particles and 'barely passing' particles. Furthermore, extrapolation of the results is discussed for implications for JET and ITER, where the effects are projected to be significant.