[Paper Review] Molecular origin of 31P-NMR chemical shifts of phosphate groups with bivalent counter ions

Leo Christanell, Karl-Jakob König+3 more|arXiv (Cornell University)|Feb 6, 2026

Chemical Reaction Mechanisms0 citations

TL;DR

The study combines enhanced sampling MD and ab initio calculations with 31P-NMR measurements to show that 31P shielding changes in dimethyl phosphate arise mainly from contact ion pair formation with Mg2+ or Ca2+, enabling quantification of CIP populations and their temperature dependence.

ABSTRACT

The electrostatic interactions of phosphate groups and counter ions critically affect the structure, function and reactivity of DNA or RNA. We present a joint experimental-theoretical investigation of dimethyl phosphate (DMP-) in aqueous solution, an established model system of the sugar-phosphate backbone. Utilizing 31P-NMR spectroscopy as probe of phosphate-ion association, variations of Mg2+ and Ca2+ content exhibit a systematic shielding of the 31P chemical shift (δiso(31P)) with moderate temperature dependence. Enhanced sampling molecular dynamics (MD) and ab initio (GIAO-DF-LMP2) level of theory are used to reveal the microscopic mechanism. Simulations are performed for a configurational ensemble of DMP-ion geometries and their first solvation shells, demonstrating (i) the spatial convergence of changes of the nuclear shielding constant σiso(31P), (ii) the intramolecular geometric origin of short-timescale σiso(31P) fluctuations and (iii) an average shift of σiso(31P) of about 3-5 ppm upon contact ion pair formation with Mg2+ or Ca2+ ions. A quantitative analysis of δiso(31P) for varying ion content and temperature allows us to extract the temperature-dependent fraction of the contact ion pair species, indicating that solvent separated or free ion pairs are the energetically preferred species. The results impose boundary conditions for improvements of phosphate ion force fields and establish the interactions underlying the changes of δiso(31P).

Motivation & Objective

• Understand how bivalent counter ions influence the 31P chemical shift of phosphate groups in a model system (DMP–).
• Identify CIP vs SSIP contributions to δiso(31P) and quantify CIP fractions as a function of ion content and temperature.
• Elucidate the microscopic origin of δiso(31P) shifts and their dependence on local hydration structure and geometry.
• Provide benchmarks to improve phosphate force fields based on NMR observations.

Proposed method

• Combine enhanced sampling molecular dynamics (Umbrella sampling) to map the PMF of DMP– with Mg2+/Ca2+ along the DMP– O1-P– Mg2+ coordination coordinate.
• Perform ab initio σiso(31P) calculations at the GIAO-DF-LMP2 level with large basis sets on ensembles of CIP and SSIP geometries extracted from MD and hydration shells.
• Investigate basis-set convergence and hydration-shell dependence to determine reliable computational protocols.
• Map σiso(31P) as a function of DMP– geometrical parameters (P–O bond lengths) to predict shifts from instantaneous geometry.
• Use a two-component model to extract CIP fractions X[CIP] from δiso(31P) measurements across ion contents and temperatures.

Experimental results

Research questions

• RQ1 What is the microscopic origin of δiso(31P) shifts in DMP– upon interaction with Mg2+ and Ca2+?
• RQ2 How do CIP and SSIP populations contribute to observed 31P-NMR shifts across varying ion content and temperature?
• RQ3 What are the spatial and geometric factors (hydration shell, DMP– bond lengths) that govern σiso(31P)?
• RQ4 Can ab initio shielding constants be reliably predicted from geometrical parameters to map experimental shifts to CIP formation?

Key findings

• δiso(31P) shifts are systematically shielded upon CIP formation with Mg2+ or Ca2+, with an average shift of about 3–5 ppm.
• σiso(31P) converges with a hydration shell extending ~4.5–5 Å around the CIP, indicating a localized electronic response within this range.
• The width of σiso(31P) distributions along DMP––Mg2+/–separation is about 10–12 ppm, driven by instantaneous DMP– geometry, especially P–O bond lengths.
• CIP formation leads to larger shielding compared to free DMP–, and the correlation between geometry and σiso(31P) is weakened when CIP is present, highlighting CIP as the dominant factor.
• A predictive geometric map using O1/O2–P and O3/O5–P bond lengths can reproduce ab initio σiso(31P) with good accuracy (R ≈ 0.92 for gas-phase references).
• Temperature dependence and ion-content dependence enable extraction of the CIP fraction X[CIP], suggesting solvent-separated or free ion pairs remain energetically favored in some regimes.