Correlating the 31P NMR chemical shielding tensor and the 2J(P,C) spin-spin coupling constants with torsion angles ζ and α in the backbone of nucleic acids

Exploring the impact of alignment media on RDC analysis of phosphorus-containing compounds: a molecular docking approach

Residual dipolar couplings (RDCs) are employed in NMR analysis when conventional methods, such as J-couplings and nuclear Overhauser effects (NOEs) fail. Low-energy (optimized) conformers are often used as input structures in RDC analysis programs. However, these low-energy structures do not necessarily resemble conformations found in anisotropic environments due to interactions with the alignment medium, especially if the analyte molecules are flexible. Considering interactions with alignment media in RDC analysis, we developed and evaluated a molecular docking-based approach to generate more accurate conformer ensembles for compounds in the presence of the poly-γ-benzyl-L-glutamate alignment medium. We designed chiral phosphorus-containing compounds that enabled us to utilize 31P NMR parameters for the stereochemical analysis. Using P3D/PALES software to evaluate diastereomer discrimination, we found that our conformer ensembles outperform moderately the standard, low-energy conformers in RDC analysis. To further improve our results, we (i) averaged the experimental values of the molecular docking-based conformers by applying the Boltzmann distribution and (ii) optimized the structures through normal mode relaxation, thereby enhancing the Pearson correlation factor R and even diastereomer discrimination in some cases. Nevertheless, we presume that significant differences between J-couplings in isotropic and in anisotropic environments may preclude RDC measurements for flexible molecules. Therefore, generating conformer ensembles based on molecular docking enhances RDC analysis for mildly flexible systems while flexible molecules may require applying more advanced approaches, in particular approaches including dynamical effects.

Probabilistic Interpretation of NMR J-Couplings Determines BI-BII State Equilibria in DNA

Interpretation of ³J_P,H3' NMR scalar spin-spin coupling constants in DNA becomes more reliable by including distinct structural states such as BI and BII, using the weighted-static or, better still, the recently implemented adiabatic-MD (Ad-MD) method. The calculation method employs an adiabatic ("Ad") dependence of ³J_P,H3' coupling on NMR-assigned torsion angle, ε, weighted by P(ε) probability distribution calculated by molecular dynamics (MD). Ad-MD calculations enable cross-validation of the bsc1, OL15, and OL21 force fields and various parametrizations of the Karplus equation describing the dependence of ³J_P,H3' coupling on ε torsion (KE). The mean absolute deviation of Ad-MD ³J_P,H3' couplings from the experimental values in Dickerson-Drew DNA is comparable to the scatter of ³J_P,H3' couplings among four separate NMR experiments. A commonly accepted assumption of homogeneity of one kind of structure-dynamic state within DNA (BI or BII) is questionable because the principal characteristics of relevant P(ε) probabilities (shapes and positioning) vary with DNA sequence. The theory outlined in the present work sets limits to future reparameterization of MD force fields, as relevant to NMR data.

31P NMR parameters may facilitate the stereochemical analysis of phosphorus-containing compounds

Conventional Nuclear Magnetic Resonance (NMR) analysis relies on H-H/C-H interactions. However, these interactions are sometimes insufficient for an accurate and precise NMR analysis. In this study, we show that ³¹P NMR parameters can provide critical structural insights into the stereochemistry of phosphorus-containing compounds. For this purpose, we prepared a set of model phosphorus-based proline derivatives, separated diastereoisomers, and determined their absolute configuration by single-crystal X-ray diffraction. After supplementing these results by electronic circular dichroism (ECD) spectroscopy, we combined experimental data and DFT calculations from our model compounds to perform a detailed conformational analysis, thereby determining their relative configuration. Overall, our findings establish an experimental paradigm for combining ³¹P NMR spectroscopy with other optical methods to facilitate the stereochemical analysis of phosphorus-containing compounds.

The Ad-MD method to calculate NMR shift including effects due to conformational dynamics: The 31 P NMR shift in DNA

A method for averaging of NMR parameters by molecular dynamics (MD) has been derived from the method of statistical averaging in MD snapshots, benchmarked and applied to structurally dynamic interpretation of the ³¹ P NMR shift (δ_31P ) in DNA phosphates. The method employs adiabatic dependence of an NMR parameter on selected geometric parameter(s) that is weighted by MD-calculated probability distribution(s) for the geometric parameter(s) (Ad-MD method). The usage of Ad-MD for polymers is computationally convenient when one pre-calculated structural dependence of an NMR parameter is employed for all chemically equivalent units differing only in dynamic behavior. The Ad-MD method is benchmarked against the statistical averaging method for δ_31P in the model phosphates featuring distinctively different structures and dynamic behavior. The applicability of Ad-MD is illustrated by calculating ³¹ P NMR spectra in the Dickerson-Drew DNA dodecamer. δ_31P was calculated with the B3LYP/IGLO-III/PCM(water) and the probability distributions for the torsion angles adjacent to the phosphorus atoms in the DNA phosphates were calculated using the OL15 force field.

Recent advances in computational 31 P NMR: Part 2. Spin-spin coupling constants

This is the second part of two closely related reviews dealing with the computation of ³¹ P nuclear magnetic resonance (NMR) parameters in a wide range of phosphorous containing compounds. The first part of this review concentrated primarily on the computation of ³¹ P NMR chemical shifts, whereas the second part concerns the calculation of spin-spin coupling constants involving phosphorus nucleus, focusing primarily on their stereochemical dependencies and stereodynamic behavior in particular classes of organophosphorus compounds. This review is dedicated to the Full Member of the Russian Academy of Sciences Professor Boris A. Trofimov in view of his invaluable contribution to the field of synthesis, NMR, and computation studies of organophosphorus compounds.

Toward the comprehensive calculations on the relationship between 1 H, 13 C, 31 P chemical shifts, 2 JPH , and the bonding structure of different phosphoryl benzamides

A comprehensive investigation was performed on ¹ H, ¹³ C, and ³¹ P nuclear magnetic resonance (NMR) chemical shifts (CSs) of phosphoryl benzamide derivatives (C₆ H₅ C(O)NHP(O)R₁ R₂ ), (R₁ , R₂  = aziridine [L₁ ], azetidine [L₂ ], pyrrolidine [L₃ ], piperidine [L₄ ], azepane [L₅ ], 4-methylpiperidine [L₆ ], propane-2-amine [L₇ ], and 2-methylpropane-2-amine [L₈ ]) by the gauge-independent atomic orbital method (GIAO) to find the most accordant level of theory with the experimental values. To achieve this goal, all the structures were optimized using the B3LYP, BP86, PBE1PBE, M06-2X, MPWB1K, and MP2 methods with 6-31+G* basis set. Computed structural parameters demonstrate that BP86 has the best agreement to the experimental values between the other methods. The def2-TZVP and aug-cc-pVDZ basis sets were also employed to inspect the effect of different types of basis sets with higher polarization and diffuse functions. The correlation between the empirical and computational values attests that 6-31+G* basis set is the optimum case regarding minimization of the costs and results. The comparison between calculated and experimental CSs at all mentioned combinations illustrated that in accordance with structural results, the best level of theory in CSs is also BP86/6-31+G*. Besides, ² J_PH values were computed with an acceptable agreement to experimental data at the optimum level of theory. The dependency between ² J_PH and the bonding structure of studied ligands was also scrutinized by the Natural Bond Orbital (NBO) analysis that interprets the relationship between the electronic properties and ² J_PH values.

One- and Two-Dimensional High-Resolution NMR from Flat Surfaces

Determining atomic-level characteristics of molecules on two-dimensional surfaces is one of the fundamental challenges in chemistry. High-resolution nuclear magnetic resonance (NMR) could deliver rich structural information, but its application to two-dimensional materials has been prevented by intrinsically low sensitivity. Here we obtain high-resolution one- and two-dimensional ³¹P NMR spectra from as little as 160 picomoles of oligonucleotide functionalities deposited onto silicate glass and sapphire wafers. This is enabled by a factor >10⁵ improvement in sensitivity compared to typical NMR approaches from combining dynamic nuclear polarization methods, multiple-echo acquisition, and optimized sample formulation. We demonstrate that, with this ultrahigh NMR sensitivity, ³¹P NMR can be used to observe DNA bound to miRNA, to sense conformational changes due to ion binding, and to follow photochemical degradation reactions.

Structural interpretation of the 31P NMR chemical shifts in thiophosphate and phosphate: key effects due to spin–orbit and explicit solvent

Structural interpretation of the ³¹P NMR shifts measured in different molecules including thiophosphate or phosphate group was obtained by means of theoretical calculations including the effects of geometry, molecular dynamics, solvent, relativistic effects and the effect of NMR reference.

The benchmark of 31P NMR parameters in phosphate: a case study on structurally constrained and flexible phosphate

A benchmark for structural interpretation of the ³¹P NMR shift and the ²J_P,C spin–spin coupling in the phosphate group was obtained by means of theoretical calculations and measurements in diethylphosphate and 5,5-dimethyl-2-hydroxy-1,3,2-dioxaphosphinane 2-oxide.

Quantum chemical calculations of31P NMR chemical shifts: scopes and limitations

High level of theory is not necessarily needed to obtain rather accurate predictions of³¹P chemical shifts by GIAO method. For example, the PBE1PBE/6-311G(2d,2p)//PBE1PBE/6-31+G(d) combination allowed to obtain good results for variety of middle-size organophosphorus compounds (M= 200–700 Da).

Toward Improved Description of DNA Backbone: Revisiting Epsilon and Zeta Torsion Force Field Parameters

We present a refinement of the backbone torsion parameters ε and ζ of the Cornell et al. AMBER force field for DNA simulations. The new parameters, denoted as εζOL1, were derived from quantum-mechanical calculations with inclusion of conformation-dependent solvation effects according to the recently reported methodology (J. Chem. Theory Comput. 2012, 7(9), 2886-2902). The performance of the refined parameters was analyzed by means of extended molecular dynamics (MD) simulations for several representative systems. The results showed that the εζOL1 refinement improves the backbone description of B-DNA double helices and G-DNA stem. In B-DNA simulations, we observed an average increase of the helical twist and narrowing of the major groove, thus achieving better agreement with X-ray and solution NMR data. The balance between populations of BI and BII backbone substates was shifted towards the BII state, in better agreement with ensemble-refined solution experimental results. Furthermore, the refined parameters decreased the backbone RMS deviations in B-DNA MD simulations. In the antiparallel guanine quadruplex (G-DNA) the εζOL1 modification improved the description of non-canonical α/γ backbone substates, which were shown to be coupled to the ε/ζ torsion potential. Thus, the refinement is suggested as a possible alternative to the current ε/ζ torsion potential, which may enable more accurate modeling of nucleic acids. However, long-term testing is recommended before its routine application in DNA simulations.

Structure, Stiffness and Substates of the Dickerson-Drew Dodecamer

The Dickerson-Drew dodecamer (DD) d-[CGCGAATTCGCG]₂ is a prototypic B-DNA molecule whose sequence-specific structure and dynamics have been investigated by many experimental and computational studies. Here, we present an analysis of DD properties based on extensive atomistic molecular dynamics (MD) simulations using different ionic conditions and water models. The 0.6-2.4-µs-long MD trajectories are compared to modern crystallographic and NMR data. In the simulations, the duplex ends can adopt an alternative base-pairing, which influences the oligomer structure. A clear relationship between the BI/BII backbone substates and the basepair step conformation has been identified, extending previous findings and exposing an interesting structural polymorphism in the helix. For a given end pairing, distributions of the basepair step coordinates can be decomposed into Gaussian-like components associated with the BI/BII backbone states. The nonlocal stiffness matrices for a rigid-base mechanical model of DD are reported for the first time, suggesting salient stiffness features of the central A-tract. The Riemann distance and Kullback-Leibler divergence are used for stiffness matrix comparison. The basic structural parameters converge very well within 300 ns, convergence of the BI/BII populations and stiffness matrices is less sharp. Our work presents new findings about the DD structural dynamics, mechanical properties, and the coupling between basepair and backbone configurations, including their statistical reliability. The results may also be useful for optimizing future force fields for DNA.