Design of Point Charge Models for Divalent Metal Cations Targeting Quantum Mechanical Ion–Water Dimer Interactions

Divalent metal cations are of vital importance in biochemistry and materials science, and their structural and thermodynamic properties in aqueous solution have often been used as targets for the development of ion models. This study presented a strategy for designing nonbonded point charge models of divalent metal cations (Mg2+ and Ca2+) and Cl− by targeting quantum mechanics (QM)-based ion–water dimer interactions. The designed models offered an accurate representation of ion–water interactions in the gas phase and showed reasonable performance for non-targeted properties in aqueous solutions, such as the ion–water oxygen distance (IOD), coordination number (CN), and density and viscosity of MgCl2 and CaCl2 solutions at low concentrations. Our metal cation models yielded considerable overestimates of the hydration free energies (HFEs) of the ions, whereas the Cl− model displayed good performance. Together with the overestimated density and viscosity of the salt solutions, these results indicated the necessity of re-optimizing ion–ion interactions and/or including polarization effects in the design of ion models. The designed Mg2+ model was capable of maintaining the crystal metal-binding networks during MD simulation of a metalloprotein, indicating great potential for biomolecular simulations. This work highlighted the potential of QM-based ion models to advance the study of metal ion interactions in biological and material systems.