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.

A novel Monte Carlo algorithm for polarizable force fields: Application to a fluctuating charge model for water

Parametrizing A Polarizable Force Field from Ab Initio Data. I. the Fluctuating Point Charge Model

A Force Consistent Method For Electrostatic Energy Calculation In Fluctuating Charge Model

A Polarizable Valence Electron Density Based Force Field for High-Energy Interactions between Atoms and Molecules

Polarizable Water Model with Ab Initio Neural Network Dynamic Charges and Spontaneous Charge Transfer

Polarizable Mean-Field Model of Water for Biological Simulations with Amber and Charmm force fields

A Transferrable Range-Separated Force Field for Water: Combining the Power of Both Physically-Motivated Models and Machine Learning Techniques

Development of an Effective Polarizable Bond Method for Biomolecular Simulation.

Convergence of stochastic-extended Lagrangian molecular dynamics method for polarizable force field simulation

Multiscale Polarizable Coarse-Graining Water Models on Cluster-Level Electrostatic Dipoles.

Pressure Calculation in Polar and Charged Systems using Ewald Summation: Results for the Extended Simple Point Charge Model of Water

Accounting for electronic polarization in nonpolarizable force fields

Efficient Formulation of Polarizable Gaussian Multipole Electrostatics for Biomolecular Simulations

Fully Polarizable Multiconfigurational Self-consistent Field/Fluctuating Charge Approach

Building Water Models Compatible with Charge Scaling Molecular Dynamics

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

High-Performance Scalable Molecular Dynamics Simulations of a Polarizable Force Field Based on Classical Drude Oscillators in NAMD

A QM/MM Derived Polarizable Water Model for Molecular Simulation

A fast, convenient, polarizable electrostatic model for molecular dynamics

Implementation of extended Lagrangian dynamics in GROMACS for polarizable simulations using the classical Drude oscillator model

Simplistic Coulomb Forces in Molecular Dynamics: Comparing the Wolf and Shifted-Force Approximations