The coupling term can be divided in two contributions, as well as in typical QM/MM frameworks: A non-electrostatic one, E LJ ML/MM , modeled at the MM level using previously assigned Lennard-Jones parameters for atoms of the ML subsystem; and an electrostatic one, E elec ML/MM , modeled similarly to the QM/MM case, but using ML-predicted properties to compute the coulombic interaction.