Our results demonstrate that this framework effectively reproduces free energy profiles, solvation structures, and vibrational spectra in aqueous environments, showing good agreement with QM/MM simulations. Additionally, the method successfully captures key features in protein-ligand interactions, with ML/MM predictions aligning closely with QM/MM calculations for electrostatic energies, particularly in the avidin complex.