Intermediates and kinetic traps in the folding of a large ribozyme revealed by circular dichroism and UV absorbance spectroscopies and catalytic activity

Watching ion-driven kinetics of ribozyme folding and misfolding caused by energetic and topological frustration one molecule at a time

Folding of ribozymes into well-defined tertiary structures usually requires divalent cations. How Mg2+ ions direct the folding kinetics has been a long-standing unsolved problem because experiments cannot detect the positions and dynamics of ions. To address this problem, we used molecular simulations to dissect the folding kinetics of the Azoarcus ribozyme by monitoring the path each molecule takes to reach the folded state. We quantitatively establish that Mg2+ binding to specific sites, coupled with counter-ion release of monovalent cations, stimulate the formation of secondary and tertiary structures, leading to diverse pathways that include direct rapid folding and trapping in misfolded structures. In some molecules, key tertiary structural elements form when Mg2+ ions bind to specific RNA sites at the earliest stages of the folding, leading to specific collapse and rapid folding. In others, the formation of non-native base pairs, whose rearrangement is needed to reach the folded state, is the rate-limiting step. Escape from energetic traps, driven by thermal fluctuations, occurs readily. In contrast, the transition to the native state from long-lived topologically trapped native-like metastable states is extremely slow. Specific collapse and formation of energetically or topologically frustrated states occur early in the assembly process.

Non-canonical nucleosides: Biomimetic triphosphorylation, incorporation into mRNA and effects on translation and structure

Recent advances in mRNA therapeutics demand efficient toolkits for the incorporation of nucleoside analogues into mRNA suitable for downstream applications. Herein, we report the application of a versatile enzyme cascade for the triphosphorylation of a broad range of nucleoside analogues, including unprotected nucleobases containing chemically labile moieties. Our biomimetic system was suitable for the preparation of nucleoside triphosphates containing adenosine, cytidine, guanosine, uridine and non-canonical core structures, as determined by capillary electrophoresis coupled to mass spectrometry. This enabled us to establish an efficient workflow for transcribing and purifying functional mRNA containing these nucleoside analogues, combined with mass spectrometric verification of analogue incorporation. Our combined methodology allows for analyses of how incorporation of nucleoside analogues that are commercially unavailable as triphosphates affect mRNA properties: The translational fidelity of the produced mRNA was demonstrated in analyses of how incorporated adenosine analogues impact translational recoding. For the SARS-CoV-2 frameshifting site, analyses of the mRNA pseudoknot structure using circular dichroism spectroscopy allowed insight into how the pharmacologically active 7-deazaadenosine destabilises RNA secondary structure, consistent with observed changes in recoding efficiency.

Cryo-EM reveals an entangled kinetic trap in the folding of a catalytic RNA

Functional RNAs fold through complex pathways that can contain misfolded "kinetic traps." A complete model of RNA folding requires understanding the formation of these misfolded states, but they are difficult to characterize because of their transient and potentially conformationally dynamic nature. We used cryo-electron microscopy (cryo-EM) to visualize a long-lived misfolded state in the folding pathway of the Tetrahymena thermophila group I intron, a paradigmatic RNA structure-function model system. The structure revealed how this state forms native-like secondary structure and tertiary contacts but contains two incorrectly crossed strands, consistent with a previous model. This incorrect topology mispositions a critical catalytic domain and cannot be resolved locally as extensive refolding is required. This work provides a structural framework for interpreting decades of biochemical and functional studies and demonstrates the power of cryo-EM for the exploration of RNA folding pathways.

Thoughts on how to think (and talk) about RNA structure

Recent events have pushed RNA research into the spotlight. Continued discoveries of RNA with unexpected diverse functions in healthy and diseased cells, such as the role of RNA as both the source and countermeasure to a severe acute respiratory syndrome coronavirus 2 infection, are igniting a new passion for understanding this functionally and structurally versatile molecule. Although RNA structure is key to function, many foundational characteristics of RNA structure are misunderstood, and the default state of RNA is often thought of and depicted as a single floppy strand. The purpose of this perspective is to help adjust mental models, equipping the community to better use the fundamental aspects of RNA structural information in new mechanistic models, enhance experimental design to test these models, and refine data interpretation. We discuss six core observations focused on the inherent nature of RNA structure and how to incorporate these characteristics to better understand RNA structure. We also offer some ideas for future efforts to make validated RNA structural information available and readily used by all researchers.

Comparative study on tertiary contacts and folding of RNase P RNAs from a psychrophilic, a mesophilic/radiation-resistant, and a thermophilic bacterium

In most bacterial type A RNase P RNAs (P RNAs), two major loop-helix tertiary contacts (L8-P4 and L18-P8) help to orient the two independently folding S- and C-domains for concerted recognition of precursor tRNA substrates. Here, we analyze the effects of mutations in these tertiary contacts in P RNAs from three different species: (i) the psychrophilic bacterium Pseudoalteromonas translucida (Ptr), (ii) the mesophilic radiation-resistant bacterium Deinococcus radiodurans (Dra), and (iii) the thermophilic bacterium Thermus thermophilus (Tth). We show by UV melting experiments that simultaneous disruption of these two interdomain contacts has a stabilizing effect on all three P RNAs. This can be inferred from reduced RNA unfolding at lower temperatures and a more concerted unfolding at higher temperatures. Thus, when the two domains tightly interact via the tertiary contacts, one domain facilitates structural transitions in the other. P RNA mutants with disrupted interdomain contacts showed severe kinetic defects that were most pronounced upon simultaneous disruption of the L8-P4 and L18-P8 contacts. At 37°C, the mildest effects were observed for the thermostable Tth RNA. A third interdomain contact, L9-P1, makes only a minor contribution to P RNA tertiary folding. Furthermore, D. radiodurans RNase P RNA forms an additional pseudoknot structure between the P9 and P12 of its S-domain. This interaction was found to be particularly crucial for RNase P holoenzyme activity at near-physiological Mg²⁺ concentrations (2 mM). We further analyzed an exceptionally stable folding trap of the G,C-rich Tth P RNA.

Proteins That Chaperone RNA Regulation

RNA-binding proteins chaperone the biological functions of noncoding RNA by reducing RNA misfolding, improving matchmaking between regulatory RNA and targets, and exerting quality control over RNP biogenesis. Recent studies of Escherichia coli CspA, HIV NCp, and E. coli Hfq are beginning to show how RNA-binding proteins remodel RNA structures. These different protein families use common strategies for disrupting or annealing RNA double helices, which can be used to understand the mechanisms by which proteins chaperone RNA-dependent regulation in bacteria.

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.

Conditioning and Robustness of RNA Boltzmann Sampling under Thermodynamic Parameter Perturbations

Understanding how RNA secondary structure prediction methods depend on the underlying nearest-neighbor thermodynamic model remains a fundamental challenge in the field. Minimum free energy (MFE) predictions are known to be "ill conditioned" in that small changes to the thermodynamic model can result in significantly different optimal structures. Hence, the best practice is now to sample from the Boltzmann distribution, which generates a set of suboptimal structures. Although the structural signal of this Boltzmann sample is known to be robust to stochastic noise, the conditioning and robustness under thermodynamic perturbations have yet to be addressed. We present here a mathematically rigorous model for conditioning inspired by numerical analysis, and also a biologically inspired definition for robustness under thermodynamic perturbation. We demonstrate the strong correlation between conditioning and robustness and use its tight relationship to define quantitative thresholds for well versus ill conditioning. These resulting thresholds demonstrate that the majority of the sequences are at least sample robust, which verifies the assumption of sampling's improved conditioning over the MFE prediction. Furthermore, because we find no correlation between conditioning and MFE accuracy, the presence of both well- and ill-conditioned sequences indicates the continued need for both thermodynamic model refinements and alternate RNA structure prediction methods beyond the physics-based ones.

Equilibrium Denaturation and Preferential Interactions of an RNA Tetraloop with Urea

Urea is an important organic cosolute with implications in maintaining osmotic stress in cells and differentially stabilizing ensembles of folded biomolecules. We report an equilibrium study of urea-induced denaturation of a hyperstable RNA tetraloop through unbiased replica exchange molecular dynamics. We find that, in addition to destabilizing the folded state, urea smooths the RNA free energy landscape by destabilizing specific configurations, and forming favorable interactions with RNA nucleobases. A linear concentration-dependence of the free energy (m-value) is observed, in agreement with the results of other RNA hairpins and proteins. Additionally, analysis of the hydrogen-bonding and stacking interactions within RNA primarily show temperature-dependence, while interactions between RNA and urea primarily show concentration-dependence. Our findings provide valuable insight into the effects of urea on RNA folding and describe the thermodynamics of a basic RNA hairpin as a function of solution chemistry.

Modulation of Conformational Equilibria in the S-Adenosylmethionine (SAM) II Riboswitch by SAM, Mg(2+), and Trimethylamine N-Oxide

The dependence of the conformation of the S-adenosylmethionine (SAM) II riboswitch on the concentration of added Mg(2+) ions and SAM, individually and in mixtures, was monitored by circular dichroism (CD) spectroscopy and by measurement of the diffusion coefficient. The results are analyzed in the context of two complementary quantitative models, both of which are consistent with a single underlying physical model. Magnesium binding sites in the open state have an affinity on average higher than the affinity of those in the compact state, but formation of the compact state is accompanied by an increase in the number of binding sites. Consequently, at low Mg(2+) concentrations, Mg(2+) binds preferentially to the open state, favoring its formation, but at high concentrations, Mg(2+) binds preferentially to the compact state. The affinity of the riboswitch for SAM increases drastically with an increased level of binding of Mg(2+) to the compact pseudoknot conformation. The effect of increasing concentrations of trimethylamine N-oxide (TMAO), a well-studied molecular crowding agent, on the conformation of the riboswitch and its affinity for SAM were also monitored by CD spectroscopy and measurement of diffusion. In the absence of added Mg(2+), high concentrations of TMAO were found to induce a conformational change compatible with the formation of the pseudoknot form but have only a small effect on the affinity of the RNA for SAM.

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

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

Tetracycline determines the conformation of its aptamer at physiological magnesium concentrations

Synthetic riboswitches are versatile tools for the study and manipulation of biological systems. Yet, the underlying mechanisms governing its structural properties and regulation under physiological conditions are poorly studied. We performed spectroscopic and calorimetric experiments to explore the folding kinetics and thermodynamics of the tetracycline-binding aptamer, which can be employed as synthetic riboswitch, in the range of physiological magnesium concentrations. The dissociation constant of the ligand-aptamer complex was found to strongly depend on the magnesium concentration. At physiological magnesium concentrations, tetracycline induces a significant conformational shift from a compact, but heterogeneous intermediate state toward the completely formed set of tertiary interactions defining the regulation-competent structure. Thus, the switching functionality of the tetracycline-binding aptamer appears to include both a conformational rearrangement toward the regulation-competent structure and its thermodynamic stabilization.

Assessing the Utility of Circular Dichroism and FTIR Spectroscopy in Monoclonal-Antibody Comparability Studies

Protein characterization is a necessary activity during development, technical transfers, and licensure. One important aspect of protein characterization is higher order structure assessment, which can be accomplished in a variety of ways. Circular dichroism (CD) and Fourier transform infrared (FTIR) spectroscopies provide global higher order structure and are routinely used to measure the overall structure for product characterization; however, their use as comparability tools is uncertain because of their insensitivity to local or small structure changes. We use a monoclonal antibody (mAb) to explore the usefulness of CD and FTIR compared with other indirect methods of structure characterization such as size-exclusion and ion-exchange chromatographies (SEC and IEC). A panel of degraded samples of a mAb was generated; their higher order structure evaluated using CD and FTIR and was found to be largely unchanged. However, the SEC and IEC chromatograms of certain degraded samples were found to have measurable changes. Based on these studies, we conclude that the application of CD and FTIR should be reserved for global higher order structure identification or product characterization only. The use of CD or FTIR comparability of mAbs should be carefully evaluated, as comparability can be sensitively determined using indirect methods based on chromatography.

Genome-wide profiling of mouse RNA secondary structures reveals key features of the mammalian transcriptome

Background

The understanding of RNA structure is a key feature toward the comprehension of RNA functions and mechanisms of action. In particular, non-coding RNAs are thought to exert their functions by specific secondary structures, but an efficient annotation on a large scale of these structures is still missing.

Results

By using a novel high-throughput method, named chemical inference of RNA structures, CIRS-seq, that uses dimethyl sulfate, and N-cyclohexyl- N'-(2-morpholinoethyl)carbodiimide metho-p-toluenesulfonate to modify RNA residues in single-stranded conformation within native deproteinized RNA secondary structures, we investigate the structural features of mouse embryonic stem cell transcripts. Our analysis reveals an unexpected higher structuring of the 5′ and 3′ untranslated regions compared to the coding regions, a reduced structuring at the Kozak sequence and stop codon, and a three-nucleotide periodicity across the coding region of messenger RNAs. We also observe that ncRNAs exhibit a higher degree of structuring with respect to protein coding transcripts. Moreover, we find that the Lin28a binding protein binds selectively to RNA motifs with a strong preference toward a single stranded conformation.

Conclusions

This work defines for the first time the complete RNA structurome of mouse embryonic stem cells, revealing an extremely distinct RNA structural landscape. These results demonstrate that CIRS-seq constitutes an important tool for the identification of native deproteinized RNA structures.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-014-0491-2) contains supplementary material, which is available to authorized users.

Using sm-FRET and denaturants to reveal folding landscapes

RNA folding studies aim to clarify the relationship among sequence, tridimensional structure, and biological function. In the last decade, the application of single-molecule fluorescence resonance energy transfer (sm-FRET) techniques to investigate RNA structure and folding has revealed the details of conformational changes and timescale of the process leading to the formation of biologically active RNA structures with subnanometer resolution on millisecond timescales. In this review, we initially summarize the first wave of single-molecule FRET-based RNA techniques that focused on analyzing the influence of mono- and divalent metal ions on RNA function, and how these studies have provided very valuable information about folding pathways and the presence of intermediate and low-populated states. Next, we describe a second generation of single-molecule techniques that combine sm-FRET with the use of chemical denaturants as an emerging powerful approach to reveal information about the dynamics and energetics of RNA folding that remains hidden using conventional sm-FRET approaches. The main advantages of using the competing interplay between folding agents such as metal ions and denaturants to observe and manipulate the dynamics of RNA folding and RNA-ligand interactions is discussed in the context of the adenine riboswitch aptamer.

Native Purification and Analysis of Long RNAs

The purification and analysis of long noncoding RNAs (lncRNAs) in vitro is a challenge, particularly if one wants to preserve elements of functional structure. Here, we describe a method for purifying lncRNAs that preserves the cotranscriptionally derived structure. The protocol avoids the misfolding that can occur during denaturation-renaturation protocols, thus facilitating the folding of long RNAs to a native-like state. This method is simple and does not require addition of tags to the RNA or the use of affinity columns. LncRNAs purified using this type of native purification protocol are amenable to biochemical and biophysical analysis. Here, we describe how to study lncRNA global compaction in the presence of divalent ions at equilibrium using sedimentation velocity analytical ultracentrifugation and analytical size-exclusion chromatography as well as how to use these uniform RNA species to determine robust lncRNA secondary structure maps by chemical probing techniques like selective 2'-hydroxyl acylation analyzed by primer extension and dimethyl sulfate probing.

Unzippers, resolvers and sensors: a structural and functional biochemistry tale of RNA helicases

The centrality of RNA within the biological world is an irrefutable fact that currently attracts increasing attention from the scientific community. The panoply of functional RNAs requires the existence of specific biological caretakers, RNA helicases, devoted to maintain the proper folding of those molecules, resolving unstable structures. However, evolution has taken advantage of the specific position and characteristics of RNA helicases to develop new functions for these proteins, which are at the interface of the basic processes for transference of information from DNA to proteins. RNA helicases are involved in many biologically relevant processes, not only as RNA chaperones, but also as signal transducers, scaffolds of molecular complexes, and regulatory elements. Structural biology studies during the last decade, founded in X-ray crystallography, have characterized in detail several RNA-helicases. This comprehensive review summarizes the structural knowledge accumulated in the last two decades within this family of proteins, with special emphasis on the structure-function relationships of the most widely-studied families of RNA helicases: the DEAD-box, RIG-I-like and viral NS3 classes.

Folding of SAM-II riboswitch explored by replica-exchange molecular dynamics simulation

Riboswitches are cis-acting RNA fragments that function via a conformational transition mechanism when a specific target molecule binds to its binding pocket, representing an inviting new class of biomolecular target for the development of antibiotics. To understand the folding mechanism of SAM-II riboswitch, occurring predominantly in proteobacteria, a 100ns replica-exchange molecular dynamics simulation in explicit solvent is performed. Our results show that this RNA pseudoknot has multiple folding pathways, and various intermediate structures. The resultant riboswitch conformational transition map is well consistent with the recent fluorescence measurement, which confirms the dynamical properties of this pseudoknot. Moreover, a novel transition pathway is predicted. The global folding dynamics is mainly coupled with the helix rather than the loop region. The potential folding pathways of the riboswitch presented here should lead to a deeper understanding of the folding mechanism of the riboswitch, as well as the conformational change of RNA pseudoknot.

Cross-platform comparison of nucleic acid hybridization: toward quantitative reference standards

Measuring interactions between biological molecules is vitally important to both basic and applied research as well as development of pharmaceuticals. Although a wide and growing range of techniques is available to measure various kinetic and thermodynamic properties of interacting biomolecules, it can be difficult to compare data across techniques of different laboratories and personnel or even across different instruments using the same technique. Here we evaluate relevant biological interactions based on complementary DNA and RNA oligonucleotides that could be used as reference standards for many experimental systems. We measured thermodynamics of duplex formation using isothermal titration calorimetry, differential scanning calorimetry, and ultraviolet-visible (UV-vis) monitored denaturation/renaturation. These standards can be used to validate results, compare data from disparate techniques, act as a teaching tool for laboratory classes, or potentially to calibrate instruments. The RNA and DNA standards have many attractive features, including low cost, high purity, easily measurable concentrations, and minimal handling concerns, making them ideal for use as a reference material.

Independent of their localization in protein the hydrophobic amino acid residues have no effect on the molten globule state of apomyoglobin and the disulfide bond on the surface of apomyoglobin stabilizes this intermediate state

At present it is unclear which interactions in proteins reveal the presence of intermediate states, their stability and formation rate. In this study, we have investigated the effect of substitutions of hydrophobic amino acid residues in the hydrophobic core of protein and on its surface on a molten globule type intermediate state of apomyoglobin. It has been found that independent of their localization in protein, substitutions of hydrophobic amino acid residues do not affect the stability of the molten globule state of apomyoglobin. It has been shown also that introduction of a disulfide bond on the protein surface can stabilize the molten globule state. However in the case of apomyoglobin, stabilization of the intermediate state leads to relative destabilization of the native state of apomyoglobin. The result obtained allows us not only to conclude which mutations can have an effect on the intermediate state of the molten globule type, but also explains why the introduction of a disulfide bond (which seems to “strengthen” the protein) can result in destabilization of the protein native state of apomyoglobin.

Fluorescence tools to investigate riboswitch structural dynamics

Riboswitches are novel regulatory elements that respond to cellular metabolites to control gene expression. They are constituted of highly conserved domains that have evolved to recognize specific metabolites. Such domains, so-called aptamers, are folded into intricate structures to enable metabolite recognition. Over the years, the development of ensemble and single-molecule fluorescence techniques has allowed to probe most of the mechanistic aspects of aptamer folding and ligand binding. In this review, we summarize the current fluorescence toolkit available to study riboswitch structural dynamics. We fist describe those methods based on fluorescent nucleotide analogues, mostly 2-aminopurine (2AP), to investigate short-range conformational changes, including some key steady-state and time-resolved examples that exemplify the versatility of fluorescent analogues as structural probes. The study of long-range structural changes by Förster resonance energy transfer (FRET) is mostly discussed in the context of single-molecule studies, including some recent developments based on the combination of single-molecule FRET techniques with controlled chemical denaturation methods. This article is part of a Special Issue entitled: Riboswitches.

Folding pathways of the Tetrahymena ribozyme

Like many structured RNAs, the Tetrahymena group I intron ribozyme folds through multiple pathways and intermediates. Under standard conditions in vitro, a small fraction reaches the native state (N) with kobs ≈ 0.6 min(-1), while the remainder forms a long-lived misfolded conformation (M) thought to differ in topology. These alternative outcomes reflect a pathway that branches late in folding, after disruption of a trapped intermediate (Itrap). Here we use catalytic activity to probe the folding transitions from Itrap to the native and misfolded states. We show that mutations predicted to weaken the core helix P3 do not increase the rate of folding from Itrap but they increase the fraction that reaches the native state rather than forming the misfolded state. Thus, P3 is disrupted during folding to the native state but not to the misfolded state, and P3 disruption occurs after the rate-limiting step. Interestingly, P3-strengthening mutants also increase native folding. Additional experiments show that these mutants are rapidly committed to folding to the native state, although they reach the native state with approximately the same rate constant as the wild-type ribozyme (~1 min(-1)). Thus, the P3-strengthening mutants populate a distinct pathway that includes at least one intermediate but avoids the M state, most likely because P3 and the correct topology are formed early. Our results highlight multiple pathways in RNA folding and illustrate how kinetic competitions between rapid events can have long-lasting effects because the "choice" is enforced by energy barriers that grow larger as folding progresses.

A chemogenetic approach to study the structural basis of protein-facilitated RNA folding

Large RNA molecules play important roles in all aspects of cellular metabolism ranging from mRNA splicing and protein biosynthesis to regulation of gene expression. In order to correctly perform its function in the cell, an RNA molecule must fold into a complex tertiary structure. Folding of many large RNAs is slow either due to formation of stable misfolded intermediates or due to high contact order or instability of obligate folding intermediates. Therefore many RNAs use protein cofactors to facilitate their folding in vivo. Folding of the yeast mitochondrial group II intron ai5γ to the native state under physiological conditions is facilitated by the protein cofactor Mss116. This chapter describes the use of Nucleotide Analog Interference Mapping (NAIM) to identify specific substructures within the intron molecule that are directly affected by the protein.

Urea-induced denaturation of preQ1-riboswitch

Urea, a polar molecule with a large dipole moment, not only destabilizes folded RNA structures but can also enhance the folding rates of large ribozymes. Unlike the mechanism of urea-induced unfolding of proteins, which is well understood, the action of urea on RNA has barely been explored. We performed extensive all-atom molecular dynamics simulations to determine the molecular underpinnings of urea-induced RNA denaturation. Urea displays its denaturing power in both secondary and tertiary motifs of the riboswitch structure. Our simulations reveal that the denaturation of RNA structures is mainly driven by the hydrogen-bonding and stacking interactions of urea with the bases. Through detailed studies of the simulation trajectories, we found that geminate pairs between urea and bases due to hydrogen bonds and stacks persist only ~0.1-1 ns, which suggests that the urea-base interaction is highly dynamic. Most importantly, the early stage of base-pair disruption is triggered by penetration of water molecules into the hydrophobic domain between the RNA bases. The infiltration of water into the narrow space between base pairs is critical in increasing the accessibility of urea to transiently disrupted bases, thus allowing urea to displace inter-base hydrogen bonds. This mechanism--water-induced disruption of base pairs resulting in the formation of a "wet" destabilized RNA followed by solvation by urea--is the exact opposite of the two-stage denaturation of proteins by urea. In the latter case, initial urea penetration creates a dry globule, which is subsequently solvated by water, leading to global protein unfolding. Our work shows that the ability to interact with both water and polar or nonpolar components of nucleotides makes urea a powerful chemical denaturant for nucleic acids.

Using VIPT-jump to distinguish between different folding mechanisms: application to BBL and a Trpzip

Protein folding involves a large number of sequential molecular steps or conformational substates. Thus, experimental characterization of the underlying folding energy landscape for any given protein is difficult. Herein, we present a new method that can be used to determine the major characteristics of the folding energy landscape in question, e.g., to distinguish between activated and barrierless downhill folding scenarios. This method is based on the idea that the conformational relaxation kinetics of different folding mechanisms at a given final condition will show different dependences on the initial condition. We show, using both simulation and experiment, that it is possible to differentiate between disparate kinetic folding models by comparing temperature jump (T-jump) relaxation traces obtained with a fixed final temperature and varied initial temperatures, which effectively varies the initial potential (VIP) of the system of interest. We apply this method (hereafter refer to as VIPT-jump) to two model systems, tryptophan zipper (Trpzip)-2c and BBL, and our results show that BBL exhibits characteristics of barrierless downhill folding, whereas Trpzip-2c folding encounters a free energy barrier. In addition, using the T-jump data of BBL we are able to provide, via Langevin dynamics simulations, a realistic estimate of its conformational diffusion coefficient.

Single-molecule chemical denaturation of riboswitches

To date, single-molecule RNA science has been developed almost exclusively around the effect of metal ions as folding promoters and stabilizers of the RNA structure. Here, we introduce a novel strategy that combines single-molecule Förster resonance energy transfer (FRET) and chemical denaturation to observe and manipulate RNA dynamics. We demonstrate that the competing interplay between metal ions and denaturant agents provides a platform to extract information that otherwise will remain hidden with current methods. Using the adenine-sensing riboswitch aptamer as a model, we provide strong evidence for a rate-limiting folding step of the aptamer domain being modulated through ligand binding, a feature that is important for regulation of the controlled gene. In the absence of ligand, the rate-determining step is dominated by the formation of long-range key tertiary contacts between peripheral stem-loop elements. In contrast, when the adenine ligand interacts with partially folded messenger RNAs, the aptamer requires specifically bound Mg²⁺ ions, as those observed in the crystal structure, to progress further towards the native form. Moreover, despite that the ligand-free and ligand-bound states are indistinguishable by FRET, their different stability against urea-induced denaturation allowed us to discriminate them, even when they coexist within a single FRET trajectory; a feature not accessible by existing methods.

Impact of urea on the three-dimensional structure, viscoelastic and thermal behavior of iota-carrageenan

Urea breaks hydrogen bonds among biopolymers leading to structural destabilization. In the case of hydrocolloids urea addition is thought to impact gelation. Detailed information about its pertinent role on influencing the structure-function relationships of hydrocolloids is still elusive, however. The present investigation is aimed at delineating hydrocolloids structural behavior in the presence of urea employing iota-carrageenan as a model system. X-ray fiber diffraction, rheological and thermal properties of two iota-carrageenan solutions with weight concentrations 4.5 and 6.0% (w/w) at two urea molar concentrations (0.5 and 2.0 M) with and without heat treatments have been analyzed. X-ray results suggest that the canonical double helical structural arrangement of iota-carrageenan is maintained even after urea addition. However, improved crystallinity, ordering and altered unit cell dimensions especially with heat treatments of the binary mixtures indicate the promotion of favorable interactions among carrageenan helices in the presence of urea. Increased elastic modulus and onset temperature of melting endotherm with the heat treatment compared to cold addition further attests the X-ray observations of enhanced structural ordering. Overall, results suggest that urea molecules synergistically aid iota-carrageenan interactions and stabilize structure of junction zones. Our findings are deemed to be helpful in the design and development of novel non-food applications of hydrocolloids.

Mss116p: a DEAD-box protein facilitates RNA folding

RNA folding is an essential aspect underlying RNA-mediated cellular processes. Many RNAs, including large, multi-domain ribozymes, are capable of folding to the native, functional state without assistance of a protein cofactor in vitro. In the cell, trans-acting factors, such as proteins, are however known to modulate the structure and thus the fate of an RNA. DEAD-box proteins, including Mss116p, were recently found to assist folding of group I and group II introns in vitro and in vivo. The underlying mechanism(s) have been studied extensively to explore the contribution of ATP hydrolysis and duplex unwinding in helicase-stimulated intron splicing. Here we summarize the ongoing efforts to understand the novel role of DEAD-box proteins in RNA folding.

Toward a molecular understanding of RNA remodeling by DEAD-box proteins

DEAD-box proteins are superfamily 2 helicases that function in all aspects of RNA metabolism. They employ ATP binding and hydrolysis to generate tight, yet regulated RNA binding, which is used to unwind short RNA helices non-processively and promote structural transitions of RNA and RNA-protein substrates. In the last few years, substantial progress has been made toward a detailed, quantitative understanding of the structural and biochemical properties of DEAD-box proteins. Concurrently, progress has been made toward a physical understanding of the RNA rearrangements and folding steps that are accelerated by DEAD-box proteins in model systems. Here, we review the recent progress on both of these fronts, focusing on the mitochondrial DEAD-box proteins Mss116 and CYT-19 and their mechanisms in promoting the splicing of group I and group II introns.

The brace for a growing scaffold: Mss116 protein promotes RNA folding by stabilizing an early assembly intermediate

The ai5γ group II intron requires a protein cofactor to facilitate native folding in the cell. Yeast protein Mss116 greatly accelerates intron folding under near-physiological conditions both in vivo and in vitro. Although the effect of Mss116 on the kinetics of ai5γ ribozyme folding and catalysis has been extensively studied, the precise structural role and interaction sites of Mss116 have been elusive. Using Nucleotide Analog Interference Mapping to study the folding of splicing precursor constructs, we have identified specific intron functional groups that participate in Mss116-facilitated folding and we have determined their role in the folding mechanism. The data indicate that Mss116 stabilizes an early, obligate folding intermediate within intron domain 1, thereby laying the foundation for productive folding to the native state. In addition, the data reveal an important role for the IBS2 exon sequence and for the terminus of domain 6, during the folding of self-splicing group IIB intron constructs.

Characterization of the kinetic and thermodynamic landscape of RNA folding using a novel application of isothermal titration calorimetry

A novel isothermal titration calorimetry (ITC) method was applied to investigate RNA helical packing driven by the GAAA tetraloop-receptor interaction in magnesium and potassium solutions. Both the kinetics and thermodynamics were obtained in individual ITC experiments, and analysis of the kinetic data over a range of temperatures provided Arrhenius activation energies (ΔH(‡)) and Eyring transition state entropies (ΔS(‡)). The resulting rich dataset reveals strongly contrasting kinetic and thermodynamic profiles for this RNA folding system when stabilized by potassium versus magnesium. In potassium, association is highly exothermic (ΔH(25°C) = -41.6 ± 1.2 kcal/mol in 150 mM KCl) and the transition state is enthalpically barrierless (ΔH(‡) = -0.6 ± 0.5). These parameters are significantly positively shifted in magnesium (ΔH(25°C) = -20.5 ± 2.1 kcal/mol, ΔH(‡) = 7.3 ± 2.2 kcal/mol in 0.5 mM MgCl(2)). Mixed salt solutions approximating physiological conditions exhibit an intermediate thermodynamic character. The cation-dependent thermodynamic landscape may reflect either a salt-dependent unbound receptor conformation, or alternatively and more generally, it may reflect a small per-cation enthalpic penalty associated with folding-coupled magnesium uptake.

Fast folding of RNA pseudoknots initiated by laser temperature-jump

RNA pseudoknots are examples of minimal structural motifs in RNA with tertiary interactions that stabilize the structures of many ribozymes. They also play an essential role in a variety of biological functions that are modulated by their structure, stability, and dynamics. Therefore, understanding the global principles that determine the thermodynamics and folding pathways of RNA pseudoknots is an important problem in biology, both for elucidating the folding mechanisms of larger ribozymes as well as addressing issues of possible kinetic control of the biological functions of pseudoknots. We report on the folding/unfolding kinetics of a hairpin-type pseudoknot obtained with microsecond time-resolution in response to a laser temperature-jump perturbation. The kinetics are monitored using UV absorbance as well as fluorescence of extrinsically attached labels as spectroscopic probes of the transiently populated RNA conformations. We measure folding times of 1-6 ms at 37 °C, which are at least 100-fold faster than previous observations of very slow folding pseudoknots that were trapped in misfolded conformations. The measured relaxation times are remarkably similar to predictions of a computational study by Thirumalai and co-workers (Cho, S. S.; Pincus, D.L.; Thirumalai, D. Proc. Natl. Acad. Sci. U. S. A. 2009, 106, 17349-17354). Thus, these studies provide the first observation of a fast-folding pseudoknot and present a benchmark against which computational models can be refined.

Effects of a protecting osmolyte on the ion atmosphere surrounding DNA duplexes

Osmolytes are small, chemically diverse, organic solutes that function as an essential component of cellular stress response. Protecting osmolytes enhance protein stability via preferential exclusion, and nonprotecting osmolytes, such as urea, destabilize protein structures. Although much is known about osmolyte effects on proteins, less is understood about osmolyte effects on nucleic acids and their counterion atmospheres. Nonprotecting osmolytes destabilize nucleic acid structures, but effects of protecting osmolytes depend on numerous factors including the type of nucleic acid and the complexity of the functional fold. To begin quantifying protecting osmolyte effects on nucleic acid interactions, we used small-angle X-ray scattering (SAXS) techniques to monitor DNA duplexes in the presence of sucrose. This protecting osmolyte is a commonly used contrast matching agent in SAXS studies of protein-nucleic acid complexes; thus, it is important to characterize interaction changes induced by sucrose. Measurements of interactions between duplexes showed no dependence on the presence of up to 30% sucrose, except under high Mg(2+) conditions where stacking interactions were disfavored. The number of excess ions associated with DNA duplexes, reported by anomalous small-angle X-ray scattering (ASAXS) experiments, was sucrose independent. Although protecting osmolytes can destabilize secondary structures, our results suggest that ion atmospheres of individual duplexes remain unperturbed by sucrose.

Solution structure of RNase P RNA

The ribonucleoprotein enzyme ribonuclease P (RNase P) processes tRNAs by cleavage of precursor-tRNAs. RNase P is a ribozyme: The RNA component catalyzes tRNA maturation in vitro without proteins. Remarkable features of RNase P include multiple turnovers in vivo and ability to process diverse substrates. Structures of the bacterial RNase P, including full-length RNAs and a ternary complex with substrate, have been determined by X-ray crystallography. However, crystal structures of free RNA are significantly different from the ternary complex, and the solution structure of the RNA is unknown. Here, we report solution structures of three phylogenetically distinct bacterial RNase P RNAs from Escherichia coli, Agrobacterium tumefaciens, and Bacillus stearothermophilus, determined using small angle X-ray scattering (SAXS) and selective 2'-hydroxyl acylation analyzed by primer extension (SHAPE) analysis. A combination of homology modeling, normal mode analysis, and molecular dynamics was used to refine the structural models against the empirical data of these RNAs in solution under the high ionic strength required for catalytic activity.

Riboswitches: structures and mechanisms

A critical feature of the hypothesized RNA world would have been the ability to control chemical processes in response to environmental cues. Riboswitches present themselves as viable candidates for a sophisticated mechanism of regulatory control in RNA-based life. These regulatory elements in the modern world are most commonly found in the 5'-untranslated regions of bacterial mRNAs, directly interacting with metabolites as a means of regulating expression of the coding region via a secondary structural switch. In this review, we focus on recent insights into how these RNAs fold into complex architectures capable of both recognizing a specific small molecule compound and exerting regulatory control over downstream sequences, with an emphasis on transcriptional regulation.

The shape-shifting quasispecies of RNA: one sequence, many functional folds

E Unus pluribum, or "Of One, Many", may be at the root of decoding the RNA sequence-structure-function relationship. RNAs embody the large majority of genes in higher eukaryotes and fold in a sequence-directed fashion into three-dimensional structures that perform functions conserved across all cellular life forms, ranging from regulating to executing gene expression. While it is the most important determinant of the RNA structure, the nucleotide sequence is generally not sufficient to specify a unique set of secondary and tertiary interactions due to the highly frustrated nature of RNA folding. This frustration results in folding heterogeneity, a common phenomenon wherein a chemically homogeneous population of RNA molecules folds into multiple stable structures. Often, these alternative conformations constitute misfolds, lacking the biological activity of the natively folded RNA. Intriguingly, a number of RNAs have recently been described as capable of adopting multiple distinct conformations that all perform, or contribute to, the same function. Characteristically, these conformations interconvert slowly on the experimental timescale, suggesting that they should be regarded as distinct native states. We discuss how rugged folding free energy landscapes give rise to multiple native states in the Tetrahymena Group I intron ribozyme, hairpin ribozyme, sarcin-ricin loop, ribosome, and an in vitro selected aptamer. We further describe the varying degrees to which folding heterogeneity impacts function in these RNAs, and compare and contrast this impact with that of heterogeneities found in protein folding. Embracing that one sequence can give rise to multiple native folds, we hypothesize that this phenomenon imparts adaptive advantages on any functionally evolving RNA quasispecies.

Removal of covalent heterogeneity reveals simple folding behavior for P4-P6 RNA

RNA folding landscapes have been described alternately as simple and as complex. The limited diversity of RNA residues and the ability of RNA to form stable secondary structures prior to adoption of a tertiary structure would appear to simplify folding relative to proteins. Nevertheless, there is considerable evidence for long-lived misfolded RNA states, and these observations have suggested rugged energy landscapes. Recently, single molecule fluorescence resonance energy transfer (smFRET) studies have exposed heterogeneity in many RNAs, consistent with deeply furrowed rugged landscapes. We turned to an RNA of intermediate complexity, the P4-P6 domain from the Tetrahymena group I intron, to address basic questions in RNA folding. P4-P6 exhibited long-lived heterogeneity in smFRET experiments, but the inability to observe exchange in the behavior of individual molecules led us to probe whether there was a non-conformational origin to this heterogeneity. We determined that routine protocols in RNA preparation and purification, including UV shadowing and heat annealing, cause covalent modifications that alter folding behavior. By taking measures to avoid these treatments and by purifying away damaged P4-P6 molecules, we obtained a population of P4-P6 that gave near-uniform behavior in single molecule studies. Thus, the folding landscape of P4-P6 lacks multiple deep furrows that would trap different P4-P6 molecules in different conformations and contrasts with the molecular heterogeneity that has been seen in many smFRET studies of structured RNAs. The simplicity of P4-P6 allowed us to reliably determine the thermodynamic and kinetic effects of metal ions on folding and to now begin to build more detailed models for RNA folding behavior.

The osmolyte TMAO stabilizes native RNA tertiary structures in the absence of Mg2+: evidence for a large barrier to folding from phosphate dehydration

The stabilization of RNA tertiary structures by ions is well known, but the neutral osmolyte trimethylamine oxide (TMAO) can also effectively stabilize RNA tertiary structure. To begin to understand the physical basis for the effects of TMAO on RNA, we have quantitated the TMAO-induced stabilization of five RNAs with known structures. So-called m values, the increment in unfolding free energy per molal of osmolyte at constant KCl activity, are ∼0 for a hairpin secondary structure and between 0.70 and 1.85 kcal mol(-1)m(-1) for four RNA tertiary structures (30-86 nt). Further analysis of two RNAs by small-angle X-ray scattering and hydroxyl radical probing shows that TMAO reduces the radius of gyration of the unfolded ensemble to the same endpoint as seen in titration with Mg(2+) and that the structures stabilized by TMAO and Mg(2+) are indistinguishable. Remarkably, TMAO induces the native conformation of a Mg(2+) ion chelation site formed in part by a buried phosphate, even though Mg(2+) is absent. TMAO interacts weakly, if at all, with KCl, ruling out the possibility that TMAO stabilizes RNA indirectly by increasing salt activity. TMAO is, however, strongly excluded from the vicinity of dimethylphosphate (unfavorable interaction free energy, +211 cal mol(-1)m(-1) for the potassium salt), an ion that mimics the RNA backbone phosphate. We suggest that formation of RNA tertiary structure is accompanied by substantial phosphate dehydration (loss of 66-173 water molecules in the RNA structures studied) and that TMAO works principally by reducing the energetic penalty associated with this dehydration. The strong parallels we find between the effects of TMAO and Mg(2+) suggest that RNA sequence is more important than specific ion interactions in specifying the native structure.

Taming free energy landscapes with RNA chaperones

Many non-coding RNAs fold into complex three-dimensional structures, yet the self-assembly of RNA structure is hampered by mispairing, weak tertiary interactions, electrostatic barriers, and the frequent requirement that the 5' and 3' ends of the transcript interact. This rugged free energy landscape for RNA folding means that some RNA molecules in a population rapidly form their native structure, while many others become kinetically trapped in misfolded conformations. Transient binding of RNA chaperone proteins destabilize misfolded intermediates and lower the transition states between conformations, producing a smoother landscape that increases the rate of folding and the probability that a molecule will find the native structure. DEAD-box proteins couple the chemical potential of ATP hydrolysis with repetitive cycles of RNA binding and release, expanding the range of conditions under which they can refold RNA structures.

Activity, folding and Z-DNA formation of the 8-17 DNAzyme in the presence of monovalent ions

The effect of monovalent ions on both the reactivity and global folding of the 8-17 DNAzyme is investigated, and the results are compared with those of the hammerhead ribozyme, which has similar size and secondary structure. In contrast to the hammerhead ribozyme, the 8-17 DNAzyme activity is not detectable in the presence of 4 M K(+), Rb(+), or Cs(+) or in the presence of 80 mM, [Co(NH(3))(6)](3+). Only 4 M Li(+), NH(4)(+) and, to a lesser extent, Na(+) conferred detectable activity. The observed rate constants (k(obs) approximately 10(-3) min(-1) for Li(+) and NH(4)(+)) are approximately 1000-fold lower than that in the presence of 10 mM Mg(2+), and approximately 200,000-fold slower than that in the presence of 100 microM Pb(2+). Since the hammerhead ribozyme displays monovalent ion-dependent activity that is often within approximately 10-fold of divalent metal ion-dependent activity, these results suggest that the 8-17 DNAzyme, obtained by in vitro selections, has evolved to have a more stringent divalent metal ion requirement for high activity as compared to the naturally occurring ribozymes, making the 8-17 DNAzyme an excellent choice as a Pb(2+) sensor with high selectivity. In contrast to the activity data, folding was observed in the presence of all the monovalent ions investigated, although those monovalent ions that do not support DNAzyme activity have weaker binding affinity (K(d) approximately 0.35 M for Rb(+) and Cs(+)), while those that confer DNAzyme activity possess stronger affinity (K(d) approximately 0.22 M for Li(+), Na(+) and NH(4)(+)). In addition, a correlation between metal ion charge density, binding affinity and enzyme activity was found among mono- and divalent metal ions except Pb(2+); higher charge density resulted in stronger affinity and higher activity, suggesting that the observed folding and activity is at least partially due to electrostatic interactions between ions and the DNAzyme. Finally, circular dichroism (CD) study has revealed Z-DNA formation with the monovalent metal ions, Zn(2+) and Mg(2+); the K(d) values obtained using CD were in the same range as those obtained from folding studies using FRET. However, Z-DNA formation was not observed with Pb(2+). These results indicate that Pb(2+)-dependent function follows a different mechanism from the monovalent metal ions and other divalent metal ions; in the presence of latter metal ions, metal-ion dependent folding and structural changes, including formation of Z-DNA, play an important role in the catalytic function of the 8-17 DNAzyme.

Compact intermediates in RNA folding

Large noncoding RNAs fold into their biologically functional structures via compact yet disordered intermediates, which couple the stable secondary structure of the RNA with the emerging tertiary fold. The specificity of the collapse transition, which coincides with the assembly of helical domains, depends on RNA sequence and counterions. It determines the specificity of the folding pathways and the magnitude of the free energy barriers to the ensuing search for the native conformation. By coupling helix assembly with nascent tertiary interactions, compact folding intermediates in RNA also play a crucial role in ligand binding and RNA-protein recognition.

Full-length hepatitis B virus core protein packages viral and heterologous RNA with similarly high levels of cooperativity

A critical feature of a viral life cycle is the ability to selectively package the viral genome. In vivo, phosphorylated hepatitis B virus (HBV) core protein specifically encapsidates a complex of pregenomic RNA (pgRNA) and viral polymerase; it has been suggested that packaging is specific for the complex. Here, we test the hypothesis that core protein has intrinsic specificity for pgRNA, independent of the polymerase. For these studies, we also evaluated the effect of core protein phosphorylation on assembly and RNA binding, using phosphorylated core protein and a phosphorylation mimic in which S155, S162, and S170 were mutated to glutamic acid. We have developed an in vitro system where capsids are disassembled and assembly-active core protein dimer is purified. With this protein, we have reassembled empty capsids and RNA-filled capsids. We found that core protein dimer bound and encapsidated both the HBV pregenomic RNA and heterologous RNA with high levels of cooperativity, irrespective of phosphorylation. In direct competition assays, no specificity for pregenomic RNA was observed. This suggests that another factor, such as the viral polymerase, is required for specific packaging. These results also beg the question of what prevents HBV core protein from assembling on nonviral RNA, preserving the protein for virus production.

Multiple unfolding events during native folding of the Tetrahymena group I ribozyme

Despite the ubiquitous nature of misfolded intermediates in RNA folding, little is known about their physical properties or the folding transitions that allow them to continue folding productively. Folding of the Tetrahymena group I ribozyme includes sequential accumulation of two intermediates, termed I(trap) and misfolded (M). Here, we probe the structure and folding transition of I(trap) and compare them to those of M. Hydroxyl radical and dimethyl sulfate footprinting show that both I(trap) and M are extensively structured and crudely resemble the native RNA. However, regions of the core P3-P8 domain are more exposed to solvent in I(trap) than in M. I(trap) rearranges to continue folding nearly 1000-fold faster than M, and urea accelerates folding of I(trap) much less than M. Thus, the rate-limiting transition from I(trap) requires a smaller increase in exposed surface. Mutations that disrupt peripheral tertiary contacts give large and nearly uniform increases in re-folding of M, whereas the same mutations give at most modest increases in folding from I(trap). Intriguingly, mutations within the peripheral element P5abc give 5- to 10-fold accelerations in escape from I(trap), whereas ablation of P13, which lies on the opposite surface in the native structure, near the P3-P8 domain, has no effect. Thus, the unfolding required from I(trap) appears to be local, whereas the unfolding of M appears to be global. Further, the modest effects from several mutations suggest that there are multiple pathways for escape from I(trap) and that escape is aided by loosening nearby native structural constraints, presumably to facilitate local movements of nucleotides or segments that have not formed native contacts. Overall, these and prior results suggest a model in which the global architecture and peripheral interactions of the RNA are achieved relatively early in folding. Multiple folding and re-folding events occur on the predominant pathway to the native state, with increasing native core interactions and cooperativity as folding progresses.

Protein-precursor tRNA contact leads to sequence-specific recognition of 5' leaders by bacterial ribonuclease P

Bacterial ribonuclease P (RNase P) catalyzes the cleavage of 5' leader sequences from precursor tRNAs (pre-tRNAs). Previously, all known substrate nucleotide specificities in this system are derived from RNA-RNA interactions with the RNase P RNA subunit. Here, we demonstrate that pre-tRNA binding affinities for Bacillus subtilis and Escherichia coli RNase P are enhanced by sequence-specific contacts between the fourth pre-tRNA nucleotide on the 5' side of the cleavage site (N(-4)) and the RNase P protein (P protein) subunit. B. subtilis RNase P has a higher affinity for pre-tRNA with adenosine at N(-4), and this binding preference is amplified at physiological divalent ion concentrations. Measurements of pre-tRNA-containing adenosine analogs at N(-4) indicate that specificity arises from a combination of hydrogen bonding to the N6 exocyclic amine of adenosine and steric exclusion of the N2 amine of guanosine. Mutagenesis of B. subtilis P protein indicates that F20 and Y34 contribute to selectivity at N(-4). The hydroxyl group of Y34 enhances selectivity, likely by forming a hydrogen bond with the N(-4) nucleotide. The sequence preference of E. coli RNase P is diminished, showing a weak preference for adenosine and cytosine at N(-4), consistent with the substitution of Leu for Y34 in the E. coli P protein. This is the first identification of a sequence-specific contact between P protein and pre-tRNA that contributes to molecular recognition of RNase P. Additionally, sequence analyses reveal that a greater-than-expected fraction of pre-tRNAs from both E. coli and B. subtilis contains a nucleotide at N(-4) that enhances RNase P affinity. This observation suggests that specificity at N(-4) contributes to substrate recognition in vivo. Furthermore, bioinformatic analyses suggest that sequence-specific contacts between the protein subunit and the leader sequences of pre-tRNAs may be common in bacterial RNase P and may lead to species-specific substrate recognition.

Ion-RNA interactions thermodynamic analysis of the effects of mono- and divalent ions on RNA conformational equilibria

RNA secondary and tertiary structures are strongly stabilized by added salts, and a quantitative thermodynamic analysis of the relevant ion-RNA interactions is an important aspect of the RNA folding problem. Because of long-range electrostatic forces, an RNA perturbs the distribution of both cations and anions throughout a large volume. Binding formalisms that require a distinction between "bound" and "free" ions become problematic in such situations. A more fundamental thermodynamic framework is developed here, based on preferential interaction coefficients; linkage equations derived from this framework provide a model-free description of the "uptake" or "release" of cations and anions that accompany an RNA conformational transition. Formulas appropriate for analyzing the dependence of RNA stability on either mono- or divalent salt concentration are presented and their application to experimental data is illustrated. Two example datasets are analyzed with respect to the monovalent salt dependence of tertiary structure formation in different RNAs, and three different experimental methods for quantitating the "uptake" of Mg(2+) ions are applied to the folding of a riboswitch RNA. Advantages and limitations of each method are discussed.

Catalytic activity as a probe of native RNA folding

As RNAs fold to functional structures, they traverse complex energy landscapes that include many partially folded and misfolded intermediates. For structured RNAs that possess catalytic activity, this activity can provide a powerful means of monitoring folding that is complementary to biophysical approaches. RNA catalysis can be used to track accumulation of the native RNA specifically and quantitatively, readily distinguishing the native structure from intermediates that resemble it and may not be differentiated by other approaches. Here, we outline how to design and interpret experiments using catalytic activity to monitor RNA folding, and we summarize adaptations of the method that have been used to probe aspects of folding well beyond determination of the folding rates.

Improvement of RNA secondary structure prediction using RNase H cleavage and randomized oligonucleotides

RNA secondary structure prediction using free energy minimization is one method to gain an approximation of structure. Constraints generated by enzymatic mapping or chemical modification can improve the accuracy of secondary structure prediction. We report a facile method that identifies single-stranded regions in RNA using short, randomized DNA oligonucleotides and RNase H cleavage. These regions are then used as constraints in secondary structure prediction. This method was used to improve the secondary structure prediction of Escherichia coli 5S rRNA. The lowest free energy structure without constraints has only 27% of the base pairs present in the phylogenetic structure. The addition of constraints from RNase H cleavage improves the prediction to 100% of base pairs. The same method was used to generate secondary structure constraints for yeast tRNA^Phe, which is accurately predicted in the absence of constraints (95%). Although RNase H mapping does not improve secondary structure prediction, it does eliminate all other suboptimal structures predicted within 10% of the lowest free energy structure. The method is advantageous over other single-stranded nucleases since RNase H is functional in physiological conditions. Moreover, it can be used for any RNA to identify accessible binding sites for oligonucleotides or small molecules.

Energy barriers, pathways, and dynamics during folding of large, multidomain RNAs

Large, multidomain RNA molecules are generally thought to fold following multiple pathways down rugged landscapes populated with intermediates and traps. A challenge to understanding RNA folding reactions is the complex relationships that exist between the structure of the RNA and its folding landscape. The identification of intermediate species that populate folding landscapes and characterization of elements of their structures are the key components to solving the RNA folding problem. This review explores recent studies that characterize the dominant pathways by which RNA folds, structural and dynamic features of intermediates that populate the folding landscape, and the energy barriers that separate the distinct steps of the folding process.

Reliable protein folding on complex energy landscapes: the free energy reaction path

A theoretical framework is developed to study the dynamics of protein folding. The key insight is that the search for the native protein conformation is influenced by the rate r at which external parameters, such as temperature, chemical denaturant, or pH, are adjusted to induce folding. A theory based on this insight predicts that 1), proteins with complex energy landscapes can fold reliably to their native state; 2), reliable folding can occur as an equilibrium or out-of-equilibrium process; and 3), reliable folding only occurs when the rate r is below a limiting value, which can be calculated from measurements of the free energy. We test these predictions against numerical simulations of model proteins with a single energy scale.