The range of the U bosonic coupling constants in neutron star matter is a very interesting but still unsolved problem which has multifaceted influences in nuclear physics, particle physics, astrophysics and cosmology. The combination of the theoretical numerical simulation and the recent observations provides a very good opportunity to solve this problem. In the present work, the range of the U bosonic coupling constants is inferred based on the three relations of the mass–radius, mass-frequency and mass-tidal deformability in neutron stars containing hyperons using the GM1, TM1 and NL3 parameter sets under the two flavor symmetries of SU(6) and SU(3) in the framework of the relativistic mean field theory. Combined with observations from PSRs J1614-2230, J0348+0432, J2215-5135, J0952-0607, J0740+6620, J0030-0451, J1748-2446ad, XTE J1739-285, GW170817 and GW190814 events, our numerical results show that the U bosonic coupling constants may tend to be within the range from 0 to 20 GeV−2 in neutron star containing hyperons. Moreover, the numerical results of the three relations obtained by the SU(3) symmetry are better in accordance with observation data than those obtained by the SU(6) symmetry. The results will help us to improve the strict constraints of the equation of state for neutron stars containing hyperons.