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.

Optimal Solution to the Torsional Coefficient Fitting Problem in Force Field Parametrization

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

Coupled Two-Dimensional Main-Chain Torsional Potential for Protein Dynamics II: Performance and Validation.

Construction of atomistic force field model for perfluorosulfonic acid polymer

Improved Peptide and Protein Torsional Energetics with the OPLSAA Force Field

Improved Peptide and Protein Torsional Energetics with the OPLS-AA Force Field

Many-body Effects Are Essential in a Physically Motivated CO2 Force Field

System-Specific Parameter Optimization for Nonpolarizable and Polarizable Force Fields

System-specific parameter optimization for non-polarizable and polarizable force fields

Open Force Field BespokeFit: Automating Bespoke Torsion Parametrization at Scale

Methodology of Parameterization of Molecular Mechanics Force Field From Quantum Chemistry Calculations using Genetic Algorithm: A case study of methanol

Fine-Tuning Property Domain Weighting Factors and the Objective Function in Force Field Parameter Optimization

Optimizations of force-field parameters for protein systems with the secondary-structure stability and instability

Practical Methods for Including Torsional Anharmonicity in Thermochemical Calculations on Complex Molecules: the Internal-Coordinate Multi-Structural Approximation.

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

Constant-pH Simulations with the Polarizable Atomic Multipole AMOEBA Force Field

Atomic permutationally invariant polynomials for fitting molecular force fields

Universal Implementation of a Residue-Specific Force Field Based on CMAP Potentials and Free Energy Decomposition

Non-Bonded Force Field for the Interaction Between Metals and Organic Molecules: A Case Study of Olefins on Aluminum

Symmetrized Drude Oscillator Force Fields Improve Numerical Performance of Polarizable Molecular Dynamics

Residue-Speci fi c Force Field Based on the Protein Coil Library