NMR methodology for the study of nucleic acids

Analyzing RNA posttranscriptional modifications to decipher the epitranscriptomic code

The discovery of RNA silencing has revealed that non-protein-coding sequences (ncRNAs) can cover essential roles in regulatory networks and their malfunction may result in severe consequences on human health. These findings have prompted a general reassessment of the significance of RNA as a key player in cellular processes. This reassessment, however, will not be complete without a greater understanding of the distribution and function of the over 170 variants of the canonical ribonucleotides, which contribute to the breathtaking structural diversity of natural RNA. This review surveys the analytical approaches employed for the identification, characterization, and detection of RNA posttranscriptional modifications (rPTMs). The merits of analyzing individual units after exhaustive hydrolysis of the initial biopolymer are outlined together with those of identifying their position in the sequence of parent strands. Approaches based on next generation sequencing and mass spectrometry technologies are covered in depth to provide a comprehensive view of their respective merits. Deciphering the epitranscriptomic code will require not only mapping the location of rPTMs in the various classes of RNAs, but also assessing the variations of expression levels under different experimental conditions. The fact that no individual platform is currently capable of meeting all such demands implies that it will be essential to capitalize on complementary approaches to obtain the desired information. For this reason, the review strived to cover the broadest possible range of techniques to provide readers with the fundamental elements necessary to make informed choices and design the most effective possible strategy to accomplish the task at hand.

A strong preference for the TA/TA dinucleotide step discovered for an acridine-based, potent antitumor dsDNA intercalator, C-1305: NMR-driven structural and sequence-specificity studies

Triazoloacridinone C-1305, a potent antitumor agent recommended for Phase I clinical trials, exhibits high activity towards a wide range of experimental colon carcinomas, in many cases associated with complete tumor regression. C-1305 is a well-established dsDNA intercalator, yet no information on its mode of binding into DNA is available to date. Herein, we present the NMR-driven and MD-refined reconstruction of the 3D structures of the d(CGATATCG)₂:C-1305 and d(CCCTAGGG)₂:C-1305 non-covalent adducts. In both cases, the ligand intercalates at the TA/TA site, forming well-defined dsDNA:drug 1:1 mol/mol complexes. Orientation of the ligand within the binding site was unambiguously established by the DNA/ligand proton-proton NOE contacts. A subsequent, NMR-driven study of the sequence-specificity of C-1305 using a series of DNA duplexes, allowed us to confirm a strong preference towards TA/TA dinucleotide steps, followed by the TG/CA steps. Interestingly, no interaction at all was observed with duplexes containing exclusively the AT/AT, GG/CC and GA/TC steps.

Structural Modification of DNA Studied by Scanning Tunneling Microscopy

In addition to the well-known double-helix structure of deoxyribonucleic acid (DNA), theoretical simulations have predicted structural variants, such as cross and T-shaped formations. The direct visualization of individual DNA molecules by scanning probe microscopy should be able to identify these structures. In this study, we examine various structures of DNA deposited on Au(111) by using an electrospray method and analyzing with a low-temperature scanning tunneling microscope. We could identify in detail several interesting structures. In particular, we found that double-stranded DNAs are partially unraveled into single-stranded DNAs and we confirmed the presence of T-shaped DNA structures.

Applications of high dimensionality experiments to biomolecular NMR

High dimensionality NMR experiments facilitate resonance assignment and precise determination of spectral parameters such as coupling constants. Sparse non-uniform sampling enables acquisition of experiments of high dimensionality with high resolution in acceptable time. In this review we present and compare some significant applications of NMR experiments of dimensionality higher than three in the field of biomolecular studies in solution.

N(6)-Methyladenosine Modification in a Long Noncoding RNA Hairpin Predisposes Its Conformation to Protein Binding

N(6)-Methyladenosine (m(6)A) is a reversible and abundant internal modification of messenger RNA (mRNA) and long noncoding RNA (lncRNA) with roles in RNA processing, transport, and stability. Although m(6)A does not preclude Watson-Crick base pairing, the N(6)-methyl group alters the stability of RNA secondary structure. Since changes in RNA structure can affect diverse cellular processes, the influence of m(6)A on mRNA and lncRNA structure has the potential to be an important mechanism for m(6)A function in the cell. Indeed, an m(6)A site in the lncRNA metastasis associated lung adenocarcinoma transcript 1 (MALAT1) was recently shown to induce a local change in structure that increases the accessibility of a U5-tract for recognition and binding by heterogeneous nuclear ribonucleoprotein C (HNRNPC). This m(6)A-dependent regulation of protein binding through a change in RNA structure, termed "m(6)A-switch", affects transcriptome-wide mRNA abundance and alternative splicing. To further characterize this first example of an m(6)A-switch in a cellular RNA, we used nuclear magnetic resonance and Förster resonance energy transfer to demonstrate the effect of m(6)A on a 32-nucleotide RNA hairpin derived from the m(6)A-switch in MALAT1. The observed imino proton nuclear magnetic resonance resonances and Förster resonance energy transfer efficiencies suggest that m(6)A selectively destabilizes the portion of the hairpin stem where the U5-tract is located, increasing the solvent accessibility of the neighboring bases while maintaining the overall hairpin structure. The m(6)A-modified hairpin has a predisposed conformation that resembles the hairpin conformation in the RNA-HNRNPC complex more closely than the unmodified hairpin. The m(6)A-induced structural changes in the MALAT1 hairpin can serve as a model for a large family of m(6)A-switches that mediate the influence of m(6)A on cellular processes.

TEMPO-derived spin labels linked to the nucleobases adenine and cytosine for probing local structural perturbations in DNA by EPR spectroscopy

Three 2´-deoxynucleosides containing semi-flexible spin labels, namely (T)A, (U)A and (U)C, were prepared and incorporated into deoxyoligonucleotides using the phosphoramidite method. All three nucleosides contain 2,2,6,6-tetramethylpiperidine-1-oxyl (TEMPO) connected to the exocyclic amino group; (T)A directly and (U)A as well as (U)C through a urea linkage. (T)A and (U)C showed a minor destabilization of a DNA duplex, as registered by a small decrease in the melting temperature, while (U)A destabilized the duplex by more than 10 °C. Circular dichroism (CD) measurements indicated that all three labels were accommodated in B-DNA duplex. The mobility of the spin label (T)A varied with different base-pairing partners in duplex DNA, with the (T)A•T pair being the least mobile. Furthermore, (T)A showed decreased mobility under acidic conditions for the sequences (T)A•C and (T)A•G, to the extent that the EPR spectrum of the latter became nearly superimposable to that of (T)A•T. The reduced mobility of the (T)A•C and (T)A•G mismatches at pH 5 is consistent with the formation of (T)AH(+)•C and (T)AH(+)•G, in which protonation of N1 of A allows the formation of an additional hydrogen bond to N3 of C and N7 of G, respectively, with G in a syn-conformation. The urea-based spin labels (U)A and (U)C were more mobile than (T)A, but still showed a minor variation in their EPR spectra when paired with A, G, C or T in a DNA duplex. (U)A and (U)C had similar mobility order for the different base pairs, with the lowest mobility when paired with C and the highest when paired with T.

G-triplex structure and formation propensity

The occurrence of a G-triplex folding intermediate of thrombin binding aptamer (TBA) has been recently predicted by metadynamics calculations, and experimentally supported by Nuclear Magnetic Resonance (NMR), Circular Dichroism (CD) and Differential Scanning Calorimetry (DSC) data collected on a 3' end TBA-truncated 11-mer oligonucleotide (11-mer-3'-t-TBA). Here we present the solution structure of 11-mer-3'-t-TBA in the presence of potassium ions. This structure is the first experimental example of a G-triplex folding, where a network of Hoogsteen-like hydrogen bonds stabilizes six guanines to form two G:G:G triad planes. The G-triplex folding of 11-mer-3'-t-TBA is stabilized by the potassium ion and destabilized by increasing the temperature. The superimposition of the experimental structure with that predicted by metadynamics shows a great similarity, with only significant differences involving two loops. These new structural data show that 11-mer-3'-t-TBA assumes a G-triplex DNA conformation as its stable form, reinforcing the idea that G-triplex folding intermediates may occur in vivo in human guanine-rich sequences. NMR and CD screening of eight different constructs obtained by removing from one to four bases at either the 3' and the 5' ends show that only the 11-mer-3'-t-TBA yields a relatively stable G-triplex.

Nucleic acid ion structures in the gas phase

Nucleic acids are diverse polymeric macromolecules that are essential for all life forms. These biomolecules possess a functional three-dimensional structure under aqueous physiological conditions. Mass spectrometry-based approaches have on the other hand opened the possibility to gain structural information on nucleic acids from gas-phase measurements. To correlate gas-phase structural probing results with solution structures, it is therefore important to grasp the extent to which nucleic acid structures are preserved, or altered, when transferred from the solution to a fully anhydrous environment. We will review here experimental and theoretical approaches available to characterize the structure of nucleic acids in the gas phase (with a focus on oligonucleotides and higher-order structures), and will summarize the structural features of nucleic acids that can be preserved in the gas phase on the experiment time scale.

Probing DNA mismatched and bulged structures by using 19F NMR spectroscopy and oligodeoxynucleotides with an 19F-labeled nucleobase

In this study, DNA local structures with bulged bases and mismatched base pairs as well as ordinary full-matched base pairs by using (19)F NMR spectroscopy with (19)F-labeled oligodeoxynucleotides (ODNs) were monitored. The chemical shift change in the (19) F NMR spectra allowed discrimination of the DNA structures. Two types of ODNs possessing the bis(trifluoromethyl)benzene unit (F-unit) at specified uridines were prepared and hybridized with their complementary or noncomplementary strands to form matched, mismatched, or bulged duplexes. By using ODN F1, in which an F-unit was connected directly to a propargyl amine-substituted uridine, three local structures, that is, full-matched, G-U mismatch, and A-bulge could be analyzed, whereas other structures could not be discriminated. A molecular modeling study revealed that the F-unit in ODN F1 interacted little with the nucleobases and sugar backbone of the opposite strand because the linker length between the F-unit and the uridine base was too short. Therefore, the capacity of ODN F1 to discriminate the DNA local structures was limited. Thus, ODN F2 was designed to improve this system; aminobenzoic acid was inserted between the F-unit and uridine base so the F-unit could interact more closely with the opposite strand. Eventually, the G-bulge and T-U mismatch and the three aforementioned local structures could be discriminated by using ODN F2. In addition, the dissociation processes of these duplexes could be monitored concurrently by (19)F NMR spectroscopy.

Automatic workflow for the classification of local DNA conformations

Background

A growing number of crystal and NMR structures reveals a considerable structural polymorphism of DNA architecture going well beyond the usual image of a double helical molecule. DNA is highly variable with dinucleotide steps exhibiting a substantial flexibility in a sequence-dependent manner. An analysis of the conformational space of the DNA backbone and the enhancement of our understanding of the conformational dependencies in DNA are therefore important for full comprehension of DNA structural polymorphism.

Results

A detailed classification of local DNA conformations based on the technique of Fourier averaging was published in our previous work. However, this procedure requires a considerable amount of manual work. To overcome this limitation we developed an automatic classification method consisting of the combination of supervised and unsupervised approaches. A proposed workflow is composed of k-NN method followed by a non-hierarchical single-pass clustering algorithm. We applied this workflow to analyze 816 X-ray and 664 NMR DNA structures released till February 2013. We identified and annotated six new conformers, and we assigned four of these conformers to two structurally important DNA families: guanine quadruplexes and Holliday (four-way) junctions. We also compared populations of the assigned conformers in the dataset of X-ray and NMR structures.

Conclusions

In the present work we developed a machine learning workflow for the automatic classification of dinucleotide conformations. Dinucleotides with unassigned conformations can be either classified into one of already known 24 classes or they can be flagged as unclassifiable. The proposed machine learning workflow permits identification of new classes among so far unclassifiable data, and we identified and annotated six new conformations in the X-ray structures released since our previous analysis. The results illustrate the utility of machine learning approaches in the classification of local DNA conformations.

Determination of (3)J((1)H3'-(31)P) couplings in a DNA oligomer with enhanced sensitivity employing a constant-time TOCSY difference experiment

A constant-time TOCSY difference experiment for the determination of (3)J((1)H3'-(31)P) coupling constants in non-isotope-labelled DNA oligonucleotides is presented. The method is tested on a DNA octamer and compared with the established constant-time NOESY difference method. Each (3)J((1)H3'-(31)P) coupling constant is determined from amplitude changes caused by phosphorous decoupling, which are observable on multiple cross-peaks, thus leading to a high accuracy of the value of the (3)J((1)H3'-(31)P) coupling constant. The new experiment delivers up to three times the sensitivity compared with previously reported methods.

Development of CHARMM polarizable force field for nucleic acid bases based on the classical Drude oscillator model

A polarizable force field for nucleic acid bases based on the classical Drude oscillator model is presented. Parameter optimization was performed to reproduce crystallographic geometries, crystal unit cell parameters, heats of sublimation, vibrational frequencies and assignments, dipole moments, molecular polarizabilities and quantum mechanical base-base and base-water interaction energies. The training and validation data included crystals of unsubstituted and alkyl-substituted adenine, guanine, cytosine, uracil, and thymine bases, hydrated crystals, and hydrogen bonded base pairs. Across all compounds, the RMSD in the calculated heats of sublimation is 4.1%. This equates to an improvement of more than 2.5 kcal/mol in accuracy compared to the nonpolarizable CHARMM27 force field. However, the level of agreement with experimental molecular volume decreased from 1.7% to 2.1% upon moving from the nonpolarizable to the polarizable model. The representation of dipole moments is significantly improved with the Drude polarizable force field. Unlike in additive force fields, there is no requirement for the gas-phase dipole moments to be overestimated, illustrating the ability of the Drude polarizable force field to treat accurately differently dielectric environments and indicating the improvements in the electrostatic model. Validation of the model was performed on the basis of the calculation of the gas-phase binding enthalpies of base pairs obtained via potential of mean force calculations; the additive and polarizable models both performed satisfactorily with average differences of 0.2 and 0.9 kcal/mol, respectively, and rms differences of 1.3 and 1.7 kcal/mol, respectively. Overall, considering the number of significant improvements versus the additive CHARMM force field, the incorporation of explicit polarizability into the force field for nucleic acid bases represents an additional step toward accurate computational modeling of biological systems.

Characterization of RNA invasion by (19)F NMR spectroscopy

(19)F NMR spectroscopy offers an efficient tool for monitoring RNA invasion. The invasion of 2'-O-methyl oligoribonucleotides into a (19)F-labeled HIV-1 TAR RNA model and the temperature-dependent behavior of the complexes obtained have been examined.

Monitoring of duplex and triplex formation by 19F NMR using oligodeoxynucleotides possessing 5-fluorodeoxyuridine unit as 19F signal transmitter

We prepared oligodeoxynucleotides (ODNs) possessing a 5-fluorodeoxyuridine (5-FU) unit as a 19F-signal transmitter, and characterized their structures including single strand, duplex, and triplex using 19F NMR. The change in chemical shift induced by incorporation of 5-FU into the ODNs and the formation of higher order structures allowed monitoring of structural changes. Data from UV melting experiments and CD spectra were consistent with the spectral changes in the NMR studies. These 19F-labeled ODNs may be promising molecular probes for the identification of DNA structures in complicated biological conditions.

Computational approaches to 3D modeling of RNA

Many exciting discoveries have recently revealed the versatility of RNA and its importance in a variety of functions within the cell. Since the structural features of RNA are of major importance to their biological function, there is much interest in predicting RNA structure, either in free form or in interaction with various ligands, including proteins, metabolites and other molecules. In recent years, an increasing number of researchers have developed novel RNA algorithms for predicting RNA secondary and tertiary structures. In this review, we describe current experimental and computational advances and discuss recent ideas that are transforming the traditional view of RNA folding. To evaluate the performance of the most recent RNA 3D folding algorithms, we provide a comparative study in order to test the performance of available 3D structure prediction algorithms for an RNA data set of 43 structures of various lengths and motifs. We find that the algorithms vary widely in terms of prediction quality across different RNA lengths and topologies; most predictions have very large root mean square deviations from the experimental structure. We conclude by outlining some suggestions for future RNA folding research.

On the parameterization of rigid base and basepair models of DNA from molecular dynamics simulations

A method is described to extract a complete set of sequence-dependent material parameters for rigid base and basepair models of DNA in solution from atomistic molecular dynamics simulations. The method is properly consistent with equilibrium statistical mechanics, leads to effective shape, stiffness and mass parameters, and employs special procedures for treating spontaneous torsion angle flips and H-bond breaks, both of which can have a significant effect on the results. The method is accompanied by various analytical consistency checks that can be used to assess the equilibration of statistical averages, and different modeling assumptions pertaining to the rigidity of the bases and basepairs and the locality of the quadratic internal energy. The practicability of the approach is verified by estimating complete parameter sets for the 16-basepair palindromic oligomer G(TA)(7)C simulated in explicit water and counterions. Our results indicate that the method is capable of resolving sequence-dependent variations in each of the material parameters. Moreover, they show that the assumptions of rigidity and locality hold rather well for the base model, but not for the basepair model. For the latter, it is shown that the non-local nature of the internal energy can be understood in terms of a certain compatibility relation involving Schur complements.

Longitudinal-relaxation-enhanced NMR experiments for the study of nucleic acids in solution

Atomic-resolution information on the structure and dynamics of nucleic acids is essential for a better understanding of the mechanistic basis of many cellular processes. NMR spectroscopy is a powerful method for studying the structure and dynamics of nucleic acids; however, solution NMR studies are currently limited to relatively small nucleic acids at high concentrations. Thus, technological and methodological improvements that increase the experimental sensitivity and spectral resolution of NMR spectroscopy are required for studies of larger nucleic acids or protein-nucleic acid complexes. Here we introduce a series of imino-proton-detected NMR experiments that yield an over 2-fold increase in sensitivity compared to conventional pulse schemes. These methods can be applied to the detection of base pair interactions, RNA-ligand titration experiments, measurement of residual dipolar (15)N-(1)H couplings, and direct measurements of conformational transitions. These NMR experiments employ longitudinal spin relaxation enhancement techniques that have proven useful in protein NMR spectroscopy. The performance of these new experiments is demonstrated for a 10 kDa TAR-TAR*(GA) RNA kissing complex and a 26 kDa tRNA.

Remarkable metal counterion effect on the internucleotide J-couplings and chemical shifts of the N-H...N hydrogen bonds in the W-C base pairs

The effects of metal ion binding on the (2h) J(NN)-coupling and delta( (1)H)/Deltadelta( (15)N) chemical shifts of N-H...N H-bond units in internucleotide base pairs were explored by a combination of density functional theory calculations and molecular dynamics (MD) simulations. Results indicate that the NMR parameters vary considerably upon cation binding to the natural GC or AT base pairs, and thus can be used to identify the status of the base pairs, if cation-perturbed. The basic trend is that cation perturbation causes (2h) J(NN) to increase, Deltadelta( (15)N) to decrease, and delta( (1)H) to shift upfield for GC, and in the opposite directions for AT. The magnitudes of variation are closely related to the Lewis acidity of the metal ions. For both base pair series (M(z+)GC and M(z+)AT), these NMR parameters are linearly correlated among themselves. Their values depend strongly on the energy gaps (Delta(ELP-->sigma*)) and the second-order interaction energies ( E(2)) between the donor N lone pair (LP(N)) and the acceptor sigma* N-H localized NBO orbitals. In addition, the (2h) J NN changes are also sensitive to the amount of sigma charge transfer from LP(N) to sigma*(N-H) NBOs or from the purine to the pyrimidine moieties. The different trends are a consequence of the different H-bond patterns combined with the polarization effect of the metal ions in the cationized M(z+)AT series, M(z+) <-- A --> T, and the cationized GC series, M(z+) <-- G <-- C. The predicted cation-induced systematic trends of (2h) J(NN) and delta( (15)N, (1)H) in N-H...N H-bond units may provide a new approach to the determination of H-bond structure and strength in Watson-Crick base pairs, and provide an alternative probe of the heterogeneity of DNA sequences.

Enhanced NMR signal detection of imino protons in RNA molecules containing 3' dangling nucleotides

We present a method for improving the quality of nuclear magnetic resonance (NMR) spectra involving exchangeable protons near the base of the stem of RNA hairpin molecules. NMR spectra of five different RNA hairpins were compared. These hairpins consisted of a native RNA structure and four molecules each having different unpaired, or dangling, nucleotides at the 3' end. NMR experiments were acquired in water for each construct and the quality of the imino proton spectral regions were examined. The imino resonances near the base of the stem of the wild type RNA structure were not observed due to breathing motions. However, a significant increase in spectral quality for molecules with dangling 3' adenosine or guanosine nucleotides was observed, with imino protons detected in these constructs that were not observed in the wild type construct. A modest improvement in spectral quality was seen for the construct with a 3' unpaired uridine, whereas no significant improvement was observed for a 3' unpaired cytidine. This improvement in NMR spectral quality mirrors the increased thermodynamic stability observed for 3' unpaired nucleotides which is dependant on the stacking interactions of these nucleotides against the base of the stem. The use of a dangling 3' adenosine nucleotide represents an easy method to significantly improve the quality of NMR spectra of RNA molecules.

RNA structure determination by NMR

This chapter reviews the methodologies for RNA structure determination by liquid-state nuclear magnetic resonance (NMR). The routine production of milligram quantities of isotopically labeled RNA remains critical to the success of NMR-based structure studies. The standard method for the preparation of isotopically labeled RNA for structural studies in solution is in vitro transcription from DNA oligonucleotide templates using T7 RNA polymerase and unlabeled or isotopically labeled nucleotide triphosphates (NTPs). The purification of the desired RNA can be performed by either denaturing polyacrylamide gel electrophoresis (PAGE) or anion-exchange chromatography. Our basic strategy for studying RNA in solution by NMR is outlined. The topics covered include RNA resonance assignment, restraint collection, and the structure calculation process. Selected examples of NMR spectra are given for a correctly folded 30 nucleotide-containing RNA.

Theoretical study of the scalar coupling constants across the noncovalent contacts in RNA base pairs: the cis- and trans-watson-crick/sugar edge base pair family

The structure and function of RNA molecules are substantially affected by non-Watson-Crick base pairs actively utilizing the 2'-hydroxyl group of ribose. Here we correlate scalar coupling constants across the noncovalent contacts calculated for the cis- and trans-WC/SE (Watson-Crick/sugar edge) RNA base pairs with the geometry of base to base and sugar to base hydrogen bond(s). 23 RNA base pairs from the 32 investigated were found in RNA crystal structures, and the calculated scalar couplings are therefore experimentally relevant with regard to the binding patterns occurring in this class of RNA base pairs. The intermolecular scalar couplings 1hJ(N,H), 2hJ(N,N), 2hJ(C,H), and 3hJ(C,N) were calculated for the N-H...N and N-H...O=C base to base contacts and various noncovalent links between the sugar hydroxyl and RNA base. Also, the intramolecular 1J(N,H) and 2J(C,H) couplings were calculated for the amino or imino group of RNA base and the ribose 2'-hydroxyl group involved in the noncovalent interactions. The calculated scalar couplings have implications for validation of local geometry, show specificity for the amino and imino groups of RNA base involved in the linkage, and can be used for discrimination between the cis- and trans-WC/SE base pairs. The RNA base pairs within an isosteric subclass of the WC/SE binding patterns can be further sorted according to the scalar couplings calculated across different local noncovalent contacts. The effect of explicit water inserted in the RNA base pairs on the magnitude of the scalar couplings was calculated, and the data for discrimination between the water-inserted and direct RNA base pairs are presented. The calculated NMR data are significant for structural interpretation of the scalar couplings in the noncanonical RNA base pairs.

Chirality as a tool in nucleic acid recognition: principles and relevance in biotechnology and in medicinal chemistry

The understanding of the interaction of chiral species with DNA or RNA is very important for the development of new tools in biology and of new drugs. Several cases in which chirality is a crucial point in determining the DNA binding mode are reviewed and discussed, with the aim of illustrating how chirality can be considered as a tool for improving the understanding of mechanisms and the effectiveness of nucleic acid recognition. The review is divided into two parts: the former describes examples of chiral species interacting with DNA: intercalators, metal complexes, and groove binders; the latter part is dedicated to chirality in DNA analogs, with discussion of phosphate stereochemistry and chirality of ribose substitutes, in particular of peptide nucleic acids (PNAs) for which a number of works have been published recently dealing with the effect of chirality in DNA recognition. The discussion is intended to show how enantiomeric recognition originates at the molecular level, by exploiting the enormous progresses recently achieved in the field of structural characterization of complexes formed by nucleic acid with their ligands by crystallographic and spectroscopic methods. Examples of application of the DNA binding molecules described and the role of chirality in DNA recognition relevant for biotechnology or medicinal chemistry are reported.

Solution structure of an rRNA substrate bound to the pseudouridylation pocket of a box H/ACA snoRNA

Base pairing between the RNA components of box H/ACA small nucleolar ribonucleoproteins (snoRNPs) and sequences in other eukaryotic RNAs target specific uridines for pseudouridylation. An RNA called HJ1 has been developed that interacts with the rRNA sequence targeted by the 5' pseudouridylation pocket of human U65 snoRNA the same way as intact U65 snoRNA. Sequences on both strands of the analog of the U65 snoRNP pseudouridylation pocket in HJ1 pair with its substrate sequence, and the resulting complex, called HJ3, is strongly stabilized by Mg(2+). The solution structure of HJ3 reveals an Omega-shaped RNA interaction motif that has not previously been described, which is likely to be common to all box H/ACA snoRNP-substrate complexes. The topology of the complex explains why the access of substrate sequences to snoRNPs is facile and how uridine selection may occur when these complexes form.

Calculation of structural behavior of indirect NMR spin-spin couplings in the backbone of nucleic acids

Calculated indirect NMR spin-spin coupling constants (J-couplings) between (31)P, (13)C, and (1)H nuclei were related to the backbone torsion angles of nucleic acids (NAs), and it was shown that J-couplings can facilitate accurate and reliable structural interpretation of NMR measurements and help to discriminate between their distinct conformational classes. A proposed stepwise procedure suggests assignment of the J-couplings to torsion angles from the sugar part to the phosphodiester link. Some J-couplings show multidimensional dependence on torsion angles, the most prominent of which is the effect of the sugar pucker. J-couplings were calculated in 16 distinct nucleic acid conformations, two principal double-helical DNAs, B- and A-, the main RNA form, A-RNA, as well as in 13 other RNA conformations. High-level quantum mechanics calculations used a baseless dinucleoside phosphate as a molecular model, and the effect of solvent was included. The predicted J-couplings correlate reliably with available experimental data from the literature.

Electron bridging dihydrogen bond in the imidazole-contained anion derivatives

The large contact distance of electron bridging dihydrogen bond (EBDB), which is over 2.4 A, is the most prominent characteristic for the imidazole-contained anion derivatives. The elongation of N-H bond and the shortening of H...H distance can be observed upon hydration and hydrogenation. Transformation from EBDB to dissociative H2 is convenient upon sequential hydrogenation. The H...H distance decreases with the enhancement of the electronegativity of the heavy atom which contacts directly with one of these two hydrogen atoms. NMR shielding of the bonding N varies significantly upon hydration and hydrogenation. The spin-spin coupling constants, 1J(H-H), is dominated predominantly by the paramagnetic spin-orbit and diamagnetic spin-orbit contributions instead of the Fermi-contact term. Enhancement of electronegativity of the heavy atom leads to the increase of 1J(H-H) coupling constants. The stabilization is enhanced upon hydration predominantly for the formation of O-H...N H bond, while it is reversed upon hydrogenation for the cleavage of big pi bond, Pi5(6). Enhancement of the stability is demonstrated by the increase of stabilization energy and vertical electron detachment energy with the electronegativity of the heavy atom. The dominant contributions for the formation of such electron bridging dihydrogen bond are the high polarity of each fragment, large electron density between two fragments, and strong bonding interaction of the bridging electron with H(N) atoms. The H...H interaction can be formed by X-Hdelta+ and Hdelta- -Y polar molecules in Hdelta+...Hdelta- and Hdelta+...e...Hdelta+ of two forms.

Dynamics-based amplification of RNA function and its characterization by using NMR spectroscopy

The ever-increasing cellular roles ascribed to RNA raise fundamental questions regarding how a biopolymer composed of only four chemically similar building-block nucleotides achieves such functional diversity. Here, I discuss how RNA achieves added mechanistic and chemical complexity by undergoing highly controlled conformational changes in response to a variety of cellular signals. I examine pathways for achieving selectivity in these conformational changes that rely to different extents on the structure and dynamics of RNA. Finally, I review solution-state NMR techniques that can be used to characterize RNA structural dynamics and its relationship to function.

NMR methods for studying the structure and dynamics of RNA

Proper functioning of RNAs requires the formation of complex three-dimensional structures combined with the ability to rapidly interconvert between multiple functional states. This review covers recent advances in isotope-labeling strategies and NMR experimental approaches that have promise for facilitating solution structure determinations and dynamics studies of biologically active RNAs. Improved methods for the production of isotopically labeled RNAs combined with new multidimensional heteronuclear NMR experiments make it possible to dramatically reduce spectral crowding and simplify resonance assignments for RNAs. Several novel applications of experiments that directly detect hydrogen-bonding interactions are discussed. These studies demonstrate how NMR spectroscopy can be used to distinguish between possible secondary structures and identify mechanisms of ligand binding in RNAs. A variety of recently developed methods for measuring base and sugar residual dipolar couplings are described. NMR residual dipolar coupling techniques provide valuable data for determining the long-range structure and orientation of helical regions in RNAs. A number of studies are also presented where residual dipolar coupling constraints are used to determine the global structure and dynamics of RNAs. NMR relaxation data can be used to probe the dynamics of macromolecules in solution. The power dependence of transverse rotating-frame relaxation rates was used here to study dynamics in the minimal hammerhead ribozyme. Improved methods for isotopically labeling RNAs combined with new types of structural data obtained from a growing repertoire of NMR experiments are facilitating structural and dynamic studies of larger RNAs.

Theoretical Calculation of the NMR Spin−Spin Coupling Constants and the NMR Shifts Allow Distinguishability between the Specific Direct and the Water-Mediated Binding of a Divalent Metal Cation to Guanine

The calculated intermolecular and intramolecular indirect NMR spin-spin coupling constants and NMR shifts were used for the discrimination between the inner-shell and the outer-shell binding motif of hydrated divalent cations Mg(2+) or Zn(2+) with a guanine base. The intermolecular coupling constants (1)J(X,O6) and (1)J(X,N7) (X = Mg(2+), Zn(2+)) can be unambiguously assigned to the specific inner-shell binding motif of the hydrated cation either with oxygen O6 or with nitrogen N7 of guanine. The calculated coupling constants (1)J(Mg,O6) and (1)J(Zn,O6) were 6.2 and -17.5 Hz, respectively, for the inner-shell complex of cation directly interacting with oxygen O6 of guanine. For the inner-shell coordination of the cation at nitrogen N7, the calculated coupling constants (1)J(Mg,N7) and (1)J(Zn,N7) were 5.6 and -36.5 Hz, respectively. When the binding of the cation is water-mediated, the coupling constant is zero. To obtain reliable shifts in NMR parameters, hydrated guanine was utilized as the reference state. The calculated change of NMR spin-spin coupling constants due to the hydration and coordination of the cation with guanine is caused mainly by the variation of Fermi-contact coupling contribution while the variation of diamagnetic spin-orbit, paramagnetic spin-orbit, and spin-dipolar coupling contributions is small. The change of s-character of guanine sigma bonding, sigma antibonding, and lone pair orbitals upon the hydration and cation coordination (calculated using the Natural Bond Orbital analysis) correlates with the variation of the Fermi-contact term. The calculated NMR shifts delta(N7) of -15.3 and -12.2 ppm upon the coordination of Mg(2+) and Zn(2+) ion are similar to the NMR shift of 19.6 ppm toward the high field measured by Tanaka for N7 of guanine upon the coordination of the Cd(2+) cation (Tanaka, Y.; Kojima, C.; Morita, E. H.; Kasai. Y.; Yamasaki, K.; Ono, A.; Kainosho, M.; Taira, K. J. Am. Chem. Soc. 2002, 124, 4595-4601). The present data indicate that measurements of NMR intermolecular coupling constants may be used to discriminate between the specific inner- and outer-shell binding of divalent cations to nucleobases in DNA and RNA.

Internal consistency of NMR data obtained in partially aligned biomacromolecules

A method of testing structure-related NMR data prior to structure calculations is presented. The test is applicable to second rank tensor interactions (dipolar coupling, anisotropic chemical shielding, and quadrupolar interaction) observed in partially aligned samples of biomacromolecules. The method utilizes the fact that only limited number of frequencies corresponding to the mentioned interactions can be measured independently in a rigid fragment of the macromolecule. Additional values can be predicted as linear combinations of the set of independent frequencies. Internal consistency of sufficiently large sets of frequencies measured in individual molecular fragments is tested by comparing the experimental data with their predicted values. The method requires only knowledge of local geometry (i.e., definition of the interaction tensors in the local coordinate frames of the fragments). No information about the alignment or shape of the molecule is required. The test is best suited for planar fragments. Application to peptide bonds and nucleic acid bases is demonstrated.

Towards structural genomics of RNA: rapid NMR resonance assignment and simultaneous RNA tertiary structure determination using residual dipolar couplings

We report a new residual dipolar couplings (RDCs) based NMR procedure for rapidly determining RNA tertiary structure demonstrated on a uniformly (15)N/(13)C-labeled 27 nt variant of the trans-activation response element (TAR) RNA from HIV-I. In this procedure, the time-consuming nuclear Overhauser enhancement (NOE)-based sequential assignment step is replaced by a fully automated RDC-based assignment strategy. This approach involves examination of all allowed sequence-specific resonance assignment permutations for best-fit agreement between measured RDCs and coordinates for sub-structures in a target RNA. Using idealized A-form geometries to model Watson-Crick helices and coordinates from a previous X-ray structure to model a hairpin loop in TAR, the best-fit RDC assignment solutions are determined very rapidly (<five minutes of computational time) and are in complete agreement with corresponding NOE-based assignments. Orientational constraints derived from RDCs are used simultaneously to assemble sub-structures into an RNA tertiary conformation. Through enhanced speeds of application and reduced reliance on chemical shift dispersion, this RDC-based approach lays the foundation for rapidly determining RNA conformations in a structural genomics context, and may increase the size limit of RNAs that can be examined by NMR.

A novel approach for uniform (13)C and (15)N labeling of DNA for NMR studies

A novel method is proposed for large-scale synthesis of (13)C- and (15)N-labeled DNA for NMR studies. In this methodology, endonuclease-sensitive repeat amplification (ESRA), a modified PCR strategy, has been used to amplify tandem repeats of the target DNA sequence. The design of the template is such that restriction enzyme (RE) sites separate repeats of the target sequence. The ESRA product is then cloned into a suitable vector. The Escherichia coli cells harboring the plasmid are grown in minimal medium containing [(13)C]glucose and (15)NH(4)Cl as the sole source of carbon and nitrogen, respectively. The target sequence is released by RE digestion of the plasmid, followed by purification using PAGE. Under optimized conditions, the yield ( approximately 5 mg/liter of culture) of (13)C/(15)N-labeled DNA prepared using this approach is found to be several times higher compared to other known enzymatic methods. Successful incorporation of the isotopes has been confirmed using 2D NMR techniques.

Concerted motions in HIV-1 TAR RNA may allow access to bound state conformations: RNA dynamics from NMR residual dipolar couplings

Ground-state dynamics in RNA is a critical precursor for structural adaptation observed ubiquitously in protein-RNA recognition. A tertiary conformational analysis of the stem-loop structural element in the transactivation response element (TAR) from human immunodeficiency virus type 1 (HIV-I) RNA is presented using recently introduced NMR methods that rely on the measurement of residual dipolar couplings (RDC) in partially oriented systems. Order matrix analysis of RDC data provides evidence for inter-helical motions that are of amplitude 46(+/-4) degrees, of random directional character, and that are executed about an average conformation with an inter-helical angle between 44 degrees and 54 degrees. The generated ensemble of TAR conformations have different organizations of functional groups responsible for interaction with the trans-activator protein Tat, including conformations similar to the previously characterized bound-state conformation. These results demonstrate the utility of RDC-NMR for simultaneously characterizing RNA tertiary dynamics and average conformation, and indicate an avenue for TAR complex formation involving tertiary structure capture.

Identifying hydrogen bond alignments in multistranded DNA architectures by NMR

NMR studies of nucleic acids have benefited tremendously from the discovery of trans-hydrogen-bond scalar coupling constants, which have enabled direct determination of N-H...N and N-H...O=C hydrogen bonds using a combination of (2h)J(NN)-, (4h)J(NN)-, and (3h)J(NC)-based spectroscopy. This is especially true of multistranded DNA architectures containing intricate hydrogen-bonded networks mediated primarily through mismatched base pairing, which often resist identification by posing serious technical, spectroscopic, and physicochemical challenges. In this Account, we present a suite of NMR pulse sequences that have been developed in our laboratory to address these issues. We demonstrate the utility of these methods for identifying hydrogen bonds in two quadruplex DNA structures, containing triad, tetrad, and hexad motifs involving Watson-Crick, G.G and sheared G.A mismatch base pairing.

Additional hydrogen bonds and base-pair kinetics in the symmetrical AMP-DNA aptamer complex

The solution structure of an adenosine monophosphate (AMP)-DNA aptamer complex has been determined previously [Lin, C. H., and Patel, D. J. (1997) Chem. Biol. 4:817-832]. On a symmetrical aptamer complex containing the same binding loop, but with better resolved spectra, we have identified two additional hydrogen bond-mediated associations in the binding loop. One of these involves a rapidly exchanging G imino proton. The phosphate group of the AMP ligand was identified as the acceptor by comparison with other aptamer complexes. Imino proton exchange measurements also yielded the dissociation constants of the stem and binding loop base pairs. This study shows that nuclear magnetic resonance-based imino proton exchange is a good probe for detection of weak hydrogen-bond associations.