Kinetic intermediates trapped by native interactions in RNA folding

Domain consolidation in Bacterial 50S assembly revealed by Anti-Sense Oligonucleotide Probing

Investigating the intricate and rapid folding kinetics of large RNA-protein complexes (RNPs), like the bacterial ribosome, remains a formidable challenge in structural biology. Previous genetic approaches to probe assembly have focused on modulating the expression of either r-proteins or assembly factors. Here, anti-sense oligonucleotides (ASOs) were used to disrupt native RNA/RNA and RNA/protein interactions, in order to generate novel folding intermediates. In an in vitro co-transcriptional assembly assay, 8 assembly inhibitor ASOs were identified. Using cryo-electron microscopy, 38 new intermediate structures were determined resulting from the specific inhibitions. In particular a novel intermediate class provided compelling evidence of independent rRNA domain folding before proper interdomain docking. Three PNAs targeting domain-I of 23S-rRNA further subdivided the previously identified assembly core into smaller blocks representing the earliest steps in assembly. The resulting assembly graph reveals template-directed RNA foldon docking and domain consolidation, which provides a hierarchical view of the RNP assembly process.

Overcoming Immunotherapeutic Resistance in PDAC: SIRPα-CD47 blockade

A daily increase in the number of new cases of pancreatic ductal adenocarcinoma remains an issue of contention in cancer research. The data revealed that a global cumulated case of about 500, 000 have been reported. This has made PDAC the fourteenth most occurring tumor case in cancer research. Furthermore, PDAC is responsible for about 466,003 deaths annually, representing the seventh prevalent type of cancer mortality. PDAC has no salient symptoms in its early stages. This has exasperated several attempts to produce a perfect therapeutic agent against PDAC. Recently, immunotherapeutic research has shifted focus to the blockade of checkpoint proteins in the management and of some cancers. Investigations have centrally focused on developing therapeutic agents that could at least to a significant extent block the SIRPα-CD47 signaling cascade (a cascade which prevent phagocytosis of tumors by dendritic cells, via the deactivation of innate immunity and subsequently resulting in tumor regression) with minimal side effects. The concept on the blockade of this interaction as a possible mechanism for inhibiting the progression of PDAC is currently being debated. This review examined the structure--function activity of SIRPα-CD47 interaction while discussing in detail the mechanism of tumor resistance in PDAC. Further, this review details how the blockade of SIRPα-CD47 interaction serve as a therapeutic option in the management of PDAC.

Iterative annealing mechanism explains the functions of the GroEL and RNA chaperones

Molecular chaperones are ATP-consuming machines, which facilitate the folding of proteins and RNA molecules that are kinetically trapped in misfolded states. Unassisted folding occurs by the kinetic partitioning mechanism according to which folding to the native state, with low probability as well as misfolding to one of the many metastable states, with high probability, occur rapidly. GroEL is an all-purpose stochastic machine that assists misfolded substrate proteins to fold. The RNA chaperones such as CYT-19, which are ATP-consuming enzymes, help the folding of ribozymes that get trapped in metastable states for long times. GroEL does not interact with the folded proteins but CYT-19 disrupts both the folded and misfolded ribozymes. The structures of GroEL and RNA chaperones are strikingly different. Despite these differences, the iterative annealing mechanism (IAM) quantitatively explains all the available experimental data for assisted folding of proteins and ribozymes. Driven by ATP binding and hydrolysis and GroES binding, GroEL undergoes a catalytic cycle during which it samples three allosteric states, T (apo), R (ATP bound), and R^″ (ADP bound). Analyses of the experimental data show that the efficiency of the GroEL-GroES machinery and mutants is determined by the resetting rate k _{R ″  → T} , which is largest for the wild-type (WT) GroEL. Generalized IAM accurately predicts the folding kinetics of Tetrahymena ribozyme and its variants. Chaperones maximize the product of the folding rate and the steady-state native state fold by driving the substrates out of equilibrium. Neither the absolute yield nor the folding rate is optimized.

Effects of flanking regions on HDV cotranscriptional folding kinetics

Hepatitis delta virus (HDV) ribozyme performs the self-cleavage activity through folding to a double pseudoknot structure. The folding of functional RNA structures is often coupled with the transcription process. In this work, we developed a new approach for predicting the cotranscriptional folding kinetics of RNA secondary structures with pseudoknots. We theoretically studied the cotranscriptional folding behavior of the 99-nucleotide (nt) HDV sequence, two upstream flanking sequences, and one downstream flanking sequence. During transcription, the 99-nt HDV can effectively avoid the trap intermediates and quickly fold to the cleavage-active state. It is different from its refolding kinetics, which folds into an intermediate trap state. For all the sequences, the ribozyme regions (from 1 to 73) all fold to the same structure during transcription. However, the existence of the 30-nt upstream flanking sequence can inhibit the ribozyme region folding into the active native state through forming an alternative helix Alt1 with the segments 70-90. The longer upstream flanking sequence of 54 nt itself forms a stable hairpin structure, which sequesters the formation of the Alt1 helix and leads to rapid formation of the cleavage-active structure. Although the 55-nt downstream flanking sequence could invade the already folded active structure during transcription by forming a more stable helix with the ribozyme region, the slow transition rate could keep the structure in the cleavage-active structure to perform the activity.

Group I Intron-Based Therapeutics Through Trans-Splicing Reaction

In 1982, the Cech group discovered that an intron structure in an rRNA precursor of Tetrahymena thermophila is sufficient to complete splicing without assistance from proteins. This was the first moment that scientists recognized RNAs can have catalytic activities derived from their own unique three-dimensional structures and thus play more various roles in biological processes than thought before. Several additional catalytic RNAs, called ribozymes, were subsequently identified in nature followed by intense studies to reveal their mechanisms of action and to engineer them for use in fields such as molecular cell biology, therapeutics, imaging, etc. Naturally occurring RNA-targeting ribozymes can be broadly classified into two categories by their abilities: Self-cleavage and self-splicing. Since ribozymes use base-pairing to recognize cleavage sites, identification of the catalytic center of naturally occurring ribozymes enables to engineer from "self" to "trans" acting ones which has accelerated to design and use ribozyme as valuable tools in gene therapy fields. Especially, group I intron-based trans-splicing ribozyme has unique property to use as a gene therapeutic agent. It can destroy and simultaneously repair (and/or reprogram) target RNAs to yield the desired therapeutic RNAs, maintaining endogenous spatial and temporal gene regulation of target RNAs. There have been progressive improvements in trans-splicing ribozymes and successful applications of these elements in gene therapy and molecular imaging approaches for various pathogenic conditions. In this chapter, current status of trans-splicing ribozyme therapeutics, focusing on Tetrahymena group I intron-based ribozymes, and their future prospects will be discussed.

Therapeutic applications of group I intron-based trans-splicing ribozymes

Since the breakthrough discovery of catalytic RNAs (ribozymes) in the early 1980s, valuable ribozyme-based gene therapies have been developed for incurable diseases ranging from genetic disorders to viral infections and cancers. Ribozymes can be engineered and used to downregulate or repair pathogenic genes via RNA cleavage mediated by trans-cleaving ribozymes or repair and reprograming mediated by trans-splicing ribozymes, respectively. Uniquely, trans-splicing ribozymes can edit target RNAs via simultaneous destruction and repair (and/or reprograming) to yield the desired therapeutic RNAs, thus selectively inducing therapeutic gene activity in cells expressing the target RNAs. In contrast to traditional gene therapy approaches, such as simple addition of therapeutic transgenes or inhibition of disease-causing genes, the selective repair and/or reprograming abilities of trans-splicing ribozymes in target RNA-expressing cells facilitates the maintenance of endogenous spatial and temporal gene regulation and reduction of disease-associated transcript expression. In molecular imaging technologies, trans-splicing ribozymes can be used to reprogram specific RNAs in living cells and organisms by the 3'-tagging of reporter RNAs. The past two decades have seen progressive improvements in trans-splicing ribozymes and the successful application of these elements in gene therapy and molecular imaging approaches for various pathogenic conditions, such as genetic, infectious, and malignant disease. This review provides an overview of the current status of trans-splicing ribozyme therapeutics, focusing on Tetrahymena group I intron-based ribozymes, and their future prospects. This article is categorized under: RNA in Disease and Development > RNA in Disease.

Pausing guides RNA folding to populate transiently stable RNA structures for riboswitch-based transcription regulation

In bacteria, the regulation of gene expression by cis-acting transcriptional riboswitches located in the 5'-untranslated regions of messenger RNA requires the temporal synchronization of RNA synthesis and ligand binding-dependent conformational refolding. Ligand binding to the aptamer domain of the riboswitch induces premature termination of the mRNA synthesis of ligand-associated genes due to the coupled formation of 3'-structural elements acting as terminators. To date, there has been no high resolution structural description of the concerted process of synthesis and ligand-induced restructuring of the regulatory RNA element. Here, we show that for the guanine-sensing xpt-pbuX riboswitch from Bacillus subtilis, the conformation of the full-length transcripts is static: it exclusively populates the functional off-state but cannot switch to the on-state, regardless of the presence or absence of ligand. We show that only the combined matching of transcription rates and ligand binding enables transcription intermediates to undergo ligand-dependent conformational refolding.

DOI:http://dx.doi.org/10.7554/eLife.21297.001

Design and Experimental Evolution of trans-Splicing Group I Intron Ribozymes

Group I intron ribozymes occur naturally as cis-splicing ribozymes, in the form of introns that do not require the spliceosome for their removal. Instead, they catalyze two consecutive trans-phosphorylation reactions to remove themselves from a primary transcript, and join the two flanking exons. Designed, trans-splicing variants of these ribozymes replace the 3′-portion of a substrate with the ribozyme’s 3′-exon, replace the 5′-portion with the ribozyme’s 5′-exon, or insert/remove an internal sequence of the substrate. Two of these designs have been evolved experimentally in cells, leading to variants of group I intron ribozymes that splice more efficiently, recruit a cellular protein to modify the substrate’s gene expression, or elucidate evolutionary pathways of ribozymes in cells. Some of the artificial, trans-splicing ribozymes are promising as tools in therapy, and as model systems for RNA evolution in cells. This review provides an overview of the different types of trans-splicing group I intron ribozymes that have been generated, and the experimental evolution systems that have been used to improve them.

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.

Bridging the gap between in vitro and in vivo RNA folding

Deciphering the folding pathways and predicting the structures of complex three-dimensional biomolecules is central to elucidating biological function. RNA is single-stranded, which gives it the freedom to fold into complex secondary and tertiary structures. These structures endow RNA with the ability to perform complex chemistries and functions ranging from enzymatic activity to gene regulation. Given that RNA is involved in many essential cellular processes, it is critical to understand how it folds and functions in vivo. Within the last few years, methods have been developed to probe RNA structures in vivo and genome-wide. These studies reveal that RNA often adopts very different structures in vivo and in vitro, and provide profound insights into RNA biology. Nonetheless, both in vitro and in vivo approaches have limitations: studies in the complex and uncontrolled cellular environment make it difficult to obtain insight into RNA folding pathways and thermodynamics, and studies in vitro often lack direct cellular relevance, leaving a gap in our knowledge of RNA folding in vivo. This gap is being bridged by biophysical and mechanistic studies of RNA structure and function under conditions that mimic the cellular environment. To date, most artificial cytoplasms have used various polymers as molecular crowding agents and a series of small molecules as cosolutes. Studies under such in vivo-like conditions are yielding fresh insights, such as cooperative folding of functional RNAs and increased activity of ribozymes. These observations are accounted for in part by molecular crowding effects and interactions with other molecules. In this review, we report milestones in RNA folding in vitro and in vivo and discuss ongoing experimental and computational efforts to bridge the gap between these two conditions in order to understand how RNA folds in the cell.

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.

Secondary structure encodes a cooperative tertiary folding funnel in the Azoarcus ribozyme

A requirement for specific RNA folding is that the free-energy landscape discriminate against non-native folds. While tertiary interactions are critical for stabilizing the native fold, they are relatively non-specific, suggesting additional mechanisms contribute to tertiary folding specificity. In this study, we use coarse-grained molecular dynamics simulations to explore how secondary structure shapes the tertiary free-energy landscape of the Azoarcus ribozyme. We show that steric and connectivity constraints posed by secondary structure strongly limit the accessible conformational space of the ribozyme, and that these so-called topological constraints in turn pose strong free-energy penalties on forming different tertiary contacts. Notably, native A-minor and base-triple interactions form with low conformational free energy, while non-native tetraloop/tetraloop–receptor interactions are penalized by high conformational free energies. Topological constraints also give rise to strong cooperativity between distal tertiary interactions, quantitatively matching prior experimental measurements. The specificity of the folding landscape is further enhanced as tertiary contacts place additional constraints on the conformational space, progressively funneling the molecule to the native state. These results indicate that secondary structure assists the ribozyme in navigating the otherwise rugged tertiary folding landscape, and further emphasize topological constraints as a key force in RNA folding.

How do metal ions direct ribozyme folding?

Ribozymes, which carry out phosphoryl-transfer reactions, often require Mg(2+) ions for catalytic activity. The correct folding of the active site and ribozyme tertiary structure is also regulated by metal ions in a manner that is not fully understood. Here we employ coarse-grained molecular simulations to show that individual structural elements of the group I ribozyme from the bacterium Azoarcus form spontaneously in the unfolded ribozyme even at very low Mg(2+) concentrations, and are transiently stabilized by the coordination of Mg(2+) ions to specific nucleotides. However, competition for scarce Mg(2+) and topological constraints that arise from chain connectivity prevent the complete folding of the ribozyme. A much higher Mg(2+) concentration is required for complete folding of the ribozyme and stabilization of the active site. When Mg(2+) is replaced by Ca(2+) the ribozyme folds, but the active site remains unstable. Our results suggest that group I ribozymes utilize the same interactions with specific metal ligands for both structural stability and chemical activity.

The Globular State of the Single-Stranded RNA: Effect of the Secondary Structure Rearrangements

The mutual influence of the slow rearrangements of secondary structure and fast collapse of the long single-stranded RNA (ssRNA) in approximation of coarse-grained model is studied with analytic calculations. It is assumed that the characteristic time of the secondary structure rearrangement is much longer than that for the formation of the tertiary structure. A nonequilibrium phase transition of the 2nd order has been observed.

Formation of Tertiary Interactions during rRNA GTPase Center Folding

The 60-nt GTPase center (GAC) of 23S rRNA has a phylogenetically conserved secondary structure with two hairpin loops and a 3-way junction. It folds into an intricate tertiary structure upon addition of Mg(2+) ions, which is stabilized by the L11 protein in cocrystal structures. Here, we monitor the kinetics of its tertiary folding and Mg(2+)-dependent intermediate states by observing selected nucleobases that contribute specific interactions to the GAC tertiary structure in the cocrystals. The fluorescent nucleobase 2-aminopurine replaced three individual adenines, two of which make long-range stacking interactions and one that also forms hydrogen bonds. Each site reveals a unique response to Mg(2+) addition and temperature, reflecting its environmental change from secondary to tertiary structure. Stopped-flow fluorescence experiments revealed that kinetics of tertiary structure formation upon addition of MgCl2 are also site specific, with local conformational changes occurring from 5 ms to 4s and with global folding from 1 to 5s. Site-specific substitution with (15)N-nucleobases allowed observation of stable hydrogen bond formation by NMR experiments. Equilibrium titration experiments indicate that a stable folding intermediate is present at stoichiometric concentrations of Mg(2+) and suggest that there are two initial sites of Mg(2+) ion association.

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.

Characterization of conformational dynamics of bistable RNA by equilibrium and non-equilibrium NMR

Unlike proteins, a given RNA sequence can adopt more than a single conformation. The two (or more) conformations are long-lived and have similar stabilities, but interconvert only slowly. Such bi- or multistability is often linked to the biological functions of the RNA. This unit describes how nuclear magnetic resonance (NMR) spectroscopy can be used to characterize the conformational dynamics of bistable RNAs.

DEAD-box helicase proteins disrupt RNA tertiary structure through helix capture

Single-molecule fluorescence experiments reveal how DEAD-box proteins unfold structured RNAs to promote conformational transitions and refolding to the native functional state.

DEAD-box helicase proteins accelerate folding and rearrangements of highly structured RNAs and RNA–protein complexes (RNPs) in many essential cellular processes. Although DEAD-box proteins have been shown to use ATP to unwind short RNA helices, it is not known how they disrupt RNA tertiary structure. Here, we use single molecule fluorescence to show that the DEAD-box protein CYT-19 disrupts tertiary structure in a group I intron using a helix capture mechanism. CYT-19 binds to a helix within the structured RNA only after the helix spontaneously loses its tertiary contacts, and then CYT-19 uses ATP to unwind the helix, liberating the product strands. Ded1, a multifunctional yeast DEAD-box protein, gives analogous results with small but reproducible differences that may reflect its in vivo roles. The requirement for spontaneous dynamics likely targets DEAD-box proteins toward less stable RNA structures, which are likely to experience greater dynamic fluctuations, and provides a satisfying explanation for previous correlations between RNA stability and CYT-19 unfolding efficiency. Biologically, the ability to sense RNA stability probably biases DEAD-box proteins to act preferentially on less stable misfolded structures and thereby to promote native folding while minimizing spurious interactions with stable, natively folded RNAs. In addition, this straightforward mechanism for RNA remodeling does not require any specific structural environment of the helicase core and is likely to be relevant for DEAD-box proteins that promote RNA rearrangements of RNP complexes including the spliceosome and ribosome.

In addition to carrying genetic information from DNA to protein, RNAs function in many essential cellular processes. This often requires the RNA to form a specific three-dimensional structure, and some functions require cycling between multiple structures. However, RNAs have a strong propensity to become trapped in nonfunctional, misfolded structures. Due to the intrinsic stability of local structure for RNA, these misfolded species can be long-lived and therefore accumulate. ATP-dependent RNA chaperone proteins called DEAD-box proteins are known to promote native RNA folding by disrupting misfolded structures. Although these proteins can unwind short RNA helices, the mechanism by which they act upon higher order tertiary contacts is unknown. Our current work shows that DEAD-box proteins capture transiently exposed RNA helices, preventing any tertiary contacts from reforming and potentially destabilizing the global RNA architecture. Helix unwinding by the DEAD-box protein then allows the product RNA strands to form new contacts. This helix capture mechanism for manipulation of RNA tertiary structure does not require a specific binding motif or structural environment and may be general for DEAD-box helicase proteins that act on structured RNAs.

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.

Mapping the kinetic barriers of a Large RNA molecule's folding landscape

The folding of linear polymers into discrete three-dimensional structures is often required for biological function. The formation of long-lived intermediates is a hallmark of the folding of large RNA molecules due to the ruggedness of their energy landscapes. The precise thermodynamic nature of the barriers (whether enthalpic or entropic) that leads to intermediate formation is still poorly characterized in large structured RNA molecules. A classic approach to analyzing kinetic barriers are temperature dependent studies analyzed with Eyring's transition state theory. We applied Eyring's theory to time-resolved hydroxyl radical (•OH) footprinting kinetics progress curves collected at eight temperature from 21.5 °C to 51 °C to characterize the thermodynamic nature of folding intermediate formation for the Mg(2+)-mediated folding of the Tetrahymena thermophila group I ribozyme. A common kinetic model configuration describes this RNA folding reaction over the entire temperature range studied consisting of primary (fast) transitions to misfolded intermediates followed by much slower secondary transitions, consistent with previous studies. Eyring analysis reveals that the primary transitions are moderate in magnitude and primarily enthalpic in nature. In contrast, the secondary transitions are daunting in magnitude and entropic in nature. The entropic character of the secondary transitions is consistent with structural rearrangement of the intermediate species to the final folded form. This segregation of kinetic control reveals distinctly different molecular mechanisms during the two stages of RNA folding and documents the importance of entropic barriers to defining rugged RNA folding landscapes.

RNA refolding studied by light-coupled NMR spectroscopy

Conformational transitions (refolding) between long-lived conformational states constitute the time-limiting step during the folding process of large RNAs. As the dynamics of these reactions dominate the regulatory and other functional behavior of RNA molecules, it is of importance to characterize them with high spatial and temporal resolution. Here, we describe a method for the investigation of RNA refolding reactions based on the photolytic generation of preselected conformations in a non-equilibrium state, followed by the observation of the folding trajectory with real-time NMR spectroscopy.

The long-range P3 helix of the Tetrahymena ribozyme is disrupted during folding between the native and misfolded conformations

RNAs are prone to misfolding, but how misfolded structures are formed and resolved remains incompletely understood. The Tetrahymena group I intron ribozyme folds in vitro to a long-lived misfolded conformation (M) that includes extensive native structure but is proposed to differ in topology from the native state (N). A leading model predicts that exchange of the topologies requires unwinding of the long-range, core helix P3, despite the presence of P3 in both conformations. To test this model, we constructed 16 mutations to strengthen or weaken P3. Catalytic activity and in-line probing showed that nearly all of the mutants form the M state before folding to N. The P3-weakening mutations accelerated refolding from M (3- to 30-fold) and the P3-strengthening mutations slowed refolding (6- to 1400-fold), suggesting that P3 indeed unwinds transiently. Upon depletion of Mg(2+), the mutations had analogous effects on unfolding from N to intermediates that subsequently fold to M. The magnitudes for the P3-weakening mutations were larger than in refolding from M, and small-angle X-ray scattering showed that the ribozyme expands rapidly to intermediates from which P3 is disrupted subsequently. These results are consistent with previous results indicating unfolding of native peripheral structure during refolding from M, which probably permits rearrangement of the core. Together, our results demonstrate that exchange of the native and misfolded conformations requires loss of a core helix in addition to peripheral structure. Further, the results strongly suggest that misfolding arises from a topological error within the ribozyme core, and a specific topology is proposed.

Understanding the role of three-dimensional topology in determining the folding intermediates of group I introns

Many RNA molecules exert their biological function only after folding to unique three-dimensional structures. For long, noncoding RNA molecules, the complexity of finding the native topology can be a major impediment to correct folding to the biologically active structure. An RNA molecule may fold to a near-native structure but not be able to continue to the correct structure due to a topological barrier such as crossed strands or incorrectly stacked helices. Achieving the native conformation thus requires unfolding and refolding, resulting in a long-lived intermediate. We investigate the role of topology in the folding of two phylogenetically related catalytic group I introns, the Twort and Azoarcus group I ribozymes. The kinetic models describing the Mg(2+)-mediated folding of these ribozymes were previously determined by time-resolved hydroxyl (∙OH) radical footprinting. Two intermediates formed by parallel intermediates were resolved for each RNA. These data and analytical ultracentrifugation compaction analyses are used herein to constrain coarse-grained models of these folding intermediates as we investigate the role of nonnative topology in dictating the lifetime of the intermediates. Starting from an ensemble of unfolded conformations, we folded the RNA molecules by progressively adding native constraints to subdomains of the RNA defined by the ∙OH time-progress curves to simulate folding through the different kinetic pathways. We find that nonnative topologies (arrangement of helices) occur frequently in the folding simulations despite using only native constraints to drive the reaction, and that the initial conformation, rather than the folding pathway, is the major determinant of whether the RNA adopts nonnative topology during folding. From these analyses we conclude that biases in the initial conformation likely determine the relative flux through parallel RNA folding pathways.

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.

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.

Kinetic mechanism of conformational switch between bistable RNA hairpins

Transitions between the different conformational states play a critical role in many RNA catalytic and regulatory functions. In this study, we use the Kinetic Monte Carlo method to investigate the kinetic mechanism for the conformational switches between bistable RNA hairpins. We find three types of conformational switch pathways for RNA hairpins: refolding after complete unfolding, folding through basepair-exchange pathways and through pseudoknot-assisted pathways, respectively. The result of the competition between the three types of pathways depends mainly on the location of the rate-limiting base stacks (such as the GC base stacks) in the structures. Depending on the structural relationships between the two bistable hairpins, the conformational switch can follow single or multiple dominant pathways. The predicted folding pathways are supported by the activation energy results derived from the Arrhenius plot as well as the NMR spectroscopy data.

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.

Cooperative tertiary interaction network guides RNA folding

Noncoding RNAs form unique 3D structures, which perform many regulatory functions. To understand how RNAs fold uniquely despite a small number of tertiary interaction motifs, we mutated the major tertiary interactions in a group I ribozyme by single-base substitutions. The resulting perturbations to the folding energy landscape were measured using SAXS, ribozyme activity, hydroxyl radical footprinting, and native PAGE. Double- and triple-mutant cycles show that most tertiary interactions have a small effect on the stability of the native state. Instead, the formation of core and peripheral structural motifs is cooperatively linked in near-native folding intermediates, and this cooperativity depends on the native helix orientation. The emergence of a cooperative interaction network at an early stage of folding suppresses nonnative structures and guides the search for the native state. We suggest that cooperativity in noncoding RNAs arose from natural selection of architectures conducive to forming a unique, stable fold.

An in vivo selection method to optimize trans-splicing ribozymes

Group I intron ribozymes can repair mutated mRNAs by replacing the 3'-terminal portion of the mRNA with their own 3'-exon. This trans-splicing reaction has the potential to treat genetic disorders and to selectively kill cancer cells or virus-infected cells. However, these ribozymes have not yet been used in therapy, partially due to a low in vivo trans-splicing efficiency. Previous strategies to improve the trans-splicing efficiencies focused on designing and testing individual ribozyme constructs. Here we describe a method that selects the most efficient ribozymes from millions of ribozyme variants. This method uses an in vivo rescue assay where the mRNA of an inactivated antibiotic resistance gene is repaired by trans-splicing group I intron ribozymes. Bacterial cells that express efficient trans-splicing ribozymes are able to grow on medium containing the antibiotic chloramphenicol. We randomized a 5'-terminal sequence of the Tetrahymena thermophila group I intron and screened a library with 9 × 10⁶ ribozyme variants for the best trans-splicing activity. The resulting ribozymes showed increased trans-splicing efficiency and help the design of efficient trans-splicing ribozymes for different sequence contexts. This in vivo selection method can now be used to optimize any sequence in trans-splicing ribozymes.

RNA folding pathways and the self-assembly of ribosomes

Many RNAs do not directly code proteins but are nonetheless indispensable to cellular function. These strands fold into intricate three-dimensional shapes that are essential structures in protein synthesis, splicing, and many other processes of gene regulation and expression. A variety of biophysical and biochemical methods are now showing, in real time, how ribosomal subunits and other ribonucleoprotein complexes assemble from their molecular components. Footprinting methods are particularly useful for studying the folding of long RNAs: they provide quantitative information about the conformational state of each residue and require little material. Data from footprinting complement the global information available from small-angle X-ray scattering or cryo-electron microscopy, as well as the dynamic information derived from single-molecule Förster resonance energy transfer (FRET) and NMR methods. In this Account, I discuss how we have used hydroxyl radical footprinting and other experimental methods to study pathways of RNA folding and 30S ribosome assembly. Hydroxyl radical footprinting probes the solvent accessibility of the RNA backbone at each residue in as little as 10 ms, providing detailed views of RNA folding pathways in real time. In conjunction with other methods such as solution scattering and single-molecule FRET, time-resolved footprinting of ribozymes showed that stable domains of RNA tertiary structure fold in less than 1 s. However, the free energy landscapes for RNA folding are rugged, and individual molecules kinetically partition into folding pathways that lead through metastable intermediates, stalling the folding or assembly process. Time-resolved footprinting was used to follow the formation of tertiary structure and protein interactions in the 16S ribosomal RNA (rRNA) during the assembly of 30S ribosomes. As previously observed in much simpler ribozymes, assembly occurs in stages, with individual molecules taking different routes to the final complex. Interactions occur concurrently in all domains of the 16S rRNA, and multistage protection of binding sites of individual proteins suggests that initial encounter complexes between the rRNA and ribosomal proteins are remodeled during assembly. Equilibrium footprinting experiments showed that one primary binding protein was sufficient to stabilize the tertiary structure of the entire 16S 5'-domain. The rich detail available from the footprinting data showed that the secondary assembly protein S16 suppresses non-native structures in the 16S 5'-domain. In doing so, S16 enables a conformational switch distant from its own binding site, which may play a role in establishing interactions with other domains of the 30S subunit. Together, the footprinting results show how protein-induced changes in RNA structure are communicated over long distances, ensuring cooperative assembly of even very large RNA-protein complexes such as the ribosome.

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.

ATP-dependent roles of the DEAD-box protein Mss116p in group II intron splicing in vitro and in vivo

The yeast DEAD-box protein Mss116p functions as a general RNA chaperone in splicing mitochondrial group I and group II introns. For most of its functions, Mss116p is thought to use ATP-dependent RNA unwinding to facilitate RNA structural transitions, but it has been suggested to assist in the folding of one group II intron (aI5γ) primarily by stabilizing a folding intermediate. Here we compare three aI5γ constructs: one with long exons, one with short exons, and a ribozyme construct lacking exons. The long exons result in slower splicing, suggesting that they misfold and/or stabilize nonnative intronic structures. Nevertheless, Mss116p acceleration of all three constructs depends on ATP and is inhibited by mutations that compromise RNA unwinding, suggesting similar mechanisms. Results of splicing assays and a new two-stage assay that separates ribozyme folding and catalysis indicate that maximal folding of all three constructs by Mss116p requires ATP-dependent RNA unwinding. ATP-independent activation is appreciable for only a subpopulation of the minimal ribozyme construct and not for constructs containing exons. As expected for a general RNA chaperone, Mss116p can also disrupt the native ribozyme, which can refold after Mss116p removal. Finally, using yeast strains with mitochondrial DNA containing only the single intron aI5γ, we show that Mss116p mutants promote splicing in vivo to degrees that correlate with their residual ATP-dependent RNA-unwinding activities. Together, our results indicate that, although DEAD-box proteins play multiple roles in RNA folding, the physiological function of Mss116p in aI5γ splicing includes a requirement for ATP-dependent local unfolding, allowing the conversion of nonfunctional RNA structure into functional RNA structure.

Nonequilibrium NMR Methods for Monitoring Protein and RNA Folding

The review introduces to time-resolved NMR spectroscopic investigations of the kinetics of protein and RNA folding. The description of the experimental investigations is discussed in the context of possible kinetic folding pathways showing the extent of information that can be gained from the various kinetic experiments. The review introduces to four different methods to initiate folding reactions in connection with time-resolved NMR experiments and discusses examples of refolding of the model proteinα-lactalbumin and of bistable RNAs.

RNA folding in living cells

RNA folding is the most essential process underlying RNA function. While significant progress has been made in understanding the forces driving RNA folding in vitro, exploring the rules governing intracellular RNA structure formation is still in its infancy. The cellular environment hosts a great diversity of factors that potentially influence RNA folding in vivo. For example, the nature of transcription and translation is known to shape the folding landscape of RNA molecules. Trans-acting factors such as proteins, RNAs and metabolites, among others, are also able to modulate the structure and thus the fate of an RNA. Here we summarize the ongoing efforts to uncover how RNA folds in living cells.

Mechanisms of StpA-mediated RNA remodeling

In bacteria, transcription, translation and gene regulation are highly coupled processes. The achievement of a certain functional structure at a distinct temporal and spatial position is therefore essential for RNA molecules. Proteins that facilitate this proper folding of RNA molecules are called RNA chaperones. Here a prominent example from E. coli is reviewed: the nucleoid associated protein StpA. Based on its various RNA remodeling functions, we propose a mechanistic model that explains how StpA promotes RNA folding. Through transient interactions via the RNA backbone, thereby shielding repelling charges in RNA, it pre-positions the RNA molecules for the successful formation of transition states from encounter complexes.

Roles of DEAD-box proteins in RNA and RNP Folding

RNAs and RNA-protein complexes (RNPs) traverse rugged energy landscapes as they fold to their native structures, and many continue to undergo conformational rearrangements as they function. Due to the inherent stability of local RNA structure, proteins are required to assist with RNA conformational transitions during initial folding and in exchange between functional structures. DEAD-box proteins are superfamily 2 RNA helicases that are ubiquitously involved in RNA-mediated processes. Some of these proteins use an ATP-dependent cycle of conformational changes to disrupt RNA structure nonprocessively, accelerating structural transitions of RNAs and RNPs in a manner that bears a strong resemblance to the activities of certain groups of protein chaperones. This review summarizes recent work using model substrates and tractable self-splicing intron RNAs, which has given new insights into how DEAD-box proteins promote RNA folding steps and conformational transitions, and it summarizes recent progress in identifying sites and mechanisms of DEAD-box protein activity within more complex cellular targets.

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.

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.

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-facilitated folding of group II intron ribozymes

Multiple studies hypothesize that DEAD-box proteins facilitate folding of the ai5gamma group II intron. However, these conclusions are generally inferred from splicing kinetics, and not from direct monitoring of DEAD-box protein-facilitated folding of the intron. Using native gel electrophoresis and dimethyl sulfate structural probing, we monitored Mss-116-facilitated folding of ai5gamma intron ribozymes and a catalytically active self-splicing RNA containing full-length intron and short exons. We found that the protein directly stimulates folding of these RNAs by accelerating formation of the compact near-native state. This process occurs in an ATP-independent manner, although ATP is required for the protein turnover. As Mss 116 binds RNA nonspecifically, most binding events do not result in the formation of the compact state, and ATP is required for the protein to dissociate from such nonproductive complexes and rebind the unfolded RNA. Results obtained from experiments at different concentrations of magnesium ions suggest that Mss 116 stimulates folding of ai5gamma ribozymes by promoting the formation of unstable folding intermediates, which is then followed by a cascade of folding events resulting in the formation of the compact near-native state. Dimethyl sulfate probing results suggest that the compact state formed in the presence of the protein is identical to the near-native state formed more slowly in its absence. Our results also indicate that Mss 116 does not stabilize the native state of the ribozyme, but that such stabilization results from binding of attached exons.

Engineering a family of synthetic splicing ribozymes

Controlling RNA splicing opens up possibilities for the synthetic biologist. The Tetrahymena ribozyme is a model group I self-splicing ribozyme that has been shown to be useful in synthetic circuits. To create additional splicing ribozymes that can function in synthetic circuits, we generated synthetic ribozyme variants by rationally mutating the Tetrahymena ribozyme. We present an alignment visualization for the ribozyme termed as structure information diagram that is similar to a sequence logo but with alignment data mapped on to secondary structure information. Using the alignment data and known biochemical information about the Tetrahymena ribozyme, we designed synthetic ribozymes with different primary sequences without altering the secondary structure. One synthetic ribozyme with 110 nt mutated retained 12% splicing efficiency in vivo. The results indicate that our biochemical understanding of the ribozyme is accurate enough to engineer a family of active splicing ribozymes with similar secondary structure but different primary sequences.

Kinetic cooperativity in Escherichia coli 30S ribosomal subunit reconstitution reveals additional complexity in the assembly landscape

The Escherichia coli 30S ribosomal subunit self-assembles in vitro in a hierarchical manner, with the RNA binding by proteins enabled by the prior binding of others under equilibrium conditions. Early 16S rRNA binding proteins also bind faster than late-binding proteins, but the specific causes for the slow binding of late proteins remain unclear. Previously, a pulse-chase monitored by quantitative mass spectrometry method was developed for monitoring 30S subunit assembly kinetics, and here a modified experimental scheme was used to probe kinetic cooperativity by including a step where subsets of ribosomal proteins bind and initiate assembly prior to the pulse-chase kinetics. In this work, 30S ribosomal subunit kinetic reconstitution experiments revealed that thermodynamic dependency does not always correlate with kinetic cooperativity. Some folding transitions that cause subsequent protein binding to be more energetically favorable do not result in faster protein binding. Although 3(') domain primary protein S7 is required for RNA binding by both proteins S9 and S19, prior binding of S7 accelerates the binding of S9, but not S19, indicating there is an additional mechanistic step required for S19 to bind. Such data on kinetic cooperativity and the presence of multiphasic assembly kinetics reveal complexity in the assembly landscape that was previously hidden.

Improving small-angle X-ray scattering data for structural analyses of the RNA world

Defining the shape, conformation, or assembly state of an RNA in solution often requires multiple investigative tools ranging from nucleotide analog interference mapping to X-ray crystallography. A key addition to this toolbox is small-angle X-ray scattering (SAXS). SAXS provides direct structural information regarding the size, shape, and flexibility of the particle in solution and has proven powerful for analyses of RNA structures with minimal requirements for sample concentration and volumes. In principle, SAXS can provide reliable data on small and large RNA molecules. In practice, SAXS investigations of RNA samples can show inconsistencies that suggest limitations in the SAXS experimental analyses or problems with the samples. Here, we show through investigations on the SAM-I riboswitch, the Group I intron P4-P6 domain, 30S ribosomal subunit from Sulfolobus solfataricus (30S), brome mosaic virus tRNA-like structure (BMV TLS), Thermotoga maritima asd lysine riboswitch, the recombinant tRNA(val), and yeast tRNA(phe) that many problems with SAXS experiments on RNA samples derive from heterogeneity of the folded RNA. Furthermore, we propose and test a general approach to reducing these sample limitations for accurate SAXS analyses of RNA. Together our method and results show that SAXS with synchrotron radiation has great potential to provide accurate RNA shapes, conformations, and assembly states in solution that inform RNA biological functions in fundamental ways.