Balancing Group 1 Monoatomic Ion–Polar Compound Interactions in the Polarizable Drude Force Field: Application in Protein and Nucleic Acid Systems
- 3 December 2024
- journal article
- research article
- Published by American Chemical Society (ACS) in The Journal of Physical Chemistry B
- Vol. 128 (49) , 12078-12091
Abstract
An accurate force field (FF) is the foundation of reliable results from molecular dynamics (MD) simulations. In our recently published work, we developed a protocol to generate atom pair-specific Lennard-Jones (known as NBFIX in CHARMM) and through-space Thole dipole screening (NBTHOLE) parameters in the context of the Drude polarizable FF based on readily accessible quantum mechanical (QM) data to fit condensed phase experimental thermodynamic benchmarks, including the osmotic pressure, diffusion coefficient, ionic conductivity, and solvation free energy, when available. In the present work, the developed protocol is applied to generate NBFIX and NBTHOLE parameters for interactions between monatomic ions (specifically Li+, Na+, K+, Rb+, Cs+, and Cl-) and common functional groups found in proteins and nucleic acids. The parameters generated for each ion-functional group pair were then applied to the corresponding functional groups within proteins or nucleic acids followed by MD simulations to analyze the distribution of ions around these biomolecules. The modified FF successfully addresses the issue of overbinding observed in a previous iteration of the Drude FF. Quantitatively, the model accurately reproduces the effective charge of proteins and demonstrates a level of charge neutralization for a double-helix B-DNA in good agreement with the counterion condensation theory. Additionally, simulations involving ion competition correlate well with experimental results, following the trend Li+ > Na+ ≈ K+ > Rb+. These results validate the refined model for group 1 ion-biomolecule interactions that will facilitate the application of the polarizable Drude FF in systems in which group 1 ions play an important role.Keywords
Funding Information
- Center for Scientific Review (R35 GM131710, R35 GM133754, T32 AI195190)
This publication has 69 references indexed in Scilit:
- Artificial Channels Based on Bottlebrush Polymers: Enhanced Ion Transport Through Polymer Topology ControlAngewandte Chemie International Edition, 2024
- Polarizable force field for RNA based on the classical drude oscillatorJournal of Computational Chemistry, 2018
- Comparative Study of Analytical Techniques for Determining Protein ChargeJournal of Pharmaceutical Sciences, 2015
- All‐atom polarizable force field for DNA based on the classical drude oscillator modelJournal of Computational Chemistry, 2014
- CHARMM36 all-atom additive protein force field: Validation based on comparison to NMR dataJournal of Computational Chemistry, 2013
- CHARMM: The biomolecular simulation programJournal of Computational Chemistry, 2009
- CHARMM fluctuating charge force field for proteins: II Protein/solvent properties from molecular dynamics simulations using a nonadditive electrostatic modelJournal of Computational Chemistry, 2004
- CHARMM fluctuating charge force field for proteins: I parameterization and application to bulk organic liquid simulationsJournal of Computational Chemistry, 2003
- An apparatus for membrane‐confined analytical electrophoresisElectrophoresis, 1998
- Water and monovalent ions in the minor groove of B-DNA oligonucleotides as seen by NMRBiopolymers, 1998