13C NMR relaxation studies of RNA base and ribose nuclei reveal a complex pattern of motions in the RNA binding site for human U1A protein

N-terminal domain of polypyrimidine-tract binding protein is a dynamic folding platform for adaptive RNA recognition

The N-terminal RNA recognition motif domain (RRM1) of polypyrimidine tract binding protein (PTB) forms an additional C-terminal helix α3, which docks to one edge of the β-sheet upon binding to a stem-loop RNA containing a UCUUU pentaloop. Importantly, α3 does not contact the RNA. The α3 helix therefore represents an allosteric means to regulate the conformation of adjacent domains in PTB upon binding structured RNAs. Here we investigate the process of dynamic adaptation by stem-loop RNA and RRM1 using NMR and MD in order to obtain mechanistic insights on how this allostery is achieved. Relaxation data and NMR structure determination of the free protein show that α3 is partially ordered and interacts with the domain transiently. Stem-loop RNA binding quenches fast time scale dynamics and α3 becomes ordered, however microsecond dynamics at the protein-RNA interface is observed. MD shows how RRM1 binding to the stem-loop RNA is coupled to the stabilization of the C-terminal helix and helps to transduce differences in RNA loop sequence into changes in α3 length and order. IRES assays of full length PTB and a mutant with altered dynamics in the α3 region show that this dynamic allostery influences PTB function in cultured HEK293T cells.

DNA base order parameter determination without influence of chemical exchange

Nuclear magnetic resonance (NMR) spectroscopy is a versatile tool used to investigate the dynamic properties of biological macromolecules and their complexes. NMR relaxation data can provide order parameters S2, which represent the mobility of bond vectors reorienting within a molecular frame. Determination of S2 parameters typically involves the use of transverse NMR relaxation rates. However, the accuracy in S2 determination can be diminished by elevation of the transverse relaxation rates through conformational or chemical exchange involving protonation/deprotonation or non-Watson-Crick base-pair states of nucleic acids. Here, we propose an approach for determination of S2 parameters without the influence of exchange processes. This approach utilizes transverse and longitudinal 13C chemical shift anisotropy (CSA) - dipole-dipole (DD) cross-correlation rates instead of 13C transverse relaxation rates. Anisotropy in rotational diffusion is taken into consideration. An application of this approach to nucleotide base CH groups of a uniformly 13C/15N-labeled DNA duplex is demonstrated.

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.

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.

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.

Access to Cα backbone dynamics of biological solids by 13C T1 relaxation and molecular dynamics simulation

We introduce a labeling scheme for magic angle spinning (MAS) solid-state NMR that is based on deuteration in combination with dilution of the carbon spin system. The labeling strategy achieves spectral editing by simplification of the HαCα and aliphatic side chain spectral region. A reduction in both proton and carbon spin density in combination with fast spinning (≥50 kHz) is essential to retrieve artifact-free (13)C-R1 relaxation data for aliphatic carbons. We obtain good agreement between the NMR experimental data and order parameters extracted from a molecular dynamics (MD) trajectory, which indicates that carbon based relaxation parameters can yield complementary information on protein backbone as well as side chain dynamics.

Dynamics of linker residues modulate the nucleic acid binding properties of the HIV-1 nucleocapsid protein zinc fingers

The HIV-1 nucleocapsid protein (NC) is a small basic protein containing two zinc fingers (ZF) separated by a short linker. It is involved in several steps of the replication cycle and acts as a nucleic acid chaperone protein in facilitating nucleic acid strand transfers occurring during reverse transcription. Recent analysis of three-dimensional structures of NC-nucleic acids complexes established a new property: the unpaired guanines targeted by NC are more often inserted in the C-terminal zinc finger (ZF2) than in the N-terminal zinc finger (ZF1). Although previous NMR dynamic studies were performed with NC, the dynamic behavior of the linker residues connecting the two ZF domains remains unclear. This prompted us to investigate the dynamic behavior of the linker residues. Here, we collected ¹⁵N NMR relaxation data and used for the first time data at several fields to probe the protein dynamics. The analysis at two fields allows us to detect a slow motion occurring between the two domains around a hinge located in the linker at the G35 position. However, the amplitude of motion appears limited in our conditions. In addition, we showed that the neighboring linker residues R29, A30, P31, R32, K33 displayed restricted motion and numerous contacts with residues of ZF1. Our results are fully consistent with a model in which the ZF1-linker contacts prevent the ZF1 domain to interact with unpaired guanines, whereas the ZF2 domain is more accessible and competent to interact with unpaired guanines. In contrast, ZF1 with its large hydrophobic plateau is able to destabilize the double-stranded regions adjacent to the guanines bound by ZF2. The linker residues and the internal dynamics of NC regulate therefore the different functions of the two zinc fingers that are required for an optimal chaperone activity.

Single cell spectroscopy: noninvasive measures of small-scale structure and function

The advancement of spectroscopy methods attained through increases in sensitivity, and often with the coupling of complementary techniques, has enabled real-time structure and function measurements of single cells. The purpose of this review is to illustrate, in light of advances, the strengths and the weaknesses of these methods. Included also is an assessment of the impact of the experimental setup and conditions of each method on cellular function and integrity. A particular emphasis is placed on noninvasive and nondestructive techniques for achieving single cell detection, including nuclear magnetic resonance, in addition to physical, optical, and vibrational methods.

Intrinsic flexibility of snRNA hairpin loops facilitates protein binding

Stem-loop II of U1 snRNA and Stem-loop IV of U2 snRNA typically have 10 or 11 nucleotides in their loops. The fluorescent nucleobase 2-aminopurine was used as a substitute for the adenines in each loop to probe the local and global structures and dynamics of these unusually long loops. Using steady-state and time-resolved fluorescence, we find that, while the bases in the loops are stacked, they are able to undergo significant local motion on the picosecond/nanosecond timescale. In addition, the loops have a global conformational change at low temperatures that occurs on the microsecond timescale, as determined using laser T-jump experiments. Nucleobase and loop motions are present at temperatures far below the melting temperature of the hairpin stem, which may facilitate the conformational change required for specific protein binding to these RNA loops.

Probing sequence-specific DNA flexibility in a-tracts and pyrimidine-purine steps by nuclear magnetic resonance (13)C relaxation and molecular dynamics simulations

Sequence-specific DNA flexibility plays a key role in a variety of cellular interactions that are critical for gene packaging, expression, and regulation, yet few studies have experimentally explored the sequence dependence of DNA dynamics that occur on biologically relevant time scales. Here, we use nuclear magnetic resonance (NMR) carbon spin relaxation combined with molecular dynamics (MD) simulations to examine the picosecond to nanosecond dynamics in a variety of dinucleotide steps as well as in varying length homopolymeric A(n)·T(n) repeats (A(n)-tracts, where n = 2, 4, or 6) that exhibit unusual structural and mechanical properties. We extend the NMR spin relaxation time scale sensitivity deeper into the nanosecond regime by using glycerol and a longer DNA duplex to slow overall tumbling. Our studies reveal a structurally unique A-tract core (for n > 3) that is uniformly rigid, flanked by junction steps that show increasing sugar flexibility with A-tract length. High sugar mobility is observed at pyrimidine residues at the A-tract junctions, which is encoded at the dinucleotide level (CA, TG, and CG steps) and increases with A-tract length. The MD simulations reproduce many of these trends, particularly the overall rigidity of A-tract base and sugar sites, and suggest that the sugar-backbone dynamics could involve transitions in sugar pucker and phosphate backbone BI ↔ BII equilibria. Our results reinforce an emerging view that sequence-specific DNA flexibility can be imprinted in dynamics occurring deep within the nanosecond time regime that is difficult to characterize experimentally at the atomic level. Such large-amplitude sequence-dependent backbone fluctuations might flag the genome for specific DNA recognition.

Intrinsic nucleic acid dynamics modulates HIV-1 nucleocapsid protein binding to its targets

HIV-1 nucleocapsid protein (NC) is involved in the rearrangement of nucleic acids occurring in key steps of reverse transcription. The protein, through its two zinc fingers, interacts preferentially with unpaired guanines in single-stranded sequences. In mini-cTAR stem-loop, which corresponds to the top half of the cDNA copy of the transactivation response element of the HIV-1 genome, NC was found to exhibit a clear preference for the TGG sequence at the bottom of mini-cTAR stem. To further understand how this site was selected among several potential binding sites containing unpaired guanines, we probed the intrinsic dynamics of mini-cTAR using (13)C relaxation measurements. Results of spin relaxation time measurements have been analyzed using the model-free formalism and completed by dispersion relaxation measurements. Our data indicate that the preferentially recognized guanine in the lower part of the stem is exempt of conformational exchange and highly mobile. In contrast, the unrecognized unpaired guanines of mini-cTAR are involved in conformational exchange, probably related to transient base-pairs. These findings support the notion that NC preferentially recognizes unpaired guanines exhibiting a high degree of mobility. The ability of NC to discriminate between close sequences through their dynamic properties contributes to understanding how NC recognizes specific sites within the HIV genome.

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.

Exploring RNA structural codes with SHAPE chemistry

RNA is the central conduit for gene expression. This role depends on an ability to encode information at two levels: in its linear sequence and in the complex structures RNA can form by folding back on itself. Understanding the global structure-function interrelationships mediated by RNA remains a great challenge in molecular and structural biology. In this Account, we discuss evolving work in our laboratory focused on creating facile, generic, quantitative, accurate, and highly informative approaches for understanding RNA structure in biologically important environments. The core innovation derives from our discovery that the nucleophilic reactivity of the ribose 2'-hydroxyl in RNA is gated by local nucleotide flexibility. The 2'-hydroxyl is reactive at conformationally flexible positions but is unreactive at nucleotides constrained by base pairing. Sites of modification in RNA can be detected efficiently either using primer extension or by protection from exoribonucleolytic degradation. This technology is now called SHAPE, for selective 2'-hydroxyl acylation analyzed by primer extension (or protection from exoribonuclease). SHAPE reactivities are largely independent of nucleotide identity but correlate closely with model-free measurements of molecular order. The simple SHAPE reaction is thus a robust, nucleotide-resolution, biophysical measurement of RNA structure. SHAPE can be used to provide an experimental correction to RNA folding algorithms and, in favorable cases, yield kilobase-scale secondary structure predictions with high accuracies. SHAPE chemistry is based on very simple reactive carbonyl centers that can be varied to yield slow- and fast-reacting reagents. Differential SHAPE reactivities can be used to detect specific RNA positions with slow local nucleotide dynamics. These positions, which are often in the C2'-endo conformation, have the potential to function as molecular timers that regulate RNA folding and function. In addition, fast-reacting SHAPE reagents can be used to visualize RNA structural biogenesis and RNA-protein assembly reactions in one second snapshots in very straightforward experiments. The application of SHAPE to challenging problems in biology has revealed surprises in well-studied systems. New regions have been identified that are likely to have critical functional roles on the basis of their high levels of RNA structure. For example, SHAPE analysis of large RNAs, such as authentic viral RNA genomes, suggests that RNA structure organizes regulatory motifs and regulates splicing, protein folding, genome recombination, and ribonucleoprotein assembly. SHAPE has also revealed limitations to the hierarchical model for RNA folding. Continued development and application of SHAPE technologies will advance our understanding of the many ways in which the genetic code is expressed through the underlying structure of RNA.

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.

Probing the folding intermediate of Bacillus subtilis RNase P protein by nuclear magnetic resonance

Protein folding intermediates are often imperative for overall folding processes and consequent biological functions. However, the low population and transient nature of the intermediate states often hinder their biochemical and biophysical characterization. Previous studies have demonstrated that Bacillus subtilis ribonuclease P protein (P protein) is conformationally heterogeneous and folds with multiphasic kinetics, indicating the presence of an equilibrium and kinetic intermediate in its folding mechanism. In this study, nuclear magnetic resonance (NMR) spectroscopy was used to study the ensemble corresponding to this intermediate (I). The results indicate that the N-terminal and C-terminal helical regions are mostly unfolded in I. 1H−15N heteronuclear single-quantum coherence NMR spectra collected as a function of pH suggest that the protonation of His 22 may play a major role in the energetics of the equilibria among the unfolded, intermediate, and folded state ensembles of P protein. NMR paramagnetic relaxation enhancement experiments were also used to locate the small anion binding sites in both the intermediate and folded ensembles. The results for the folded protein are consistent with the previously modeled binding regions. These structural insights suggest a possible role for I in the RNase P holoenzyme assembly process.

Multistep kinetics of the U1A-SL2 RNA complex dissociation

The U1A-SL2 RNA complex is a model system for studying interactions between RNA and the RNA recognition motif (RRM), which is one of the most common RNA binding domains. We report here kinetic studies of dissociation of the U1A-SL2 RNA complex, using laser temperature jump and stopped-flow fluorescence methods with U1A proteins labeled with the intrinsic chromophore tryptophan. An analysis of the kinetic data suggests three phases of dissociation with time scales of ∼100 μs, ∼50 ms, and ∼2 s. We propose that the first step of dissociation is a fast rearrangement of the complex to form a loosely bound complex. The intermediate step is assigned to be the dissociation of the U1A-SL2 RNA complex, and the final step is assigned to a reorganization of the U1A protein structure into the conformation of the free protein. These assignments are consistent with previous proposals based on thermodynamic, NMR, and surface plasmon resonance experiments and molecular dynamics simulations. Together, these results begin to build a comprehensive model of the complex dynamic processes involved in the formation and dissociation of an RRM-RNA complex.

Slow exchange model of nonrigid rotational motion in RNA for combined solid-state and solution NMR studies

Functional RNA molecules are conformationally dynamic and sample a multitude of dynamic modes over a wide range of frequencies. Thus, a comprehensive description of RNA dynamics requires the inclusion of a broad range of motions across multiple dynamic rates which must be derived from multiple spectroscopies. Here we describe a slow conformational exchange theoretical approach to combining the description of local motions in RNA that occur in the nanosecond to microsecond window and are detected by solid-state NMR with nonrigid rotational motion of the HIV-1 transactivation response element (TAR) RNA in solution as observed by solution NMR. This theoretical model unifies the experimental results generated by solution and solid-state NMR and provides a comprehensive view of the dynamics of HIV-1 TAR RNA, a well-known paradigm of an RNA where function requires extensive conformational rearrangements. This methodology provides a quantitative atomic level view of the amplitudes and rates of the local and collective displacements of the TAR RNA molecule and provides directly motional parameters for the conformational capture hypothesis of this classical RNA-ligand interaction.

The unstable part of the apical stem of duck hepatitis B virus epsilon shows enhanced base pair opening but not pico- to nanosecond dynamics and is essential for reverse transcriptase binding

Hepatitis B virus (HBV) replication starts with binding of reverse transcriptase (RT) to the apical stem-loop region of epsilon, a conserved element of the RNA pregenome. For duck HBV, an in vitro replication system has provided molecular details of this interaction. Further insights can be obtained from the structure and dynamics of the duck and human apical stem-loops. Previously, we reported these for the human apical stem-loop. Here, we present the same for the duck counterpart. Unlike its human counterpart, the duck apical stem is unstable in its middle/upper part and contains noncanonical base pairs. This dynamics study is the first of an unstable RNA-DNA stem. Similar to the human stem, the duck apical stem comprises two helical segments with a bend angle of ca. 10 degrees , separated by a nonpaired mobile U residue. It is capped by a well-structured conserved UGUU loop with two residues mobile on the pico- to nanosecond time scale, one of which is involved in RT binding. Remarkably, the unstable middle/upper part of the stem does not show enhanced pico- to nanosecond time scale dynamics. Instead, adenine dispersion relaxation studies indicate enhanced millisecond time scale dynamics involving base pair opening. It can then be concluded that base pair opening is essential for epsilon-RT binding, because stabilization of the stem abolishes binding. We hypothesize that binding occurs by conformational capture of bases in the base pair open state. The unstable secondary structure of the apical stem-loop makes duck epsilon-RT binding unusual in light of recent classifications of RNA target interactions that assume stable secondary structures.

RNA phosphodiester backbone dynamics of a perdeuterated cUUCGg tetraloop RNA from phosphorus-31 NMR relaxation analysis

We have analyzed the relaxation properties of all (31)P nuclei in an RNA cUUCGg tetraloop model hairpin at proton magnetic field strengths of 300, 600 and 900 MHz in solution. Significant H, P dipolar contributions to R (1) and R (2) relaxation are observed in a protonated RNA sample at 600 MHz. These contributions can be suppressed using a perdeuterated RNA sample. In order to interpret the (31)P relaxation data (R (1), R (2)), we measured the (31)P chemical shift anisotropy (CSA) by solid-state NMR spectroscopy under various salt and hydration conditions. A value of 178.5 ppm for the (31)P CSA in the static state (S (2) = 1) could be determined. In order to obtain information about fast time scale dynamics we performed a modelfree analysis on the basis of our relaxation data. The results show that subnanosecond dynamics detected around the phosphodiester backbone are more pronounced than the dynamics detected for the ribofuranosyl and nucleobase moieties of the individual nucleotides (Duchardt and Schwalbe, J Biomol NMR 32:295-308, 2005; Ferner et al., Nucleic Acids Res 36:1928-1940, 2008). Furthermore, the dynamics of the individual phosphate groups seem to be correlated to the 5' neighbouring nucleobases.

Hydration dependent dynamics in RNA

The essential role played by local and collective motions in RNA function has led to a growing interest in the characterization of RNA dynamics. Recent investigations have revealed that even relatively simple RNAs experience complex motions over multiple time scales covering the entire ms-ps motional range. In this work, we use deuterium solid-state NMR to systematically investigate motions in HIV-1 TAR RNA as a function of hydration. We probe dynamics at three uridine residues in different structural environments ranging from helical to completely unrestrained. We observe distinct and substantial changes in (2)H solid-state relaxation times and lineshapes at each site as hydration levels increase. By comparing solid-state and solution state (13)C relaxation measurements, we establish that ns-micros motions that may be indicative of collective dynamics suddenly arise in the RNA as hydration reaches a critical point coincident with the onset of bulk hydration. Beyond that point, we observe smaller changes in relaxation rates and lineshapes in these highly hydrated solid samples, compared to the dramatic activation of motion occurring at moderate hydration.

Preparation, resonance assignment, and preliminary dynamics characterization of residue specific 13C/15N-labeled elongated DNA for the study of sequence-directed dynamics by NMR

DNA is a highly flexible molecule that undergoes functionally important structural transitions in response to external cellular stimuli. Atomic level spin relaxation NMR studies of DNA dynamics have been limited to short duplexes in which sensitivity to biologically relevant fluctuations occurring at nanosecond timescales is often inadequate. Here, we introduce a method for preparing residue-specific (13)C/(15)N-labeled elongated DNA along with a strategy for establishing resonance assignments and apply the approach to probe fast inter-helical bending motions induced by an adenine tract. Preliminary results suggest the presence of elevated A-tract independent end-fraying internal motions occurring at nanosecond timescales, which evade detection in short DNA constructs and that penetrate deep (7 bp) within the DNA helix and gradually fade away towards the helix interior.

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

Conformational dynamics play a key role in the properties and functions of proteins and nucleic acids. Heteronuclear NMR spin relaxation is a uniquely powerful site-specific probe of dynamics in proteins and has found increasing applications to nucleotide base side chains and anomeric sites in RNA. Applications to the nucleic acid ribose backbone, however, have been hampered by strong magnetic coupling among ring carbons in uniformly 13C-labeled samples. In this work, we apply a recently developed, metabolically directed isotope labeling scheme that places 13C with high efficiency and specificity at the nucleotide ribose C2' and C4' sites. We take advantage of this scheme to explore backbone dynamics in the well-studied GCAA RNA tetraloop. Using a combination of CPMG (Carr-Purcell-Meiboom-Gill) and R(1rho) relaxation dispersion spectroscopy to explore exchange processes on the microsecond to millisecond time scale, we find an extensive pattern of dynamic transitions connecting a set of relatively well-defined conformations. In many cases, the observed transitions appear to be linked to C3'-endo/C2'-endo sugar pucker transitions of the corresponding nucleotides, and may also be correlated across multiple nucleotides within the tetraloop. These results demonstrate the power of NMR spin relaxation based on alternate-site isotope labeling to open a new window into the dynamic properties of ribose backbone groups in RNA.

Stabilization of RNA oligomers through reverse micelle encapsulation

The cellular milieu is rich in diversity of both simple and complex molecules and is also quite crowded. By contrast, typical sample concentrations employed for in vitro investigation of biophysics and structural biology make use of purified macromolecules in simple buffer systems at concentrations that range from micromolar to millimolar. Although this formulation has proven to be compatible with a wide range of biological and structural studies, it is quite different from the relatively crowded conditions typically found within cells. The importance of these crowding effects for proteins has been recognized for some time, but the equivalent analysis is underexplored in nucleic acids. Encapsulation with surfactant-based reverse micelles has emerged as an effective biophysical tool, allowing study of the influence of ionic strength, pH, hydration, and crowding on biologically active macromolecules over a wide range of conditions. We have encapsulated an oligonucleotide model of TAR RNA from HIV and the 5' stem loop oligonucleotide of the U4 snRNA. Observation of imino (1)H resonances is an established method for evaluating the stability of nucleic acid oligonucleotides, implying the presence of stacked, hydrogen bonded base pairs. Inspection of (1)H NMR spectra of the RNA molecules reveals that the intensity of several of the imino resonances increases upon encapsulation. Additional resonances not observed in spectra of the oligonucleotides free in solution support the suggestion that the molecules have gained stability. These results indicate that RNA oligonucleotides may acquire significant stability in the presence of cellular levels of crowding.

How binding of small molecule and peptide ligands to HIV-1 TAR alters the RNA motional landscape

The HIV-1 TAR RNA represents a well-known paradigm to study the role of dynamics and conformational change in RNA function. This regulatory RNA changes conformation in response to binding of Tat protein and of a variety of peptidic and small molecule ligands, indicating that its conformational flexibility and intrinsic dynamics play important roles in molecular recognition. We have used ¹³C NMR relaxation experiments to examine changes in the motional landscape of HIV-1 TAR in the presence of three ligands of different affinity and specificity. The ligands are argininamide, a linear peptide mimic of the Tat basic domain and a cyclic peptide that potently inhibits Tat-dependent activation of transcription. All three molecules induce the same motional characteristics within the three nucleotides bulge that represents the Tat-binding site. However, the cyclic peptide has a unique motional signature in the apical loop, which represents a binding site for the essential host co-factor cyclin T1. These results suggest that all peptidic mimics of Tat induce the same dynamics in TAR within this protein binding site. However, the new cyclic peptide mimic of Tat represents a new class of ligands with a unique effect on the dynamics and the structure of the apical loop.

Strong correlation between SHAPE chemistry and the generalized NMR order parameter (S2) in RNA

The functions of most RNA molecules are critically dependent on the distinct local dynamics that characterize secondary structure and tertiary interactions and on structural changes that occur upon binding by proteins and small molecule ligands. Measurements of RNA dynamics at nucleotide resolution set the foundation for understanding the roles of individual residues in folding, catalysis, and ligand recognition. In favorable cases, local order in small RNAs can be quantitatively analyzed by NMR in terms of a generalized order parameter, S2. Alternatively, SHAPE (selective 2'-hydroxyl acylation analyzed by primer extension) chemistry measures local nucleotide flexibility in RNAs of any size using structure-sensitive reagents that acylate the 2'-hydroxyl position. In this work, we compare per-residue RNA dynamics, analyzed by both S2 and SHAPE, for three RNAs: the HIV-1 TAR element, the U1A protein binding site, and the Tetrahymena telomerase stem loop 4. We find a very strong correlation between the two measurements: nucleotides with high SHAPE reactivities consistently have low S2 values. We conclude that SHAPE chemistry quantitatively reports local nucleotide dynamics and can be used with confidence to analyze dynamics in large RNAs, RNA-protein complexes, and RNAs in vivo.

Conformational dynamics of an intact virus: order parameters for the coat protein of Pf1 bacteriophage

This study has examined the atomic-level dynamics of the protein in the capsid of filamentous phage Pf1. This capsid consists of approximately 7,300 small subunits of only 46 aa in a helical array around a highly extended, circular single-stranded DNA molecule of 7,349 nt. Measurements were made of site-specific, solid-state NMR order parameters, S, the values which are dimensionless quantities between 0 (mobile) and 1 (static) that characterize the amplitudes of molecular bond angular motions that are faster than microseconds. It was found that the protein subunit backbone is very static, and of particular interest, it appears to be static at residues glycine 15 and glutamine 16 where it had been previously thought to be mobile. In contrast to the backbone, several side chains display large-amplitude angular motions. Side chains on the virion exterior that interact with solvent are highly mobile, but surprisingly, the side chains of residues arginine 44 and lysine 45 near the DNA deep in the interior of the virion are also highly mobile. The large-amplitude dynamic motion of these positively charged side chains in their interactions with the DNA were not previously expected. The results reveal a highly dynamic aspect of a DNA-protein interface within a virus.

13C relaxation studies of the DNA target sequence for hhai methyltransferase reveal unique motional properties

The goal of this work was to examine if sequence-dependent conformational flexibility in DNA plays a role in base extrusion, a common conformational change induced by many DNA-modifying enzymes. We studied the dynamics of the double-stranded DNA target of the HhaI methyltransferase by recording an extensive set of (13)C NMR relaxation parameters. We observe that the cytidine furanose rings experience fast (picosecond to nanosecond) motions that are not present in other nucleotides; the methylation site experiences particularly high mobility. We also observe that the bases of guanosine and cytidine residues within the HhaI recognition sequence GCGC experience motions on a much slower (1-100 micros) time scale. We compare these observations with previous solution and solid-state NMR studies of the EcoRI nuclease target sequence, and solid-state NMR studies of a similar HhaI target construct. While an increased mobility of cytidine furanose rings compared to those of other nucleotides is observed for both sequences, the slower motions are only observed in the HhaI target DNA. We propose that this inherent flexibility lowers the energetic barriers that must occur when the DNA binds to the HhaI methyltransferase and for extrusion of the cytidine prior to its methylation.

Cytosine ribose flexibility in DNA: a combined NMR 13C spin relaxation and molecular dynamics simulation study

Using (13)C spin relaxation NMR in combination with molecular dynamic (MD) simulations, we characterized internal motions within double-stranded DNA on the pico- to nano-second time scale. We found that the C-H vectors in all cytosine ribose moieties within the Dickerson-Drew dodecamer (5'-CGCGAATTCGCG-3') are subject to high amplitude motions, while the other nucleotides are essentially rigid. MD simulations showed that repuckering is a likely motional model for the cytosine ribose moiety. Repuckering occurs with a time constant of around 100 ps. Knowledge of DNA dynamics will contribute to our understanding of the recognition specificity of DNA-binding proteins such as cytosine methyltransferase.

Analyzing the flexibility of RNA structures by constraint counting

RNA requires conformational dynamics to undergo its diverse functional roles. Here, a new topological network representation of RNA structures is presented that allows analyzing RNA flexibility/rigidity based on constraint counting. The method extends the FIRST approach, which identifies flexible and rigid regions in atomic detail in a single, static, three-dimensional molecular framework. Initially, the network rigidity of a canonical A-form RNA is analyzed by counting on constraints of network elements of increasing size. These considerations demonstrate that it is the inclusion of hydrophobic contacts into the RNA topological network that is crucial for an accurate flexibility prediction. The counting also explains why a protein-based parameterization results in overly rigid RNA structures. The new network representation is then validated on a tRNA(ASP) structure and all NMR-derived ensembles of RNA structures currently available in the Protein Data Bank (with chain length >/=40). The flexibility predictions demonstrate good agreement with experimental mobility data, and the results are superior compared to predictions based on two previously used network representations. Encouragingly, this holds for flexibility predictions as well as mobility predictions obtained by constrained geometric simulations on these networks. Potential applications of the approach to analyzing the flexibility of DNA and RNA/protein complexes are discussed.

NMR and MD studies of the temperature-dependent dynamics of RNA YNMG-tetraloops

In a combined NMR/MD study, the temperature-dependent changes in the conformation of two members of the RNA YNMG-tetraloop motif (cUUCGg and uCACGg) have been investigated at temperatures of 298, 317 and 325 K. The two members have considerable different thermal stability and biological functions. In order to address these differences, the combined NMR/MD study was performed. The large temperature range represents a challenge for both, NMR relaxation analysis (consistent choice of effective bond length and CSA parameter) and all-atom MD simulation with explicit solvent (necessity to rescale the temperature). A convincing agreement of experiment and theory is found. Employing a principle component analysis of the MD trajectories, the conformational distribution of both hairpins at various temperatures is investigated. The ground state conformation and dynamics of the two tetraloops are indeed found to be very similar. Furthermore, both systems are initially destabilized by a loss of the stacking interactions between the first and the third nucleobase in the loop region. While the global fold is still preserved, this initiation of unfolding is already observed at 317 K for the uCACGg hairpin but at a significantly higher temperature for the cUUCGg hairpin.

Solution structure of stem-loop alpha of the hepatitis B virus post-transcriptional regulatory element

Chronic hepatitis B virus (HBV) infections may lead to severe diseases like liver cirrhosis or hepatocellular carcinoma (HCC). The HBV post-transcriptional regulatory element (HPRE) facilitates the nuclear export of unspliced viral mRNAs, contains a splicing regulatory element and resides in the 3′-region of all viral transcripts. The HPRE consists of three sub-elements α (nucleotides 1151–1346), β1 (nucleotides 1347–1457) and β2 (nucleotides 1458–1582), which confer together full export competence. Here, we present the NMR solution structure (pdb 2JYM) of the stem-loop α (SLα, nucleotides 1292–1321) located in the sub-element α. The SLα contains a CAGGC pentaloop highly conserved in hepatoviruses, which essentially adopts a CUNG-like tetraloop conformation. Furthermore, the SLα harbours a single bulged G residue flanked by A-helical regions. The structure is highly suggestive of serving two functions in the context of export of unspliced viral RNA: binding sterile alpha motif (SAM-) domain containing proteins and/or preventing the utilization of a 3′-splice site contained within SLα.

DNA mimicry by a high-affinity anti-NF-kappaB RNA aptamer

The binding of RNA molecules to proteins or other ligands can require extensive RNA folding to create an induced fit. Understanding the generality of this principle involves comparing structures of RNA before and after complex formation. Here we report the NMR solution structure of a 29-nt RNA aptamer whose crystal structure had previously been determined in complex with its transcription factor target, the p50₂ form of NF-κB. The RNA aptamer internal loop structure has pre-organized features that are also found in the complex, including non-canonical base pairing and cross-strand base stacking. Remarkably, the free RNA aptamer structure possesses a major groove that more closely resembles B-form DNA than RNA. Upon protein binding, changes in RNA structure include the kinking of the internal loop and distortion of the terminal tetraloop. Thus, complex formation involves both pre-formed and induced fit binding interactions. The high affinity of the NF-κB transcription factor for this RNA aptamer may largely be due to the structural pre-organization of the RNA that results in its ability to mimic DNA.

Dynamics of large elongated RNA by NMR carbon relaxation

We present an NMR strategy for characterizing picosecond-to-nanosecond internal motions in uniformly 13C/15N-labeled RNAs that combines measurements of R1, R1rho, and heteronuclear 13C{1H} NOEs for protonated base (C2, C5, C6, and C8) and sugar (C1') carbons with a domain elongation strategy for decoupling internal from overall motions and residual dipolar coupling (RDC) measurements for determining the average RNA global conformation and orientation of the principal axis of the axially symmetric rotational diffusion. TROSY-detected pulse sequences are presented for the accurate measurement of nucleobase carbon R1 and R1rho rates in large RNAs. The relaxation data is analyzed using a model free formalism which takes into account the very high anisotropy of overall rotational diffusion (Dratio approximately 4.7), asymmetry of the nucleobase CSAs and noncollinearity of C-C, C-H dipolar and CSA interactions under the assumption that all interaction tensors for a given carbon experience identical isotropic internal motions. The approach is demonstrated and validated on an elongated HIV-1 TAR RNA (taum approximately 18 ns) both in free form and bound to the ligand argininamide (ARG). Results show that, while ARG binding reduces the amplitude of collective helix motions and local mobility at the binding pocket, it leads to a drastic increase in the local mobility of "spacer" bulge residues linking the two helices which undergo virtually unrestricted internal motions (S2 approximately 0.2) in the ARG bound state. Our results establish the ability to quantitatively study the dynamics of RNAs which are significantly larger and more anisotropic than customarily studied by NMR carbon relaxation.

Conserved nucleotides in an RNA essential for hepatitis B virus replication show distinct mobility patterns

The number of regulatory RNAs with identified non-canonical structures is increasing, and structural transitions often play a role in their biological function. This stimulates interest in internal motions of RNA, which can underlie structural transitions. Heteronuclear NMR relaxation measurements, which are commonly used to study internal motion, only report on local motions of few sites within the molecule. Here we have studied a 27-nt segment of the human hepatitis B virus (HBV) pregenomic RNA, which is essential for viral replication. We combined heteronuclear relaxation with the new off-resonance ROESY technique, which reports on internal motions of H,H contacts. Using off-resonance ROESY, we could for the first time detect motion of through-space H,H contacts, such as in intra-residue base-ribose contacts or inter-nucleotide contacts, both essential for NMR structure determination. Motions in non-canonical structure elements were found primarily on the sub-nanosecond timescale. Different patterns of mobility were observed among several mobile nucleotides. The most mobile nucleotides are highly conserved among different HBV strains, suggesting that their mobility patterns may be necessary for the RNA’s biological function.

NMR studies of dynamics in RNA and DNA by 13C relaxation

RNA and DNA molecules experience motions on a wide range of time scales, ranging from rapid localized motions to much slower collective motions of entire helical domains. The many functions of RNA in biology very often require this molecule to change its conformation in response to biological signals in the form of small molecules, proteins or other nucleic acids, whereas local motions in DNA may facilitate protein recognition and allow enzymes acting on DNA to access functional groups on the bases that would otherwise be buried in Watson-Crick base pairs. Although these statements make a compelling case to study the sequence dependent dynamics in nucleic acids, there are few residue-specific studies of nucleic acid dynamics. Fortunately, NMR studies of dynamics of nucleic acids and nucleic acids-protein complexes are gaining increased attention. The aim of this review is to provide an update of the recent progress in studies of nucleic acid dynamics by NMR based on the application of solution relaxation techniques.

Binding of U1A protein changes RNA dynamics as observed by 13C NMR relaxation studies

Recognition of RNA by proteins and small molecules often involves large changes in RNA structure and dynamics, yet very few studies have so far characterized these motional changes. Here we extend to the protein-bound RNA recent 13C relaxation studies of motions in the RNA recognized by human U1A protein, a well-known model for protein-RNA recognition. Changes in relaxation observed upon complex formation demonstrate that the protein-binding site becomes rigid in the complex, but the upper stem-loop that defines the secondary structure of this RNA experiences unexpected motional freedom. By using a helix elongation strategy, we observe that the upper stem-loop moves independently of the remainder of the structure also in the absence of U1A. Surprisingly, RNA residues making important intermolecular contacts in the structure of the complex exhibit increased flexibility in the presence of the protein. Both of these results support the hypothesis that RNA-binding proteins select a structure that optimizes intermolecular contacts in the manifold of conformations sampled by the free RNA and that protein binding quenches these motions. Together with previous studies of the RNA-bound protein, they also demonstrate that protein-RNA interfaces experience complex motions that modulate the strength of individual interactions.

Structural plasticity and Mg2+ binding properties of RNase P P4 from combined analysis of NMR residual dipolar couplings and motionally decoupled spin relaxation

The P4 helix is an essential element of ribonuclease P (RNase P) that is believed to bind catalytically important metals. Here, we applied a combination of NMR residual dipolar couplings (RDCs) and a recently introduced domain-elongation strategy for measuring "motionally decoupled" relaxation data to characterize the structural dynamics of the P4 helix from Bacillus subtilis RNase P. In the absence of divalent ions, the two P4 helical domains undergo small amplitude (approximately 13 degrees) collective motions about an average interhelical angle of 10 degrees. The highly conserved U7 bulge and helical residue C8, which are proposed to be important for substrate recognition and metal binding, are locally mobile at pico- to nanosecond timescales and together form the pivot point for the collective domain motions. Chemical shift mapping reveals significant association of Mg2+ ions at the P4 major groove near the flexible pivot point at residues (A5, G22, G23) previously identified to bind catalytically important metals. The Mg2+ ions do not, however, significantly alter the structure or dynamics of P4. Analysis of results in the context of available X-ray structures of the RNA component of RNase P and structural models that include the pre-tRNA substrate suggest that the internal motions observed in P4 likely facilitate adaptive changes in conformation that take place during folding and substrate recognition, possibly aided by interactions with Mg2+ ions. Our results add to a growing view supporting the existence of functionally important internal motions in RNA occurring at nanosecond timescales.

Internal Bulge and Tetraloop of the Catalytic Domain 5 of a Group II Intron Ribozyme Are Flexible: Implications for Catalysis

RNA molecules have an inherent flexibility that enables recognition of other interacting partners through potential disorder-order transitions, yet studies to quantify such motional dynamics remain few. With an increasing database of three-dimensional structures of biologically important RNA molecules, quantifying such motions becomes important to link structural deformations with function. One such system studied intensely is domain 5 (D5) from the self-splicing group II introns, which is at the heart of its catalytic machinery. We report the dynamics of a 36 nucleotide D5 from the Pylaiella littoralis group II intron in the presence and absence of magnesium ions, and at a range of temperatures (298K-318 K). Using high-resolution NMR experiments of heteronuclear nuclear Overhauser enhancement (NOE), spin-lattice (R(1)), and spin-spin (R(2)) (13)C relaxation rates, we determined the rotational diffusion tensor of D5 using the ROTDIF program modified for RNA dynamic analysis (ROTDIF_RNA). The D5 rotational diffusion tensor has an axial symmetric ratio (D(||)/D(perpendicular)) of 1.7+/-0.3, consistent with an estimated overall rotational correlation time of tau(m)=(2D(||)+4D(perpendicular))(-1) of 6.1(+/-0.3) ns at 298 K and 4.1(+/-0.2) ns at 318 K. The measured relaxation data were analyzed with the reduced spectral density mapping formalism using assumed values of the chemical shift anisotropy of the (13)C spins. Both the relaxation data and the values of the spectral density function reveal that the functional groups in D5 implicated in magnesium ion binding and catalysis (catalytic triad, internal bulge, and tetraloop regions) exhibit thermally induced motion on a wide variety of timescales. Because these motions parallel those observed in the intramolecular stem-loop of the U6 element within the spliceosome, we hypothesize that such extensive dynamic disorder likely facilitates D5 engaging both binding and catalytic regions of the ribozyme, and these may be a conserved feature of the catalytic machinery essential for catalysis.

Impact of static and dynamic A-form heterogeneity on the determination of RNA global structural dynamics using NMR residual dipolar couplings

We examined how static and dynamic deviations from the idealized A-form helix propagate into errors in the principal order tensor parameters determined using residual dipolar couplings (rdcs). A 20-ns molecular dynamics (MD) simulation of the HIV-1 transactivation response element (TAR) RNA together with a survey of spin relaxation studies of RNA dynamics reveals that pico-to-nanosecond local motions in non-terminal Watson-Crick base-pairs will uniformly attenuate base and sugar one bond rdcs by approximately 7%. Gaussian distributions were generated for base and sugar torsion angles through statistical comparison of 40 RNA X-ray structures solved to <3.0 A resolution. For a typical number (>or=11) of one bond C-H base and sugar rdcs, these structural deviations together with rdc uncertainty (1.5 Hz) lead to average errors in the magnitude and orientation of the principal axis of order that are <9% and <4 degrees, respectively. The errors decrease to <5% and <4 degrees for >or=17 rdcs. A protocol that allows for estimation of error in A-form order tensors due to both angular deviations and rdc uncertainty (Aform-RDC) is validated using theoretical simulations and used to analyze rdcs measured previously in TAR in the free state and bound to four distinct ligands. Results confirm earlier findings that the two TAR helices undergo large changes in both their mean relative orientation and dynamics upon binding to different targets.

Alternate-site isotopic labeling of ribonucleotides for NMR studies of ribose conformational dynamics in RNA

Heteronuclear NMR spin relaxation studies of conformational dynamics are coming into increasing use to help understand the functions of ribozymes and other RNAs. Due to strong 13C-13C magnetic interactions within the ribose ring, however, these studies have thus far largely been limited to (13)C and (15)N resonances on the nucleotide base side chains. We report here the application of the alternate-site (13)C isotopic labeling scheme, pioneered by LeMaster for relaxation studies of amino acid side chains, to nucleic acid systems. We have used different strains of E. coli to prepare mononucleotides containing (13)C label in one of two patterns: Either C1' or C2' in addition to C4', termed (1'/2',4') labeling, or nearly complete labeling at the C2' and C4' sites only, termed (2',4') labeling. These patterns provide isolated 13C-1H spin systems on the labeled carbon atoms and thus allow spin relaxation studies without interference from 13C-13C scalar or dipolar coupling. Using relaxation studies of AMP dissolved in glycerol at varying temperature to produce systems with correlation times characteristic of different size RNAs, we demonstrate the removal of errors due to 13C-13C interaction in T (1) measurements of larger nucleic acids and in T (1rho) measurements in RNA molecules. By extending the applicability of spin relaxation measurements to backbone ribose groups, this technology should greatly improve the flexibility and completeness of NMR analyses of conformational dynamics in RNA.

Decoding RNA motional codes

When proteins and small molecules bind to RNA, they often alter its conformation. These structural changes are an essential aspect of the ability of RNA to sense signaling molecules and modulate gene expression. Thus far, few studies have been dedicated to understanding how RNA moves at a residue level and how these motions change upon complex formation. A recent report highlights how intrinsic motions in RNA correlate with its ability to bind to cognate ligands.

New RNA Labeling Methods Offer Dramatic Sensitivity Enhancements in 2H NMR Relaxation Spectra

A new labeling strategy is presented that greatly facilitates the measurement of 2H spin relaxation rates in RNA molecules as a probe of pico- to nanosecond time scale dynamics. In this labeling scheme the sugar positions are uniformly 13C-labeled, with position 2' protonated and all other sites on the sugar deuterated. Pulse sequences are presented for measurement of 2H R1 and R2 relaxation rates at positions 1', 3', and 4' with sensitivity gains that are on the order of 5-fold relative to previous methods that employed random fractional deuteration. The improved sensitivity is transformative and facilitates the study of motion in moderately sized RNA molecules with good sensitivity. The utility of the approach is demonstrated with an application to HIV-2 TAR, where the site-specific measures of molecular dynamics at sugar positions obtained here complement previous studies of dynamics at aromatic sites in the molecule.

Structure of stem-loop IV of Tetrahymena telomerase RNA

Conserved domains within the RNA component of telomerase provide the template for reverse transcription, recruit protein components to the holoenzyme and are required for enzymatic activity. Among the functionally essential domains in ciliate telomerase RNA is stem-loop IV, which strongly stimulates telomerase activity and processivity even when provided in trans. The NMR structure of Tetrahymena thermophila stem-loop IV shows a highly structured distal stem-loop linked to a conformationally flexible template-proximal region by a bulge that severely kinks the entire RNA. Through extensive structure-function studies, we identify residues that contribute to both these structural features and to enzymatic activity, with no apparent effect on the binding of TERT protein. We propose that the bending induced by the GA bulge and the flexibility of the template-proximal region allow positioning of the prestructured apical loop during the catalytic cycle.

Insight into the CSA tensors of nucleobase carbons in RNA polynucleotides from solution measurements of residual CSA: towards new long-range orientational constraints

Using residual chemical shift anisotropies (RCSAs) measured in a weakly aligned stem-loop RNA, we examined the carbon chemical shift anisotropy (CSA) tensors of nucleobase adenine C2, pyrimidine C5 and C6, and purine C8. The differences between the measured RCSAs and the values back-calculated using three nucleobase carbon CSA sets [D. Stueber, D.M. Grant, 13C and 15N chemical shift tensors in adenosine, guanosine dihydrate, 2'-deoxythymidine, and cytidine, J. Am. Chem. Soc. 124 (2002) 10539-10551; D. Sitkoff, D.A. Case, Theories of chemical shift anisotropies in proteins and nucleic acids, Prog. NMR Spectrosc. 32 (1998) 165-190; R. Fiala, J. Czernek, V. Sklenar, Transverse relaxation optimized triple-resonance NMR experiments for nucleic acids, J. Biomol. NMR 16 (2000) 291-302] reported previously for mononucleotides (1.4 Hz) is significantly smaller than the predicted RCSA range (-10-10 Hz) but remains larger than the RCSA measurement uncertainty (0.8 Hz). Fitting of the traceless principal CSA values to the measured RCSAs using a grid search procedure yields a cytosine C5 CSA magnitude (CSAa=(3/2.(delta11(2)+delta22(2)+delta33(2)))1/2=173+/-21 ppm), which is significantly higher than the reported mononucleotide values (131-138 ppm) and a guanine C8 CSAa (148+/-13 ppm) that is in very good agreement with the mononucleotide value reported by solid-state NMR [134 ppm, D. Stueber, D.M. Grant, 13C and (15)N chemical shift tensors in adenosine, guanosine dihydrate, 2'-deoxythymidine, and cytidine, J. Am. Chem. Soc. 124 (2002) 10539-10551]. Owing to a unique sensitivity to directions normal to the base plane, the RCSAs can be translated into useful long-range orientational constraints for RNA structure determination even after allowing for substantial uncertainty in the nucleobase carbon CSA tensors.

Carbon-13 chemical shift anisotropy in DNA bases from field dependence of solution NMR relaxation rates

Knowledge of (13)C chemical shift anisotropy (CSA) in nucleotide bases is important for the interpretation of solution-state NMR relaxation data in terms of local dynamic properties of DNA and RNA. Accurate knowledge of the CSA becomes particularly important at high magnetic fields, prerequisite for adequate spectral resolution in larger oligonucleotides. Measurement of (13)C relaxation rates of protonated carbons in the bases of the so-called Dickerson dodecamer, d(CGCGAATTCGCG)(2), at 500 and 800 MHz (1)H frequency, together with the previously characterized structure and diffusion tensor yields CSA values for C5 in C, C6 in C and T, C8 in A and G, and C2 in A that are closest to values previously reported on the basis of solid-state FIREMAT NMR measurements, and mostly larger than values obtained by in vacuo DFT calculations. Owing to the noncollinearity of dipolar and CSA interactions, interpretation of the NMR relaxation rates is particularly sensitive to anisotropy of rotational diffusion, and use of isotropic diffusion models can result in considerable errors.