Abstract:A large diversity in the targeted hydration free energies (HFEs) during model parameterization of metal ions was reported in the literature with a difference by dozens of kcal/mol. Here, we developed a series of nonbonded dummy models of the Mg2+ ion targeting different HFE references in TIP3P water, followed by assessments of the designed models in the simulations of MgCl2 solution and biological systems. Together with the comparison of existing models, we conclude that the difference in the targeted HFEs has a limited influence on the model performance, while the usability of these models differs from case to case. The feasibility of reproducing more properties of Mg2+ such as diffusion constants and water exchange rates using a nonbonded dummy model is demonstrated. Underestimated activity derivative and osmotic coefficient of MgCl2 solutions in high concentration reveal a necessity for further optimization of ion-pair interactions. The developed dummy models are applicable to metal coordination with Asp, Glu, and His residues in metalloenzymes, and the performance in predicting monodentate or bidentate binding modes of Asp/Glu residues depends on the complexity of metal centers and the choice of protein force fields. When both the binding modes coexist, the nonbonded dummy models outperform point charge models, probably in need of considering polarization of metal-binding residues by, for instance, charge calibration in classical force fields. This work is valuable for the use and further development of magnesium ion models for simulations of metal-containing systems with good accuracy.The Supporting Information is available free of charge at <a class="ext-link" href="/doi/10.1021/acs.jcim.1c00281?goto=supporting-info">https://pubs.acs.org/doi/10.1021/acs.jcim.1c00281</a>.Description of simulated systems (Table S1); calculated real HFEs (Table S2); IOD and CN in the first hydration shell (Table S3); IOD and CN in the second hydration shell (Table S4); diffusion coefficient (Table S5); free energy of activation for water exchange (Table S6) in the LJ parameter space, ion–water, and ion–ion structural characterization of the MgCl2 solution (Table S7); calculated activity derivative (Table S8); osmotic coefficient (Table S9) of MgCl2 solutions, simulated interatomic distances between Mg2+ and ligands for protein <a href="https://www.rcsb.org/pdb/search/structidSearch.do?structureId=1M0W">1M0W</a> (Table S10), <a href="https://www.rcsb.org/pdb/search/structidSearch.do?structureId=1PYX">1PYX</a> (Table S11), 4WH3 (Table S12), and <a href="https://www.rcsb.org/pdb/search/structidSearch.do?structureId=5TKA">5TKA</a> (Table S13); protocol validation of osmotic pressure calculation (Figure S1); comparison of simulated MgCl2 density with experiment (Figure S2); and the RMSDs of RNA <a href="https://www.rcsb.org/pdb/search/structidSearch.do?structureId=1Y26">1Y26</a> (Figure S3) and proteins <a href="https://www.rcsb.org/pdb/search/structidSearch.do?structureId=1M0W">1M0W</a> (Figure S4), <a href="https://www.rcsb.org/pdb/search/structidSearch.do?structureId=1PYX">1PYX</a> (Figure S5), <a href="https://www.rcsb.org/pdb/search/structidSearch.do?structureId=4WH3">4WH3</a> (Figure S6), and <a href="https://www.rcsb.org/pdb/search/structidSearch.do?structureId=5TKA">5TKA</a> (Figure S7); and the LJ parameters of Mg2+ models (<a class="ext-link" href="/doi/suppl/10.1021/acs.jcim.1c00281/suppl_file/ci1c00281_si_001.pdf">PDF</a>)Simulation parameters (<a class="ext-link" href="/doi/suppl/10.1021/acs.jcim.1c00281/suppl_file/ci1c00281_si_002.zip">ZIP</a>)This article has not yet been cited by other publications.

Optimized Magnesium Force Field Parameters for Biomolecular Simulations with Accurate Solvation, Ion-Binding, and Water-Exchange Properties

Optimized Magnesium Force Field Parameters for Biomolecular Simulations with Accurate Solvation, Ion-Binding, and Water-Exchange Properties in SPC/E, TIP3P-fb, TIP4P/2005, TIP4P-Ew, and TIP4P-D

Extended magnesium and calcium force field parameters for accurate ion-nucleic acid interactions in biomolecular simulations.

Magnesium Force Fields for OPC Water with Accurate Solvation, Ion-Binding, and Water-Exchange Properties: Successful Transfer from SPC/E

Including Explicit Water Molecules as Part of the Protein Structure in MM/PBSA Calculations.

Force fields for monovalent and divalent metal cations in TIP3P water based on thermodynamic and kinetic properties

Developing and Assessing Nonbonded Dummy Models of Magnesium Ion with Different Hydration Free Energy References

Balancing Group I Monoatomic Ion-Polar Compound Interactions for Condensed Phase Simulation in the Polarizable Drude Force Field

Taking into Account the Ion-Induced Dipole Interaction in the Nonbonded Model of Ions

Polarizable AMOEBA Model for Simulating Mg2+·Protein·Nucleotide Complexes

Force fields for divalent cations based on single-ion and ion-pair properties

Improved Parameterization of Amine-Carboxyate and Amine-Phosphate Interactions for Molecular Dynamics Simulations Using the CHARMM and AMBER Force Fields

Developing Force Fields for Accurate Mg2+ and Triphosphate Interactions in ATP∙Mg2+ Complexes

Improved model of hydrated calcium ion for molecular dynamics simulations using classical biomolecular force fields

Development of Accurate Force Fields for Mg2+ and Triphosphate Interactions in ATP·Mg2+ and GTP·Mg2+ Complexes.

Developing Force Fields for Accurate Mg2+ and Triphosphate Interactions in ATP∙Mg2+ and GTP∙Mg2+ Complexes

Ionic force field optimization based on single-ion and ion-pair solvation properties

Coarse-Grained Parameters for Divalent Cations within the SIRAH Force Field

Optimized Lennard-Jones Parameters for Druglike Small Molecules.

Protein Simulation Using Coarse-Grained Two-Bead Multipole Force Field with Polarizable Water Models.

AUTOMATED FORCE FIELD PARAMETERIZATION FOR NON-POLARIZABLE AND POLARIZABLE ATOMIC MODELS BASED ON AB INITIO TARGET DATA.