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.

Molecular Dynamics Simulations Of Local Field Factors

Modeling of Dynamic Molecular Solvent Properties Using Local and Cavity Field Approaches

Learning the Quantum Centroid Force Correction in Molecular Systems: A Localized Approach

Hybrid Particle-Continuum Simulations Coupling Brownian Dynamics and Local Dynamic Density Functional Theory

New Formulation for Non-Equilibrium Solvation: Spectral Shifts and Cavity Radii of 6-Propanoyl-2(n,n-dimethylamino) Naphthalene and 4-(N,n-dimethylamino) Benzonitrile

Modelling nonlinear optical effects in guest-host systems

Sph Simulations of 2d Transient Viscoelastic Flows Using Brownian Configuration Fields

Molecular dynamics with quantum forces: Vibrational spectra of localized systems

Molecular Dynamics Simulation for Dmf Aqueous Solution

Auxiliary field simulation and Coulomb's Law

Efficient anharmonic vibrational spectroscopy for large molecules using local-mode coordinates.

Molecular dynamics simulation of a chemical reaction in solution

Local geometry of electromagnetic fields and its role in molecular multipole transitions

FIELD IONIZATION OF MOLECULES IN AN INTENSE LASER FIELD

Simulations of the molecular dynamics of nucleic acids

Franck-Condon Factors Perturbed by Damped Harmonic Oscillators: Solvent Enhanced X 1ag ↔ A1B1u Absorption and Fluorescence Spectra of Perylene

Microscopic Origin of the Deviation from Stokes-Einstein Behavior Observed in Dynamics of the Kscn Aqueous Solutions: A Md Simulation Study

Accelerating Anharmonic Spectroscopy Simulations via Local-Mode, Multilevel Methods.

Molecular Dynamics Simulations of Ion Distribution in Nanochannels

Computation of Polarizability Anisotropies for Molecules in Dilute Solution: an Approach to the Local Field Correction

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