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.

Polarizable Coarse-Grained Force Field for Noble Metal Nanoparticles: Theory and Applications

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

Polarizable Potentials For Metals: The Density Readjusting Embedded Atom Method (DR-EAM)

Balancing Group 1 Monoatomic Ion-Polar Compound Interactions in the Polarizable Drude Force Field: Application in Protein and Nucleic Acid Systems

Ionic structure around polarizable metal nanoparticles in aqueous electrolytes

Physically motivated, robust, ab initio force fields for CO2 and N2.

Phase Field Simulations of Ferroelectric Nanoparticles with Different Long-Range-electrostatic and -Elastic Interactions

Polarisable force fields: what do they add in biomolecular simulations?

Highly Heterogeneous Polarization and Solvation of Gold Nanoparticles in Aqueous Electrolytes

Metallicious : Automated Force-Field Parameterization of Covalently Bound Metals for Supramolecular Structures

Repulsive Forces Between Neutral Surfaces Induced by Adatoms

The electrostatic potential of atomic nanostructures on a metal surface

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

Polarizable Force Fields Based On Physical Models And Quantum Chemical Calculations

Understanding the behavior of exchange bias and coercive fields and connections with anisotropic interfaces in nanoparticle system

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

Interaction between charge-regulated metal nanoparticles in an electrolyte solution

metallicious: Automated force-field parametrization of covalently bound metals for supramolecular structures

DFT-based force field development for noble gas adsorption in metal organic frameworks

Evaluation of Force Fields for Molecular Dynamics Simulations of Platinum in Bulk and Nanoparticle Forms

An Empirical Polarizable Force Field Based on the Classical Drude Oscillator Model: Development History and Recent Applications