This work promotes the room temperature energy storage properties of the multiferroics. In this approach, impacts of PrFeO₃ doping on PT-based solid solutions (Pb${}_{1-}$${}_x$Pr${}_x$Ti${}_{1-}$${}_x$Fe${}_x$O₃, $x$ = 0.21, 0.22, 0.23, 0.24, 0.25 and 0.26) have been explored. X-ray diffraction (XRD) patterns were used to estimate the crystallographic parameters, confirming the single phase tetragonal structure. The ferroelectric Curie temperature ($T_c$${}^{\mathrm{\text{FE}}}$) is observed to drop from 410 K to below room temperature as the Pr concentration increases. The ferroelectric P-E loops were used to determine the energy storage values at room temperature. The sample $x$ = 0.24 achieved the maximum value of energy storage density of 362.25 mJ/cm³ with the efficiency of 40.5%. The ferroelectric P-E loops were used to determine the energy storage values at room temperature. The validity of magnetoelectric coupling in all samples was confirmed by magneto-dielectric studies and found that the sample $x$ = 0.24 shows the maximum response with the coupling coefficient ($\gamma$) = 15.54 g²/emu².This work promotes the room temperature energy storage properties of the multiferroics. In this approach, impacts of PrFeO₃ doping on PT-based solid solutions (Pb${}_{1-}$${}_x$Pr${}_x$Ti${}_{1-}$${}_x$Fe${}_x$O₃, $x$ = 0.21, 0.22, 0.23, 0.24, 0.25 and 0.26) have been explored. X-ray diffraction (XRD) patterns were used to estimate the crystallographic parameters, confirming the single phase tetragonal structure. The ferroelectric Curie temperature ($T_c$${}^{\mathrm{\text{FE}}}$) is observed to drop from 410 K to below room temperature as the Pr concentration increases. The ferroelectric P-E loops were used to determine the energy storage values at room temperature. The sample $x$ = 0.24 achieved the maximum value of energy storage density of 362.25 mJ/cm³ with the efficiency of 40.5%. The ferroelectric P-E loops were used to determine the energy storage values at room temperature. The validity of magnetoelectric coupling in all samples was confirmed by magneto-dielectric studies and found that the sample $x$ = 0.24 shows the maximum response with the coupling coefficient ($\gamma$) = 15.54 g²/emu².