Abstract:Molecular dynamics (MD) simulations are a commonly used method for investigating molecular behavior at the atomic level. Achieving reliable MD simulation results necessitates the use of an accurate force field. In the present work, we present a new protocol to enhance the quality of group 1 monoatomic ions (specifically Li+, Na+, K+, Rb+, and Cs+) with respect to their interactions with common polar model compounds in biomolecules in condensed phases in the context of the Drude polarizable force field. Instead of adjusting pre-existing individual parameters for ions, model compounds, and water, we employ atom-pair specific LJ (known as NBFIX in CHARMM) and through-space Thole dipole screening (NBTHOLE) terms to fine-tune the balance of ion-model compound, ion-water, and model compound-water interactions. This involved establishing a protocol for the optimization of NBFIX and NBTHOLE parameters targeting the difference between molecular mechanical (MM) and quantum mechanical (QM) potential energy scans (PES). It is shown that targeting PES involving complexes that include multiple model compounds and/or ions as trimers and tetramers yields parameters that produce condensed phase properties in agreement with experimental data. Validation of this protocol involved the reproduction of experimental thermodynamic benchmarks, including solvation free energies of ions in methanol and N-methyl acetamide, osmotic pressures, ionic conductivities, and diffusion coefficients within the condensed phase. These results show the importance of including more complex ion-model compound complexes beyond dimers in the QM target data to account for many-body effects during parameter fitting. The presented parameters represent a significant refinement of the Drude polarizable force field, which will lead to improved accuracy for modeling ion-biomolecular interactions.

Coupling all-atom molecular dynamics simulations of ions in water with Brownian dynamics

Mode coupling theory of electrolyte dynamics: Time dependent diffusion, dynamic structure factor and solvation dynamics

Quantifying the Structure of Water and Hydrated Monovalent Ions by Density Functional Theory-Based Molecular Dynamics.

Molecular Dynamics in Physiological Solutions: Force Fields, Alkali Metal Ions, and Ionic Strength

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

Molecular dynamics simulation of ionic concentration effects on ionic hydration in aqueous KCl solutions

Understanding the Anomalous Diffusion of Water in Aqueous Electrolytes Using Machine Learned Potentials

Grand-canonical molecular dynamics simulations powered by a hybrid 4D nonequilibrium MD/MC method: Implementation in LAMMPS and applications to electrolyte solutions

A Coupled Molecular Dynamics/Kinetic Monte Carlo Approach for Protonation Dynamics in Extended Systems.

Ion clusters and ion-water potentials in MD-simulations

Study on the transport properties of aqueous electrolyte solution by brownian dynamics simulation

A Coarse-Grained SPICA Makeover for Solvated and Bare Sodium and Chloride Ions

Bridging atomistic simulations and thermodynamic hydration models of aqueous electrolyte solutions

Scalable molecular simulation of electrolyte solutions with quantum chemical accuracy

Consistent Density Functional Theory-Based Description of Ion Hydration Through Density-Corrected Many-Body Representations

Hydration structure and diffusion of ions in supercooled water: Ion size effects

Molecular dynamics study on ionic hydration

Solvation Structures and Ion Dynamics of CaCl$_2$ Aqueous Electrolytes Using Metadynamics and Machine Learning Molecular Dynamics Simulations

Fluctuating hydrodynamics for multiscale simulation of inhomogeneous fluids: mapping all-atom molecular dynamics to capillary waves

Molecular Theories and Simulation of Ions and Polar Molecules in Water

Hydration Structure of Sodium and Potassium Ions with DFT-MD