• A novel attempt has been put together to design a Model Predictive Control based on the leader Harris Hawks Optimization approach for a PV integrated LFC system. Further, a combined AVR-LFC loop has been designed for a PV, wind and hydro integrated hybrid power system. • The validation of the controller's better performance has been obtained through a thorough comparative analysis where the transient response attributes from the proposed controller have been compared with the exiting controllers available in the literature. Also, the convergence mobility of the proposed LHHO optimizer has been compared with HHO, WOA, and PSO algorithms. • To further strengthen the frequency stability of the network, the authors have proposed and implemented the virtual inertia and capacitive energy storage devices in association with the LHHO-MPC controller. • The robust nature of the proposed controller algorithm has been assessed under different test case scenarios where the step load perturbation (SLP) values have been varied randomly, and the stochastic load model has been considered. Also, various nonlinearities and reduced inertia cases have been presented to show the efficiency of the proposed LHHO-MPC controller working in coordination with CES and VI units. The work in this manuscript aims at designing a novel control strategy based on the Model Predictive Controller aided with Leader Harris Hawks Optimization (MPC-LHHO) algorithm for the regulation of frequency and voltage in renewable penetrated power systems. Herein, two cases have been considered where the first case discusses the frequency stabilization in a thermal power system linked with the solar photovoltaic (PV) plant. Subsequently, the second case features the simultaneous voltage and frequency regulation in a hybrid interconnected power system including thermal, wind, diesel, photovoltaic, and hydro generation units. Additionally, for both cases, the proposed MPC-LHHO algorithm has been evaluated in coordination with the capacitive energy storage and virtual inertia units. Furthermore, the stability analysis of the controlled PV-thermal system has been conducted using eigenvalues and pole-zero plots. In both cases, the effectiveness of the proposed algorithm has been tested after considering diverse scenarios including random loading, nonlinearities, time delays, reduced inertia and stochastic variations in load. For all these cases the presented algorithm has given a stable response and has ensured proper regulation of voltage and frequency in the renewable integrated power system network.