Using the Multiconfiguration Dirac-Hartree-Fock method as implemented in the General Relativistic Atomic Structure Package, the magnetic dipole and electric quadrupole hyperfine structure constants were determined for the ground and first excited levels of 135,137Ba II isotopes, as well as for 137Ba I and 87Sr II, to assess the robustness of the developed model. This study builds upon and extends previous investigations by examining the levels involved in resonance lines, with the aim of resolving persistent discrepancies in the hyperfine structure of 137Ba II and 87Sr II. New code developments such as the use of natural orbitals, as well as the addition of polarization effects and Configuration State Function Generators, as implemented in GRASPG, were tested for these heavy elements. The developed strategy allowed us to achieve encouraging results that satisfactorily agree with experiments for all studied levels but D5/22 in the 137Ba II isotope. This disagreement was also observed in 135Ba II isotope as well as in 87Sr II. With two valence electrons, 137Ba I is definitely more complex, requiring a multireference approach. Even with the latter, the theory-observation disagreement observed for the hyperfine structure of the low-lying levels remains large in comparison with the alkali-like systems. Possible ongoing developments to remediate this issue are discussed in the conclusions.