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.

Theoretical study of phase behaviour of DLVO model for lysozyme and $γ$-crystalline aqueous electrolyte solutions

Molecular Dynamics Simulation Studies of Cryoprotective Agent Solutions: the Relation Between Melting Temperature and the Ratio of Hydrogen Bonding Acceptor to Donor Number

A simple patchy colloid model for the phase behavior of lysozyme dispersions

Ion specific effects on phase transitions in protein solutions

Molecular Theories and Simulation of Ions and Polar Molecules in Water

Mean-field theory of the interaction of the magnesium ion with biopolymers: the case of lysozyme

Insights from Molecular Dynamics Simulations on Structural Organization and Diffusive Dynamics of an Ionic Liquid at Solid and Vacuum Interfaces

A Modified Landau-de Gennes Theory for Smectic Liquid Crystals: Phase Transitions and Structural Transitions

Machine learning force fields for molecular liquids: Ethylene Carbonate/Ethyl Methyl Carbonate binary solvent

Study on Osmotic Pressures for Aqueous Lysozyme and Α-Chymotrypsin—electrolyte Solutions with Two Yukawa Potentials

Balancing Local Order and Long-Ranged Interactions in the Molecular Theory of Liquid Water

Study of the Precrystallization Solution of Lysozyme by Accelerated Molecular Dynamics Simulation

Computational investigation of structure, dynamics and nucleation kinetics of a family of modified Stillinger-Weber model fluids in bulk and free-standing thin films

Bridging atomistic simulations and thermodynamic hydration models of aqueous electrolyte solutions

Phase behaviour of primitive models of molecular ionic liquids in porous media: effects of cation shape, ion association and disordered confinement

Condensed-phase molecular representation to link structure and thermodynamics in molecular dynamics

Explicit Solvent Effects on Macromolecular Interactions From a Solvent-Augmented Contact Value Theorem

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

A comparative study of aqueous DMSO mixtures by computer simulations and integral equation theories

Interactions in protein solutions close to liquid-liquid phase separation: Ethanol reduces attractions via changes of the dielectric solution properties

Inelastic X-ray Scattering Studies of the Short-Time Collective Vibrational Motions in Hydrated Lysozyme Powders and Their Possible Relation to Enzymatic Function.