Tetrahymena Ribozyme Research Articles

Goldilocks and RNA: where Mg2+ concentration is just right.

Magnesium, the most abundant divalent cation in cells, catalyzes RNA cleavage but also promotes RNA folding. Because folding can protect RNA from cleavage, we predicted a 'Goldilocks landscape', with local maximum in RNA lifetime at Mg2+ concentrations required for folding. Here, we use simulation and experiment to discover an innate and sophisticated mechanism of control of RNA lifetime. By simulation we characterized RNA Goldilocks landscapes and their dependence on cleavage and folding parameters. Experiments with yeast tRNAPhe and the Tetrahymena ribozyme P4-P6 domain show that structured RNAs can inhabit Goldilocks peaks. The Goldilocks peaks are tunable by differences in folded and unfolded cleavage rate constants, Mg2+ binding cooperativity, and Mg2+ affinity. Different folding and cleavage parameters produce Goldilocks landscapes with a variety of features. Goldilocks behavior allows ultrafine control of RNA chemical lifetime, whereas non-folding RNAs do not display Goldilocks peaks of protection. In sum, the effects of Mg2+ on RNA persistence are expected to be pleomorphic, both protecting and degrading RNA. In evolutionary context, Goldilocks behavior may have been a selectable trait of RNA in an early Earth environment containing Mg2+ and other metals.

Snapshots of the second-step self-splicing of Tetrahymena ribozyme revealed by cryo-EM

Group I introns are catalytic RNAs that coordinate two consecutive transesterification reactions for self-splicing. To understand how the group I intron promotes catalysis and coordinates self-splicing reactions, we determine the structures of L-16 Tetrahymena ribozyme in complex with a 5′-splice site analog product and a 3′-splice site analog substrate using cryo-EM. We solve six conformations from a single specimen, corresponding to different splicing intermediates after the first ester-transfer reaction. The structures reveal dynamics during self-splicing, including large conformational changes of the internal guide sequence and the J5/4 junction as well as subtle rearrangements of active-site metals and the hydrogen bond formed between the 2′-OH group of A261 and the N2 group of guanosine substrate. These results help complete a detailed structural and mechanistic view of this paradigmatic group I intron undergoing the second step of self-splicing.

Snapshots of the first-step self-splicing of Tetrahymena ribozyme revealed by cryo-EM.

Tetrahymena ribozyme is a group I intron, whose self-splicing is the result of two sequential ester-transfer reactions. To understand how it facilitates catalysis in the first self-splicing reaction, we used cryogenic electron microscopy (cryo-EM) to resolve the structures of L-16 Tetrahymena ribozyme complexed with a 11-nucleotide 5'-splice site analog substrate. Four conformations were achieved to 4.14, 3.18, 3.09 and 2.98 Å resolutions, respectively, corresponding to different splicing intermediates during the first enzymatic reaction. Comparison of these structures reveals structural alterations, including large conformational changes in IGS/IGSext (P1-P1ext duplex) and J5/4, as well as subtle local rearrangements in the G-binding site. These structural changes are required for the enzymatic activity of the Tetrahymena ribozyme. Our study demonstrates the ability of cryo-EM to capture dynamic RNA structural changes, ushering in a new era in the analysis of RNA structure-function by cryo-EM.

Topological crossing in the misfolded Tetrahymena ribozyme resolved by cryo-EM

The Tetrahymena group I intron has been a key system in the understanding of RNA folding and misfolding. The molecule folds into a long-lived misfolded intermediate (M) invitro, which has been known to form extensive native-like secondary and tertiary structures but is separated by an unknown kinetic barrier from the native state (N). Here, we used cryogenic electron microscopy (cryo-EM) to resolve misfolded structures of the Tetrahymena L-21 ScaI ribozyme. Maps of three M substates (M1, M2, M3) and one N state were achieved from a single specimen with overall resolutions of 3.5 Å, 3.8 Å, 4.0 Å, and 3.0 Å, respectively. Comparisons of the structures reveal that all the M substates are highly similar to N, except for rotation of a core helix P7 that harbors the ribozyme's guanosine binding site and the crossing of the strands J7/3 and J8/7 that connect P7 to the other elements in the ribozyme core. This topological difference between the M substates and N state explains the failure of 5'-splice site substrate docking in M, supports a topological isomer model for the slow refolding of M to N due to a trapped strand crossing, and suggests pathways for M-to-N refolding.

RiboDraw: semiautomated two-dimensional drawing of RNA tertiary structure diagrams

Publishing, discussing, envisioning, modeling, designing and experimentally determining RNA three-dimensional (3D) structures involve preparation of two-dimensional (2D) drawings that depict critical functional features of the subject molecules, such as noncanonical base pairs and protein contacts. Here, we describe RiboDraw, new software for crafting these drawings. We illustrate the features of RiboDraw by applying it to several RNAs, including the Escherichia coli tRNA-Phe, the P4–P6 domain of Tetrahymena ribozyme, a −1 ribosomal frameshift stimulation element from beet western yellows virus and the 5′ untranslated region of SARS-CoV-2. We show secondary structure diagrams of the 23S and 16S subunits of the E. coli ribosome that reflect noncanonical base pairs, ribosomal proteins and structural motifs, and that convey the relative positions of these critical features in 3D space. This software is a MATLAB package freely available at https://github.com/DasLab/RiboDraw.

Cryo-EM structures of full-length Tetrahymena ribozyme at 3.1Å resolution.

Single-particle cryogenic electron microscopy (cryo-EM) has become a standard technique for determining protein structures at atomic resolution1-3. However, cryo-EM studies of protein-free RNA are in their early days. The Tetrahymena thermophila group I self-splicing intron was the first ribozyme to be discovered and has been a prominent model system for the study of RNA catalysis and structure-function relationships4, but its full structure remains unknown. Here we report cryo-EM structures of the full-length Tetrahymena ribozyme in substrate-free and bound states at a resolution of 3.1Å. Newly resolved peripheral regions form two coaxially stacked helices; these are interconnected by two kissing loop pseudoknots that wrap around the catalytic core and include two previously unforeseen (to our knowledge) tertiary interactions. The global architecture is nearly identical in both states; only the internal guide sequence and guanosine binding site undergo a large conformational change and a localized shift, respectively, upon binding of RNA substrates. These results provide a long-sought structural view of a paradigmatic RNA enzyme and signal a new era for the cryo-EM-based study of structure-function relationships in ribozymes.

Effects of chain length of polyethylene glycol molecular crowders on a mutant Tetrahymena group I ribozyme lacking large peripheral module

While current group I ribozymes use several distinct strategies to function under conditions of low Mg2+ concentration (≤ 3 mM), a deletion mutant of the Tetrahymena ribozyme (ΔP5 ribozyme) is virtually inactive with 3 mM Mg2+ due to removal of the large peripheral module, P5abc, supporting the active conformation of the core module. We investigated the molecular crowding effects of synthetic polyethylene glycols (PEGs) on the activity of the ΔP5 ribozyme. Among PEG molecules with different chain lengths, PEG600 improved the activity of the ΔP5 ribozyme most effectively in the presence of 3 mM Mg2+.

Flexible Assembly of Engineered Tetrahymena Ribozymes Forming Polygonal RNA Nanostructures with Catalytic Ability.

Ribozymes with modular architecture constitute an attractive class of structural platforms for design and construction of nucleic acid nanostructures with biological functions. Through modular engineering of the Tetrahymena ribozyme, we have designed unit RNAs (L-RNAs), assembly of which formed ribozyme-based closed trimers and closed tetramers. Their catalytic activity was dependent on oligomer formation. In this study, the structural variety of L-RNA oligomers was extended by tuning their structural elements, yielding closed pentamers and closed hexamers. Their assembly properties were analyzed by electrophoretic mobility shift assay (EMSA) and atomic force microscopy (AFM).

Accelerated cryo-EM-guided determination of three-dimensional RNA-only structures.

The discovery and design of biologically important RNA molecules is outpacing three-dimensional structural characterization. Here, we demonstrate that cryo-EM can routinely resolve maps of RNA-only systems and that these maps enable sub-nanometer resolution coordinate estimation when complemented with multidimensional chemical mapping and Rosetta DRRAFTER computational modeling. This hybrid ‘Ribosolve’ pipeline detects and falsifies homologies and conformational rearrangements in eleven previously unknown 119- to 338-nucleotide protein-free RNA structures: full-length Tetrahymena ribozyme, hc16 ligase with and without substrate, full-length V. cholerae and F. nucleatum glycine riboswitch aptamers with and without glycine, Mycobacterium SAM-IV riboswitch with and without S-adenosylmethionine, and computer-designed ATP-TTR-3 aptamer with and without AMP. Simulation benchmarks, blind challenges, compensatory mutagenesis, cross-RNA homologies, and internal controls demonstrate that Ribosolve can accurately resolve the global architectures of RNA molecules, but does not resolve atomic details. These tests offer guidelines for making inferences in future RNA structural studies with similarly accelerated throughput.

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.

How the Conformations of an Internal Junction Contribute to Fold an RNA Domain.

Like proteins, some RNAs fold to compact structures. We can model functional RNAs as a series of short, rigid, base-paired elements, connected by non-base-paired nucleotides that serve as junctions. These connecting regions bend and twist, facilitating the formation of tertiary contacts that stabilize compact states. Here, we explore the roles of salt and junction sequence in determining the structures of a ubiquitous connector: an asymmetric internal loop. We focus on the J5/5a junction from the widely studied P4-P6 domain of the Tetrahymena ribozyme. Following the addition of magnesium ions to fold P4-P6, this junction bends dramatically, bringing the two halves of the RNA domain together for tertiary contact engagement. Using single-molecule fluorescence resonance energy transfer (smFRET), we examine the role of sequence and salt on model RNA constructs that contain these junction regions. We explore the wild-type J5/5a junction as well as two sequence variants. These junctions display distinct, salt-dependent conformations. Small-angle X-ray scattering (SAXS) measurements verify that these effects persist in the full-length P4-P6 domain. These measurements underscore the importance of junction sequence and interactions with ions in facilitating RNA folding.

In Crystallo Selection to Establish New RNA Crystal Contacts

Crystallography is a major technique for determining large RNA structures. Obtaining diffraction-quality crystals has been the bottleneck. Although several RNA crystallization methods have been developed, the field strongly needs additional approaches. Here we invented an in crystallo selection strategy for identifying mutations that enhance a target RNA's crystallizability. The strategy includes constructing an RNA pool containing random mutations, obtaining crystals, and amplifying the sequences enriched by crystallization. We demonstrated a proof-of-principle application to the P4-P6 domain from the Tetrahymena ribozyme. We further determined the structures of four selected mutants. All four establish new crystal lattice contacts while maintaining the native structure. Three mutants achieve this by relocating bulges and one by making a helix more flexible. In crystallo selection provides opportunities to improve crystals of RNAs or RNA-ligand complexes. Our results also suggest that mutants may be rationally designed for crystallization by "walking" a bulge along the RNA chain.

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

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

Differential Assembly of Catalytic Interactions within the Conserved Active Sites of Two Ribozymes.

Molecular recognition is central to biology and a critical aspect of RNA function. Yet structured RNAs typically lack the preorganization needed for strong binding and precise positioning. A striking example is the group I ribozyme from Tetrahymena, which binds its guanosine substrate (G) orders of magnitude slower than diffusion. Binding of G is also thermodynamically coupled to binding of the oligonucleotide substrate (S) and further work has shown that the transition from E•G to E•S•G accompanies a conformational change that allows G to make the active site interactions required for catalysis. The group I ribozyme from Azoarcus has a similarly slow association rate but lacks the coupled binding observed for the Tetrahymena ribozyme. Here we test, using G analogs and metal ion rescue experiments, whether this absence of coupling arises from a higher degree of preorganization within the Azoarcus active site. Our results suggest that the Azoarcus ribozyme forms cognate catalytic metal ion interactions with G in the E•G complex, interactions that are absent in the Tetrahymena E•G complex. Thus, RNAs that share highly similar active site architectures and catalyze the same reactions can differ in the assembly of transition state interactions. More generally, an ability to readily access distinct local conformational states may have facilitated the evolutionary exploration needed to attain RNA machines that carry out complex, multi-step processes.

Probing the kinetic and thermodynamic consequences of the tetraloop/tetraloop receptor monovalent ion-binding site in P4-P6 RNA by smFRET.

Structured RNA molecules play roles in central biological processes and understanding the basic forces and features that govern RNA folding kinetics and thermodynamics can help elucidate principles that underlie biological function. Here we investigate one such feature, the specific interaction of monovalent cations with a structured RNA, the P4-P6 domain of the Tetrahymena ribozyme. We employ single molecule FRET (smFRET) approaches as these allow determination of folding equilibrium and rate constants over a wide range of stabilities and thus allow direct comparisons without the need for extrapolation. These experiments provide additional evidence for specific binding of monovalent cations, Na+ and K+, to the RNA tetraloop-tetraloop receptor (TL-TLR) tertiary motif. These ions facilitate both folding and unfolding, consistent with an ability to help order the TLR for binding and further stabilize the tertiary contact subsequent to attainment of the folding transition state.

Twenty years of RNA: reflections from the RNP world.

Twenty years of RNA: reflections from the RNP world.

Empirical demonstration of environmental sensing in catalytic RNA: evolution of interpretive behavior at the origins of life.

BackgroundThe origins of life on the Earth required chemical entities to interact with their environments in ways that could respond to natural selection. The concept of interpretation, where biotic entities use signs in their environment as proxy for the existence of other items of selective value in their environment, has been proposed on theoretical grounds to be relevant to the origins and early evolution of life. However this concept has not been demonstrated empirically.ResultsHere, we present data that certain catalytic RNA sequences have properties that would enable interpretation of divalent cation levels in their environment. By assaying the responsiveness of two variants of the Tetrahymena ribozyme to the Ca2+ ion as a sign for the more catalytically useful Mg2+ ion, we show an empirical proof-of-principle that interpretation can be an evolvable trait in RNA, often suggested as a model system for early life. In particular we demonstrate that in vitro, the wild-type version of the Tetrahymena ribozyme is not interpretive, in that it cannot use Ca2+ as a sign for Mg2+. Yet a variant of this sequence containing five mutations that alter its ability to utilize the Ca2+ ion engenders a strong interpretive characteristic in this RNA.ConclusionsWe have shown that RNA molecules in a test tube can meet the minimum criteria for the evolution of interpretive behaviour in regards to their responses to divalent metal ion concentrations in their environment. Interpretation in RNA molecules provides a property entirely dependent on natural physico-chemical interactions, but capable of shaping the evolutionary trajectory of macromolecules, especially in the earliest stages of life’s history.

A kinetic and thermodynamic framework for the Azoarcus group I ribozyme reaction.

Determination of quantitative thermodynamic and kinetic frameworks for ribozymes derived from the Azoarcus group I intron and comparisons to their well-studied analogs from the Tetrahymena group I intron reveal similarities and differences between these RNAs. The guanosine (G) substrate binds to the Azoarcus and Tetrahymena ribozymes with similar equilibrium binding constants and similar very slow association rate constants. These and additional literature observations support a model in which the free ribozyme is not conformationally competent to bind G and in which the probability of assuming the binding-competent state is determined by tertiary interactions of peripheral elements. As proposed previously, the slow binding of guanosine may play a role in the specificity of group I intron self-splicing, and slow binding may be used analogously in other biological processes. The internal equilibrium between ribozyme-bound substrates and products is similar for these ribozymes, but the Azoarcus ribozyme does not display the coupling in the binding of substrates that is observed with the Tetrahymena ribozyme, suggesting that local preorganization of the active site and rearrangements within the active site upon substrate binding are different for these ribozymes. Our results also confirm the much greater tertiary binding energy of the 5′-splice site analog with the Azoarcus ribozyme, binding energy that presumably compensates for the fewer base-pairing interactions to allow the 5′-exon intermediate in self splicing to remain bound subsequent to 5′-exon cleavage and prior to exon ligation. Most generally, these frameworks provide a foundation for design and interpretation of experiments investigating fundamental properties of these and other structured RNAs.

Development and Applications of Coarse‐Grained Models for RNA

AbstractIn contrast to proteins, much less attention has been focused on the development of computational models for describing RNA molecules, which are being recognized as playing key roles in many cellular functions. Current atomically detailed force fields are not accurate enough to capture the properties of even simple nucleic acid constructs. In this article, we review our efforts to develop coarse‐grained (CG) models that capture the underlying physics for the particular length scale of interest. Two models are discussed. One of them is the three interaction site (TIS) model, in which each nucleotide is represented by three beads corresponding to sugar, phosphate, and base. The other is the self‐organized polymer (SOP) model, in which each nucleotide is represented as a single interaction center. Applications of the TIS model to study the complexity of hairpin formation and the effects of crowding in a shifting equilibrium between two conformations in human telomerase pseudoknot are described. The work on crowding illustrates a direct link to the activity of telomerase. We use the SOP model to describe the response of the Tetrahymena ribozyme to force. The simulated unfolding pathways agree well with single molecule pulling experiments. We also review predictions for the unfolding pathways for the Azoarcus ribozyme. The success of the CG applications to describe dynamics in RNA gives hope that more complex processes involving RNA‐protein interactions can be tackled using variants of the proposed models.

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.6min−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 (~1min−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.