Chemical Shift and Structure Relationship in Nucleic Acids: Correlation of Backbone Torsion Angles γ and α with 13C Chemical Shifts

Structural analysis of ATP bound to the F1-ATPase β-subunit monomer by solid-state NMR- insight into the hydrolysis mechanism in F1

ATP-hydrolysis-associated conformational change of the β-subunit during the rotation of F1-ATPase (F1) has been discussed using cryo-electron microscopy (cryo-EM). Since it is worthwhile to further investigate the conformation of ATP at the catalytic subunit through an alternative approach, the structure of ATP bound to the F1β-subunit monomer (β) was analyzed by solid-state NMR. The adenosine conformation of ATP-β was similar to that of ATP analog in F1 crystal structures. 31P chemical shift analysis showed that the Pα and Pβ conformations of ATP-β are gauche-trans and trans-trans, respectively. The triphosphate chain is more extended in ATP-β than in ATP analog in F1 crystals. This appears to be in the state just before ATP hydrolysis. Furthermore, the ATP-β conformation is known to be more closed than the closed form in F1 crystal structures. In view of the cryo-EM results, ATP-β would be a model of the most closed β-subunit with ATP ready for hydrolysis in the hydrolysis stroke of the F1 rotation.

Why are Hoogsteen base pairs energetically disfavored in A-RNA compared to B-DNA?

A(syn)-U/T and G(syn)-C⁺ Hoogsteen (HG) base pairs (bps) are energetically more disfavored relative to Watson–Crick (WC) bps in A-RNA as compared to B-DNA by >1 kcal/mol for reasons that are not fully understood. Here, we used NMR spectroscopy, optical melting experiments, molecular dynamics simulations and modified nucleotides to identify factors that contribute to this destabilization of HG bps in A-RNA. Removing the 2′-hydroxyl at single purine nucleotides in A-RNA duplexes did not stabilize HG bps relative to WC. In contrast, loosening the A-form geometry using a bulge in A-RNA reduced the energy cost of forming HG bps at the flanking sites to B-DNA levels. A structural and thermodynamic analysis of purine-purine HG mismatches reveals that compared to B-DNA, the A-form geometry disfavors syn purines by 1.5–4 kcal/mol due to sugar-backbone rearrangements needed to sterically accommodate the syn base. Based on MD simulations, an additional penalty of 3–4 kcal/mol applies for purine-pyrimidine HG bps due to the higher energetic cost associated with moving the bases to form hydrogen bonds in A-RNA versus B-DNA. These results provide insights into a fundamental difference between A-RNA and B-DNA duplexes with important implications for how they respond to damage and post-transcriptional modifications.

Atomic structures of excited state A-T Hoogsteen base pairs in duplex DNA by combining NMR relaxation dispersion, mutagenesis, and chemical shift calculations

NMR relaxation dispersion studies indicate that in canonical duplex DNA, Watson-Crick base pairs (bps) exist in dynamic equilibrium with short-lived low abundance excited state Hoogsteen bps. N1-methylated adenine (m¹A) and guanine (m¹G) are naturally occurring forms of damage that stabilize Hoogsteen bps in duplex DNA. NMR dynamic ensembles of DNA duplexes with m¹A-T Hoogsteen bps reveal significant changes in sugar pucker and backbone angles in and around the Hoogsteen bp, as well as kinking of the duplex towards the major groove. Whether these structural changes also occur upon forming excited state Hoogsteen bps in unmodified duplexes remains to be established because prior relaxation dispersion probes provided limited information regarding the sugar-backbone conformation. Here, we demonstrate measurements of C3' and C4' spin relaxation in the rotating frame (R1ρ) in uniformly ¹³C/¹⁵N labeled DNA as sensitive probes of the sugar-backbone conformation in DNA excited states. The chemical shifts, combined with structure-based predictions using an automated fragmentation quantum mechanics/molecular mechanics method, show that the dynamic ensemble of DNA duplexes containing m¹A-T Hoogsteen bps accurately model the excited state Hoogsteen conformation in two different sequence contexts. Formation of excited state A-T Hoogsteen bps is accompanied by changes in sugar-backbone conformation that allow the flipped syn adenine to form hydrogen-bonds with its partner thymine and this in turn results in overall kinking of the DNA toward the major groove. Results support the assignment of Hoogsteen bps as the excited state observed in canonical duplex DNA, provide an atomic view of DNA dynamics linked to formation of Hoogsteen bps, and lay the groundwork for a potentially general strategy for solving structures of nucleic acid excited states.

Solid-state NMR studies of nucleic acid components

Recent applications of solid-state NMR spectroscopy to studies of nucleic acids and their components.

Nucleotide-type chemical shift assignment of the encapsulated 40 kbp dsDNA in intact bacteriophage T7 by MAS solid-state NMR

The icosahedral bacteriophage T7 is a 50 MDa double-stranded DNA (dsDNA) virus that infects Escherichia coli. Although there is substantial information on the physical and morphological properties of T7, structural information, based mostly on Raman spectroscopy and cryo-electron microscopy, is limited. Here, we apply the magic-angle spinning (MAS) solid-state NMR (SSNMR) technique to study a uniformly (13)C and (15)N labeled wild-type T7 phage. We describe the details of the large-scale preparation and purification of an isotopically enriched phage sample under fully hydrated conditions, and show a complete (13)C and a near-complete (15)N nucleotide-type specific assignment of the sugar and base moieties in the 40 kbp dsDNA of T7 using two-dimensional (13)C-(13)C and (15)N-(13)C correlation experiments. The chemical shifts are interpreted as reporters of a B-form conformation of the encapsulated dsDNA. While MAS SSNMR was found to be extremely useful in determining the structures of proteins in native-like environments, its application to nucleic acids has lagged behind, leaving a missing (13)C and (15)N chemical shift database. This work therefore expands the (13)C and (15)N database of real B-form DNA systems, and opens routes to characterize more complex nucleic acid systems by SSNMR.

A geometrical parametrization of C1'-C5' RNA ribose chemical shifts calculated by density functional theory

It has been recently shown that NMR chemical shifts can be used to determine the structures of proteins. In order to begin to extend this type of approach to nucleic acids, we present an equation that relates the structural parameters and the (13)C chemical shifts of the ribose group. The parameters in the equation were determined by maximizing the agreement between the DFT-derived chemical shifts and those predicted through the equation for a database of ribose structures. Our results indicate that this type of approach represents a promising way of establishing quantitative and computationally efficient analytical relationships between chemical shifts and structural parameters in nucleic acids.

Complete chemical shift assignment of the ssDNA in the filamentous bacteriophage fd reports on its conformation and on its interface with the capsid shell

The fd bacteriophage is a filamentous virus consisting of a circular single-stranded DNA (ssDNA) wrapped by thousands of copies of a major coat protein subunit (the capsid). The coat protein subunits are mostly α-helical and curved, and are arranged in the capsid in consecutive pentamers related by a translation along the main viral axis and a rotation of ~36° (C5S2 symmetry). The DNA is right-handed and helical, but information on its structure and on its interface with the capsid is incomplete. We present here an approach for assigning the DNA nucleotides and studying its interactions with the capsid by magic-angle spinning solid-state NMR. Capsid contacts with the ssDNA are obtained using a two-dimensional (13)C-(13)C correlation experiment and a proton-mediated (31)P-(13)C polarization transfer experiment, both acquired on an aromatic-unlabeled phage sample. Our results allow us to map the residues that face the interior of the capsid and to show that the ssDNA-capsid interactions are sustained mainly by electrostatic interactions between the positively charged lysine side chains and the phosphate backbone. The use of natural abundance aromatic amino acids in the growth media facilitated the complete assignment of the four nucleotides and the observation of internucleotide contacts. Using chemical shift analysis, our study shows that structural features of the deoxyribose carbons reporting on the sugar pucker are strikingly similar to those observed recently for the Pf1 phage. However, the ssDNA-protein interface is different, and chemical shift markers of base pairing are different. This experimental approach can be utilized in other filamentous and icosahedral bacteriophages, and also in other biomolecular complexes involving structurally and functionally important DNA-protein interactions.

Chemical shifts for the unusual DNA structure in Pf1 bacteriophage from dynamic-nuclear-polarization-enhanced solid-state NMR spectroscopy

Solid-state NMR spectra, including dynamic nuclear polarization enhanced 400 MHz spectra acquired at 100 K, as well as non-DNP spectra at a variety of field strengths and at temperatures in the range 213-243 K, have allowed the assignment of the (13)C and (15)N resonances of the unusual DNA structure in the Pf1 virion. The (13)C chemical shifts of C3' and C5', considered to be key reporters of deoxyribose conformation, fall near or beyond the edges of their respective ranges in available databases. The (13)C and (15)N chemical shifts of the DNA bases have above-average values for AC4, AC5, CC5, TC2, and TC5, and below average values for AC8, GC8, and GN2, pointing to an absence of Watson-Crick hydrogen bonding, yet the presence of some type of aromatic ring interaction. Crosspeaks between Tyr40 of the coat protein and several DNA atoms suggest that Tyr40 is involved in this ring interaction. In addition, these crosspeak resonances and several deoxyribose resonances are multiply split, presumably through the effects of ordered but differing interactions between capsid protein subunits and each type of nucleotide in each of the two DNA strands. Overall, these observations characterize and support the DNA model proposed by Liu and Day and refined by Tsuboi et al., which calls for the most highly stretched and twisted naturally occurring DNA yet encountered.

DNA structures from phosphate chemical shifts

For B-DNA, the strong linear correlation observed by nuclear magnetic resonance (NMR) between the (31)P chemical shifts (deltaP) and three recurrent internucleotide distances demonstrates the tight coupling between phosphate motions and helicoidal parameters. It allows to translate deltaP into distance restraints directly exploitable in structural refinement. It even provides a new method for refining DNA oligomers with restraints exclusively inferred from deltaP. Combined with molecular dynamics in explicit solvent, these restraints lead to a structural and dynamical view of the DNA as detailed as that obtained with conventional and more extensive restraints. Tests with the Jun-Fos oligomer show that this deltaP-based strategy can provide a simple and straightforward method to capture DNA properties in solution, from routine NMR experiments on unlabeled samples.

Conformational signatures of 13C chemical shifts in RNA ribose

The conformational dependence of (13)C chemical shift values of RNA riboses determined by liquid-state NMR spectroscopy was evaluated using data deposited for RNA structures in the RCSD and BMRB data bases. Results derived support the applicability of the canonical coordinates approach of Rossi and Harbison (J Magn Reson 151:1-8, 2001) in liquid-state NMR to assess the sugar pucker of ribose units in RNA.

Recent advances in RNA structure determination by NMR

Despite recent advances in the solution of NMR structures of RNA and RNA-ligand complexes, the rate limiting step remains the gathering of a large number of NOE and torsion restraints. Additional sources of information for structure determination of larger RNA molecules have recently become available, and it is possible to supplement NOE and J-coupling data with the measurement of dipolar couplings and cross-correlated relaxation rates in high-resolution NMR spectroscopy.

Molecular dynamics simulations of the d(CCAACGTTGG)(2) decamer in crystal environment: comparison of atomic point-charge, extra-point, and polarizable force fields

Molecular dynamics simulations of the DNA duplex d(CCAACGTTGG)(2) were used to study the relationship between DNA sequence and structure in a crystal environment. Three different force fields were used: a traditional description based on atomic point charges, a polarizable force field, and an "extra-point" force field (with additional charges on extranuclear sites). It is found that all the force fields reproduce fairly well the sequence-dependent features of the experimental structure. The polarizable force field, however, provides the most accurate representation of the crystal structure and the sequence-dependent effects observed in the experiment. These results point out to the need of the inclusion of polarization for accurate descriptions of DNA.

Predicting the efficiency of UAG translational stop signal through studies of physicochemical properties of its composite mono- and dinucleotides

In this study, we explored the problem of predicting the UAG stop-codon read-through efficiency. The reported nucleotide sequences were first converted into physicochemical property vectors before being presented to a machine learning algorithm. Two sets of physicochemical properties were applied: one for mononucleosides (in terms of steric bulk, hydrophobicity and electronics) and another for dinucleotides. To the best of our knowledge, this is the first report of how dinucleotides are converted into principle components derived from NMR chemical shift data. A few efficiency prediction models were then derived and a comparison between mononucleoside and dinucleotide-based models was shown. In the derived models, the coefficients of these property based predictors lend themselves to bio-physical interpretations, an advantage which is demonstrated in this study via a prediction model based on the steric bulk factor. Although it is quite simple, the steric bulk factor model explained well the effect of sequence variations surrounding the amber stop codon and the tRNA bearing UCCU anticodon. We further proposed new alternatives at position -1 and +4 of a UAG stop codon sequence to enhance the readthrough efficiency. This research may contribute to a better understanding of the readthrough mechanisms and may also help to study the normal translation termination process.

Changes in drug 13C NMR chemical shifts as a tool for monitoring interactions with DNA

The antibiotic drug, netropsin, was complexed with the DNA oligonucleotide duplex [d(GGTATACC)]2 to monitor drug 13C NMR chemical shifts changes. The binding mode of netropsin to the minor groove of DNA is well-known, and served as a good model for evaluating the relative sensitivity of 13C chemical shifts to hydrogen bonding. Large downfield shifts were observed for four resonances of carbons that neighbor sites which are known to form hydrogen bond interactions with the DNA minor groove. Many of the remaining resonances of netropsin exhibit shielding or relatively smaller deshielding changes. Based on the model system presented here, large deshielding NMR shift changes of a ligand upon macromolecule binding can likely be attributed to hydrogen bond formation at nearby sites.

Random coil carbon chemical shifts of deoxyribonucleic acids

The sequence and temperature effects on random coil DNA carbon chemical shifts have been investigated using sixteen 17-nucleotide sequences. Temperature effect correction parameters have been determined for the aromatic C6/C8 carbons and the deoxyribose C1', C2', and C3' carbons. The carbon chemical shifts of a specific nucleotide in a random coil sequence have been shown to depend mainly on the type of its nearest neighbors. A carbon chemical shift database containing all 64 different types of triplets has been established for predicting random coil DNA carbon chemical shifts. The use of this triplet database for carbon chemical shift predictions shows good accuracy with experimental data, with root-mean-square deviations of 0.09, 0.10, 0.10, and 0.10 ppm and correlation coefficients of 0.999, 0.996, 0.978, and 0.974 for C6/C8, C1', C2', and C3', respectively.

Solid-state NMR determination of sugar ring pucker in (13)C-labeled 2'-deoxynucleosides

The H3'-C3'-C4'-H4' torsional angles of two microcrystalline 2'-deoxynucleosides, thymidine and 2'-deoxycytidine.HCl, doubly (13)C-labeled at the C3' and C4' positions of the sugar ring, have been measured by solid-state magic-angle-spinning nuclear magnetic resonance (NMR). A double-quantum heteronuclear local field experiment with frequency-switched Lee-Goldberg homonuclear decoupling was used. The H3'-C3'-C4'-H4' torsional angles were obtained by comparing the experimental curves with numerical simulations, including the two (13)C nuclei, the directly bonded (1)H nuclei, and five remote protons. The H3'-C3'-C4'-H4' angles were converted into sugar pucker angles and compared with crystallographic data. The delta torsional angles determined by solid-state NMR and x-ray crystallography agree within experimental error. Evidence is also obtained that the proton positions may be unreliable in the x-ray structures. This work confirms that double-quantum solid-state NMR is a feasible tool for studying sugar pucker conformations in macromolecular complexes that are unsuitable for solution NMR or crystallography.

Calculation of 13C chemical shifts in rna nucleosides: structure-13C chemical shift relationships

Isotropic 13C chemical shifts of the ribose sugar in model RNA nucleosides are calculated using SCF and DFT-GIAO ab initio methods for different combinations of ribose sugar pucker, exocyclic torsion angle, and glycosidic torsion angle. Idealized conformations were obtained using structures that were fully optimized by ab initio DFT methods starting with averaged parameters from a collection of crystallographic data. Solid-state coordinates of accurate crystal or neutron diffraction structures were also examined directly without optimization. The resulting 13C chemical shifts for the two sets of calculations are then compared. The GIAO-DFT method overestimates the shifts by an average of 5 ppm while the GIAO-SCF underestimates the shifts by the same amount. However, in the majority of cases the errors appear to be systematic, as the slope of a plot of calculated vs experimental shifts is very close to unity, with minimal scatter. The values of the 13C NMR shifts of the ribose sugar are therefore sufficiently precise to allow for statistical separation of sugar puckering modes and exocyclic torsion angle conformers, based on the canonical equation model formulated in a previous paper.

Changes in 13C NMR chemical shifts of DNA as a tool for monitoring drug interactions

The antibiotic drug, netropsin, was complexed with the DNA oligonucleotide duplex [d(GGTATACC)]2 to explore the effects of ligand binding on the 13C NMR chemical shifts of the DNA base and sugar carbons. The binding mode of netrospin to TA-rich tracts of DNA has been well documented and served as an attractive model system. For the base carbons, four large changes in resonance chemical shifts were observed upon complex formation: -0.64 ppm for carbon 4 of either Ado4 or Ado6, 1.36 ppm for carbon 2 of Thd5, 1.33 ppm for carbon 5 of Thd5 and 0.94 for carbon 6 of Thd5. AdoC4 is covalently bonded to a heteroatom that is hydrogen bonded to netropsin; this relatively large deshielding is consistent with the known hydrogen bond formed at AdoN3. The three large shielding increases are consistent with hydrogen bonds to water in the minor groove being disrupted upon netropsin binding. For the DNA sugar resonances, large changes in chemical shifts were observed upon netropsin complexation. The 2', 3' and 5' 13C resonances of Thd3 and Thd5 were shielded whereas those of Ado4 and Ado6 were deshielded; the 13C resonances of 1' and 4' could not be assigned. These changes are consistent with alteration of the dynamic pseudorotational states occupied by the DNA sugars. A significant alteration in the pseudorotational states of Ado4 or Ado6 must occur as suggested by the large change in chemical shift of -1.65 ppm of the C3' carbon. In conclusion, 13C NMR may serve as a practical tool for analyzing structural changes in DNA-ligand complexes.

Dependence of 13C NMR chemical shifts on conformations of rna nucleosides and nucleotides

Cross-polarization magic-angle spinning solid-state NMR spectroscopy has been used to investigate the dependence of (13)C sugar chemical shifts on specific conformational parameters of a variety of ribonucleotides and ribonucleosides. Solid-state NMR is a valuable tool for nucleoside and nucleotide structural studies since it provides the means to acquire spectra that correspond to single conformations, as opposed to (13)C solution NMR methods. The distinct effects of sugar puckering on the C1', C4', and C5' resonances of C2' endo (S type) and C3' endo (N type) furanoid conformations allow us to separate them into two groups. Further analysis of each group reveals an additional dependence of the C1' and C5' resonances on the glycosidic and C4'-C5' exocyclic torsion angles, respectively. However, it is found that the glycosidic conformation cannot independently be determined from sugar (13)C chemical shift data. The statistical methods of exploratory data analysis and discriminant analysis are used to construct two canonical coordinates-linear combinations of chemical shifts which give the statistically optimal determination of the conformation from the NMR data.

Interpretation of chemical shifts and coupling constants in macromolecules

Recent developments in NMR spectroscopy, along with advances in computational techniques, have produced new approaches to the interpretation of chemical shifts and spin-spin coupling constants in biomolecules. Quantum chemical studies of useful accuracy are now becoming more routine and are increasingly being used in conjunction with experimental studies to map out expected structural patterns for peptides and oligonucleotides. Topics of recent special interest include spin couplings across hydrogen bonds and patterns of chemical shift anisotropies, in both diamagnetic and paramagnetic proteins.

Molecular dynamics simulations of the d(CCAACGTTGG)(2) decamer: influence of the crystal environment

Molecular dynamics (MD) simulations of the DNA duplex d(CCAACGTTGG)(2) were used to study the relationship between DNA sequence and structure. Two crystal simulations were carried out; one consisted of one unit cell containing two duplexes, and the other of two unit cells containing four duplexes. Two solution simulations were also carried out, one starting from canonical B-DNA and the other starting from the crystal structure. For many helicoidal parameters, the results from the crystal and solution simulations were essentially identical. However, for other parameters, in particular, alpha, gamma, delta, (epsilon - zeta), phase, and helical twist, differences between crystal and solution simulations were apparent. Notably, during crystal simulations, values of helical twist remained comparable to those in the crystal structure, to include the sequence-dependent differences among base steps, in which values ranged from 20 degrees to 50 degrees per base step. However, in the solution simulations, not only did the average values of helical twist decrease to approximately 30 degrees per base step, but every base step was approximately 30 degrees, suggesting that the sequence-dependent information may be lost. This study reveals that MD simulations of the crystal environment complement solution simulations in validating the applicability of MD to the analysis of DNA structure.

Nitrogen-15 chemical shifts in AT (adenine-thymine) and CG (cytosine-guanine) nucleic acid base pairs

This paper presents ab initio (DFT) calculations of the 15N chemical shifts in AT (Adenine-Thymine) and CG (Cytosine-Guanine) nucleic acid base pairs. Calculations were done on 14 AT and 18 CG base pairs using experimental (X-ray) geometries obtained from several DNA decamers. The calculated chemical shifts are compared with the experimental values in the pure bases and subjected to statistical analysis to explore their sensitivity to the local geometry and pair helix parameters. The results indicate that the 15N chemical shifts, isotropic and principal components are quite sensitive to small changes in the geometry of the pairs, but they do not correlate well with the helix pair parameters. From the statistical analysis, several linear correlations between structural parameters and chemical shifts emerge. These relationships may serve as a foundation to extract information on molecular structure from 15N chemical shift measurements.

The use of chemical shifts and their anisotropies in biomolecular structure determination

The existence of chemical shift dispersion is crucial for the application of NMR spectroscopy to biomolecules, but the direct interpretation of shift tensors in terms of structure and dynamics is often difficult. Proton shifts reflect environmental influences from nearby aromatic groups, metal sites or hydrogen-bonding partners. These effects can be reasonably modeled with empirical equations, but multiple contributions to shifts can be difficult to disentangle. Shifts for carbon and nitrogen generally reflect local bonding interactions, often in ways that allow the local structure to be inferred. The anisotropy of the shielding tensor is also of interest. It influences the resonance position in partially-ordered samples and has consequences for spin relaxation, even in isotropic systems. There has been recent progress in measuring and interpreting these anisotropies.