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.

Modeling solvent effects and convergence of 31P‐NMR shielding calculations with COBRAMM

Accelerating NMR Shielding Calculations Through Machine Learning Methods: Application to Magnesium Sodium Silicate Glasses

Revving up 13C NMR shielding predictions across chemical space: Benchmarks for atoms-in-molecules kernel machine learning with new data for 134 kilo molecules

An in-silico NMR laboratory for nuclear magnetic shieldings computed via finite fields: Exploring nucleus-specific renormalizations of MP2 and MP3

Modeling solvation effects in real-space and real-time within Density Functional Approaches

Quantum Calculation of Protein Nmr Chemical Shifts Based on the Automated Fragmentation Method

Protein NMR chemical shift calculations based on the automated fragmentation QM/MM approach.

Systematic studies on the computation of nuclear magnetic resonance shielding constants and chemical shifts: the density functional models.

AFNMR: Automated Fragmentation Quantum Mechanical Calculation of NMR Chemical Shifts for Biomolecules

The importance of nuclear quantum effects for NMR crystallography

Automated Fragmentation QM/MM Calculation of NMR Chemical Shifts for Protein-Ligand Complexes.

Automated Fragmentation QM/MM Calculation of Amide Proton Chemical Shifts in Proteins with Explicit Solvent Model.

Investigating solvent effects on the magnetic properties of molybdate ions (MoO$_{4}^{2-}$) with relativistic embedding

Nonequilibrium Solvent Polarization Effects in Real-Time Electronic Dynamics of Solute Molecules Subject to Time-Dependent Electric Fields: A New Feature of the Polarizable Continuum Model

Tackling solvent effect by coupling electronic and molecular Density Functional Theory

Advanced Modelling of Polyoxometalate Counter-Cations: Bridging Explicit Counter-Cation Inclusion and Solvent Effects for Electronic Applications

The Pseudo-Resonant-Isotope-Model for Predicting the Resonance-Interference Effect in the Self-Shielding Calculation

Constrained Nuclear-Electronic Orbital QM/MM Approach for Simulating Complex Systems with Quantum Nuclear Delocalization Effects Incorporated

Improving Quantum Chemical Solvation Models by Dynamic Radii Adjustment for Continuum Solvation (DRACO)

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

Accurate Prediction of Nuclear Magnetic Resonance Shielding Constants: Towards the Accuracy of CCSD(T) Complete Basis Set Limit