文摘
Dependence of NMR 31P shielding tensor and 2J(P,C) coupling constants on solvation of nucleic acid phosphate by Mg2+ and water was studied using methods of bioinformatic structural analyses of crystallographic data and DFT B3LYP calculations of NMR parameters. The effect of solvent dynamics on NMR parameters was calculated using molecular dynamic. The NMR calculations for representative solvation patterns determined in crystals of B-DNA and A-RNA molecules pointed out the crucial importance of local Mg2+ coordination geometry, including hydration by explicit water molecules and necessity of dynamical averaging over the solvent reorientation. The dynamically averaged 31P chemical shift decreased by 2鈭?.5 ppm upon Mg2+ coordination, the chemical shielding anisotropy increased by 0鈭?0 ppm, and the 2J(P,C5鈥? coupling magnitude decreased by 0.2鈭?.8 Hz upon Mg2+ coordination. The calculated decrease of the 31P chemical shift is in excellent agreement with the 1.5鈭?0 ppm decrease of the phosphorothioate 31P chemical shift upon Cd2+ coordination probed experimentally in hammerhead ribozyme (Suzumura; et al. J. Am. Chem. Soc. 2002, 124, 8230鈭?236; Osborne; et al., Biochemistry 2009, 48, 10654鈭?0664). None of the dynamically averaged NMR parameters unequivocally distinguishes the site-specific Mg2+ coordination to one of the two nonesterified phosphate oxygen atoms of the phosphate determined by bioinformatic analyses. By comparing the limit cases of static and dynamically averaged solvation, we propose that mobility of the solvent has a dramatic impact on NMR parameters of nucleic acid phosphate and must be taken into account for their accurate modeling.