Structural Interpretation of J Coupling Constants in Guanosine and Deoxyguanosine: Modeling the Effects of Sugar Pucker, Backbone Conformation, and Base Pairing

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.

Impact of modified ribose sugars on nucleic acid conformation and function

The modification of the ribofuranose in nucleic acids is a widespread method of manipulating the activity of nucleic acids. These alterations, however, impact the local conformation and chemical reactivity of the sugar. Changes in the conformation and dynamics of the sugar moiety alter the local and potentially global structure and plasticity of nucleic acids, which in turn contributes to recognition, binding of ligands and enzymatic activity of proteins. This review article introduces the conformational properties of the (deoxy)ribofuranose ring and then explores sugar modifications and how they impact local and global structure and dynamics in nucleic acids.

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.

Empirical Corrections to the Amber RNA Force Field with Target Metadynamics

The computational study of conformational transitions in nucleic acids still faces many challenges. For example, in the case of single stranded RNA tetranucleotides, agreement between simulations and experiments is not satisfactory due to inaccuracies in the force fields commonly used in molecular dynamics simulations. We here use experimental data collected from high-resolution X-ray structures to attempt an improvement of the latest version of the AMBER force field. A modified metadynamics algorithm is used to calculate correcting potentials designed to enforce experimental distributions of backbone torsion angles. Replica-exchange simulations of tetranucleotides including these correcting potentials show significantly better agreement with independent solution experiments for the oligonucleotides containing pyrimidine bases. Although the proposed corrections do not seem to be portable to generic RNA systems, the simulations revealed the importance of the α and ζ backbone angles for the modulation of the RNA conformational ensemble. The correction protocol presented here suggests a systematic procedure for force-field refinement.

Impact of nucleic acid self-alignment in a strong magnetic field on the interpretation of indirect spin-spin interactions

Heteronuclear and homonuclear direct (D) and indirect (J) spin-spin interactions are important sources of structural information about nucleic acids (NAs). The Hamiltonians for the D and J interactions have the same functional form; thus, the experimentally measured apparent spin-spin coupling constant corresponds to a sum of J and D. In biomolecular NMR studies, it is commonly presumed that the dipolar contributions to Js are effectively canceled due to random molecular tumbling. However, in strong magnetic fields, such as those employed for NMR analysis, the tumbling of NA fragments is anisotropic because the inherent magnetic susceptibility of NAs causes an interaction with the external magnetic field. This motional anisotropy is responsible for non-zero D contributions to Js. Here, we calculated the field-induced D contributions to 33 structurally relevant scalar coupling constants as a function of magnetic field strength, temperature and NA fragment size. We identified two classes of Js, namely (1)JCH and (3)JHH couplings, whose quantitative interpretation is notably biased by NA motional anisotropy. For these couplings, the magnetic field-induced dipolar contributions were found to exceed the typical experimental error in J-coupling determinations by a factor of two or more and to produce considerable over- or under-estimations of the J coupling-related torsion angles, especially at magnetic field strengths >12 T and for NA fragments longer than 12 bp. We show that if the non-zero D contributions to J are not properly accounted for, they might cause structural artifacts/bias in NA studies that use solution NMR spectroscopy.

Gas-Phase Conformations and Energetics of Protonated 2'-Deoxyguanosine and Guanosine: IRMPD Action Spectroscopy and Theoretical Studies

The gas-phase structures of protonated 2'-deoxyguanosine, [dGuo+H](+), and its RNA analogue protonated guanosine, [Guo+H](+), are investigated by infrared multiple photon dissociation (IRMPD) action spectroscopy and theoretical electronic structure calculations. IRMPD action spectra are measured over the range extending from ∼550 to 1900 cm(-1) using the FELIX free electron laser and from ∼2800 to 3800 cm(-1) using an optical parametric oscillator/amplifier (OPO/OPA) laser system. The measured IRMPD spectra of [dGuo+H](+) and [Guo+H](+) are compared to each other and to B3LYP/6-311+G(d,p) linear IR spectra predicted for the stable low-energy conformations computed for these species to determine the most favorable site of protonation, identify the structures accessed in the experiments, and elucidate the influence of the 2'-hydroxyl substituent on the structures and the IRMPD spectral features. Theoretical energetics and the measured IRMPD spectra find that N7 protonation is preferred for both [dGuo+H](+) and [Guo+H](+), whereas O6 and N3 protonated conformers are found to be much less stable. The 2'-hydroxyl substituent does not exert a significant influence on the structures and relative stabilities of the stable low-energy conformations of [dGuo+H](+) versus [Guo+H](+) but does provide additional opportunities for hydrogen bonding such that more low-energy structures are found for [Guo+H](+). [dGuo+H](+) and [Guo+H](+) share very parallel IRMPD spectral features in the FELIX and OPO regions, whereas the effect of the 2'-hydroxyl substituent is primarily seen in the relative intensities of the measured IR bands. The measured OPO/OPA spectral signatures, primarily reflecting the IR features associated with the O-H and N-H stretches, provide complementary information to that of the FELIX region and enable the conformers that arise from different protonation sites to be more readily distinguished. Insight gained from this and parallel studies of other DNA and RNA nucleosides and nucleotides should help better elucidate the effects of the local environment on the overall structures of DNA and RNA.

Restricted puckering of mineralized RNA-like riboses

The pseudorotational motions of highly hydroxylated pentafuranose sugars in the free state and tethered to hydroxyapatite have been compared. The conformation pentafuranose ring remains restricted at the North region of the pseudorotational wheel, which is the one typically observed for nucleosides and nucleotides in the double helix A-RNA, when the phosphate-bearing sugar is anchored to the mineral surface. Results indicate that the severe restrictions imposed by the mineral are responsible of the double helix preservation when DNA and RNA are encapsulated in crystalline nanorods.

The DNA and RNA sugar-phosphate backbone emerges as the key player. An overview of quantum-chemical, structural biology and simulation studies

Knowledge of geometrical and physico-chemical properties of the sugar-phosphate backbone substantially contributes to the comprehension of the structural dynamics, function and evolution of nucleic acids. We provide a side by side overview of structural biology/bioinformatics, quantum chemical and molecular mechanical/simulation studies of the nucleic acids backbone. We highlight main features, advantages and limitations of these techniques, with a special emphasis given to their synergy. The present status of the research is then illustrated by selected examples which include classification of DNA and RNA backbone families, benchmark structure-energy quantum chemical calculations, parameterization of the dihedral space of simulation force fields, incorporation of arsenate into DNA, sugar-phosphate backbone self-cleavage in small RNA enzymes, and intricate geometries of the backbone in recurrent RNA building blocks. Although not apparent from the current literature showing limited overlaps between the QM, simulation and bioinformatics studies of the nucleic acids backbone, there in fact should be a major cooperative interaction between these three approaches in studies of the sugar-phosphate backbone.

Chemical shifts in nucleic acids studied by density functional theory calculations and comparison with experiment

NMR chemical shifts are highly sensitive probes of local molecular conformation and environment and form an important source of structural information. In this study, the relationship between the NMR chemical shifts of nucleic acids and the glycosidic torsion angle, χ, has been investigated for the two commonly occurring sugar conformations. We have calculated by means of DFT the chemical shifts of all atoms in the eight DNA and RNA mono-nucleosides as a function of these two variables. From the DFT calculations, structures and potential energy surfaces were determined by using constrained geometry optimizations at the BP86/TZ2P level of theory. The NMR parameters were subsequently calculated by single-point calculations at the SAOP/TZ2P level of theory. Comparison of the (1)H and (13)C NMR shifts calculated for the mono-nucleosides with the shifts determined by NMR spectroscopy for nucleic acids demonstrates that the theoretical shifts are valuable for the characterization of nucleic acid conformation. For example, a clear distinction can be made between χ angles in the anti and syn domains. Furthermore, a quantitative determination of the χ angle in the syn domain is possible, in particular when (13)C and (1)H chemical shift data are combined. The approximate linear dependence of the C1' shift on the χ angle in the anti domain provides a good estimate of the angle in this region. It is also possible to derive the sugar conformation from the chemical shift information. The DFT calculations reported herein were performed on mono-nucleosides, but examples are also provided to estimate intramolecularly induced shifts as a result of hydrogen bonding, polarization effects, or ring-current effects.