In this study, we investigate the impact of modified gravity (MG) on the merger rate of compact binaries within dark matter spikes surrounding supermassive black holes (SMBHs). Specifically, we calculate the binary merger rates involving primordial black holes (PBHs) and/or neutron stars (NSs) in Hu–Sawicki f(R) gravity and the normal branch of Dvali–Gabadadze–Porrati (nDGP) gravity, with three SMBH mass functions: Benson, Vika, and Shankar. The results show consistently higher merger rates predicted for PBH–PBH and PBH–NS binaries in these gravity models compared to general relativity (GR), in particular at lower SMBH masses and for steeper dark matter spike density profiles. The predicted merger rates are compared to the LIGO–Virgo–KAGRA observations in constraining the parameters of the theory. In particular, we find steeper dark matter spike density profiles in the MG scenarios compared to GR. When compared to current observational constraints on PBH abundance, the mass ranges allowed by Hu–Sawicki f(R) models are found to be wider than those allowed by nDGP models, for given merger rates. The results are highly dependent on the choice of SMBH mass function, with the Vika and Shankar mass functions predicting lower abundances. The considerable sensitivity of the results to the assumed gravity scenario and SMBH mass function demonstrates the necessity of incorporating the corresponding theoretical uncertainties when making relatively robust predictions on compact binary merger rates and, as a result, on PBH properties.