In this Chapter, we have studied a spatially homogeneous and anisotropic Kantowski-Sachs universe in the presence of Barrow holographic dark energy (BHDE) in the background of Saez-Ballester (SB) scalar-tensor theory of gravitation. To find the exact solution of the SB field equations, we have assumed that the shear scalar is directly proportional to expansion scalar. This assumption leads to relation between metric potentials of the models. We have discussed non-interacting and interacting cosmological models. Moreover, we have discussed several cosmological parameters such as energy densities of DM and DE ($\rho_m$ \& $\rho_b$), deceleration $(q)$, equation of state ($\omega_b$) and skewness ($\alpha$) parameters, squared sound speed ($v_s^2$), $\omega_b$-$\omega_b^{'}$ plane statefinders and $Om$-diagnostics parameters through graphical representation for both the interacting models. Also, we have observed that the current values of deceleration and EoS parameters of our constructed models are coincide with the recent observational data.