Once the force field impact on the active site characterized, we generated 4 long 500 ns NPT simulations starting from the structure obtained at the end of the equilibration protocol, with restraints on the distance between Mg 2+ ions and the non-water oxygen atoms in their hydration sphere. We used here the pair-specific Lennard Jones force field for Mg 2+