Abstract
Magnesium plays a vital role in a large variety of biological processes. To model such processes by molecular dynamics simulations, researchers rely on accurate force field parameters for Mg2+ and water. OPC is one of the most promising water models yielding an improved description of biomolecules in water. The aim of this work is to provide force field parameters for Mg2+ that lead to accurate simulation results in combination with OPC water. Using twelve different Mg2+ parameter sets, that were previously optimized with different water models, we systematically assess the transferability to OPC based on a large variety of experimental properties. The results show that the Mg2+ parameters for SPC/E are transferable to OPC and closely reproduce the experimental solvation free energy, radius of the first hydration shell, coordination number, activity derivative, and binding affinity toward the phosphate oxygens on RNA. Two optimal parameter sets are presented: MicroMg yields water exchange in OPC on the microsecond timescale in agreement with experiments. NanoMg yields accelerated exchange on the nanosecond timescale and facilitates the direct observation of ion binding events for enhanced sampling purposes.
Publisher
Cold Spring Harbor Laboratory
Reference68 articles.
1. 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-D1;J. Chem. Theory Comput.,2022
2. On the transferability of ion parameters to the TIP4P/2005 water model using molecular dynamics simulations;J. Chem. Phys.,2020
3. Kinetic pathways of water exchange in the first hydration shell of magnesium: Influence of water model and ionic force field;J Chem. Phys.,2021
4. Simulating water with rigid non-polarizable models: a general perspective
5. A systematic comparison of the structural and dynamic properties of commonly used water models for molecular dynamics simulations;J. Chem. Inf. Model,2021