Title: Binding of Phosphate Species to Ca2+ and Mg2+ in Aqueous Solution

Journal Article · · Journal of Chemical Theory and Computation

Phosphate derivatives and their interaction with metal cations are involved in many important biological phenomena, so an accurate characterization of the phosphate-metal interaction is necessary to properly understand the role of phosphate-metal contacts in mediating biological function. Herein, we improved the standard 12-6 Lennard-Jones (LJ) potential via the usage of the 12-6-4 LJ model, which incorporates ion-induced dipole interactions. Via parameter scanning, we fine-tuned the 12-6-4 LJ polarizability values to obtain accurate absolute binding free energies for the phosphate anions H2PO4–, HPO42–, PO43– coordinating with Ca2+ and Mg2+. First, we modified the phosphate 12-6-4 LJ parameters to reproduce the solvation free energies of the series of phosphate anions using the thermodynamic integration (TI) method. Then, using the potential mean force (PMF) method, the polarizability of the metal-phosphate interaction was obtained. We show that the free energy profiles of phosphate ions coordinated to Ca2+ and Mg2+ generally show similar trends at longer metal–phosphate distances, while the absolute binding energy values increased with deprotonation. The resulting parameters demonstrate the flexibility of the 12-6-4 LJ-type nonbonded model and its usefulness in accurately describing cation–anion interactions.

Research Organization:
Lawrence Livermore National Laboratory (LLNL), Livermore, CA (United States)
Sponsoring Organization:
National Institutes of Health (NIH); USDOE
Grant/Contract Number:
AC52-07NA27344
OSTI ID:
2471165
Journal Information:
Journal of Chemical Theory and Computation, Journal Name: Journal of Chemical Theory and Computation Journal Issue: 10 Vol. 20; ISSN 1549-9618
Publisher:
American Chemical SocietyCopyright Statement
Country of Publication:
United States
Language:
English

References (44)

Raman spectroscopic study of the effects of Ca2+, Mg2+, Zn2+, and Cd2+ ions on calf thymus DNA: Binding sites and conformational changes journal January 1990
The M06 suite of density functionals for main group thermochemistry, thermochemical kinetics, noncovalent interactions, excited states, and transition elements: two new functionals and systematic testing of four M06-class functionals and 12 other functionals journal July 2007
Avoiding the integral storage bottleneck in LCAO calculations of electron correlation journal January 1989
A direct MP2 gradient method journal February 1990
Semi-direct algorithms for the MP2 energy and gradient journal February 1990
Analytic MP2 frequencies without fifth-order storage. Theory and application to bifurcated hydrogen bonds in the water hexamer journal March 1994
VMD: Visual molecular dynamics journal February 1996
The effect of metal cations on the phase behavior and hydration characteristics of phospholipid membranes journal May 2002
Multiplicity of metal ion binding patterns to nucleobases journal May 2000
Molecular dynamics simulations of water and biomolecules with a Monte Carlo constant pressure algorithm journal January 2004
Automatic atom type and bond type perception in molecular mechanical calculations journal October 2006
Metal Ion Modeling Using Classical Mechanics journal January 2017
Parameterization of Monovalent Ions for the OPC3, OPC, TIP3P-FB, and TIP4P-FB Water Models journal February 2021
Rational Design of Nonbonded Point Charge Models for Divalent Metal Cations with Lennard-Jones 12-6 Potential journal July 2021
Systematic Parametrization of Divalent Metal Ions for the OPC3, OPC, TIP3P-FB, and TIP4P-FB Water Models journal June 2020
Optimized Magnesium Force Field Parameters for Biomolecular Simulations with Accurate Solvation, Ion-Binding, and Water-Exchange Properties journal March 2021
Parametrization of Trivalent and Tetravalent Metal Ions for the OPC3, OPC, TIP3P-FB, and TIP4P-FB Water Models journal April 2021
Development of Site-Specific Mg 2+ –RNA Force Field Parameters: A Dream or Reality? Guidelines from Combined Molecular Dynamics and Quantum Mechanics Simulations journal December 2016
Thermodynamics of Metal–Acetate Interactions journal January 2024
Balancing the Interactions of Mg 2+ in Aqueous Solution and with Nucleic Acid Moieties For a Polarizable Force Field Based on the Classical Drude Oscillator Model journal October 2016
A Polarizable Cationic Dummy Metal Ion Model journal June 2022
Extending the Nonbonded Cationic Dummy Model to Account for Ion-Induced Dipole Interactions journal October 2017
Charge Model 5: An Extension of Hirshfeld Population Analysis for the Accurate Description of Molecular Interactions in Gaseous and Condensed Phases
  • Marenich, Aleksandr V.; Jerome, Steven V.; Cramer, Christopher J.
  • Journal of Chemical Theory and Computation, Vol. 8, Issue 2, p. 527-541 https://doi.org/10.1021/ct200866d
journal February 2012
Revised AMBER Parameters for Bioorganic Phosphates journal October 2012
Rational Design of Particle Mesh Ewald Compatible Lennard-Jones Parameters for +2 Metal Cations in Explicit Solvent journal May 2013
Taking into Account the Ion-Induced Dipole Interaction in the Nonbonded Model of Ions journal November 2013
Systematic Parameterization of Monovalent Ions Employing the Nonbonded Model journal March 2015
Self-Consistent Reaction Field Model for Aqueous and Nonaqueous Solutions Based on Accurate Polarized Partial Charges journal October 2007
Solubility of magnesium hydrogen phosphate trihydrate and ion-pair formation in the system magnesium hydroxide-phosphoric acid-water at 25.degree.C journal June 1984
Complexes in calcium phosphate solutions journal January 1968
Influence of the Dynamic Positions of Cations on the Structure of the DNA Minor Groove:  Sequence-Dependent Effects journal August 2001
Thermodynamics of Transition Metal Ion Binding to Proteins journal March 2020
Simulating the Chelate Effect journal November 2018
Parameterization of Highly Charged Metal Ions Using the 12-6-4 LJ-Type Nonbonded Model in Explicit Water journal August 2014
Universal Solvation Model Based on Solute Electron Density and on a Continuum Model of the Solvent Defined by the Bulk Dielectric Constant and Atomic Surface Tensions
  • Marenich, Aleksandr V.; Cramer, Christopher J.; Truhlar, Donald G.
  • The Journal of Physical Chemistry B, Vol. 113, Issue 18, p. 6378-6396 https://doi.org/10.1021/jp810292n
journal May 2009
Metal complex interactions with DNA journal January 2015
Thermodynamics of solvation of ions. Part 5.—Gibbs free energy of hydration at 298.15 K journal January 1991
Comparison of simple potential functions for simulating liquid water journal July 1983
A fast and high-quality charge model for the next generation general AMBER force field journal September 2020
Extended magnesium and calcium force field parameters for accurate ion–nucleic acid interactions in biomolecular simulations journal May 2021
Software for the frontiers of quantum chemistry: An overview of developments in the Q-Chem 5 package journal August 2021
Why nature chose phosphates journal March 1987
Unraveling Mg 2+ –RNA binding with atomistic molecular dynamics journal February 2017
Chemical and Biochemical Basis of Cell-Bone Matrix Interaction in Health and Disease journal May 2009