High-frequency gravitational waves, particularly in the range f ≳ 10¹⁰ Hz, represent a compelling probe of physics beyond the Standard Model. Due to the absence of direct detection methods in this frequency regime, alternative strategies may be pursued. One promising approach involves the conversion of gravitons into photons in the presence of magnetic fields, a process known as the inverse Gertsenshtein effect. In this study, we explore such graviton-to-photon conversions occurring within the magnetic field environment of the M87 galaxy, utilizing realistic models for the galactic magnetic field and plasma density structure. We use the broadband electromagnetic spectrum of M87, ranging from millimeter to TeV gamma rays, to search for hidden contributions from graviton-photon conversions. In the well-constrained frequency range 10¹⁰-10²⁷ Hz, the lack of excess emission allows us to place improved bounds on the gravitational wave strain amplitude h_c or on spectral energy density Ω_gw h². We find that our results from M87 yield substantially stronger constraints compared to existing bounds derived from Milky Way magnetic field considerations, with improvements ranging from one to five orders of magnitude depending on the frequency band, thereby enhancing the prospects for probing high-frequency gravitational wave backgrounds through indirect electromagnetic signatures.