Extensive backbone dynamics in the GCAA RNA tetraloop analyzed using 13C NMR spin relaxation and specific isotope labeling

Investigating the Intramolecular Competition of Different RNA Binding Motifs for Neomycin B by Native Top-Down Mass Spectrometry

The ongoing search for small molecule drugs that target ribonucleic acids (RNA) is complicated by a limited understanding of the principles that govern RNA-small molecule interactions. Here we have used stoichiometry-resolved native top-down mass spectrometry (MS) to study the binding of neomycin B to small model hairpin RNAs, an unstructured RNA, and a viral RNA construct. For 15-22 nt model RNAs with hairpin structure, we found that neomycin B binding to hairpin loops relies on interactions with both the nucleobases and the 2'-OH groups, and that a simple 5' or 3' overhang can introduce an additional binding motif. For a 47 nt RNA construct derived from stem IA of the human immunodeficiency virus 1 (HIV-1) rev response element (RRE) RNA, native top-down MS identified four different binding motifs, of which the purine-rich internal loop showed the highest affinity for neomycin B. Stoichiometry-resolved binding site mapping by native top-down MS allows for a new perspective on binding specificity, and has the potential to reveal unexpected principles of small molecule binding to RNA.

Sequence-Dependent Melting and Refolding Dynamics of RNA UNCG Tetraloops Using Temperature-Jump/Drop Infrared Spectroscopy

Time-resolved temperature-jump/drop infrared (IR) spectroscopy has been used to measure the impact of stem base sequence on the melting and refolding dynamics of ribonucleic acid (RNA) tetraloops. A series of three 12-nucleotide RNA hairpin sequences were studied, each featuring a UACG tetraloop motif and a double-stranded stem containing four base pairs. In each case, the stem comprised three GC pairs plus a single AU base pair inserted at the closing point of the loop (RNA_loop), in the middle of the stem (RNA_mid), or at the stem terminus (RNA_end). Results from analogous DNA tetraloop (TACG) sequences were also obtained. Inclusion of AU or AT base pairs in the stem leads to faster melting of the stem-loop structure compared to a stem sequence featuring four GC base pairs while refolding times were found to be slower, consistent with a general reduction in stem-loop stability caused by the AU/AT pair. Independent measurement of the dynamic timescales for melting and refolding of ring vibrational modes of guanine (G_R) and adenine (A_R) provided position-specific insight into hairpin dynamics. The G_R-derived data showed that DNA sequences melted more quickly (0.5 ± 0.1 to 0.7 ± 0.1 μs at 70 °C) than analogous RNA sequences (4.3 ± 0.4 to 4.4 ± 0.3 μs at 70 °C). Position-sensitive data from the A_R modes suggests that DNA hairpins begin melting from the terminal end of the stem toward the loop while RNA sequences begin melting from the loop. Refolding timescales for both RNA and DNA hairpins were found to be similar (250 ± 50 μs at 70 °C) except for RNA_end and DNA_loop which refolded much more slowly (746 ± 36 and 430 ± 31 μs, respectively), showing that the refolding pathway is significantly impaired by the placement of AU/AT pairs at different points in the stem. We conclude that conformational changes of analogous pairs of RNA and DNA tetraloops proceed by different mechanisms.

Characterizing the Conformational Free-Energy Landscape of RNA Stem-Loops Using Single-Molecule Field-Effect Transistors

We have developed and used single-molecule field-effect transistors (smFETs) to characterize the conformational free-energy landscape of RNA stem-loops. Stem-loops are one of the most common RNA structural motifs and serve as building blocks for the formation of complex RNA structures. Given their prevalence and integral role in RNA folding, the kinetics of stem-loop (un)folding has been extensively characterized using both experimental and computational approaches. Interestingly, these studies have reported vastly disparate timescales of (un)folding, which has been interpreted as evidence that (un)folding of even simple stem-loops occurs on a highly rugged conformational energy landscape. Because smFETs do not rely on fluorophore reporters of conformation or mechanical (un)folding forces, they provide a unique approach that has allowed us to directly monitor tens of thousands of (un)folding events of individual stem-loops at a 200 μs time resolution. Our results show that under our experimental conditions, stem-loops (un)fold over a 1-200 ms timescale during which they transition between ensembles of unfolded and folded conformations, the latter of which is composed of at least two sub-populations. The 1-200 ms timescale of (un)folding we observe here indicates that smFETs report on complete (un)folding trajectories in which unfolded conformations of the RNA spend long periods of time wandering the free-energy landscape before sampling one of several misfolded conformations or the natively folded conformation. Our findings highlight the extremely rugged landscape on which even the simplest RNA structural elements fold and demonstrate that smFETs are a unique and powerful approach for characterizing the conformational free-energy of RNA.

Isotope Labels Combined with Solution NMR Spectroscopy Make Visible the Invisible Conformations of Small-to-Large RNAs

RNA is central to the proper function of cellular processes important for life on earth and implicated in various medical dysfunctions. Yet, RNA structural biology lags significantly behind that of proteins, limiting mechanistic understanding of RNA chemical biology. Fortunately, solution NMR spectroscopy can probe the structural dynamics of RNA in solution at atomic resolution, opening the door to their functional understanding. However, NMR analysis of RNA, with only four unique ribonucleotide building blocks, suffers from spectral crowding and broad linewidths, especially as RNAs grow in size. One effective strategy to overcome these challenges is to introduce NMR-active stable isotopes into RNA. However, traditional uniform labeling methods introduce scalar and dipolar couplings that complicate the implementation and analysis of NMR measurements. This challenge can be circumvented with selective isotope labeling. In this review, we outline the development of labeling technologies and their application to study biologically relevant RNAs and their complexes ranging in size from 5 to 300 kDa by NMR spectroscopy.

Dynamics-Function Analysis in Catalytic RNA Using NMR Spin Relaxation and Conformationally Restricted Nucleotides

A full understanding of biomolecular function requires an analysis of both the dynamic properties of the system of interest and the identification of those dynamics that are required for function. We describe NMR methods based on metabolically directed specific isotope labeling for the identification of molecular disorder and/or conformational transitions on the RNA backbone ribose groups. These analyses are complemented by the use of synthetic covalently modified nucleotides constrained to a single sugar pucker, which allow functional assessment of dynamics by selectively removing a minor conformer identified by NMR from the structural ensemble.

2'-O-Trifluoromethylated RNA - a powerful modification for RNA chemistry and NMR spectroscopy

New RNA modifications are needed to advance our toolbox for targeted manipulation of RNA. In particular, the development of high-performance reporter groups facilitating spectroscopic analysis of RNA structure and dynamics, and of RNA-ligand interactions has attracted considerable interest. To this end, fluorine labeling in conjunction with 19F-NMR spectroscopy has emerged as a powerful strategy. Appropriate probes for RNA previously focused on single fluorine atoms attached to the 5-position of pyrimidine nucleobases or at the ribose 2'-position. To increase NMR sensitivity, trifluoromethyl labeling approaches have been developed, with the ribose 2'-SCF3 modification being the most prominent one. A major drawback of the 2'-SCF3 group, however, is its strong impact on RNA base pairing stability. Interestingly, RNA containing the structurally related 2'-OCF3 modification has not yet been reported. Therefore, we set out to overcome the synthetic challenges toward 2'-OCF3 labeled RNA and to investigate the impact of this modification. We present the syntheses of 2'-OCF3 adenosine and cytidine phosphoramidites and their incorporation into oligoribonucleotides by solid-phase synthesis. Importantly, it turns out that the 2'-OCF3 group has only a slight destabilizing effect when located in double helical regions which is consistent with the preferential C3'-endo conformation of the 2'-OCF3 ribose as reflected in the 3 J (H1'-H2') coupling constants. Furthermore, we demonstrate the exceptionally high sensitivity of the new label in 19F-NMR analysis of RNA structure equilibria and of RNA-small molecule interactions. The study is complemented by a crystal structure at 0.9 Å resolution of a 27 nt hairpin RNA containing a single 2'-OCF3 group that well integrates into the minor groove. The new label carries high potential to outcompete currently applied fluorine labels for nucleic acid NMR spectroscopy because of its significantly advanced performance.

Deleterious effects of carbon-carbon dipolar coupling on RNA NMR dynamics

Many regulatory RNAs undergo dynamic exchanges that are crucial for their biological functions and NMR spectroscopy is a versatile tool for monitoring dynamic motions of biomolecules. Meaningful information on biomolecular dynamics requires an accurate measurement of relaxation parameters such as longitudinal (R₁) rates, transverse (R₂) rates and heteronuclear Overhauser effect (hNOE). However, earlier studies have shown that the large ¹³C-¹³C interactions complicate analysis of the carbon relaxation parameters. To investigate the effect of ¹³C-¹³C interactions on RNA dynamic studies, we performed relaxation measurements on various RNA samples with different labeling patterns and compared these measurements with the computational simulations. For uniformly labeled samples, contributions of the neighboring carbon to R₁ measurements were observed. These contributions increased with increasing magnetic field and overall correlation time ([Formula: see text]) for R₁ rates, necessitating more careful analysis for uniformly labeled large RNAs. In addition, the hNOE measurements were also affected by the adjacent carbon nuclei. Unlike R₁ rates, R_1ρ rates showed relatively good agreement between uniformly- and site-selectively labeled samples, suggesting no dramatic effect from their attached carbon, in agreement with previous observations. Overall, having more accurate rate measurements avoids complex analysis and will be a key for interpreting ¹³C relaxation rates for molecular motion that can provide valuable insights into cellular molecular recognition events.

Benefits of stable isotope labeling in RNA analysis

RNAs are key players in life as they connect the genetic code (DNA) with all cellular processes dominated by proteins. They contain a variety of chemical modifications and many RNAs fold into complex structures. Here, we review recent progress in the analysis of RNA modification and structure on the basis of stable isotope labeling techniques. Mass spectrometry (MS) and nuclear magnetic resonance (NMR) spectroscopy are the key tools and many breakthrough developments were made possible by the analysis of stable isotope labeled RNA. Therefore, we discuss current stable isotope labeling techniques such as metabolic labeling, enzymatic labeling and chemical synthesis. RNA structure analysis by NMR is challenging due to two major problems that become even more salient when the size of the RNA increases, namely chemical shift overlaps and line broadening leading to complete signal loss. Several isotope labeling strategies have been developed to provide solutions to these major issues, such as deuteration, segmental isotope labeling or site-specific labeling. Quantification of modified nucleosides in RNA by MS is only possible through the application of stable isotope labeled internal standards. With nucleic acid isotope labeling coupled mass spectrometry (NAIL-MS), it is now possible to analyze the dynamic processes of post-transcriptional RNA modification and demodification. The trend, in both NMR and MS RNA analytics, is without doubt shifting from the analysis of snapshot moments towards the development and application of tools capable of analyzing the dynamics of RNA structure and modification profiles.

Probing excited conformational states of nucleic acids by nitrogen CEST NMR spectroscopy

Characterizing low-populated and short-lived excited conformational states has become increasingly important for understanding mechanisms of RNA function. Interconversion between RNA ground and excited conformational states often involves base pairing rearrangements that lead to changes in the hydrogen-bond network. Here, we present two 15N chemical exchange saturation transfer (CEST) NMR experiments that utilize protonated and non-protonated nitrogens, which are key hydrogen-bond donors and acceptors, for characterizing excited conformational states in RNA. We demonstrated these approaches on the B. Cereus fluoride riboswitch, where 15N CEST profiles complement 13C CEST profiles in depicting a potential pathway for ligand-dependent allosteric regulation of the excited conformational state of the fluoride riboswitch.

RNA Dynamics by NMR Spectroscopy

An ever-increasing number of functional RNAs require a mechanistic understanding. RNA function relies on changes in its structure, so-called dynamics. To reveal dynamic processes and higher energy structures, new NMR methods have been developed to elucidate these dynamics in RNA with atomic resolution. In this Review, we provide an introduction to dynamics novices and an overview of methods that access most dynamic timescales, from picoseconds to hours. Examples are provided as well as insight into theory, data acquisition and analysis for these different methods. Using this broad spectrum of methodology, unprecedented detail and invisible structures have been obtained and are reviewed here. RNA, though often more complicated and therefore neglected, also provides a great system to study structural changes, as these RNA structural changes are more easily defined-Lego like-than in proteins, hence the numerous revelations of RNA excited states.

Characterizing micro-to-millisecond chemical exchange in nucleic acids using off-resonance R1ρ relaxation dispersion

This review describes off-resonance R1ρ relaxation dispersion NMR methods for characterizing microsecond-to-millisecond chemical exchange in uniformly 13C/15N labeled nucleic acids in solution. The review opens with a historical account of key developments that formed the basis for modern R1ρ techniques used to study chemical exchange in biomolecules. A vector model is then used to describe the R1ρ relaxation dispersion experiment, and how the exchange contribution to relaxation varies with the amplitude and frequency offset of an applied spin-locking field, as well as the population, exchange rate, and differences in chemical shifts of two exchanging species. Mathematical treatment of chemical exchange based on the Bloch-McConnell equations is then presented and used to examine relaxation dispersion profiles for more complex exchange scenarios including three-state exchange. Pulse sequences that employ selective Hartmann-Hahn cross-polarization transfers to excite individual 13C or 15N spins are then described for measuring off-resonance R1ρ(13C) and R1ρ(15N) in uniformly 13C/15N labeled DNA and RNA samples prepared using commercially available 13C/15N labeled nucleotide triphosphates. Approaches for analyzing R1ρ data measured at a single static magnetic field to extract a full set of exchange parameters are then presented that rely on numerical integration of the Bloch-McConnell equations or the use of algebraic expressions. Methods for determining structures of nucleic acid excited states are then reviewed that rely on mutations and chemical modifications to bias conformational equilibria, as well as structure-based approaches to calculate chemical shifts. Applications of the methodology to the study of DNA and RNA conformational dynamics are reviewed and the biological significance of the exchange processes is briefly discussed.

NMR characterization of RNA small molecule interactions

Exciting discoveries of naturally occurring ligand-sensing and disease-linked noncoding RNAs have promoted significant interests in understanding RNA-small molecule interactions. NMR spectroscopy is a powerful tool for characterizing intermolecular interactions. In this review, we describe protocols and approaches for applying NMR spectroscopy to investigate interactions between RNA and small molecules. We review protocols for RNA sample preparation, methods for identifying RNA-binding small molecules, approaches for mapping RNA-small molecule interactions, determining complex structures, and characterizing binding kinetics. We hope this review will provide a guideline to streamline NMR applications in studying RNA-small molecule interactions, facilitating both basic mechanistic understandings of RNA functions and translational efforts in developing RNA-targeted therapeutics.

Advances that facilitate the study of large RNA structure and dynamics by nuclear magnetic resonance spectroscopy

The characterization of functional yet nonprotein coding (nc) RNAs has expanded the role of RNA in the cell from a passive player in the central dogma of molecular biology to an active regulator of gene expression. The misregulation of ncRNA function has been linked with a variety of diseases and disorders ranging from cancers to neurodegeneration. However, a detailed molecular understanding of how ncRNAs function has been limited; due, in part, to the difficulties associated with obtaining high-resolution structures of large RNAs. Tertiary structure determination of RNA as a whole is hampered by various technical challenges, all of which are exacerbated as the size of the RNA increases. Namely, RNAs tend to be highly flexible and dynamic molecules, which are difficult to crystallize. Biomolecular nuclear magnetic resonance (NMR) spectroscopy offers a viable alternative to determining the structure of large RNA molecules that do not readily crystallize, but is itself hindered by some technical limitations. Recently, a series of advancements have allowed the biomolecular NMR field to overcome, at least in part, some of these limitations. These advances include improvements in sample preparation strategies as well as methodological improvements. Together, these innovations pave the way for the study of ever larger RNA molecules that have important biological function. This article is categorized under: RNA Structure and Dynamics > RNA Structure, Dynamics, and Chemistry Regulatory RNAs/RNAi/Riboswitches > Regulatory RNAs RNA Structure and Dynamics > Influence of RNA Structure in Biological Systems.

Coupling between conformational dynamics and catalytic function at the active site of the lead-dependent ribozyme

In common with other self-cleaving RNAs, the lead-dependent ribozyme (leadzyme) undergoes dynamic fluctuations to a chemically activated conformation. We explored the connection between conformational dynamics and self-cleavage function in the leadzyme using a combination of NMR spin-relaxation analysis of ribose groups and conformational restriction via chemical modification. The functional studies were performed with a North-methanocarbacytidine modification that prevents fluctuations to C2'-endo conformations while maintaining an intact 2'-hydroxyl nucleophile. Spin-relaxation data demonstrate that the active-site Cyt-6 undergoes conformational exchange attributed to sampling of a minor C2'-endo state with an exchange lifetime on the order of microseconds to tens of microseconds. A conformationally restricted species in which the fluctuations to the minor species are interrupted shows a drastic decrease in self-cleavage activity. Taken together, these data indicate that dynamic sampling of a minor species at the active site of this ribozyme, and likely of related naturally occurring motifs, is strongly coupled to catalytic function. The combination of NMR dynamics analysis with functional probing via conformational restriction is a general methodology for dissecting dynamics-function relationships in RNA.

NMR probing of invisible excited states using selectively labeled RNAs

Carr-Purcell-Meiboom-Gill (CPMG) relaxation dispersion NMR experiments are invaluable for probing sparsely and transiently populated biomolecular states that cannot be directly detected by traditional NMR experiments and that are invisible by other biophysical approaches. A notable gap for RNA is the absence of CPMG experiments for measurement of methine base ¹H and methylene C5' chemical shifts of ribose moieties in the excited state, partly because of complications from homonuclear ¹³C-¹³C scalar couplings. Here we present site-specific ¹³C labeling that makes possible the design of pulse sequences for recording accurate ¹H-¹³C MQ and SQ CPMG experiments for ribose methine H1'-C1' and H2'-C2', base and ribose ¹H CPMG, as well as a new ¹H-¹³C TROSY-detected methylene (CH₂) C5' CPMG relaxation pulse schemes. We demonstrate the utility of these experiments for two RNAs, the A-Site RNA known to undergo exchange and the IRE RNA suspected of undergoing exchange on microseconds to millisecond time-scale. We anticipate the new labeling approaches will facilitate obtaining structures of invisible states and provide insights into the relevance of such states for RNA-drug interactions.

Studying sparsely populated conformational states in RNA combining chemical synthesis and solution NMR spectroscopy

Using chemical synthesis and solution NMR spectroscopy, RNA structural ensembles including a major ground state and minor populated excited states can be studied at atomic resolution. In this work, atom-specific ¹³C labeled RNA building blocks - a 5-¹³C-uridine and a 2,8-¹³C₂-adenosine building block - are used to introduce isolated ¹³C-¹H-spin topologies into a target RNA to probe such structural ensembles via NMR spectroscopy. First, the 5-¹³C-uridine 2'-O-TBDMS-phosphoramidite building block was introduced into a 21 nucleotide (nt) tP5c stem construct of the tP5abc subdomain of the Tetrahymena group I ribozyme. Then, the 2,8-¹³C₂-adenosine 2'-O-TBDMS-phosphoramidite building block was incorporated into a 9 kDa and a 15 kD construct derived from the epsilon (ε) RNA element of the duck Hepatitis B virus. The 2,8-¹³C₂-adenosine resonances of the 9 kDa 28 nt sequence could be mapped to the full-length 53 nt construct. The isolated NMR active nuclei pairs were used to probe for low populated excited states (<10%) via ¹³C-Carr-Purcell-Meiboom-Gill (CPMG)-relaxation dispersion NMR spectroscopy. The ¹³C-CPMG relaxation dispersion experiment recapitulated a secondary structure switching event in the P5c hairpin of the group I intron construct previously revealed by ¹⁵N relaxation dispersion experiments. In the ε-HBV RNA an unfolding event occurring on the millisecond time scale was found in the upper stem in-line with earlier observations. This unpaired conformational state is presumed to be important for the binding of the epsilon reverse transcriptase (RT) enzyme. Thus, a full description of an RNA's folding landscape helps to obtain a deeper understanding of its function, as these high energy conformational states often represent functionally important intermediates involved in (un)folding or ribozyme catalysis.

RNA Structural Dynamics As Captured by Molecular Simulations: A Comprehensive Overview

With both catalytic and genetic functions, ribonucleic acid (RNA) is perhaps the most pluripotent chemical species in molecular biology, and its functions are intimately linked to its structure and dynamics. Computer simulations, and in particular atomistic molecular dynamics (MD), allow structural dynamics of biomolecular systems to be investigated with unprecedented temporal and spatial resolution. We here provide a comprehensive overview of the fast-developing field of MD simulations of RNA molecules. We begin with an in-depth, evaluatory coverage of the most fundamental methodological challenges that set the basis for the future development of the field, in particular, the current developments and inherent physical limitations of the atomistic force fields and the recent advances in a broad spectrum of enhanced sampling methods. We also survey the closely related field of coarse-grained modeling of RNA systems. After dealing with the methodological aspects, we provide an exhaustive overview of the available RNA simulation literature, ranging from studies of the smallest RNA oligonucleotides to investigations of the entire ribosome. Our review encompasses tetranucleotides, tetraloops, a number of small RNA motifs, A-helix RNA, kissing-loop complexes, the TAR RNA element, the decoding center and other important regions of the ribosome, as well as assorted others systems. Extended sections are devoted to RNA-ion interactions, ribozymes, riboswitches, and protein/RNA complexes. Our overview is written for as broad of an audience as possible, aiming to provide a much-needed interdisciplinary bridge between computation and experiment, together with a perspective on the future of the field.

Resolving sugar puckers in RNA excited states exposes slow modes of repuckering dynamics

Recent studies have shown that RNAs exist in dynamic equilibrium with short-lived low-abundance 'excited states' that form by reshuffling base pairs in and around non-canonical motifs. These conformational states are proposed to be rich in non-canonical motifs and to play roles in the folding and regulatory functions of non-coding RNAs but their structure proves difficult to characterize given their transient nature. Here, we describe an approach for determining sugar pucker conformation in RNA excited states through nuclear magnetic resonance measurements of C1΄ and C4΄ rotating frame spin relaxation (R1ρ) in uniformly 13C/15N labeled RNA samples. Application to HIV-1 TAR exposed slow modes of sugar repuckering dynamics at the μs and ms timescale accompanying transitions between non-helical (C2΄-endo) to helical (C3΄-endo) conformations during formation of two distinct excited states. In contrast, we did not obtain any evidence for slow sugar repuckering dynamics for nucleotides in a variety of structural contexts that do not undergo non-helical to helical transitions. Our results outline a route for significantly improving the conformational characterization of RNA excited states and suggest that slow modes of repuckering dynamics gated by transient changes in secondary structure are quite common in RNA.

Intrinsic Base-Pair Rearrangement in the Hairpin Ribozyme Directs RNA Conformational Sampling and Tertiary Interface Formation

Dynamic fluctuations in RNA structure enable conformational changes that are required for catalysis and recognition. In the hairpin ribozyme, the catalytically active structure is formed as an intricate tertiary interface between two RNA internal loops. Substantial alterations in the structure of each loop are observed upon interface formation, or docking. The very slow on-rate for this relatively tight interaction has led us to hypothesize a double conformational capture mechanism for RNA-RNA recognition. We used extensive molecular dynamics simulations to assess conformational sampling in the undocked form of the loop domain containing the scissile phosphate (loop A). We observed several major accessible conformations with distinctive patterns of hydrogen bonding and base stacking interactions in the active-site internal loop. Several important conformational features characteristic of the docked state were observed in well-populated substates, consistent with the kinetic sampling of docking-competent states by isolated loop A. Our observations suggest a hybrid or multistage binding mechanism, in which initial conformational selection of a docking-competent state is followed by induced-fit adjustment to an in-line, chemically reactive state only after formation of the initial complex with loop B.

Visualizing the formation of an RNA folding intermediate through a fast highly modular secondary structure switch

Intermediates play important roles in RNA folding but can be difficult to characterize when short-lived or not significantly populated. By combining (15)N relaxation dispersion NMR with chemical probing, we visualized a fast (kex=k1+k-1≈423 s(-1)) secondary structural switch directed towards a low-populated (∼3%) partially folded intermediate in tertiary folding of the P5abc subdomain of the 'Tetrahymena' group I intron ribozyme. The secondary structure switch changes the base-pairing register across the P5c hairpin, creating a native-like structure, and occurs at rates of more than two orders of magnitude faster than tertiary folding. The switch occurs robustly in the absence of tertiary interactions, Mg(2+) or even when the hairpin is excised from the three-way junction. Fast, highly modular secondary structural switches may be quite common during RNA tertiary folding where they may help smoothen the folding landscape by allowing folding to proceed efficiently via additional pathways.

Free-energy landscape of a hyperstable RNA tetraloop

We report the characterization of the energy landscape and the folding/unfolding thermodynamics of a hyperstable RNA tetraloop obtained through high-performance molecular dynamics simulations at microsecond timescales. Sampling of the configurational landscape is conducted using temperature replica exchange molecular dynamics over three isochores at high, ambient, and negative pressures to determine the thermodynamic stability and the free-energy landscape of the tetraloop. The simulations reveal reversible folding/unfolding transitions of the tetraloop into the canonical A-RNA conformation and the presence of two alternative configurations, including a left-handed Z-RNA conformation and a compact purine Triplet. Increasing hydrostatic pressure shows a stabilizing effect on the A-RNA conformation and a destabilization of the left-handed Z-RNA. Our results provide a comprehensive description of the folded free-energy landscape of a hyperstable RNA tetraloop and highlight the significant advances of all-atom molecular dynamics in describing the unbiased folding of a simple RNA secondary structure motif.

Chemo-enzymatic synthesis of site-specific isotopically labeled nucleotides for use in NMR resonance assignment, dynamics and structural characterizations

Stable isotope labeling is central to NMR studies of nucleic acids. Development of methods that incorporate labels at specific atomic positions within each nucleotide promises to expand the size range of RNAs that can be studied by NMR. Using recombinantly expressed enzymes and chemically synthesized ribose and nucleobase, we have developed an inexpensive, rapid chemo-enzymatic method to label ATP and GTP site specifically and in high yields of up to 90%. We incorporated these nucleotides into RNAs with sizes ranging from 27 to 59 nucleotides using in vitro transcription: A-Site (27 nt), the iron responsive elements (29 nt), a fluoride riboswitch from Bacillus anthracis (48 nt), and a frame-shifting element from a human corona virus (59 nt). Finally, we showcase the improvement in spectral quality arising from reduced crowding and narrowed linewidths, and accurate analysis of NMR relaxation dispersion (CPMG) and TROSY-based CEST experiments to measure μs-ms time scale motions, and an improved NOESY strategy for resonance assignment. Applications of this selective labeling technology promises to reduce difficulties associated with chemical shift overlap and rapid signal decay that have made it challenging to study the structure and dynamics of large RNAs beyond the 50 nt median size found in the PDB.

Measuring Residual Dipolar Couplings in Excited Conformational States of Nucleic Acids by CEST NMR Spectroscopy

Nucleic acids undergo structural transitions to access sparsely populated and transiently lived conformational states--or excited conformational states--that play important roles in diverse biological processes. Despite ever-increasing detection of these functionally essential states, 3D structure determination of excited states (ESs) of RNA remains elusive. This is largely due to challenges in obtaining high-resolution structural constraints in these ESs by conventional structural biology approaches. Here, we present nucleic-acid-optimized chemical exchange saturation transfer (CEST) NMR spectroscopy for measuring residual dipolar couplings (RDCs), which provide unique long-range angular constraints in ESs of nucleic acids. We demonstrate these approaches on a fluoride riboswitch, where one-bond (13)C-(1)H RDCs from both base and sugar moieties provide direct structural probes into an ES of the ligand-free riboswitch.

Stable isotope-labeled RNA phosphoramidites to facilitate dynamics by NMR

Given that Ribonucleic acids (RNAs) are a central hub of various cellular processes, methods to synthesize these RNAs for biophysical studies are much needed. Here, we showcase the applicability of 6-(13)C-pyrimidine phosphoramidites to introduce isolated (13)C-(1)H spin pairs into RNAs up to 40 nucleotides long. The method allows the incorporation of 6-(13)C-uridine and -cytidine residues at any desired position within a target RNA. By site-specific positioning of the (13)C-label using RNA solid phase synthesis, these stable isotope-labeling patterns are especially well suited to resolve resonance assignment ambiguities. Of even greater importance, the labeling pattern affords accurate quantification of important functional transitions of biologically relevant RNAs (e.g., riboswitch aptamer domains, viral RNAs, or ribozymes) in the μs- to ms time regime and beyond without complications of one bond carbon scalar couplings. We outline the chemical synthesis of the 6-(13)C-pyrimidine building blocks and their use in RNA solid phase synthesis and demonstrate their utility in Carr Purcell Meiboom Gill relaxation dispersion, ZZ exchange, and chemical exchange saturation transfer NMR experiments.

Assessment of chemical exchange in tryptophan-albumin solution through (19)F multicomponent transverse relaxation dispersion analysis

A number of NMR methods possess the capability of probing chemical exchange dynamics in solution. However, certain drawbacks limit the applications of these NMR approaches, particularly, to a complex system. Here, we propose a procedure that integrates the regularized nonnegative least squares (NNLS) analysis of multiexponential T2 relaxation into Carr-Purcell-Meiboom-Gill (CPMG) relaxation dispersion experiments to probe chemical exchange in a multicompartmental system. The proposed procedure was validated through analysis of (19)F T2 relaxation data of 6-fluoro-DL-tryptophan in a two-compartment solution with and without bovine serum albumin. Given the regularized NNLS analysis of a T2 relaxation curve acquired, for example, at the CPMG frequency υ CPMG  = 125, the nature of two distinct peaks in the associated T2 distribution spectrum indicated 6-fluoro-DL-tryptophan either retaining the free state, with geometric mean */multiplicative standard deviation (MSD) = 1851.2 ms */1.51, or undergoing free/albumin-bound interconversion, with geometric mean */MSD = 236.8 ms */1.54, in the two-compartment system. Quantities of the individual tryptophan species were accurately reflected by the associated T2 peak areas, with an interconversion state-to-free state ratio of 0.45 ± 0.11. Furthermore, the CPMG relaxation dispersion analysis estimated the exchange rate between the free and albumin-bound states in this fluorinated tryptophan analog and the corresponding dissociation constant of the fluorinated tryptophan-albumin complex in the chemical-exchanging, two-compartment system.

Characterizing RNA Excited States Using NMR Relaxation Dispersion

Changes in RNA secondary structure play fundamental roles in the cellular functions of a growing number of noncoding RNAs. This chapter describes NMR-based approaches for characterizing microsecond-to-millisecond changes in RNA secondary structure that are directed toward short-lived and low-populated species often referred to as "excited states." Compared to larger scale changes in RNA secondary structure, transitions toward excited states do not require assistance from chaperones, are often orders of magnitude faster, and are localized to a small number of nearby base pairs in and around noncanonical motifs. Here, we describe a procedure for characterizing RNA excited states using off-resonance R1ρ NMR relaxation dispersion utilizing low-to-high spin-lock fields (25-3000 Hz). R1ρ NMR relaxation dispersion experiments are used to measure carbon and nitrogen chemical shifts in base and sugar moieties of the excited state. The chemical shift data are then interpreted with the aid of secondary structure prediction to infer potential excited states that feature alternative secondary structures. Candidate structures are then tested by using mutations, single-atom substitutions, or by changing physiochemical conditions, such as pH and temperature, to either stabilize or destabilize the candidate excited state. The resulting chemical shifts of the mutants or under different physiochemical conditions are then compared to those of the ground and excited states. Application is illustrated with a focus on the transactivation response element from the human immune deficiency virus type 1, which exists in dynamic equilibrium with at least two distinct excited states.

Characterizing excited conformational states of RNA by NMR spectroscopy

Conformational dynamics is a hallmark of diverse non-coding RNA functions. During these functional processes, RNA molecules almost ubiquitously undergo conformational transitions that are tuned to meet distinct structural and kinetic requirements for proper function. A complete mechanistic understanding of RNA function requires comprehensive structural and dynamic knowledge of these complex transitions, which often involve alternative higher-energy conformational states that pose a major challenge for high-resolution structural study by conventional methods. In this review, we describe recent progress in RNA NMR that has started to unveil detailed structural, thermodynamic and kinetic insights into some of these excited conformational states of RNA and their functional roles in biology.

Dynamic motions of the HIV-1 frameshift site RNA

The HIV-1 frameshift site (FS) plays a critical role in viral replication. During translation, the HIV-1 FS transitions from a 3-helix to a 2-helix junction RNA secondary structure. The 2-helix junction structure contains a GGA bulge, and purine-rich bulges are common motifs in RNA secondary structure. Here, we investigate the dynamics of the HIV-1 FS 2-helix junction RNA. Interhelical motions were studied under different ionic conditions using NMR order tensor analysis of residual dipolar couplings. In 150 mM potassium, the RNA adopts a 43°(±4°) interhelical bend angle (β) and displays large amplitude, anisotropic interhelical motions characterized by a 0.52(±0.04) internal generalized degree of order (GDOint) and distinct order tensor asymmetries for its two helices (η = 0.26(±0.04) and 0.5(±0.1)). These motions are effectively quenched by addition of 2 mM magnesium (GDOint = 0.87(±0.06)), which promotes a near-coaxial conformation (β = 15°(±6°)) of the two helices. Base stacking in the bulge was investigated using the fluorescent purine analog 2-aminopurine. These results indicate that magnesium stabilizes extrahelical conformations of the bulge nucleotides, thereby promoting coaxial stacking of helices. These results are highly similar to previous studies of the HIV transactivation response RNA, despite a complete lack of sequence similarity between the two RNAs. Thus, the conformational space of these RNAs is largely determined by the topology of their interhelical junctions.

Chemical exchange in biomacromolecules: past, present, and future

The perspective reviews quantitative investigations of chemical exchange phenomena in proteins and other biological macromolecules using NMR spectroscopy, particularly relaxation dispersion methods. The emphasis is on techniques and applications that quantify the populations, interconversion kinetics, and structural features of sparsely populated conformational states in equilibrium with a highly populated ground state. Applications to folding, molecular recognition, catalysis, and allostery by proteins and nucleic acids are highlighted.

Hierarchy of RNA functional dynamics

RNA dynamics play a fundamental role in many cellular functions. However, there is no general framework to describe these complex processes, which typically consist of many structural maneuvers that occur over timescales ranging from picoseconds to seconds. Here, we classify RNA dynamics into distinct modes representing transitions between basins on a hierarchical free-energy landscape. These transitions include large-scale secondary-structural transitions at >0.1-s timescales, base-pair/tertiary dynamics at microsecond-to-millisecond timescales, stacking dynamics at timescales ranging from nanoseconds to microseconds, and other "jittering" motions at timescales ranging from picoseconds to nanoseconds. We review various modes within these three different tiers, the different mechanisms by which they are used to regulate function, and how they can be coupled together to achieve greater functional complexity.

The free energy landscape of pseudorotation in 3'-5' and 2'-5' linked nucleic acids

The five-membered furanose ring is a central component of the chemical structure of biological nucleic acids. The conformations of the furanose ring can be analytically described using the concept of pseudorotation, and for RNA and DNA they are dominated by the C2'-endo and C3'-endo conformers. While the free energy difference between these two conformers can be inferred from NMR measurements, a free energy landscape of the complete pseudorotation cycle of nucleic acids in solution has remained elusive. Here, we describe a new free energy calculation method for molecular dynamics (MD) simulations using the two pseudorotation parameters directly as the collective variables. To validate our approach, we calculated the free energy surface of ribose pseudorotation in guanosine and 2'-deoxyguanosine. The calculated free energy landscape reveals not only the relative stability of the different pseudorotation conformers, but also the main transition path between the stable conformations. Applying this method to a standard A-form RNA duplex uncovered the expected minimum at the C3'-endo state. However, at a 2'-5' linkage, the minimum shifts to the C2'-endo conformation. The free energy of the C3'-endo conformation is 3 kcal/mol higher due to a weaker hydrogen bond and a reduced base stacking interaction. Unrestrained MD simulations suggest that the conversion from C3'-endo to C2'-endo and vice versa is on the nanosecond and microsecond time scale, respectively. These calculations suggest that 2'-5' linkages may enable folded RNAs to sample a wider spectrum of their pseudorotation conformations.

Role of helical constraints of the EBS1-IBS1 duplex of a group II intron on demarcation of the 5' splice site

Recognition of the 5' splice site by group II introns involves pairing between an exon binding sequence (EBS) 1 within the ID3 stem-loop of domain 1 and a complementary sequence at the 3' end of exon 1 (IBS1). To identify the molecular basis for splice site definition of a group IIB ai5γ intron, we probed the solution structure of the ID3 stem-loop alone and upon binding of its IBS1 target by solution NMR. The ID3 stem was structured. The base of the ID3 loop was stacked but displayed a highly flexible EBS1 region. The flexibility of EBS1 appears to be a general feature of the ai5γ and the smaller Oceanobacillus iheyensis (O.i.) intron and may help in effective search of conformational space and prevent errors in splicing as a result of fortuitous base-pairing. Binding of IBS1 results in formation of a structured seven base pair duplex that terminates at the 5' splice site in spite of the potential for additional A-U and G•U pairs. Comparison of these data with conformational features of EBS1-IBS1 duplexes extracted from published structures suggests that termination of the duplex and definition of the splice site are governed by constraints of the helical geometry within the ID3 loop. This feature and flexibility of the uncomplexed ID3 loop appear to be common for both the ai5γ and O.i. introns and may help to fine-tune elements of recognition in group II introns.

Characterizing slow chemical exchange in nucleic acids by carbon CEST and low spin-lock field R(1ρ) NMR spectroscopy

Quantitative characterization of dynamic exchange between various conformational states provides essential insights into the molecular basis of many regulatory RNA functions. Here, we present an application of nucleic-acid-optimized carbon chemical exchange saturation transfer (CEST) and low spin-lock field R(1ρ) relaxation dispersion (RD) NMR experiments in characterizing slow chemical exchange in nucleic acids that is otherwise difficult if not impossible to be quantified by the ZZ-exchange NMR experiment. We demonstrated the application on a 47-nucleotide fluoride riboswitch in the ligand-free state, for which CEST and R(1ρ) RD profiles of base and sugar carbons revealed slow exchange dynamics involving a sparsely populated (p ~ 10%) and shortly lived (τ ~ 10 ms) NMR "invisible" state. The utility of CEST and low spin-lock field R(1ρ) RD experiments in studying slow exchange was further validated in characterizing an exchange as slow as ~60 s(-1).

NMR studies of nucleic acid dynamics

Nucleic acid structures have to satisfy two diametrically opposite requirements; on one hand they have to adopt well-defined 3D structures that can be specifically recognized by proteins; on the other hand, their structures must be sufficiently flexible to undergo very large conformational changes that are required during key biochemical processes, including replication, transcription, and translation. How do nucleic acids introduce flexibility into their 3D structure without losing biological specificity? Here, I describe the development and application of NMR spectroscopic techniques in my laboratory for characterizing the dynamic properties of nucleic acids that tightly integrate a broad set of NMR measurements, including residual dipolar couplings, spin relaxation, and relaxation dispersion with sample engineering and computational approaches. This approach allowed us to obtain fundamental new insights into directional flexibility in nucleic acids that enable their structures to change in a very specific functional manner.

A personal perspective on chemistry-driven RNA research

In this mini review, we discuss how our understanding of ribonucleic acid (RNA) properties becomes significantly deepened when a broad range of modern chemical and biophysical methods is applied. We span our perspective from RNA solid-phase synthesis and site-specific labeling to single-molecule fluorescence-resonance-energy-transfer imaging and NMR spectroscopy approaches to explore the dynamics of RNA over a broad timescale. We then move on to Fourier-transform-ion-cyclotron-resonance mass spectrometry (FT-ICR-MS) as a powerful technique for RNA sequencing and modification analysis. The novel methodological developments are discussed for selected biological systems that include the thiamine-pyrophosphate riboswitch, HIV and ribosomal A-site RNA, and transfer RNA.

High-resolution reversible folding of hyperstable RNA tetraloops using molecular dynamics simulations

We report the de novo folding of three hyperstable RNA tetraloops to 1-3 Å rmsd from their experimentally determined structures using molecular dynamics simulations initialized in the unfolded state. RNA tetraloops with loop sequences UUCG, GCAA, or CUUG are hyperstable because of the formation of noncanonical loop-stabilizing interactions, and they are all faithfully reproduced to angstrom-level accuracy in replica exchange molecular dynamics simulations, including explicit solvent and ion molecules. This accuracy is accomplished using unique RNA parameters, in which biases that favor rigid, highly stacked conformations are corrected to accurately capture the inherent flexibility of ssRNA loops, accurate base stacking energetics, and purine syn-anti interconversions. In a departure from traditional quantum chemistrycentric approaches to force field optimization, our parameters are calibrated directly from thermodynamic and kinetic measurements of intra- and internucleotide structural transitions. The ability to recapitulate the signature noncanonical interactions of the three most abundant hyperstable stem loop motifs represents a significant milestone to the accurate prediction of RNA tertiary structure using unbiased all-atom molecular dynamics simulations.

Visualizing transient low-populated structures of RNA

The visualization of RNA conformational changes has provided fundamental insights into how regulatory RNAs carry out their biological functions. The RNA structural transitions that have been characterized to date involve long-lived species that can be captured by structure characterization techniques. Here, we report the Nuclear Magnetic Resonance visualization of RNA transitions towards invisible ‘excited states’ (ES), which exist in too little abundance (2–13%) and for too short periods of time (45–250 μs) to allow structural characterization by conventional techniques. Transitions towards ESs result in localized rearrangements in base-pairing that alter building block elements of RNA architecture, including helix-junction-helix motifs and apical loops. The ES can inhibit function by sequestering residues involved in recognition and signaling or promote ATP-independent strand exchange. Thus, RNAs do not adopt a single conformation, but rather exist in rapid equilibrium with alternative ESs, which can be stabilized by cellular cues to affect functional outcomes.

Csy4 relies on an unusual catalytic dyad to position and cleave CRISPR RNA

CRISPR-Cas adaptive immune systems protect prokaryotes against foreign genetic elements. crRNAs derived from CRISPR loci base pair with complementary nucleic acids, leading to their destruction. In Pseudomonas aeruginosa, crRNA biogenesis requires the endoribonuclease Csy4, which binds and cleaves the repetitive sequence of the CRISPR transcript. Biochemical assays and three co-crystal structures of wild-type and mutant Csy4/RNA complexes reveal a substrate positioning and cleavage mechanism in which a histidine deprotonates the ribosyl 2'-hydroxyl pinned in place by a serine, leading to nucleophilic attack on the scissile phosphate. The active site catalytic dyad lacks a general acid to protonate the leaving group and positively charged residues to stabilize the transition state, explaining why the observed catalytic rate constant is ∼10(4)-fold slower than that of RNase A. We show that this RNA cleavage step is essential for assembly of the Csy protein-crRNA complex that facilitates target recognition. Considering that Csy4 recognizes a single cellular substrate and sequesters the cleavage product, evolutionary pressure has likely selected for substrate specificity and high-affinity crRNA interactions at the expense of rapid cleavage kinetics.

Synthesis of (6-(13)C)pyrimidine nucleotides as spin-labels for RNA dynamics

We present a (13)C-based isotope labeling protocol for RNA. Using (6-(13)C)pyrimidine phosphoramidite building blocks, site-specific labels can be incorporated into a target RNA via chemical oligonucleotide solid-phase synthesis. This labeling scheme is particularly useful for studying milli- to microsecond dynamics via NMR spectroscopy, as an isolated spin system is a crucial prerequisite to apply Carr-Purcell-Meiboom-Gill (CPMG) relaxation dispersion type experiments. We demonstrate the applicability for the characterization and detection of functional dynamics on various time scales by incorporating the (6-(13)C)uridine and -cytidine labels into biologically relevant RNAs. The refolding kinetics of a bistable terminator antiterminator segment involved in the gene regulation process controlled by the preQ(1) riboswitch class I was investigated. Using (13)C CPMG relaxation dispersion NMR spectroscopy, the milli- to microsecond dynamics of the HIV-1 transactivation response element RNA and the Varkud satellite stem loop V motif was addressed.

Biomass production of site selective 13C/15N nucleotides using wild type and a transketolase E. coli mutant for labeling RNA for high resolution NMR

Characterization of the structure and dynamics of nucleic acids by NMR benefits significantly from position specifically labeled nucleotides. Here an E. coli strain deficient in the transketolase gene (tktA) and grown on glucose that is labeled at different carbon sites is shown to facilitate cost-effective and large scale production of useful nucleotides. These nucleotides are site specifically labeled in C1' and C5' with minimal scrambling within the ribose ring. To demonstrate the utility of this labeling approach, the new site-specific labeled and the uniformly labeled nucleotides were used to synthesize a 36-nt RNA containing the catalytically essential domain 5 (D5) of the brown algae group II intron self-splicing ribozyme. The D5 RNA was used in binding and relaxation studies probed by NMR spectroscopy. Key nucleotides in the D5 RNA that are implicated in binding Mg(2+) ions are well resolved. As a result, spectra obtained using selectively labeled nucleotides have higher signal-to-noise ratio compared to those obtained using uniformly labeled nucleotides. Thus, compared to the uniformly (13)C/(15)N-labeled nucleotides, these specifically labeled nucleotides eliminate the extensive (13)C-(13)C coupling within the nitrogenous base and ribose ring, give rise to less crowded and more resolved NMR spectra, and accurate relaxation rates without the need for constant-time or band-selective decoupled NMR experiments. These position selective labeled nucleotides should, therefore, find wide use in NMR analysis of biologically interesting RNA molecules.

Asymmetry of 13C labeled 3-pyruvate affords improved site specific labeling of RNA for NMR spectroscopy

Selective isotopic labeling provides an unparalleled window within which to study the structure and dynamics of RNAs by high resolution NMR spectroscopy. Unlike commonly used carbon sources, the asymmetry of (13)C-labeled pyruvate provides selective labeling in both the ribose and base moieties of nucleotides using Escherichia coli variants, that until now were not feasible. Here we show that an E. coli mutant strain that lacks succinate and malate dehydrogenases (DL323) and grown on [3-(13)C]-pyruvate affords ribonucleotides with site specific labeling at C5' (~95%) and C1' (~42%) and minimal enrichment elsewhere in the ribose ring. Enrichment is also achieved at purine C2 and C8 (~95%) and pyrimidine C5 (~100%) positions with minimal labeling at pyrimidine C6 and purine C5 positions. These labeling patterns contrast with those obtained with DL323 E. coli grown on [1, 3-(13)C]-glycerol for which the ribose ring is labeled in all but the C4' carbon position, leading to multiplet splitting of the C1', C2' and C3' carbon atoms. The usefulness of these labeling patterns is demonstrated with a 27-nt RNA fragment derived from the 30S ribosomal subunit. Removal of the strong magnetic coupling within the ribose and base leads to increased sensitivity, substantial simplification of NMR spectra, and more precise and accurate dynamic parameters derived from NMR relaxation measurements. Thus these new labels offer valuable probes for characterizing the structure and dynamics of RNA that were previously limited by the constraint of uniformly labeled nucleotides.

GUARDD: user-friendly MATLAB software for rigorous analysis of CPMG RD NMR data

Molecular dynamics are essential for life, and nuclear magnetic resonance (NMR) spectroscopy has been used extensively to characterize these phenomena since the 1950s. For the past 15 years, the Carr-Purcell Meiboom-Gill relaxation dispersion (CPMG RD) NMR experiment has afforded advanced NMR labs access to kinetic, thermodynamic, and structural details of protein and RNA dynamics in the crucial μs-ms time window. However, analysis of RD data is challenging because datasets are often large and require many non-linear fitting parameters, thereby confounding assessment of accuracy. Moreover, novice CPMG experimentalists face an additional barrier because current software options lack an intuitive user interface and extensive documentation. Hence, we present the open-source software package GUARDD (Graphical User-friendly Analysis of Relaxation Dispersion Data), which is designed to organize, automate, and enhance the analytical procedures which operate on CPMG RD data ( http://code.google.com/p/guardd/). This MATLAB-based program includes a graphical user interface, permits global fitting to multi-field, multi-temperature, multi-coherence data, and implements χ (2)-mapping procedures, via grid-search and Monte Carlo methods, to enhance and assess fitting accuracy. The presentation features allow users to seamlessly traverse the large amount of results, and the RD Simulator feature can help design future experiments as well as serve as a teaching tool for those unfamiliar with RD phenomena. Based on these innovative features, we expect that GUARDD will fill a well-defined gap in service of the RD NMR community.

Mapping the landscape of RNA dynamics with NMR spectroscopy

Among the three major classes of biomacromolecules (DNA, RNA, and proteins) RNA's pronounced dynamics are the most explicitly linked to its wide variety of functions, which include catalysis and the regulation of transcription, translation, and splicing. These functions are mediated by a range of RNA biomachinery, including such varied examples as macromolecular noncoding RNAs, microRNAs, small interfering RNAs, riboswitch RNAs, and RNA thermometers. In each case, the functional dynamics of an interconversion is characterized by an associated rate constant. In this Account, we provide an introduction to NMR spectroscopic characterization of the landscape of RNA dynamics. We introduce strategies for measuring NMR parameters at various time scales as well as the underlying models for describing the corresponding rate constants. RNA exhibits significant dynamic motion, which can be modulated by (i) intermolecular interactions, including specific and nonspecific binding of ions (such as Mg(2+) and tertiary amines), (ii) metabolites in riboswitches or RNA aptamers, and (iii) macromolecular interactions within ribonucleic protein particles, including the ribosome and the spliceosome. Our understanding of the nature of these dynamic changes in RNA targets is now being incorporated into RNA-specific approaches in the design of RNA inhibitors. Interactions of RNA with proteins, other RNAs, or small molecules often occur through binding mechanisms that follow an induced fit mechanism or a conformational selection mechanism, in which one of several populated RNA conformations is selected through ligand binding. The extent of functional dynamics, including the kinetic formation of a specific RNA tertiary fold, is dependent on the messenger RNA (mRNA) chain length. Thus, during de novo synthesis of mRNA, both in prokaryotes and eukaryotes, nascent mRNA of various lengths will adopt different secondary and tertiary structures. The speed of transcription has a critical influence on the functional dynamics of the RNA being synthesized. In addition to modulating the local dynamics of a conformational RNA ensemble, a given RNA sequence may adopt more than one global, three-dimensional structure. RNA modification is one way to select among these alternative structures, which are often characterized by nearly equal stability, but with high energy barriers for conformational interconversion. The refolding of different secondary and tertiary structures has been found to be a major regulatory mechanism for transcription and translation. These conformational transitions can be characterized with NMR spectroscopy, for any given RNA sequence, in response to external stimuli.

Asymmetry of (13)C labeled 3-pyruvate affords improved site specific labeling of RNA for NMR spectroscopy

Selective isotopic labeling provides an unparalleled window within which to study the structure and dynamics of RNAs by high resolution NMR spectroscopy. Unlike commonly used carbon sources, the asymmetry of (13)C-labeled pyruvate provides selective labeling in both the ribose and base moieties of nucleotides using E. coli variants, that until now were not feasible. Here we show that an E. coli mutant strain that lacks succinate and malate dehydrogenases (DL323) and grown on [3-(13)C]-pyruvate affords ribonucleotides with site specific labeling at C5' (~95%) and C1' (~42%) and minimal enrichment elsewhere in the ribose ring. Enrichment is also achieved at purine C2 and C8 (~95%) and pyrimidine C5 (~100%) positions with minimal labeling at pyrimidine C6 and purine C5 positions. These labeling patterns contrast with those obtained with DL323 E. coli grown on [1, 3-(13)C]-glycerol for which the ribose ring is labeled in all but the C4' carbon position, leading to multiplet splitting of the C1', C2' and C3' carbon atoms. The usefulness of these labeling patterns is demonstrated with a 27-nt RNA fragment derived from the 30S ribosomal subunit. Removal of the strong magnetic coupling within the ribose and base leads to increased sensitivity, substantial simplification of NMR spectra, and more precise and accurate dynamic parameters derived from NMR relaxation measurements. Thus these new labels offer valuable probes for characterizing the structure and dynamics of RNA that were previously limited by the constraint of uniformly labeled nucleotides.

Characterizing RNA dynamics at atomic resolution using solution-state NMR spectroscopy

Many recently discovered noncoding RNAs do not fold into a single native conformation but sample many different conformations along their free-energy landscape to carry out their biological function. Here we review solution-state NMR techniques that measure the structural, kinetic and thermodynamic characteristics of RNA motions spanning picosecond to second timescales at atomic resolution, allowing unprecedented insights into the RNA dynamic structure landscape. From these studies a basic description of the RNA dynamic structure landscape is emerging, bringing new insights into how RNA structures change to carry out their function as well as applications in RNA-targeted drug discovery and RNA bioengineering.

Examining the relationship between RNA function and motion using nuclear magnetic resonance

The biological function of proteins and nucleic acids relies on their complex structures, yet dynamics provides an additional layer of functional adaptability. Numerous studies have demonstrated that RNA is only able to perform the multitude of functions for which it is responsible by readily changing its conformation in response to binding of proteins or small molecules. Examination of RNA dynamics is therefore essential to understanding its biological function. Nuclear magnetic resonance (NMR) has emerged as a leading technique for the examination of RNA motion and conformational transitions. It can examine domain motions as well as motion with atomic level resolution over a wide range of time scales. This review examines how NMR spectroscopy can be applied to examine the relationship between function and dynamics in RNA.

Transient Hoogsteen base pairs in canonical duplex DNA

Sequence-directed variations in the canonical DNA double helix structure that retain Watson-Crick base-pairing have important roles in DNA recognition, topology and nucleosome positioning. By using nuclear magnetic resonance relaxation dispersion spectroscopy in concert with steered molecular dynamics simulations, we have observed transient sequence-specific excursions away from Watson-Crick base-pairing at CA and TA steps inside canonical duplex DNA towards low-populated and short-lived A•T and G•C Hoogsteen base pairs. The observation of Hoogsteen base pairs in DNA duplexes specifically bound to transcription factors and in damaged DNA sites implies that the DNA double helix intrinsically codes for excited state Hoogsteen base pairs as a means of expanding its structural complexity beyond that which can be achieved based on Watson-Crick base-pairing. The methods presented here provide a new route for characterizing transient low-populated nucleic acid structures, which we predict will be abundant in the genome and constitute a second transient layer of the genetic code.