[论文解读] 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 Mechanisms被引用 0

一句话总结

该研究将增强采样MD与从头计算结合，并结合31P-NMR测量，显示二甲基磷酸盐中31P屏蔽的变化主要来自与Mg2+或Ca2+形成的接触离子对(CIP)，从而实现CIP人口及其温度依赖性的量化。

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).

研究动机与目标

• 理解双价对离子如何影响模型体系中磷酸基团的31P化学位移（DMP–）。
• 识别δiso(31P)中CIP与SSIP的贡献，并随离子含量和温度量化CIP分数。
• 阐明δiso(31P)位移的微观起源及其对局部水合结构和几何的依赖。
• 提供基于NMR观测改进磷酸盐力场的基准。

提出的方法

• 结合增强采样分子动力学（Umbrella sampling）以绘制DMP–在Mg2+/Ca2+作用下的PMF，沿DMP– O1-P– Mg2+配位坐标映射。
• 在MD提取的CIP和SSIP几何构型及水合层上，对大基组基于GIAO-DF-LMP2水平进行σiso(31P)从头计算。
• 研究基组收敛性和水合层依赖性，以确定可靠的计算协议。
• 将σiso(31P)作为DMP–几何参数（P–O键长）的函数进行映射，以预测来自瞬时几何的位移。
• 使用两分量模型从δiso(31P)在不同离子含量和温度下的测量中提取CIP分数X[CIP]。

实验结果

研究问题

• RQ1 DMP–在与Mg2+和Ca2+相互作用时，δiso(31P)位移的微观起源是什么？
• RQ2 CIP与SSIP群体在不同离子含量和温度下对31P-NMR位移有何贡献？
• RQ3 决定σiso(31P)的空间和几何因素（水合层、DMP–键长）有哪些？
• RQ4 能否从几何参数可靠预测从头计算的屏蔽常数，以将实验位移映射到CIP形成？

主要发现

• δiso(31P)在与Mg2+或Ca2+形成CIP时系统性产生屏蔽，平均位移约为3–5 ppm。
• σiso(31P)在围绕CIP延展约4.5–5 Å的水合壳时趋于收敛，表明在该半径内存在局部电子响应。
• DMP––Mg2+/–距离沿σiso(31P)分布的宽度约为10–12 ppm，由瞬时DMP–几何，特别是P–O键长驱动。
• CIP形成导致的屏蔽强于游离DMP–，且当存在CIP时几何与σiso(31P)之间的相关性减弱，强调CIP为主导因素。
• 使用O1/O2–P和O3/O5–P键长的预测几何图可以以较高精度再现从头计算的σiso(31P)（气相参考的相关性R约0.92）。
• 温度依赖性和离子含量依赖性使得可提取CIP分数X[CIP]，表明在某些条件下仍偏好溶剂分离或自由离子对的能量状态。