We compared our results with calculations using the CANDLE solvation model, which belongs to the class of non-local linearized Poisson-Boltzmann models. We used the ORR as a case study to assess the performance of the models for catalytic reactions. Our results indicate that current RISM parametrizations do not enhance formation energies when compared to explicit water models. Both RISM water models predict destabilization of ORR intermediates, while explicit models suggest stabilization of most species. In contrast, the CANDLE model, in line with similar implicit models, predicts small stabilization. This suggests room for improvement in the parametrization of O and H DFT solutes within the RISM framework, potentially leading to more accurate predictions of ORR intermediates formation energies.