Solution structure of a GAAA tetraloop receptor RNA

High-throughput determination of RNA tertiary contact thermodynamics by quantitative DMS chemical mapping

Structured RNAs often contain long-range tertiary contacts that are critical to their function. Despite the importance of tertiary contacts, methods to measure their thermodynamics are low throughput or require specialized instruments. Here, we introduce a new quantitative chemical mapping method (qMaPseq) to measure Mg2+-induced formation of tertiary contact thermodynamics in a high-throughput manner using standard biochemistry equipment. With qMaPseq, we measured the ΔG of 98 unique tetraloop/tetraloop receptor (TL/TLR) variants in a one-pot reaction. These results agree well with measurements from specialized instruments (R2= 0.64). Furthermore, the DMS reactivity of the TL directly correlates to the stability of the contact (R2= 0.68), the first direct evidence that a single DMS reactivity measurement reports on thermodynamics. Combined with structure prediction, DMS reactivity allowed the development of experimentally accurate 3D models of TLR mutants. These results demonstrate that qMaPseq is broadly accessible, high-throughput and directly links DMS reactivity to thermodynamics.

Acylation probing of "generic" RNA libraries reveals critical influence of loop constraints on reactivity

The reactivity of RNA 2'-OH acylation is broadly useful both in probing structure and in preparing conjugates. To date, this reactivity has been analyzed in limited sets of biological RNA sequences, leaving open questions of how reactivity varies inherently without regard to sequence in structured contexts. We constructed and probed "generic" structured RNA libraries using homogeneous loop sequences, employing deep sequencing to carry out a systematic survey of reactivity. We find a wide range of RNA reactivities among single-stranded sequences, with nearest neighbors playing substantial roles. Remarkably, certain small loops are found to be far more reactive on average (up to 4,000-fold) than single-stranded RNAs, due to conformational constraints that enhance reactivity. Among loops, we observe large variations in reactivity based on size, type, and position. The results lend insights into RNA designs for achieving high-efficiency local conjugation and provide new opportunities to refine structure analysis.

Mg2+ Sensing by an RNA Fragment: Role of Mg2+-Coordinated Water Molecules

RNA molecules selectively bind to specific metal ions to populate their functional active states, making it important to understand their source of ion selectivity. In large RNA systems, metal ions interact with the RNA at multiple locations, making it difficult to decipher the precise role of ions in folding. To overcome this complexity, we studied the role of different metal ions (Mg²⁺, Ca²⁺, and K⁺) in the folding of a small RNA hairpin motif (5'-ucCAAAga-3') using unbiased all-atom molecular dynamics simulations. The advantage of studying this system is that it requires specific binding of a single metal ion to fold to its native state. We find that even for this small RNA, the folding free energy surface (FES) is multidimensional as different metal ions present in the solution can simultaneously facilitate folding. The FES shows that specific binding of a metal ion is indispensable for its folding. We further show that in addition to the negatively charged phosphate groups, the spatial organization of electronegative nucleobase atoms drives the site-specific binding of the metal ions. Even though the binding site cannot discriminate between different metal ions, RNA folds efficiently only in a Mg²⁺ solution. We show that the rigid network of Mg²⁺-coordinated water molecules facilitates the formation of important interactions in the transition state. The other metal ions such as K⁺ and Ca²⁺ cannot facilitate the formation of such interactions. These results allow us to hypothesize possible metal-sensing mechanisms in large metalloriboswitches and also provide useful insights into the design of appropriate collective variables for studying large RNA molecules using enhanced sampling methods.

Responsive self-assembly of tectoRNAs with loop-receptor interactions from the tetrahydrofolate (THF) riboswitch

Naturally occurring RNAs are known to exhibit a high degree of modularity, whereby specific structural modules (or motifs) can be mixed and matched to create new molecular architectures. The modular nature of RNA also affords researchers the ability to characterize individual structural elements in controlled synthetic contexts in order to gain new and critical insights into their particular structural features and overall performance. Here, we characterized the binding affinity of a unique loop–receptor interaction found in the tetrahydrofolate (THF) riboswitch using rationally designed self-assembling tectoRNAs. Our work suggests that the THF loop–receptor interaction has been fine-tuned for its particular role as a riboswitch component. We also demonstrate that the thermodynamic stability of this interaction can be modulated by the presence of folinic acid, which induces a local structural change at the level of the loop–receptor. This corroborates the existence of a THF binding site within this tertiary module and paves the way for its potential use as a THF responsive module for RNA nanotechnology and synthetic biology.

Deducing putative ancestral forms of GNRA/receptor interactions from the ribosome

Stable RNAs rely on a vast repertoire of long-range interactions to assist in the folding of complex cellular machineries such as the ribosome. The universally conserved L39/H89 interaction is a long-range GNRA-like/receptor interaction localized in proximity to the peptidyl transferase center of the large subunit of the ribosome. Because of its central location, L39/H89 likely originated at an early evolutionary stage of the ribosome and played a significant role in its early function. However, L39/H89 self-assembly is impaired outside the ribosomal context. Herein, we demonstrate that structural modularity principles can be used to re-engineer L39/H89 to self-assemble in vitro. The new versions of L39/H89 improve affinity and loop selectivity by several orders of magnitude and retain the structural and functional features of their natural counterparts. These versions of L39/H89 are proposed to be ancestral forms of L39/H89 that were capable of assembling and folding independently from proteins and post-transcriptional modifications. This work demonstrates that novel RNA modules can be rationally designed by taking advantage of the modular syntax of RNA. It offers the prospect of creating new biochemical models of the ancestral ribosome and increases the tool kit for RNA nanotechnology and synthetic biology.

The Story of RNA Folding, as Told in Epochs

The past decades have witnessed tremendous developments in our understanding of RNA biology. At the core of these advances have been studies aimed at discerning RNA structure and at understanding the forces that influence the RNA folding process. It is easy to take the present state of understanding for granted, but there is much to be learned by considering the path to our current understanding, which has been tortuous, with the birth and death of models, the adaptation of experimental tools originally developed for characterization of protein structure and catalysis, and the development of novel tools for probing RNA. In this review we tour the stages of RNA folding studies, considering them as "epochs" that can be generalized across scientific disciplines. These epochs span from the discovery of catalytic RNA, through biophysical insights into the putative primordial RNA World, to characterization of structured RNAs, the building and testing of models, and, finally, to the development of models with the potential to yield generalizable predictive and quantitative models for RNA conformational, thermodynamic, and kinetic behavior. We hope that this accounting will aid others as they navigate the many fascinating questions about RNA and its roles in biology, in the past, present, and future.

Single-Molecule Fluorescence Reveals Commonalities and Distinctions among Natural and in Vitro-Selected RNA Tertiary Motifs in a Multistep Folding Pathway

Decades of study of the RNA folding problem have revealed that diverse and complex structured RNAs are built from a common set of recurring structural motifs, leading to the perspective that a generalizable model of RNA folding may be developed from understanding of the folding properties of individual structural motifs. We used single-molecule fluorescence to dissect the kinetic and thermodynamic properties of a set of variants of a common tertiary structural motif, the tetraloop/tetraloop-receptor (TL/TLR). Our results revealed a multistep TL/TLR folding pathway in which preorganization of the ubiquitous AA-platform submotif precedes the formation of the docking transition state and tertiary A-minor hydrogen bond interactions form after the docking transition state. Differences in ion dependences between TL/TLR variants indicated the occurrence of sequence-dependent conformational rearrangements prior to and after the formation of the docking transition state. Nevertheless, varying the junction connecting the TL/TLR produced a common kinetic and ionic effect for all variants, suggesting that the global conformational search and compaction electrostatics are energetically independent from the formation of the tertiary motif contacts. We also found that in vitro-selected variants, despite their similar stability at high Mg²⁺ concentrations, are considerably less stable than natural variants under near-physiological ionic conditions, and the occurrence of the TL/TLR sequence variants in Nature correlates with their thermodynamic stability in isolation. Overall, our findings are consistent with modular but complex energetic properties of RNA structural motifs and will aid in the eventual quantitative description of RNA folding from its secondary and tertiary structural elements.

Combining asymmetric 13C-labeling and isotopic filter/edit NOESY: a novel strategy for rapid and logical RNA resonance assignment

Although ∼98% of the human genomic output is transcribed as non-protein coding RNA, <2% of the protein data bank structures comprise RNA. This huge structural disparity stems from combined difficulties of crystallizing RNA for X-ray crystallography along with extensive chemical shift overlap and broadened linewidths associated with NMR of RNA. While half of the deposited RNA structures in the PDB were solved by NMR methods, the usefulness of NMR is still limited by the high cost of sample preparation and challenges of resonance assignment. Here we propose a novel strategy for resonance assignment that combines new strategic 13C labeling technologies with filter/edit type NOESY experiments to greatly reduce spectral complexity and crowding. This new strategy allowed us to assign important non-exchangeable resonances of proton and carbon (1', 2', 2, 5, 6 and 8) nuclei using only one sample and <24 h of NMR instrument time for a 27 nt model RNA. The method was further extended to assigning a 6 nt bulge from a 61 nt viral RNA element justifying its use for a wide range RNA chemical shift resonance assignment problems.

Quantitative tests of a reconstitution model for RNA folding thermodynamics and kinetics

Decades of study of the architecture and function of structured RNAs have led to the perspective that RNA tertiary structure is modular, made of locally stable domains that retain their structure across RNAs. We formalize a hypothesis inspired by this modularity-that RNA folding thermodynamics and kinetics can be quantitatively predicted from separable energetic contributions of the individual components of a complex RNA. This reconstitution hypothesis considers RNA tertiary folding in terms of ΔG_align, the probability of aligning tertiary contact partners, and ΔG_tert, the favorable energetic contribution from the formation of tertiary contacts in an aligned state. This hypothesis predicts that changes in the alignment of tertiary contacts from different connecting helices and junctions (ΔG_HJH) or from changes in the electrostatic environment (ΔG_+/-) will not affect the energetic perturbation from a mutation in a tertiary contact (ΔΔG_tert). Consistent with these predictions, single-molecule FRET measurements of folding of model RNAs revealed constant ΔΔG_tert values for mutations in a tertiary contact embedded in different structural contexts and under different electrostatic conditions. The kinetic effects of these mutations provide further support for modular behavior of RNA elements and suggest that tertiary mutations may be used to identify rate-limiting steps and dissect folding and assembly pathways for complex RNAs. Overall, our model and results are foundational for a predictive understanding of RNA folding that will allow manipulation of RNA folding thermodynamics and kinetics. Conversely, the approaches herein can identify cases where an independent, additive model cannot be applied and so require additional investigation.

Revised RNA Dihedral Parameters for the Amber Force Field Improve RNA Molecular Dynamics

The backbone dihedral parameters of the Amber RNA force field were improved by fitting using multiple linear regression to potential energies determined by quantum chemistry calculations. Five backbone and four glycosidic dihedral parameters were fit simultaneously to reproduce the potential energies determined by a high-level density functional theory calculation (B97D3 functional with the AUG-CC-PVTZ basis set). Umbrella sampling was used to determine conformational free energies along the dihedral angles, and these better agree with the population of conformations observed in the protein data bank for the new parameters than for the conventional parameters. Molecular dynamics simulations performed on a set of hairpin loops, duplexes and tetramers with the new parameter set show improved modeling for the structures of tetramers CCCC, CAAU, and GACC, and an RNA internal loop of noncanonical pairs, as compared to the conventional parameters. For the tetramers, the new parameters largely avoid the incorrect intercalated structures that dominate the conformational samples from the conventional parameters. For the internal loop, the major conformation solved by NMR is stable with the new parameters, but not with the conventional parameters. The new force field performs similarly to the conventional parameters for the UUCG and GCAA hairpin loops and the [U(UA)6A]2 duplex.

Amino Acid Specific Effects on RNA Tertiary Interactions: Single-Molecule Kinetic and Thermodynamic Studies

In light of the current models for an early RNA-based universe, the potential influence of simple amino acids on tertiary folding of ribozymal RNA into biochemically competent structures is speculated to be of significant evolutionary importance. In the present work, the folding-unfolding kinetics of a ubiquitous tertiary interaction motif, the GAAA tetraloop-tetraloop receptor (TL-TLR), is investigated by single-molecule fluorescence resonance energy transfer spectroscopy in the presence of natural amino acids both with (e.g., lysine, arginine) and without (e.g., glycine) protonated side chain residues. By way of control, we also investigate the effects of a special amino acid (e.g., proline) and amino acid mimetic (e.g., betaine) that contain secondary or quaternary amine groups rather than a primary amine group. This combination permits systematic study of amino acid induced (or amino acid like) RNA folding dynamics as a function of side chain complexity, pK_a, charge state, and amine group content. Most importantly, each of the naturally occurring amino acids is found to destabilize the TL-TLR tertiary folding equilibrium, the kinetic origin of which is dominated by a decrease in the folding rate constant (k_dock), also affected by a strongly amino acid selective increase in the unfolding rate constant (k_undock). To further elucidate the underlying thermodynamics, single-molecule equilibrium constants (K_eq) for TL-TLR folding have been probed as a function of temperature, which reveal an amino acid dependent decrease in both overall exothermicity (ΔΔH° > 0) and entropic cost (-TΔΔS° < 0) for the overall folding process. Temperature-dependent studies on the folding/unfolding kinetic rate constants reveal analogous amino acid specific changes in both enthalpy (ΔΔH^⧧) and entropy (ΔΔS^⧧) for accessing the transition state barrier. The maximum destabilization of the TL-TLR tertiary interaction is observed for arginine, which is consistent with early studies of arginine and guanidine ion-inhibited self-splicing kinetics for the full Tetrahymena ribozyme [ Yarus , M. ; Christian , E. L. Nature 1989 , 342 , 349 - 350 ; Yarus , M. Science 1988 , 240 , 1751 - 1758 ].

Kinetic and thermodynamic framework for P4-P6 RNA reveals tertiary motif modularity and modulation of the folding preferred pathway

The past decade has seen a wealth of 3D structural information about complex structured RNAs and identification of functional intermediates. Nevertheless, developing a complete and predictive understanding of the folding and function of these RNAs in biology will require connection of individual rate and equilibrium constants to structural changes that occur in individual folding steps and further relating these steps to the properties and behavior of isolated, simplified systems. To accomplish these goals we used the considerable structural knowledge of the folded, unfolded, and intermediate states of P4-P6 RNA. We enumerated structural states and possible folding transitions and determined rate and equilibrium constants for the transitions between these states using single-molecule FRET with a series of mutant P4-P6 variants. Comparisons with simplified constructs containing an isolated tertiary contact suggest that a given tertiary interaction has a stereotyped rate for breaking that may help identify structural transitions within complex RNAs and simplify the prediction of folding kinetics and thermodynamics for structured RNAs from their parts. The preferred folding pathway involves initial formation of the proximal tertiary contact. However, this preference was only ∼10 fold and could be reversed by a single point mutation, indicating that a model akin to a protein-folding contact order model will not suffice to describe RNA folding. Instead, our results suggest a strong analogy with a modified RNA diffusion-collision model in which tertiary elements within preformed secondary structures collide, with the success of these collisions dependent on whether the tertiary elements are in their rare binding-competent conformations.

The RNA 3D Motif Atlas: Computational methods for extraction, organization and evaluation of RNA motifs

RNA 3D motifs occupy places in structured RNA molecules that correspond to the hairpin, internal and multi-helix junction "loops" of their secondary structure representations. As many as 40% of the nucleotides of an RNA molecule can belong to these structural elements, which are distinct from the regular double helical regions formed by contiguous AU, GC, and GU Watson-Crick basepairs. With the large number of atomic- or near atomic-resolution 3D structures appearing in a steady stream in the PDB/NDB structure databases, the automated identification, extraction, comparison, clustering and visualization of these structural elements presents an opportunity to enhance RNA science. Three broad applications are: (1) identification of modular, autonomous structural units for RNA nanotechnology, nanobiology and synthetic biology applications; (2) bioinformatic analysis to improve RNA 3D structure prediction from sequence; and (3) creation of searchable databases for exploring the binding specificities, structural flexibility, and dynamics of these RNA elements. In this contribution, we review methods developed for computational extraction of hairpin and internal loop motifs from a non-redundant set of high-quality RNA 3D structures. We provide a statistical summary of the extracted hairpin and internal loop motifs in the most recent version of the RNA 3D Motif Atlas. We also explore the reliability and accuracy of the extraction process by examining its performance in clustering recurrent motifs from homologous ribosomal RNA (rRNA) structures. We conclude with a summary of remaining challenges, especially with regard to extraction of multi-helix junction motifs.

From static to dynamic: the need for structural ensembles and a predictive model of RNA folding and function

To understand RNA, it is necessary to move beyond a descriptive categorization towards quantitative predictions of its molecular conformations and functional behavior. An incisive approach to understanding the function and folding of biological RNA systems involves characterizing small, simple components that are largely responsible for the behavior of complex systems including helix-junction-helix elements and tertiary motifs. State-of-the-art methods have permitted unprecedented insight into the conformational ensembles of these elements revealing, for example, that conformations of helix-junction-helix elements are confined to a small region of the ensemble, that this region is highly dependent on the junction's topology, and that the correct alignment of tertiary motifs may be a rare conformation on the overall folding landscape. Further characterization of RNA components and continued development of experimental and computational methods with the goal of quantitatively predicting RNA folding and functional behavior will be critical to understanding biological RNA systems.

NMR characterization of an oligonucleotide model of the miR-21 pre-element

We have used NMR spectroscopy to characterize an oligonucleotide stem loop structure based on the pre-element of an oncogenic microRNA, miR-21. This predicted stem-loop structure is cleaved from the precursor of miR-21 (pre-miR-21) by the nuclease Dicer. It is also a critical feature recognized by the protein complex that converts the primary transcript (pri-miR-21) into the pre-miRNA. The secondary structure of the native sequence is poorly defined by NMR due to rapid exchange of imino protons with solvent; however, replacement of two adjacent putative G•U base pairs with G•C base pairs retains the conformation of the hairpin observed by chemical probing and stabilizes it sufficiently to observe most of the imino proton resonances of the molecule. The observed resonances are consistent with the predicted secondary structure. In addition, a peak due to a loop uridine suggests an interaction between it and a bulged uridine in the stem. Assignment of non-exchangeable proton resonances and characterization of NOEs and coupling constants allows inference of the following features of the structure: extrahelicity of a bulged adenosine, deviation from A-form geometry in a base-paired stem, and consecutive stacking of the adenosines in the 5′ side of the loop, the guanosine of the closing base pair, and a cross-strand adenosine. Modeling of the structure by restrained molecular dynamics suggests a basis for the interaction between the loop uridine, the bulged uridine in the stem, and an A•U base pair in the stem.

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.

Exploring the electrostatic energy landscape for tetraloop-receptor docking

It has long been appreciated that Mg(2+) is essential for the stabilization of RNA tertiary structure. However, the problem of quantitative prediction for the ion effect in tertiary structure folding remains. By using the virtual bond RNA folding model (Vfold) to generate RNA conformations and the newly improved tightly bound ion model (TBI) to treat ion-RNA interactions, we investigate Mg(2+)-facilitated tetraloop-receptor docking. For the specific construct of the tetraloop-receptor system, the theoretical analysis shows that the Mg(2+)-induced stabilizing force for the docked state is predominantly entropic and the major contribution comes from the entropy of the diffusive ions. Furthermore, our results show that Mg(2+) ions promote tetraloop-receptor docking mainly through the entropy of the diffusive ions. The theoretical prediction agrees with experimental analysis. The method developed in this paper, which combines the theory for the (Mg(2+)) ion effects in RNA folding and RNA conformational sampling, may provide a useful framework for studying the ion effect in the folding of more complex RNA structures.

An RNA folding motif: GNRA tetraloop-receptor interactions

Nearly two decades after Westhof and Michel first proposed that RNA tetraloops may interact with distal helices, tetraloop–receptor interactions have been recognized as ubiquitous elements of RNA tertiary structure. The unique architecture of GNRA tetraloops (N=any nucleotide, R=purine) enables interaction with a variety of receptors, e.g., helical minor grooves and asymmetric internal loops. The most common example of the latter is the GAAA tetraloop–11 nt tetraloop receptor motif. Biophysical characterization of this motif provided evidence for the modularity of RNA structure, with applications spanning improved crystallization methods to RNA tectonics. In this review, we identify and compare types of GNRA tetraloop–receptor interactions. Then we explore the abundance of structural, kinetic, and thermodynamic information on the frequently occurring and most widely studied GAAA tetraloop–11 nt receptor motif. Studies of this interaction have revealed powerful paradigms for structural assembly of RNA, as well as providing new insights into the roles of cations, transition states and protein chaperones in RNA folding pathways. However, further research will clearly be necessary to characterize other tetraloop–receptor and long-range tertiary binding interactions in detail – an important milestone in the quantitative prediction of free energy landscapes for RNA folding.

An in vitro-selected RNA receptor for the GAAC loop: modular receptor for non-GNRA-type tetraloop

Although artificial RNA motifs that can functionally replace the GNRA/receptor interaction, a class of RNA-RNA interacting motifs, were isolated from RNA libraries and used to generate designer RNA structures, receptors for non-GNRA tetraloops have not been found in nature or selected from RNA libraries. In this study, we report successful isolation of a receptor motif interacting with GAAC, a non-GNRA tetraloop, from randomized sequences embedded in a catalytic RNA. Biochemical characterization of the GAAC/receptor interacting motif within three structural contexts showed its binding affinity, selectivity and structural autonomy. The motif has binding affinity comparable with that of a GNRA/receptor, selectivity orthogonal to GNRA/receptors and structural autonomy even in a large RNA context. These features would be advantageous for usage of the motif as a building block for designer RNAs. The isolated motif can also be used as a query sequence to search for unidentified naturally occurring GANC receptor motifs.

Molecular paleontology: a biochemical model of the ancestral ribosome

Ancient components of the ribosome, inferred from a consensus of previous work, were constructed in silico, in vitro and in vivo. The resulting model of the ancestral ribosome presented here incorporates ∼20% of the extant 23S rRNA and fragments of five ribosomal proteins. We test hypotheses that ancestral rRNA can: (i) assume canonical 23S rRNA-like secondary structure, (ii) assume canonical tertiary structure and (iii) form native complexes with ribosomal protein fragments. Footprinting experiments support formation of predicted secondary and tertiary structure. Gel shift, spectroscopic and yeast three-hybrid assays show specific interactions between ancestral rRNA and ribosomal protein fragments, independent of other, more recent, components of the ribosome. This robustness suggests that the catalytic core of the ribosome is an ancient construct that has survived billions of years of evolution without major changes in structure. Collectively, the data here support a model in which ancestors of the large and small subunits originated and evolved independently of each other, with autonomous functionalities.

The structural stabilization of the κ three-way junction by Mg(II) represents the first step in the folding of a group II intron

Folding of group II introns is characterized by a first slow compaction of domain 1 (D1) followed by the rapid docking of other domains to this scaffold. D1 compaction initiates in a small subregion encompassing the κ and ζ elements. These two tertiary elements are also the major interaction sites with domain 5 to form the catalytic core. Here, we provide the first characterization of the structure adopted at an early folding step and show that the folding control element can be narrowed down to the three-way junction with the κ motif. In our nuclear magnetic resonance studies of this substructure derived from the yeast mitochondrial group II intron Sc.ai5γ, we show that a high affinity Mg(II) ion stabilizes the κ element and enables coaxial stacking between helices d′ and d′′, favoring a rigid duplex across the three-way junction. The κ-element folds into a stable GAAA-tetraloop motif and engages in A-minor interactions with helix d′. The addition of cobalt(III)hexammine reveals three distinct binding sites. The Mg(II)-promoted structural rearrangement and rigidification of the D1 core can be identified as the first micro-step of D1 folding.

Denaturation of RNA secondary and tertiary structure by urea: simple unfolded state models and free energy parameters account for measured m-values

To investigate the mechanism by which urea destabilizes RNA structure, urea-induced unfolding of four different RNA secondary and tertiary structures was quantified in terms of an m-value, the rate at which the free energy of unfolding changes with urea molality. From literature data and our osmometric study of a backbone analogue, we derived average interaction potentials (per square angstrom of solvent accessible surface) between urea and three kinds of RNA surfaces: phosphate, ribose, and base. Estimates of the increases in solvent accessible surface areas upon RNA denaturation were based on a simple model of unfolded RNA as a combination of helical and single-strand segments. These estimates, combined with the three interaction potentials and a term to account for interactions of urea with released ions, yield calculated m-values that are in good agreement with experimental values (200 mm monovalent salt). Agreement was obtained only if single-stranded RNAs were modeled in a highly stacked, A-form conformation. The primary driving force for urea-induced denaturation is the strong interaction of urea with the large surface areas of bases that become exposed upon denaturation of either RNA secondary or tertiary structure, though interactions of urea with backbone and released ions may account for up to a third of the m-value. Urea m-values for all four RNAs are salt-dependent, which we attribute to an increased extension (or decreased charge density) of unfolded RNAs with an increased urea concentration. The sensitivity of the urea m-value to base surface exposure makes it a potentially useful probe of the conformations of RNA unfolded states.

The role of counterion valence and size in GAAA tetraloop-receptor docking/undocking kinetics

For RNA to fold into compact, ordered structures, it must overcome electrostatic repulsion between negatively charged phosphate groups by counterion recruitment. A physical understanding of the counterion-assisted folding process requires addressing how cations kinetically and thermodynamically control the folding equilibrium for each tertiary interaction in a full-length RNA. In this work, single-molecule FRET (fluorescence resonance energy transfer) techniques are exploited to isolate and explore the cation-concentration-dependent kinetics for formation of a ubiquitous RNA tertiary interaction, that is, the docking/undocking of a GAAA tetraloop with its 11-nt receptor. Rate constants for docking (k(dock)) and undocking (k(undock)) are obtained as a function of cation concentration, size, and valence, specifically for the series Na(+), K(+), Mg(2+), Ca(2+), Co(NH(3))(6)(3+), and spermidine(3+). Increasing cation concentration acceleratesk(dock)dramatically but achieves only a slight decrease in k(undock). These results can be kinetically modeled using parallel cation-dependent and cation-independent docking pathways, which allows for isolation of the folding kinetics from the interaction energetics of the cations with the undocked and docked states, respectively. This analysis reveals a preferential interaction of the cations with the transition state and docked state as compared to the undocked RNA, with the ion-RNA interaction strength growing with cation valence. However, the corresponding number of cations that are taken up by the RNA upon folding decreases with charge density of the cation. The only exception to these behaviors is spermidine(3+), whose weaker influence on the docking equilibria with respect to Co(NH(3))(6)(3+) can be ascribed to steric effects preventing complete neutralization of the RNA phosphate groups.

Folding path of P5abc RNA involves direct coupling of secondary and tertiary structures

Folding mechanisms in which secondary structures are stabilized through the formation of tertiary interactions are well documented in protein folding but challenge the folding hierarchy normally assumed for RNA. However, it is increasingly clear that RNA could fold by a similar mechanism. P5abc, a small independently folding tertiary domain of the Tetrahymena thermophila group I ribozyme, is known to fold by a secondary structure rearrangement involving helix P5c. However, the extent of this rearrangement and the precise stage of folding that triggers it are unknown. We use experiments and simulations to show that the P5c helix switches to the native secondary structure late in the folding pathway and is directly coupled to the formation of tertiary interactions in the A-rich bulge. P5c mutations show that the switch in P5c is not rate-determining and suggest that non-native interactions in P5c aid folding rather than impede it. Our study illustrates that despite significant differences in the building blocks of proteins and RNA, there may be common ways in which they self-assemble.

Thermodynamic origins of monovalent facilitated RNA folding

Cations have long been associated with formation of native RNA structure and are commonly thought to stabilize the formation of tertiary contacts by favorably interacting with the electrostatic potential of the RNA, giving rise to an "ion atmosphere". A significant amount of information regarding the thermodynamics of structural transitions in the presence of an ion atmosphere has accumulated and suggests stabilization is dominated by entropic terms. This work provides an analysis of how RNA-cation interactions affect the entropy and enthalpy associated with an RNA tertiary transition. Specifically, temperature-dependent single-molecule fluorescence resonance energy transfer studies have been exploited to determine the free energy (ΔG°), enthalpy (ΔH°), and entropy (ΔS°) of folding for an isolated tetraloop-receptor tertiary interaction as a function of Na(+) concentration. Somewhat unexpectedly, increasing the Na(+) concentration changes the folding enthalpy from a strongly exothermic process [e.g., ΔH° = -26(2) kcal/mol at 180 mM] to a weakly exothermic process [e.g., ΔH° = -4(1) kcal/mol at 630 mM]. As a direct corollary, it is the strong increase in folding entropy [Δ(ΔS°) > 0] that compensates for this loss of exothermicity for the achievement of more favorable folding [Δ(ΔG°) < 0] at higher Na(+) concentrations. In conjunction with corresponding measurements of the thermodynamics of the transition state barrier, these data provide a detailed description of the folding pathway associated with the GAAA tetraloop-receptor interaction as a function of Na(+) concentration. The results support a potentially universal mechanism for monovalent facilitated RNA folding, whereby an increasing monovalent concentration stabilizes tertiary structure by reducing the entropic penalty for folding.

Entropic origin of Mg2+-facilitated RNA folding

Mg(2+) is essential for the proper folding and function of RNA, though the effect of Mg(2+) concentration on the free energy, enthalpy, and entropy landscapes of RNA folding is unknown. This work exploits temperature-controlled single-molecule FRET methods to address the thermodynamics of RNA folding pathways by probing the intramolecular docking/undocking kinetics of the ubiquitous GAAA tetraloop-receptor tertiary interaction as a function of [Mg(2+)]. These measurements yield the barrier and standard state enthalpies, entropies, and free energies for an RNA tertiary transition, in particular, revealing the thermodynamic origin of [Mg(2+)]-facilitated folding. Surprisingly, these studies reveal that increasing [Mg(2+)] promotes tetraloop-receptor interaction by reducing the entropic barrier (-TΔS(++)(dock)) and the overall entropic penalty (-TΔS(+) (dock)) for docking, with essentially negligible effects on both the activation enthalpy (ΔH(++)(dock)) and overall exothermicity (ΔH(+)(dock)). These observations contrast with the conventional notion that increasing [Mg(2+)] facilitates folding by minimizing electrostatic repulsion of opposing RNA helices, which would incorrectly predict a decrease in ΔH(++)(dock)) and ΔH(+)(dock)) with [Mg(2+)]. Instead we propose that higher [Mg(2+)] can aid RNA folding by decreasing the entropic penalty of counterion uptake and by reducing disorder of the unfolded conformational ensemble.

Functional roles of a tetraloop/receptor interacting module in a cyclic di-GMP riboswitch

Riboswitches are a class of structural RNAs that regulate transcription and translation through specific recognition of small molecules. Riboswitches are attractive not only as drug targets for novel antibiotics but also as modular tools for controlling gene expression. Sequence comparison of a class of riboswitches that sense cyclic di-GMP (type-I c-di-GMP riboswitches) revealed that this type of riboswitch frequently shows a GAAA loop/receptor interaction between P1 and P3 elements. In the crystal structures of a type-I c-di-GMP riboswitch from Vibrio cholerae (the Vc2 riboswitch), the GNRA loop/receptor interaction assembled P2 and P3 stems to organize a ligand-binding pocket. In this study, the functional importance of the GAAA loop-receptor interaction in the Vc2 riboswitch was examined. A series of variant Vc2 riboswitches with mutations in the GAAA loop/receptor interaction were assayed for their switching abilities. In mutants with mutations in the P2 GAAA loop, expression of the reporter gene was reduced to approximately 40% - 60% of that in the wild-type. However, mutants in which the P3 receptor motif was substituted with base pairs were as active as the wild-type. These results suggested that the GAAA loop/receptor interaction does not simply establish the RNA 3D structure but docking of P2 GAAA loop reduces the flexibility of the GAAA receptor motif in the P3 element. This mechanism was supported by a variant riboswitch bearing a theophylline aptamer module in P3 the structural rigidity of which could be modulated by the small molecule theophylline.

New insights into the fundamental role of topological constraints as a determinant of two-way junction conformation

Recent studies have shown that topological constraints encoded at the RNA secondary structure level involving basic steric and stereochemical forces can significantly restrict the orientations sampled by helices across two-way RNA junctions. Here, we formulate these topological constraints in greater quantitative detail and use this topological framework to rationalize long-standing but poorly understood observations regarding the basic behavior of RNA two-way junctions. Notably, we show that the asymmetric nature of the A-form helix and the finite length of a bulge provide a physical basis for the experimentally observed directionality and bulge-length amplitude dependence of bulge induced inter-helical bends. We also find that the topologically allowed space can be modulated by variations in sequence, particularly with the addition of non-canonical GU base pairs at the junction, and, surprisingly, by the length of the 5′ and 3′ helices. A survey of two-way RNA junctions in the protein data bank confirms that junction residues have a strong preference to adopt looped-in, non-canonically base-paired conformations, providing a route for extending our bulge-directed framework to internal loop motifs and implying a simplified link between secondary and tertiary structure. Finally, our results uncover a new simple mechanism for coupling junction-induced topological constraints with tertiary interactions.

Structural model of the p14/SF3b155 · branch duplex complex

Human p14 (SF3b14), a component of the spliceosomal U2 snRNP, interacts directly with the pre-mRNA branch adenosine within the context of the bulged duplex formed between the pre-mRNA branch region and U2 snRNA. This association occurs early in spliceosome assembly and persists within the fully assembled spliceosome. Analysis of the crystal structure of a complex containing p14 and a peptide derived from p14-associated SF3b155 combined with the results of cross-linking studies has suggested that the branch nucleotide interacts with a pocket on a non-canonical RNA binding surface formed by the complex. Here we report a structural model of the p14 · bulged duplex interaction based on a combination of X-ray crystallography of an adenine p14/SF3b155 peptide complex, biochemical comparison of a panel of disulfide cross-linked protein-RNA complexes, and small-angle X-ray scattering (SAXS). These studies reveal specific recognition of the branch adenosine within the p14 pocket and establish the orientation of the bulged duplex RNA bound on the protein surface. The intimate association of one surface of the bulged duplex with the p14/SF3b155 peptide complex described by this model buries the branch nucleotide at the interface and suggests that p14 · duplex interaction must be disrupted before the first step of splicing.

Enthalpy-driven RNA folding: single-molecule thermodynamics of tetraloop-receptor tertiary interaction

RNA folding thermodynamics are crucial for structure prediction, which requires characterization of both enthalpic and entropic contributions of tertiary motifs to conformational stability. We explore the temperature dependence of RNA folding due to the ubiquitous GAAA tetraloop-receptor docking interaction, exploiting immobilized and freely diffusing single-molecule fluorescence resonance energy transfer (smFRET) methods. The equilibrium constant for intramolecular docking is obtained as a function of temperature (T = 21-47 degrees C), from which a van't Hoff analysis yields the enthalpy (DeltaH degrees) and entropy (DeltaS degrees) of docking. Tetraloop-receptor docking is significantly exothermic and entropically unfavorable in 1 mM MgCl(2) and 100 mM NaCl, with excellent agreement between immobilized (DeltaH degrees = -17.4 +/- 1.6 kcal/mol, and DeltaS degrees = -56.2 +/- 5.4 cal mol(-1) K(-1)) and freely diffusing (DeltaH degrees = -17.2 +/- 1.6 kcal/mol, and DeltaS degrees = -55.9 +/- 5.2 cal mol(-1) K(-1)) species. Kinetic heterogeneity in the tetraloop-receptor construct is unaffected over the temperature range investigated, indicating a large energy barrier for interconversion between the actively docking and nondocking subpopulations. Formation of the tetraloop-receptor interaction can account for approximately 60% of the DeltaH degrees and DeltaS degrees of P4-P6 domain folding in the Tetrahymena ribozyme, suggesting that it may act as a thermodynamic clamp for the domain. Comparison of the isolated tetraloop-receptor and other tertiary folding thermodynamics supports a theme that enthalpy- versus entropy-driven folding is determined by the number of hydrogen bonding and base stacking interactions.

Solution structure of the cap-independent translational enhancer and ribosome-binding element in the 3' UTR of turnip crinkle virus

The 3(') untranslated region (3(') UTR) of turnip crinkle virus (TCV) genomic RNA contains a cap-independent translation element (CITE), which includes a ribosome-binding structural element (RBSE) that participates in recruitment of the large ribosomal subunit. In addition, a large symmetric loop in the RBSE plays a key role in coordinating the incompatible processes of viral translation and replication, which require enzyme progression in opposite directions on the viral template. To understand the structural basis for the large ribosomal subunit recruitment and the intricate interplay among different parts of the molecule, we determined the global structure of the 102-nt RBSE RNA using solution NMR and small-angle x-ray scattering. This RNA has many structural features that resemble those of a tRNA in solution. The hairpins H1 and H2, linked by a 7-nucleotide linker, form the upper part of RBSE and hairpin H3 is relatively independent from the rest of the structure and is accessible to interactions. This global structure provides insights into the three-dimensional layout for ribosome binding, which may serve as a structural basis for its involvement in recruitment of the large ribosomal subunit and the switch between viral translation and replication. The experimentally determined three-dimensional structure of a functional element in the 3(') UTR of an RNA from any organism has not been previously reported. The RBSE structure represents a prototype structure of a new class of RNA structural elements involved in viral translation/replication processes.

Identification of catalytic metal ion ligands in ribozymes

Site-bound metal ions participate in the catalytic mechanisms of many ribozymes. Understanding these mechanisms therefore requires knowledge of the specific ligands on both substrate and ribozyme that coordinate these catalytic metal ions. A number of different structural and biochemical strategies have been developed and refined for identifying metal ion binding sites within ribozymes, and for assessing the catalytic contributions of the metal ions bound at those sites. We review these approaches and provide examples of their application, focusing in particular on metal ion rescue experiments and their roles in the construction of the transition state models for the Tetrahymena group I and RNase P ribozymes.

The influence of monovalent cation size on the stability of RNA tertiary structures

Many RNA tertiary structures are stable in the presence of monovalent ions alone. To evaluate the degree to which ions at or near the surfaces of such RNAs contribute to stability, the salt-dependent stability of a variety of RNA structures was measured with each of the five group I cations. The stability of hairpin secondary structures and a pseudoknot tertiary structure are insensitive to the ion identity, but the tertiary structures of two other RNAs, an adenine riboswitch and a kissing loop complex, become more stable by 2-3 kcal/mol as ion size decreases. This "default" trend is attributed to the ability of smaller ions to approach the RNA surface more closely. The degree of cation accumulation around the kissing loop complex was also inversely proportional to ion radius, perhaps because of the presence of sterically restricted pockets that can be accessed only by smaller ions. An RNA containing the tetraloop-receptor motif shows a strong (up to approximately 3 kcal/mol) preference for Na(+) or K(+) over other group I ions, consistent with the chelation of K(+) by this motif in some crystal structures. This RNA reverts to the default dependence on ion size when a base forming part of the chelation site is mutated. Lastly, an RNA aptamer for cobinamide, which was originally selected in the presence of high concentrations of LiCl, binds ligand more strongly in the presence of Li(+) than other monovalent ions. On the basis of these trends in RNA stability with group I ion size, it is argued that two features of RNA tertiary structures may promote strong interactions with ions at or near the RNA surface: negative charge densities that are higher than that in secondary structures, and the occasional presence of chelation sites, which are electronegative pockets that selectively bind ions of an optimum size.

Trypanosoma brucei ATPase subunit 6 mRNA bound to gA6-14 forms a conserved three-helical structure

T. brucei survival relies on the expression of mitochondrial genes, most of which require RNA editing to become translatable. In trypanosomes, RNA editing involves the insertion and deletion of uridylates, a developmentally regulated process directed by guide RNAs (gRNAs) and catalyzed by the editosome, a complex of proteins. The pathway for mRNA/gRNA complex formation and assembly with the editosome is still unknown. Work from our laboratory has suggested that distinct mRNA/gRNA complexes anneal to form a conserved core structure that may be important for editosome assembly. The secondary structure for the apocytochrome b (CYb) pair has been previously determined and is consistent with our model of a three-helical structure. Here, we used cross-linking and solution structure probing experiments to determine the structure of the ATPase subunit 6 (A6) mRNA hybridized to its cognate gA6-14 gRNA in different stages of editing. Our results indicate that both unedited and partially edited A6/gA6-14 pairs fold into a three-helical structure similar to the previously characterized CYb/gCYb-558 pair. These results lead us to conclude that at least two mRNA/gRNA pairs with distinct editing sites and distinct primary sequences fold to a three-helical secondary configuration that persists through the first few editing events.

Thermodynamics and folding pathway of tetraloop receptor-mediated RNA helical packing

Little is known about the thermodynamic forces that drive the folding pathways of higher-order RNA structure. In this study, we employ calorimetric [isothermal titration calorimetry (ITC) and differential scanning calorimetry (DSC)] and spectroscopic (NMR and UV) methods to characterize the thermodynamics of the GAAA tetraloop-receptor interaction, utilizing a previously described bivalent construct. ITC studies indicate that the bivalent interaction is enthalpy driven and highly stable, with a binding constant (K(obs)) of 5.5x10(6) M(-1) and enthalpy (DeltaH(obs)(o)) of -33.8 kcal/mol at 45 degrees C in 20 mM KCl and 2 mM MgCl(2). Thus, we derive the DeltaH(obs)(o) for a single tetraloop-receptor interaction to be -16.9 kcal/mol at these conditions. UV absorbance data indicate that an increase in base stacking quality contributes to the enthalpy of complex formation. These highly favorable thermodynamics are consistent with the known critical role for the tetraloop-receptor motif in the folding of large RNAs. Additionally, a significant heat capacity change (DeltaC(p,obs)(o)) of -0.24 kcal mol(-1) K(-1) was determined by ITC. DSC and UV-monitored thermal denaturation experiments indicate that the bivalent tetraloop-receptor construct follows a minimally five-state unfolding pathway and suggest the observed DeltaC(p,obs)(o) for the interaction results from a temperature-dependent unbound receptor RNA structure.

Monovalent and divalent promoted GAAA tetraloop-receptor tertiary interactions from freely diffusing single-molecule studies

Proper assembly of RNA into catalytically active three-dimensional structures requires multiple tertiary binding interactions, individual characterization of which is crucial to a detailed understanding of global RNA folding. This work focuses on single-molecule fluorescence studies of freely diffusing RNA constructs that isolate the GAAA tetraloop-receptor tertiary interaction. Freely diffusing conformational dynamics are explored as a function of Mg(2+) and Na(+) concentration, both of which promote facile docking, but with 500-fold different affinities. Systematic shifts in mean fluorescence resonance energy transfer efficiency values and line widths with increasing [Na(+)] are observed for the undocked species and can be interpreted with a Debye model in terms of electrostatic relaxation and increased flexibility in the RNA. Furthermore, we identify a 34 +/- 2% fraction of freely diffusing RNA constructs remaining undocked even at saturating [Mg(2+)] levels, which agrees quantitatively with the 32 +/- 1% fraction previously reported for immobilized constructs. This verifies that the kinetic heterogeneity observed in the docking rates is not the result of surface tethering. Finally, the K(D) value and Hill coefficient for [Mg(2+)]-dependent docking decrease significantly for [Na(+)] = 25 mM vs. 125 mM, indicating Mg(2+) and Na(+) synergy in the RNA folding process.

Biophysical analysis of nucleic acids

This overview unit provides a thorough overview of biophysical methods used for structure analysis, including X-ray diffraction, nuclear magnetic resonance, optical spectroscopy, theoretical and computational methods, and single-molecule methods. Advantages and disadvantages of the methods are compared.

Site-directed spin labeling studies on nucleic acid structure and dynamics

Site-directed spin labeling (SDSL) uses electron paramagnetic resonance (EPR) spectroscopy to monitor the behavior of a stable nitroxide radical attached at specific locations within a macromolecule such as protein, DNA, or RNA. Parameters obtained from EPR measurements, such as internitroxide distances and descriptions of the rotational motion of a nitroxide, provide unique information on features near the labeling site. With recent advances in solid-phase synthesis of nucleic acids and developments in EPR methodologies, particularly pulsed EPR technologies, SDSL has been increasingly used to study the structure and dynamics of DNA and RNA at the level of the individual nucleotides. This chapter summarizes the current SDSL studies on nucleic acids, with discussions focusing on literature from the last decade.

Comprehensive features of natural and in vitro selected GNRA tetraloop-binding receptors

Specific recognitions of GNRA tetraloops by small helical receptors are among the most widespread long-range packing interactions in large ribozymes. However, in contrast to GYRA and GAAA tetraloops, very few GNRA/receptor interactions have yet been identified to involve GGAA tetraloops in nature. A novel in vitro selection scheme based on a rigid self-assembling tectoRNA scaffold designed for isolation of intermolecular interactions with A-minor motifs has yielded new GGAA tetraloop-binding receptors with affinity in the nanomolar range. One of the selected receptors is a novel 12 nt RNA motif, (CCUGUG … AUCUGG), that recognizes GGAA tetraloop hairpin with a remarkable specificity and affinity. Its physical and chemical characteristics are comparable to those of the well-studied ‘11nt’ GAAA tetraloop receptor motif. A second less specific motif (CCCAGCCC … GAUAGGG) binds GGRA tetraloops and appears to be related to group IC3 tetraloop receptors. Mutational, thermodynamic and comparative structural analysis suggests that natural and in vitro selected GNRA receptors can essentially be grouped in two major classes of GNRA binders. New insights about the evolution, recognition and structural modularity of GNRA and A-minor RNA–RNA interactions are proposed.

NMR structures of (rGCUGAGGCU)2 and (rGCGGAUGCU)2: probing the structural features that shape the thermodynamic stability of GA pairs

The NMR structures of [see text] and [see text] are reported. The internal loop, [see text], is about 2 kcal/mol more stable than [see text] at 37 degrees C. The duplexes assemble into similar global folds characterized by the formation of tandem sheared GA pairs. The different stabilities of the loops are accompanied by differences in the local structure of the closing GU pairs. In the [see text] internal loop, the GU pairs form canonical wobble configurations with two hydrogen bonds, whereas in [see text], the GU pairs form a single hydrogen bond involving the amino group, GH22, and the carbonyl group, UO4. This pairing is similar to the GU closing pair of the 690 hairpin loop found in E. coli 16S rRNA. The [see text] and [see text] structures reveal how the subtle interplay between stacking and hydrogen bonding determines sequence dependent conformation and thermodynamic stability. Thus, this work provides structural and thermodynamic benchmarks for theoreticians in the ongoing effort to understand the sequence dependence of RNA physicochemical properties.

A folding control element for tertiary collapse of a group II intron ribozyme

Ribozymes derived from the group II intron ai5gamma collapse to a compact intermediate, folding to the native state through a slow, direct pathway that is unperturbed by kinetic traps. Molecular collapse of ribozyme D135 requires high magnesium concentrations and is thought to involve a structural element in domain 1 (D1). We used nucleotide analog interference mapping, in combination with nondenaturing gel electrophoresis, to identify RNA substructures and functional groups that are essential for D135 tertiary collapse. This revealed that the most crucial atoms for compaction are located within a small section of D1 that includes the kappa and zeta elements. This small substructure controls specific collapse of the molecule and, in later steps of the folding pathway, it forms the docking site for catalytic D5. In this way, the stage is set for proper active site formation during the earliest steps of ribozyme folding.

An alternating sheared AA pair and elements of stability for a single sheared purine-purine pair flanked by sheared GA pairs in RNA

A previous NMR structure of the duplex 5'GGU GGA GGCU/PCCG AAG CCG5' revealed an unusually stable RNA internal loop with three consecutive sheared GA pairs. Here, we report NMR studies of two duplexes, 5'GGU GGA GGCU/PCCA AAG CCG5' (replacing the UG pair with a UA closing pair) and 5'GGU GAA GGCU/PCCG AAG CCG5' (replacing the middle GA pair with an AA pair). An unusually stable loop with three consecutive sheared GA pairs forms in the duplex 5'GGU GGA GGCU/PCCA AAG CCG5'. The structure contrasts with that reported for this loop in the crystal structure of the large ribosomal subunit of Deinococcus radiodurans [Harms, J., Schluenzen, F., Zarivach, R., Bashan, A., Gat, S., Agmon, I., Bartels, H., Franceschi, F., and Yonath, A. (2001) Cell 107, 679-688]. The middle AA pair in the duplex 5'GGU GAA GGCU/PCCG AAG CCG5' rapidly exchanges orientations, resulting in alternative base stacking and pseudosymmetry with exclusively sheared pairs. The U GAA G/G AAG C internal loop is 2.1 kcal/mol less stable than the U GGA G/G AAG C internal loop at 37 degrees C. Structural, energetic, and dynamic consequences upon functional group substitutions within related 3 x 3 and 3 x 6 internal loops are also reported.

Metal ion dependence, thermodynamics, and kinetics for intramolecular docking of a GAAA tetraloop and receptor connected by a flexible linker

The GAAA tetraloop-receptor motif is a commonly occurring tertiary interaction in RNA. This motif usually occurs in combination with other tertiary interactions in complex RNA structures. Thus, it is difficult to measure directly the contribution that a single GAAA tetraloop-receptor interaction makes to the folding properties of a RNA. To investigate the kinetics and thermodynamics for the isolated interaction, a GAAA tetraloop domain and receptor domain were connected by a single-stranded A(7) linker. Fluorescence resonance energy transfer (FRET) experiments were used to probe intramolecular docking of the GAAA tetraloop and receptor. Docking was induced using a variety of metal ions, where the charge of the ion was the most important factor in determining the concentration of the ion required to promote docking {[Co(NH(3))(6)(3+)] << [Ca(2+)], [Mg(2+)], [Mn(2+)] << [Na(+)], [K(+)]}. Analysis of metal ion cooperativity yielded Hill coefficients of approximately 2 for Na(+)- or K(+)-dependent docking versus approximately 1 for the divalent ions and Co(NH(3))(6)(3+). Ensemble stopped-flow FRET kinetic measurements yielded an apparent activation energy of 12.7 kcal/mol for GAAA tetraloop-receptor docking. RNA constructs with U(7) and A(14) single-stranded linkers were investigated by single-molecule and ensemble FRET techniques to determine how linker length and composition affect docking. These studies showed that the single-stranded region functions primarily as a flexible tether. Inhibition of docking by oligonucleotides complementary to the linker was also investigated. The influence of flexible versus rigid linkers on GAAA tetraloop-receptor docking is discussed.

Single nucleotide RNA choreography

New structural analysis methods, and a tree formalism re-define and expand the RNA motif concept, unifying what previously appeared to be disparate groups of structures. We find RNA tetraloops at high frequencies, in new contexts, with unexpected lengths, and in novel topologies. The results, with broad implications for RNA structure in general, show that even at this most elementary level of organization, RNA tolerates astounding variation in conformation, length, sequence and context. However the variation is not random; it is well-described by four distinct modes, which are 3-2 switches (backbone topology variations), insertions, deletions and strand clips.

Statistical analysis of RNA backbone

Local conformation is an important determinant of RNA catalysis and binding. The analysis of RNA conformation is particularly difficult due to the large number of degrees of freedom (torsion angles) per residue. Proteins, by comparison, have many fewer degrees of freedom per residue. In this work, we use and extend classical tools from statistics and signal processing to search for clusters in RNA conformational space. Results are reported both for scalar analysis, where each torsion angle is separately studied, and for vectorial analysis, where several angles are simultaneously clustered. Adapting techniques from vector quantization and clustering to the RNA structure, we find torsion angle clusters and RNA conformational motifs. We validate the technique using well-known conformational motifs, showing that the simultaneous study of the total torsion angle space leads to results consistent with known motifs reported in the literature and also to the finding of new ones.

Structure and dynamics of an RNA tetraloop: a joint molecular dynamics and NMR study

Molecular dynamics simulations of the RNA tetraloop 5'-CGCUUUUGCG-3' with high melting temperature and significant conformational heterogeneity in explicit water solvent are presented and compared to NMR studies. The NMR data allow for a detailed test of the theoretical model, including the quality of the force field and the conformational sampling. Due to the conformational heterogeneity of the tetraloop, high temperature (350 K) and locally enhanced sampling simulations need to be invoked. The Amber98 force field leads to a good overall agreement with experimental data. Based on NMR data and a principal component analysis of the 350 K trajectory, the dynamic structure of the tetraloop is revealed. The principal component free energy surface exhibits four minima, which correspond to well-defined conformational structures that differ mainly by their base stacking in the loop region. No correlation between the motion of the sugar rings and the stacking dynamics of the loop bases is found.