A hybrid simulation of non-resonant fishbone (NRF) instabilities with reversed safety factor profile is investigated with a global kinetic-magnetohydrodynamic (MHD) code (M3D-K). With EAST parameters, NRF instability can be driven by energetic particles as minimum safety factor is a little greater than unity. With analysis on energy resonance, physical mechanism of NRF excited by energetic particles at different angles is studied. Parameters such as injection energy E, P_hot/P_total, are scanned to find parameter thresholds and change of fishbone instability. In addition, we analyze slowing-down distribution of energetic particles, evolution of mode structure and fishbone mode sweep in nonlinear process. A flattening trend of nonlinear slowing-down distribution at Λ = 0.6, 0.7, 1.0 is consistent with that of energetic particles in experiments. Mode structure evolves from that of NRF mode to those of other high-frequency modes at Λ = 0.6. The upward sweep frequency of fishbone is consistent with the sweep frequency of classical fishbone at Λ = 0.6.