[논문 리뷰] Molecular origin of 31P-NMR chemical shifts of phosphate groups with bivalent counter ions

Leo Christanell, Karl-Jakob König외 3명|arXiv (Cornell University)|2026. 02. 06.

Chemical Reaction Mechanisms인용 수 0

한 줄 요약

본 연구는 강화 샘플링 MD와 ab initio 계산을 31P-NMR 측정과 결합하여 디메틸 포스페이트의 31P 방위가 Mg2+ 또는 Ca2+와의 접촉 이온쌍 형성에서 주로 기인함을 보여주고, 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).

연구 동기 및 목표

• 모형 시스템(DMP–)에서 이원당 이온이 인산기의 31P 화학적 시프트에 어떤 영향을 미치는지 이해한다.
• δiso(31P)에 대한 CIP 대 SSIP 기여를 식별하고 이온 함량과 온도에 따른 CIP 분율을 정량화한다.
• 인접 수화 구조 및 기하학에 대한 의존성과 함께 δiso(31P) 시프트의 미시적 기원을 설명한다.
• NMR 관찰에 근거하여 포스페이트 포스필드를 개선하기 위한 벤치마크를 제공한다.

제안 방법

• Mg2+/Ca2+를 따른 DMP–의 PMF를 매핑하기 위해 강화 샘플링 분자동역학(Umbrella sampling)을 결합한다.
• MD 및 수화 껍질에서 추출한 CIP 및 SSIP 기하에서 대규모 기저집합으로 GIAO-DF-LMP2 수준에서 기초 σiso(31P) 계산을 수행한다.
• 신뢰할 수 있는 계산 프로토콜을 결정하기 위해 기저집합 수렴성과 수화 껍질 의존성을 조사한다.
• DMP– 기하학적 매개변수(P–O 결합 길이)의 함수로 σiso(31P)를 매핑하여 순간 기하학으로부터 시프트를 예측한다.
• 이온 함량과 온도에 따라 δiso(31P) 측정에서 CIP 분율 X[CIP]를 추출하기 위해 이성분 모델을 사용한다.

실험 결과

연구 질문

• RQ1 Mg2+ 및 Ca2+와의 상호작용 시 DMP–에서 δiso(31P) 시프트의 미시적 기원은 무엇인가?
• RQ2 다양한 이온 함량 및 온도에서 관찰된 31P-NMR 시프트에 CIP와 SSIP 인구가 어떻게 기여하는가?
• RQ3 σiso(31P)를 지배하는 공간적 및 기하학적 요인(수화 껍질, DMP– 결합 길이)은 무엇인가?
• RQ4 기하학적 매개변수로부터 ab initio 차폐 상수를 신뢰하게 예측하여 실험 시프트를 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 결합 길이를 이용한 예측 기하지도는 좋은 정확도로 ab initio σiso(31P)을 재현할 수 있으며(기상상 reference에서 R ≈ 0.92).
• 온도 의존성과 이온 함량 의존성은 CIP 분율 X[CIP]를 추출하게 하며, 어떤 영역에서는 용매로 분리된 이온쌍 또는 자유 이온쌍이 에너지적으로 우세한 것으로 보인다.