Solvation effects can play a substantial role in influencing reaction rates, mechanisms, and selectivity of chemical reactions . Incorporating solvation effects into RMLPs is challenging, especially for explicit solvation effects. One possible approach is to utilize RMLPs in conjunction with QM/MM methods by accelerating QM calculations with RMLPs. The main challenge lies in describing the electrostatic coupling between QM and MM parts (E elec QM/MM term in eq1 and eq2), which is usually handled with mechanical (ME) or electrostatic embedding (EE). The EE scheme is more accurate since the QM region is polarized by the MM region, while the ME scheme omits the polarization between QM and MM regions and simply considers the E elec QM/MM term at a classical level. In this section, we will introduce different ways where ML is hybridized with QM/MM with consideration of how the electrostatic coupling is treated. Meanwhile, we will also include several cases on pure RMLPs that have been trained with solvated systems against expensive AIMD data. . In addition, other systems such as thiol disulfide exchange in proteins and the chorismate mutase enzyme reaction have been simulated with the ∆ML corrected SE-QM/MM approach.