Crystal structure of a phage Twort group I ribozyme–product complex

Exploring the energetic and conformational properties of the sequence space connecting naturally occurring RNA tetraloop receptor motifs

Folded RNAs contain tertiary contact motifs whose structures and energetics are conserved across different RNAs. The transferable properties of RNA motifs simplify the RNA folding problem, but measuring energetic and conformational properties of many motifs remains a challenge. Here, we use a high-throughput thermodynamic approach to investigate how sequence changes alter the binding properties of naturally occurring motifs, the GAAA tetraloop • tetraloop receptor (TLR) interactions. We measured the binding energies and conformational preferences of TLR sequences that span mutational pathways from the canonical 11ntR to two other natural TLRs, the IC3R and Vc2R. While the IC3R and Vc2R share highly similar energetic and conformational properties, the landscapes that map the sequence changes for their conversion from the 11ntR to changes in these properties differ dramatically. Differences in the energetic landscapes stem from the mutations needed to convert the 11ntR to the IC3R and Vc2R rather than a difference in the intrinsic energetic architectures of these TLRs. The conformational landscapes feature several nonnative TLR variants with conformational preferences that differ from both the initial and final TLRs; these species represent potential branching points along the multidimensional sequence space to sequences with greater fitness in other RNA contexts with alternative conformational preferences. Our high-throughput, quantitative approach reveals the complex nature of sequence-fitness landscapes and leads to models for their molecular origins. Systematic and quantitative molecular approaches provide critical insights into understanding the evolution of natural RNAs as they traverse complex landscapes in response to selective pressures.

Antagonistic conflict between transposon-encoded introns and guide RNAs

TnpB nucleases represent the evolutionary precursors to CRISPR-Cas12 and are widespread in all domains of life. IS605-family TnpB homologs function as programmable RNA-guided homing endonucleases in bacteria, driving transposon maintenance through DNA double-strand break-stimulated homologous recombination. In this work, we uncovered molecular mechanisms of the transposition life cycle of IS607-family elements that, notably, also encode group I introns. We identified specific features for a candidate "IStron" from Clostridium botulinum that allow the element to carefully control the relative levels of spliced products versus functional guide RNAs. Our results suggest that IStron transcripts evolved an ability to balance competing and mutually exclusive activities that promote selfish transposon spread while limiting adverse fitness costs on the host. Collectively, this work highlights molecular innovation in the multifunctional utility of transposon-encoded noncoding RNAs.

Antagonistic conflict between transposon-encoded introns and guide RNAs

TnpB nucleases represent the evolutionary precursors to CRISPR-Cas12 and are widespread in all domains of life, presumably due to the critical roles they play in transposon proliferation. IS605family TnpB homologs function in bacteria as programmable homing endonucleases by exploiting transposon-encoded guide RNAs to cleave vacant genomic sites, thereby driving transposon maintenance through DSB-stimulated homologous recombination. Whether this pathway is conserved in other genetic contexts, and in association with other transposases, is unknown. Here we uncover molecular mechanisms of transposition and RNA-guided DNA cleavage by IS607-family elements that, remarkably, also encode catalytic, self-splicing group I introns. After reconstituting and systematically investigating each of these biochemical activities for a candidate 'IStron' derived from Clostridium botulinum, we discovered sequence and structural features of the transposon-encoded RNA that satisfy molecular requirements of a group I intron and TnpB guide RNA, while still retaining the ability to be faithfully mobilized at the DNA level by the TnpA transposase. Strikingly, intron splicing was strongly repressed not only by TnpB, but also by the secondary structure of ωRNA alone, allowing the element to carefully control the relative levels of spliced products versus functional guide RNAs. Our results suggest that IStron transcripts have evolved a sensitive equilibrium to balance competing and mutually exclusive activities that promote transposon maintenance while limiting adverse fitness costs on the host. Collectively, this work explains how diverse enzymatic activities emerged during the selfish spread of IS607-family elements and highlights molecular innovation in the multi-functional utility of transposon-encoded noncoding RNAs.

RNA as a Major-Groove Ligand: RNA-RNA and RNA-DNA Triplexes Formed by GAA and UUC or TTC Sequences

Friedreich's ataxia is associated with noncanonical nucleic acid structures that emerge when GAA:TTC repeats in the first intron of the FXN gene expand beyond a critical number of repeats. Specifically, the noncanonical repeats are associated with both triplexes and R-loops. Here, we present an in silico investigation of all possible triplexes that form by attaching a third RNA strand to an RNA:RNA or DNA:DNA duplex, complementing previous DNA-based triplex studies. For both new triplexes results are similar. For a pyridimine UUC⁺ third strand, the parallel orientation is stable while its antiparallel counterpart is unstable. For a neutral GAA third strand, the parallel conformation is stable. A protonated GA⁺A third strand is stable in both parallel and antiparallel orientations. We have also investigated Na⁺ and Mg²⁺ ion distributions around the triplexes. The presence of Mg²⁺ ions helps stabilize neutral, antiparallel GAA triplexes. These results (along with previous DNA-based studies) allow for the emergence of a complete picture of the stability and structural characteristics of triplexes based on the GAA and TTC/UUC sequences, thereby contributing to the field of trinucleotide repeats and the associated unusual structures that trigger expansion.

Sub-3-Å cryo-EM structure of RNA enabled by engineered homomeric self-assembly

High-resolution structural studies are essential for understanding the folding and function of diverse RNAs. Herein, we present a nanoarchitectural engineering strategy for efficient structural determination of RNA-only structures using single-particle cryogenic electron microscopy (cryo-EM). This strategy-ROCK (RNA oligomerization-enabled cryo-EM via installing kissing loops)-involves installing kissing-loop sequences onto the functionally nonessential stems of RNAs for homomeric self-assembly into closed rings with multiplied molecular weights and mitigated structural flexibility. ROCK enables cryo-EM reconstruction of the Tetrahymena group I intron at 2.98-Å resolution overall (2.85 Å for the core), allowing de novo model building of the complete RNA, including the previously unknown peripheral domains. ROCK is further applied to two smaller RNAs-the Azoarcus group I intron and the FMN riboswitch, revealing the conformational change of the former and the bound ligand in the latter. ROCK holds promise to greatly facilitate the use of cryo-EM in RNA structural studies.

Integral equation models for solvent in macromolecular crystals

The solvent can occupy up to ∼70% of macromolecular crystals, and hence, having models that predict solvent distributions in periodic systems could improve the interpretation of crystallographic data. Yet, there are few implicit solvent models applicable to periodic solutes, and crystallographic structures are commonly solved assuming a flat solvent model. Here, we present a newly developed periodic version of the 3D-reference interaction site model (RISM) integral equation method that is able to solve efficiently and describe accurately water and ion distributions in periodic systems; the code can compute accurate gradients that can be used in minimizations or molecular dynamics simulations. The new method includes an extension of the Ornstein–Zernike equation needed to yield charge neutrality for charged solutes, which requires an additional contribution to the excess chemical potential that has not been previously identified; this is an important consideration for nucleic acids or any other charged system where most or all the counter- and co-ions are part of the “disordered” solvent. We present several calculations of proteins, RNAs, and small molecule crystals to show that x-ray scattering intensities and the solvent structure predicted by the periodic 3D-RISM solvent model are in closer agreement with the experiment than are intensities computed using the default flat solvent model in the refmac5 or phenix refinement programs, with the greatest improvement in the 2 to 4 Å range. Prospects for incorporating integral equation models into crystallographic refinement are discussed.

Naturally occurring three-way junctions can be repurposed as genetically encoded RNA-based sensors

Small molecules can be imaged in living cells using biosensors composed of RNA. However, RNA-based devices are difficult to design. Here, we describe a versatile platform for designing RNA-based fluorescent small-molecule sensors using naturally occurring highly stable three-way junction RNAs. We show that ligand-binding aptamers and fluorogenic aptamers can be inserted into three-way junctions and connected in a way that enables the three-way junction to function as a small-molecule-regulated fluorescent sensor in vitro and in cells. The sensors are designed so that the interhelical stabilizing interactions in the three-way junction are only induced upon ligand binding. We use these RNA-based devices to measure the dynamics of S-adenosylmethionine levels in mammalian cells in real time. We show that this strategy is compatible with diverse metabolite-binding RNA aptamers, fluorogenic aptamers, and three-way junctions. Overall, these data demonstrate a versatile method for readily generating RNA devices that function in living cells.

Refinement of RNA Structures Using Amber Force Fields

Atomic models for nucleic acids derived from X-ray diffraction data at low resolution provide much useful information, but the observed scattering intensities can be fit with models that can differ in structural detail. Tradtional geometric restraints favor models that have bond length and angle terms derived from small molecule crystal structures. Here we explore replacing these restraints with energy gradients derived from force fields, including recently developed integral equation models to account for the effects of water molecules and ions that are not part of the explicit model. We compare conventional and force-field based refinements for 22 RNA crystals, ranging in resolution from 1.1 to 3.6 Å. As expected, it can be important to account for solvent screening of charge–charge interactions, especially in the crowded environment of a nucleic acid crystal. The newly refined models can show improvements in torsion angles and hydrogen-bonding interactions, and can significantly reduce unfavorable atomic clashes, while maintaining or improving agreement with observed scattering intensities.

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.

Multipolar electrostatics for hairpin and pseudoknots in RNA: Improving the accuracy of force field potential energy function

Molecular dynamics (MD) simulations that rely on force field methods has been widely used to explore the structure and function of RNAs. However, the current commonly used force fields are limited by the electrostatic description offered by atomic charge, dipole and at most quadrupole moments, failing to capture the anisotropic picture of electronic features. Actually, the distribution of electrons around atomic nuclei is not spherically symmetric but is geometry dependent. A multipolar electrostatic model based on high rank multipole moments is described in this work, which allows us to combine polarizability and anisotropy of electron density. RNA secondary structure was taken as a research system, and its substructures including stem, loops (hairpin loop, bulge loop, internal loop, and multi-branch loop), and pseudoknots (H-type and K-type) were investigated, respectively, as well as the hairpin. First, the atom-atom electrostatic properties derived from one chain of a duplex RNA 2MVY in our previous work (Ref. 58) were measured by the pilot RNA systems of hairpin, hairpin loop, stem, and H-type pseudoknot, respectively. The prediction results were not satisfactory. Consequently, to obtain a general set of electrostatic parameters for RNA force fields, the convergence behavior of the atom-atom electrostatic interactions in the pilot RNA systems was explored using high rank atomic multipole moments. The pilot RNA systems were cut into four types of different-sized molecular fragments, and the single nucleotide fragment and nucleotide-paired fragment proved to be the most reasonable systems for base-unpairing regions and base-pairing regions to investigate the convergence behavior of all types of atom-atom electrostatic interactions, respectively. Transferability of the electrostatic properties drawn from the pilot RNA systems to the corresponding test systems was also investigated. Furthermore, the convergence behavior of atomic electrostatic interactions in other substructures including bulge loop, internal loop, multi-branch loop, and K-type pseudoknot was expected to be modeled via the hairpin.

Affinity and Structural Analysis of the U1A RNA Recognition Motif with Engineered Methionines to Improve Experimental Phasing

RNA plays a central role in all organisms and can fold into complex structures to orchestrate function. Visualization of such structures often requires crystallization, which can be a bottleneck in the structure-determination process. To promote crystallization, an RNA-recognition motif (RRM) of the U1A spliceosomal protein has been co-opted as a crystallization module. Specifically, the U1-snRNA hairpin II (hpII) single-stranded loop recognized by U1A can be transplanted into an RNA target to promote crystal contacts and to attain phase information via molecular replacement or anomalous diffraction methods using selenomethionine. Herein, we produced the F37M/F77M mutant of U1A to augment the phasing capability of this powerful crystallization module. Selenomethionine-substituted U1A(F37M/F77M) retains high affinity for hpII (K _D of 59.7 ± 11.4 nM). The 2.20 Å resolution crystal structure reveals that the mutated sidechains make new S-π interactions in the hydrophobic core and are useful for single-wavelength anomalous diffraction. Crystals were also attained of U1A(F37M/F77M) in complex with a bacterial preQ₁-II riboswitch. The F34M/F37M/F77M mutant was introduced similarly into a lab-evolved U1A variant (TBP6.9) that recognizes the internal bulged loop of HIV-1 TAR RNA. We envision that this short RNA sequence can be placed into non-essential duplex regions to promote crystallization and phasing of target RNAs. We show that selenomethionine-substituted TBP6.9(F34M/F37M/F77M) binds a TAR variant wherein the apical loop was replaced with a GNRA tetraloop (K _D of 69.8 ± 2.9 nM), laying the groundwork for use of TBP6.9(F34M/F37M/F77M) as a crystallization module. These new tools are available to the research community.

Base-pair Opening Dynamics of Nucleic Acids in Relation to Their Biological Function

Base-pair opening is a conformational transition that is required for proper biological function of nucleic acids. Hydrogen exchange, observed by NMR spectroscopic experiments, is a widely used method to study the thermodynamics and kinetics of base-pair opening in nucleic acids. The hydrogen exchange data of imino protons are analyzed based on a two-state (open/closed) model for the base-pair, where hydrogen exchange only occurs from the open state. In this review, we discuss examples of how hydrogen exchange data provide insight into several interesting biological processes involving functional interactions of nucleic acids: i) selective recognition of DNA by proteins; ii) regulation of RNA cleavage by site-specific mutations; iii) intermolecular interaction of proteins with their target DNA or RNA; iv) formation of PNA:DNA hybrid duplexes.

•

This review systematically summarizes hydrogen exchange theory for base-paired imino protons of nucleic acids.

•

Base-pair opening kinetics explain how the DNA can be selectively recognized by its target proteins.

•

Base-pair opening kinetics explain the mechanisms by which site-specific mutations regulate RNA cleavage.

•

Hydrogen exchange studies can elucidate the intermolecular interaction of proteins with their target DNA or RNA.

Integrated evolution of ribosomal RNAs, introns, and intron nurseries

The initial components of ribosomes first appeared more than 3.8 billion years ago during a time when many types of RNAs were evolving. While modern ribosomes are complex molecular machines consisting of rRNAs and proteins, they were assembled during early evolution by the association and joining of small functional RNA units. Introns may have provided the means to ligate many of these pieces together. All four classes of introns (group I, group II, spliceosomal, and archaeal) are present in many rRNA gene loci over a broad phylogenetic range. A survey of rRNA intron sequences across the three major life domains suggests that some of the classes of introns may have diverged from one another within rRNA gene loci. Analyses of rRNA sequences revealed self-splicing group I and group II introns are present in ancestral regions of the SSU (small subunit) and LSU (large subunit), whereas spliceosomal and archaeal introns appeared in sections of the rRNA that evolved later. Most classes of introns increased in number for approximately 1 billion years. However, their frequencies are low in the most recently evolved regions added to the SSU and LSU rRNAs. Furthermore, many of the introns appear to have been in the same locations for billions of years, suggesting an ancient origin for these sequences. In this Perspectives paper, I reviewed and analyzed rRNA intron sequences, locations, structural characteristics, and splicing mechanisms; and suggest that rRNA gene loci may have served as evolutionary nurseries for intron formation and diversification.

Revisiting the Closed-Loop Model and the Nature of mRNA 5'-3' Communication

Communication between the 5' and 3' ends of mature eukaryotic mRNAs lies at the heart of gene regulation, likely arising at the same time as the eukaryotic lineage itself. Our view of how and why it occurs has been shaped by elegant experiments that led to nearly universal acceptance of the "closed-loop model." However, new observations suggest that this classic model needs to be reexamined, revised, and expanded. Here, we address fundamental questions about the closed-loop model and discuss how a growing understanding of mRNA structure, dynamics, and intermolecular interactions presents new experimental opportunities. We anticipate that the application of emerging methods will lead to expanded models that include the role of intrinsic mRNA structure and quantitative dynamic descriptions of 5'-3' proximity linked to the functional status of an mRNA and will better reflect the messy realities of the crowded and rapidly changing cellular environment.

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.

Contact Statistics Highlight Distinct Organizing Principles of Proteins and RNA

Although both RNA and proteins have densely packed native structures, chain organizations of these two biopolymers are fundamentally different. Motivated by the recent discoveries in chromatin folding that interphase chromosomes have territorial organization with signatures pointing to metastability, we analyzed the biomolecular structures deposited in the Protein Data Bank and found that the intrachain contact probabilities, P(s) as a function of the arc length s, decay in power-law ∼s(-γ) over the intermediate range of s, 10 ≲ s ≲ 110. We found that the contact probability scaling exponent is γ ≈ 1.11 for large RNA (N > 110), γ ≈ 1.41 for small-sized RNA (N < 110), and γ ≈ 1.65 for proteins. Given that Gaussian statistics is expected for a fully equilibrated chain in polymer melts, the deviation of γ-value from γ = 1.5 for the subchains of large RNA in the native state suggests that the chain configuration of RNA is not fully equilibrated. It is visually clear that folded structures of large-sized RNA (N ≳ 110) adopt crumpled structures, partitioned into modular multidomains assembled by proximal sequences along the chain, whereas the polypeptide chain of folded proteins looks better mixed with the rest of the structure. Our finding of γ ≈ 1 for large RNA might be an ineluctable consequence of the hierarchical ordering of the secondary to tertiary elements in the folding process.

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.

Use of a Fluorescent Aptamer RNA as an Exonic Sequence to Analyze Self-Splicing Ability of aGroup I Intron from Structured RNAs

Group I self-splicing intron constitutes an important class of functional RNA molecules that can promote chemical transformation. Although the fundamental mechanism of the auto-excision from its precursor RNA has been established, convenient assay systems for its splicing activity are still useful for a further understanding of its detailed mechanism and of its application. Because some host RNA sequences, to which group I introns inserted form stable three-dimensional (3D) structures, the effects of the 3D structures of exonic elements on the splicing efficiency of group I introns are important but not a fully investigated issue. We developed an assay system for group I intron self-splicing by employing a fluorescent aptamer RNA (spinach RNA) as a model exonic sequence inserted by the Tetrahymena group I intron. We investigated self-splicing of the intron from spinach RNA, serving as a model exonic sequence with a 3D structure.

Accumulation of Stable Full-Length Circular Group I Intron RNAs during Heat-Shock

Group I introns in nuclear ribosomal RNA of eukaryotic microorganisms are processed by splicing or circularization. The latter results in formation of full-length circular introns without ligation of the exons and has been proposed to be active in intron mobility. We applied qRT-PCR to estimate the copy number of circular intron RNA from the myxomycete Didymium iridis. In exponentially growing amoebae, the circular introns are nuclear and found in 70 copies per cell. During heat-shock, the circular form is up-regulated to more than 500 copies per cell. The intron harbours two ribozymes that have the potential to linearize the circle. To understand the structural features that maintain circle integrity, we performed chemical and enzymatic probing of the splicing ribozyme combined with molecular modeling to arrive at models of the inactive circular form and its active linear counterpart. We show that the two forms have the same overall structure but differ in key parts, including the catalytic core element P7 and the junctions at which reactions take place. These differences explain the relative stability of the circular species, demonstrate how it is prone to react with a target molecule for circle integration and thus supports the notion that the circular form is a biologically significant molecule possibly with a role in intron mobility.

RNA Structure Refinement Using the ERRASER-Phenix Pipeline

The final step of RNA crystallography involves the fitting of coordinates into electron density maps. The large number of backbone atoms in RNA presents a difficult and tedious challenge, particularly when experimental density is poor. The ERRASER-Phenix pipeline can improve an initial set of RNA coordinates automatically based on a physically realistic model of atomic-level RNA interactions. The pipeline couples diffraction-based refinement in Phenix with the Rosetta-based real-space refinement protocol ERRASER (Enumerative Real-Space Refinement ASsisted by Electron density under Rosetta). The combination of ERRASER and Phenix can improve the geometrical quality of RNA crystallographic models while maintaining or improving the fit to the diffraction data (as measured by R free). Here we present a complete tutorial for running ERRASER-Phenix through the Phenix GUI, from the command-line, and via an application in the Rosetta On-line Server that Includes Everyone (ROSIE).

An active site rearrangement within the Tetrahymena group I ribozyme releases nonproductive interactions and allows formation of catalytic interactions

Biological catalysis hinges on the precise structural integrity of an active site that binds and transforms its substrates and meeting this requirement presents a unique challenge for RNA enzymes. Functional RNAs, including ribozymes, fold into their active conformations within rugged energy landscapes that often contain misfolded conformers. Here we uncover and characterize one such "off-pathway" species within an active site after overall folding of the ribozyme is complete. The Tetrahymena group I ribozyme (E) catalyzes cleavage of an oligonucleotide substrate (S) by an exogenous guanosine (G) cofactor. We tested whether specific catalytic interactions with G are present in the preceding E•S•G and E•G ground-state complexes. We monitored interactions with G via the effects of 2'- and 3'-deoxy (-H) and -amino (-NH(2)) substitutions on G binding. These and prior results reveal that G is bound in an inactive configuration within E•G, with the nucleophilic 3'-OH making a nonproductive interaction with an active site metal ion termed MA and with the adjacent 2'-OH making no interaction. Upon S binding, a rearrangement occurs that allows both -OH groups to contact a different active site metal ion, termed M(C), to make what are likely to be their catalytic interactions. The reactive phosphoryl group on S promotes this change, presumably by repositioning the metal ions with respect to G. This conformational transition demonstrates local rearrangements within an otherwise folded RNA, underscoring RNA's difficulty in specifying a unique conformation and highlighting Nature's potential to use local transitions of RNA in complex function.

Increased efficiency of evolved group I intron spliceozymes by decreased side product formation

The group I intron ribozyme from Tetrahymena was recently reengineered into a trans-splicing variant that is able to remove 100-nt introns from pre-mRNA, analogous to the spliceosome. These spliceozymes were improved in this study by 10 rounds of evolution in Escherichia coli cells. One clone with increased activity in E. coli cells was analyzed in detail. Three of its 10 necessary mutations extended the substrate binding duplexes, which led to increased product formation and reduced cleavage at the 5'-splice site. One mutation in the conserved core of the spliceozyme led to a further reduction of cleavage at the 5'-splice site but an increase in cleavage side products at the 3'-splice site. The latter was partially reduced by six additional mutations. Together, the mutations increased product formation while reducing activity at the 5'-splice site and increasing activity at the 3'-splice site. These results show the adaptation of a ribozyme that evolved in nature for cis-splicing to trans-splicing, and they highlight the interdependent function of nucleotides within group I intron ribozymes. Implications for the possible use of spliceozymes as tools in research and therapy, and as a model for the evolution of the spliceosome, are discussed.

Evolution of RNA-protein interactions: non-specific binding led to RNA splicing activity of fungal mitochondrial tyrosyl-tRNA synthetases

Studies of tRNA synthetases that adapted to assist the splicing of group I introns provide insight into how proteins can evolve new RNA-binding functions.

The Neurospora crassa mitochondrial tyrosyl-tRNA synthetase (mtTyrRS; CYT-18 protein) evolved a new function as a group I intron splicing factor by acquiring the ability to bind group I intron RNAs and stabilize their catalytically active RNA structure. Previous studies showed: (i) CYT-18 binds group I introns by using both its N-terminal catalytic domain and flexibly attached C-terminal anticodon-binding domain (CTD); and (ii) the catalytic domain binds group I introns specifically via multiple structural adaptations that occurred during or after the divergence of Peziomycotina and Saccharomycotina. However, the function of the CTD and how it contributed to the evolution of splicing activity have been unclear. Here, small angle X-ray scattering analysis of CYT-18 shows that both CTDs of the homodimeric protein extend outward from the catalytic domain, but move inward to bind opposite ends of a group I intron RNA. Biochemical assays show that the isolated CTD of CYT-18 binds RNAs non-specifically, possibly contributing to its interaction with the structurally different ends of the intron RNA. Finally, we find that the yeast mtTyrRS, which diverged from Pezizomycotina fungal mtTyrRSs prior to the evolution of splicing activity, binds group I intron and other RNAs non-specifically via its CTD, but lacks further adaptations needed for group I intron splicing. Our results suggest a scenario of constructive neutral (i.e., pre-adaptive) evolution in which an initial non-specific interaction between the CTD of an ancestral fungal mtTyrRS and a self-splicing group I intron was “fixed” by an intron RNA mutation that resulted in protein-dependent splicing. Once fixed, this interaction could be elaborated by further adaptive mutations in both the catalytic domain and CTD that enabled specific binding of group I introns. Our results highlight a role for non-specific RNA binding in the evolution of RNA-binding proteins.

The acquisition of new modes of post-transcriptional gene regulation played an important role in the evolution of eukaryotes and was achieved by an increase in the number of RNA-binding proteins with new functions. RNA-binding proteins bind directly to double- or single-stranded RNA and regulate many cellular processes. Here, we address how proteins evolve new RNA-binding functions by using as a model system a fungal mitochondrial tyrosyl-tRNA synthetase that evolved to acquire a novel function in splicing group I introns. Group I introns are RNA enzymes (or “ribozymes”) that catalyze their own removal from transcripts, but can become dependent upon proteins to stabilize their active structure. We show that the C-terminal domain of the synthetase is flexibly attached and has high non-specific RNA-binding activity that likely pre-dated the evolution of splicing activity. Our findings suggest an evolutionary scenario in which an initial non-specific interaction between an ancestral synthetase and a self-splicing group I intron was fixed by an intron RNA mutation, thereby making it dependent upon the protein for structural stabilization. The interaction then evolved by the acquisition of adaptive mutations throughout the protein and RNA that increased both the splicing efficiency and its protein-dependence. Our results suggest a general mechanism by which non-specific binding interactions can lead to the evolution of new RNA-binding functions and provide novel insights into splicing and synthetase mechanisms.

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.

Nuclear magnetic resonance structure of the III-IV-V three-way junction from the Varkud satellite ribozyme and identification of magnesium-binding sites using paramagnetic relaxation enhancement

The VS ribozyme is a catalytic RNA found within some natural isolates of Neurospora that is being used as a model system to improve our understanding of RNA structure, catalysis, and engineering. The catalytic domain contains five helical domains (SLII-SLVI) that are organized by two three-way junctions. The III-IV-V junction is required for high-affinity binding of the substrate domain (SLI) through formation of a kissing loop interaction with SLV. Here, we determine the high-resolution nuclear magnetic resonance (NMR) structure of a 47-nucleotide RNA containing the III-IV-V junction (J345). The J345 RNA adopts a Y-shaped fold typical of the family C three-way junctions, with coaxial stacking between stems III and IV and an acute angle between stems III and V. The NMR structure reveals that the core of the III-IV-V junction contains four stacked base triples, a U-turn motif, a cross-strand stacking interaction, an A-minor interaction, and a ribose zipper. In addition, the NMR structure shows that the cCUUGg tetraloop used to stabilize stem IV adopts a novel RNA tetraloop fold, different from the known gCUUGc tetraloop structure. Using Mn(2+)-induced paramagnetic relaxation enhancement, we identify six Mg(2+)-binding sites within J345, including one associated with the cCUUGg tetraloop and two with the junction core. The NMR structure of J345 likely represents the conformation of the III-IV-V junction in the context of the active VS ribozyme and suggests that this junction functions as a dynamic hinge that contributes to substrate recognition and catalysis. Moreover, this study highlights a new role for family C three-way junctions in long-range tertiary interactions.

Spliceozymes: ribozymes that remove introns from pre-mRNAs in trans

Group I introns are pre-mRNA introns that do not require the spliceosome for their removal. Instead, they fold into complex three-dimensional structures and catalyze two transesterification reactions, thereby excising themselves and joining the flanking exons. These catalytic RNAs (ribozymes) have been modified previously to work in trans, whereby the ribozymes can recognize a splice site on a substrate RNA and replace the 5′- or 3′-portion of the substrate. Here we describe a new variant of the group I intron ribozyme from Tetrahymena that recognizes two splice sites on a substrate RNA, removes the intron sequences between the splice sites, and joins the flanking exons, analogous to the action of the spliceosome. This ‘group I spliceozyme’ functions in vitro and in vivo, and it is able to mediate a growth phenotype in E. coli cells. The intron sequences of the target pre-mRNAs are constrained near the splice sites but can carry a wide range of sequences in their interior. Because the splice site recognition sequences can be adjusted to different splice sites, the spliceozyme may have the potential for wide applications as tool in research and therapy.

Speciation of a group I intron into a lariat capping ribozyme

The lariat-capping (LC) ribozyme is a natural ribozyme isolated from eukaryotic microorganisms. Despite apparent structural similarity to group I introns, the LC ribozyme catalyzes cleavage by a 2',5' branching reaction, leaving the 3' product with a 3-nt lariat cap that functionally substitutes for a conventional mRNA cap in the downstream pre-mRNA encoding a homing endonuclease. We describe the crystal structures of the precleavage and postcleavage LC ribozymes, which suggest that structural features inherited from group I ribozymes have undergone speciation due to profound changes in molecular selection pressure, ultimately giving rise to an original branching ribozyme family. The structures elucidate the role of key elements that regulate the activity of the LC ribozyme by conformational switching and suggest a mechanism by which the signal for branching is transmitted to the catalytic core. The structures also show how conserved interactions twist residues, forming the lariat to join chemical groups involved in branching.

Detecting RNA tertiary folding by sedimentation velocity analytical ultracentrifugation

Analytical Ultracentrifugation (AUC) is a highly sensitive technique for detecting global conformational features of biological molecules and molecular interactions in solution. When operated in a sedimentation velocity (SV) recording mode, it reports precisely on the hydrodynamic properties of a molecule, including its sedimentation and diffusion coefficients, which can be used to calculate its hydrated radius, as well as, to estimate its global shape. This chapter describes the application of SV-AUC to the detection of global conformational changes accompanying equilibrium counterion induced tertiary folding of structured RNA molecules. A brief theoretical background is provided at the beginning, aimed at familiarizing the readers with the operational principle of the technique; then, a detailed set of instructions is provided on how to design, conduct, and analyze the data from an equilibrium RNA folding experiment, using SV-AUC.

Bacterial group I introns: mobile RNA catalysts

Group I introns are intervening sequences that have invaded tRNA, rRNA and protein coding genes in bacteria and their phages. The ability of group I introns to self-splice from their host transcripts, by acting as ribozymes, potentially renders their insertion into genes phenotypically neutral. Some group I introns are mobile genetic elements due to encoded homing endonuclease genes that function in DNA-based mobility pathways to promote spread to intronless alleles. Group I introns have a limited distribution among bacteria and the current assumption is that they are benign selfish elements, although some introns and homing endonucleases are a source of genetic novelty as they have been co-opted by host genomes to provide regulatory functions. Questions regarding the origin and maintenance of group I introns among the bacteria and phages are also addressed.

RNA synthesis and purification for structural studies

RNAs play pivotal roles in the cell, ranging from catalysis (e.g., RNase P), acting as adaptor molecule (tRNA) to regulation (e.g., riboswitches). Precise understanding of its three-dimensional structures has given unprecedented insight into the molecular basis for all of these processes. Nevertheless, structural studies on RNA are still limited by the very special nature of this polymer. The most common methods for the determination of 3D RNA structures are NMR and X-ray crystallography. Both methods have their own set of requirements and give different amounts of information about the target RNA. For structural studies, the major bottleneck is usually obtaining large amounts of highly pure and homogeneously folded RNA. Especially for X-ray crystallography it can be necessary to screen a large number of variants to obtain well-ordered single crystals. In this mini-review we give an overview about strategies for the design, in vitro production, and purification of RNA for structural studies.

On the propagation of errors

The policy of the Protein Data Bank (PDB) that the first deposition of a small-molecule ligand, even with erroneous atom numbering, sets a precedent over accepted nomenclature rules is disputed. Recommendations regarding ligand molecules in the PDB are suggested.

Structure and function of the T-loop structural motif in noncoding RNAs

The T-loop is a frequently occurring five-nucleotide motif found in the structure of noncoding RNAs where it is commonly assumed to play an important role in stabilizing the tertiary RNA structure by facilitating long-range interactions between different regions of the molecule. T-loops were first identified in tRNA(Phe) and a formal consensus sequence for this motif was formulated and later revised based on analyses of the crystal structures of prokaryotic ribosomal RNAs and RNase P and the corresponding primary sequence of their orthologues. In the past decade, several new structures of large RNA molecules have been added to the RCSB Protein Data Bank, including the eukaryotic ribosome, a self-splicing group II intron, numerous synthetases in complex with their cognate transfer RNAs (tRNAs), transfer-messenger RNA (tmRNA) in complex with SmpB, several riboswitches, and a complex of bacterial RNase P bound to its tRNA substrate. In this review, the search for T-loops is extended to these new RNA molecules based on the previously established structure-based criteria. The review highlights and discusses the function and additional roles of T-loops in four broad categories of RNA molecules, namely tRNAs, ribosomal RNAs (rRNAs), P RNAs, and RNA genetic elements. Additionally, the potential application for T-loops as interaction modules is also discussed.

Structure of an A-form RNA duplex obtained by degradation of 6S RNA in a crystallization droplet

In the course of a crystallographic study of a 132 nt variant of Aquifex aeolicus 6S RNA, a crystal structure of an A-form RNA duplex containing 12 base pairs was solved at a resolution of 2.6 Å. In fact, the RNA duplex is part of the 6S RNA and was obtained by accidental but precise degradation of the 6S RNA in a crystallization droplet. 6S RNA degradation was confirmed by microscopic observation of crystals and gel electrophoresis of crystallization droplets. The RNA oligomers obtained form regular A-form duplexes containing three GoU wobble-type base pairs, one of which engages in intermolecular contacts through a ribose-zipper motif at the crystal-packing interface.

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.

Metal-ion rescue revisited: biochemical detection of site-bound metal ions important for RNA folding

Within the three-dimensional architectures of RNA molecules, divalent metal ions populate specific locations, shedding their water molecules to form chelates. These interactions help the RNA adopt and maintain specific conformations and frequently make essential contributions to function. Defining the locations of these site-bound metal ions remains challenging despite the growing database of RNA structures. Metal-ion rescue experiments have provided a powerful approach to identify and distinguish catalytic metal ions within RNA active sites, but the ability of such experiments to identify metal ions that contribute to tertiary structure acquisition and structural stability is less developed and has been challenged. Herein, we use the well-defined P4-P6 RNA domain of the Tetrahymena group I intron to reevaluate prior evidence against the discriminatory power of metal-ion rescue experiments and to advance thermodynamic descriptions necessary for interpreting these experiments. The approach successfully identifies ligands within the RNA that occupy the inner coordination sphere of divalent metal ions and distinguishes them from ligands that occupy the outer coordination sphere. Our results underscore the importance of obtaining complete folding isotherms and establishing and evaluating thermodynamic models in order to draw conclusions from metal-ion rescue experiments. These results establish metal-ion rescue as a rigorous tool for identifying and dissecting energetically important metal-ion interactions in RNAs that are noncatalytic but critical for RNA tertiary structure.

Exploring purine N7 interactions via atomic mutagenesis: the group I ribozyme as a case study

Atomic mutagenesis has emerged as a powerful tool to unravel specific interactions in complex RNA molecules. An early extensive study of analogs of the exogenous guanosine nucleophile in group I intron self-splicing by Bass and Cech demonstrated structure-function relationships analogous to those seen for protein ligands and provided strong evidence for a well-formed substrate binding site made of RNA. Subsequent functional and structural studies have confirmed these interacting sites and extended our understanding of them, with one notable exception. Whereas 7-methyl guanosine did not affect reactivity in the original study, a subsequent study revealed a deleterious effect of the seemingly more conservative 7-deaza substitution. Here we investigate this paradox, studying these and other analogs with the more thoroughly characterized ribozyme derived from the Tetrahymena group I intron. We found that the 7-deaza substitution lowers binding by ~20-fold, relative to the cognate exogenous guanosine nucleophile, whereas binding and reaction with 7-methyl and 8-aza-7-deaza substitutions have no effect. These and additional results suggest that there is no functionally important contact between the N7 atom of the exogenous guanosine and the ribozyme. Rather, they are consistent with indirect effects introduced by the N7 substitution on stacking interactions and/or solvation that are important for binding. The set of analogs used herein should be valuable in deciphering nucleic acid interactions and how they change through reaction cycles for other RNAs and RNA/protein complexes.

Computational prediction of efficient splice sites for trans-splicing ribozymes

Group I introns have been engineered into trans-splicing ribozymes capable of replacing the 3'-terminal portion of an external mRNA with their own 3'-exon. Although this design makes trans-splicing ribozymes potentially useful for therapeutic application, their trans-splicing efficiency is usually too low for medical use. One factor that strongly influences trans-splicing efficiency is the position of the target splice site on the mRNA substrate. Viable splice sites are currently determined using a biochemical trans-tagging assay. Here, we propose a rapid and inexpensive alternative approach to identify efficient splice sites. This approach involves the computation of the binding free energies between ribozyme and mRNA substrate. We found that the computed binding free energies correlate well with the trans-splicing efficiency experimentally determined at 18 different splice sites on the mRNA of chloramphenicol acetyl transferase. In contrast, our results from the trans-tagging assay correlate less well with measured trans-splicing efficiency. The computed free energy components suggest that splice site efficiency depends on the following secondary structure rearrangements: hybridization of the ribozyme's internal guide sequence (IGS) with mRNA substrate (most important), unfolding of substrate proximal to the splice site, and release of the IGS from the 3'-exon (least important). The proposed computational approach can also be extended to fulfill additional design requirements of efficient trans-splicing ribozymes, such as the optimization of 3'-exon and extended guide sequences.

Synthesis, properties, and applications of oligonucleotides containing an RNA dinucleotide phosphorothiolate linkage

RNA represents a prominent class of biomolecules. Present in all living systems, RNA plays many essential roles in gene expression, regulation, and development. Accordingly, many biological processes depend on the accurate enzymatic processing, modification, and cleavage of RNA. Understanding the catalytic mechanisms of these enzymes therefore represents an important goal in defining living systems at the molecular level. In this context, RNA molecules bearing 3'- or 5'-S-phosphorothiolate linkages comprise what are arguably among the most incisive mechanistic probes available. They have been instrumental in showing that RNA splicing systems are metalloenzymes and in mapping the ligands that reside within RNA active sites. The resulting models have in turn verified the functional relevance of crystal structures. In other cases, phosphorothiolates have offered an experimental strategy to circumvent the classic problem of kinetic ambiguity; mechanistic enzymologists have used this tool to assign precise roles to catalytic groups as general acids or bases. These insights into macromolecular function are enabled by the synthesis of nucleic acids bearing phosphorothiolate linkages and the unique chemical properties they impart. In this Account, we review the synthesis, properties, and applications of oligonucleotides and oligodeoxynucleotides containing an RNA dinucleotide phosphorothiolate linkage. Phosphorothioate linkages are structurally very similar to phosphorothiolate linkages, as reflected in the single letter of difference in nomenclature. Phosphorothioate substitutions, in which sulfur replaces one or both nonbridging oxygens within a phosphodiester linkage, are now widely available and are used routinely in numerous biochemical and medicinal applications. Indeed, synthetic phosphorothioate linkages can be introduced readily via a sulfurization step programmed into automated solid-phase oligonucleotide synthesizers. In contrast, phosphorothiolate oligonucleotides, in which sulfur replaces a specific 3'- or 5'-bridging oxygen, have presented a more difficult synthetic challenge, requiring chemical alterations to the attached sugar moiety. Here we begin by outlining the synthetic strategies used to access these phosphorothiolate RNA analogues. The Arbuzov reaction and phosphoramidite chemistry are often brought to bear in creating either 3'- or 5'-S-phosphorothiolate dinucleotides. We then summarize the responses of the phosphorothiolate derivatives to chemical and enzymatic cleavage agents, as well as mechanistic insights their use has engendered. They demonstrate particular utility as probes of metal-ion-dependent phosphotransesterification, general acid-base-catalyzed phosphotransesterification, and rate-limiting chemistry. The 3'- and 5'-S-phosphorothiolates have proven invaluable in elucidating the mechanisms of enzymatic and nonenzymatic phosphoryl transfer reactions. Considering that RNA cleavage represents a fundamental step in the maturation, degradation, and regulation of this important macromolecule, the significant synthetic challenges that remain offer rich research opportunities.

2'-Fluoro substituents can mimic native 2'-hydroxyls within structured RNA

The ability of fluorine in a C-F bond to act as a hydrogen bond acceptor is controversial. To test such ability in complex RNA macromolecules, we have replaced native 2'-OH groups with 2'-F and 2'-H groups in two related systems, the Tetrahymena group I ribozyme and the ΔC209 P4-P6 RNA domain. In three cases the introduced 2'-F mimics the native 2'-OH group, suggesting that the fluorine atom can accept a hydrogen bond. In each of these cases the native hydroxyl group interacts with a purine exocyclic amine. Our results give insight about the properties of organofluorine and suggest a possible general biochemical signature for tertiary interactions between 2'-hydroxyl groups and exocyclic amino groups within RNA.

Helix-length compensation studies reveal the adaptability of the VS ribozyme architecture

Compensatory mutations in RNA are generally regarded as those that maintain base pairing, and their identification forms the basis of phylogenetic predictions of RNA secondary structure. However, other types of compensatory mutations can provide higher-order structural and evolutionary information. Here, we present a helix-length compensation study for investigating structure–function relationships in RNA. The approach is demonstrated for stem-loop I and stem-loop V of the Neurospora VS ribozyme, which form a kissing–loop interaction important for substrate recognition. To rapidly characterize the substrate specificity (k_cat/K_M) of several substrate/ribozyme pairs, a procedure was established for simultaneous kinetic characterization of multiple substrates. Several active substrate/ribozyme pairs were identified, indicating the presence of limited substrate promiscuity for stem Ib variants and helix-length compensation between stems Ib and V. 3D models of the I/V interaction were generated that are compatible with the kinetic data. These models further illustrate the adaptability of the VS ribozyme architecture for substrate cleavage and provide global structural information on the I/V kissing–loop interaction. By exploring higher-order compensatory mutations in RNA our approach brings a deeper understanding of the adaptability of RNA structure, while opening new avenues for RNA research.

Two distinct catalytic strategies in the hepatitis δ virus ribozyme cleavage reaction

The hepatitis delta virus (HDV) ribozyme and related RNAs are widely dispersed in nature. This RNA is a small nucleolytic ribozyme that self-cleaves to generate products with a 2',3'-cyclic phosphate and a free 5'-hydroxyl. Although small ribozymes are dependent on divalent metal ions under biologically relevant buffer conditions, they function in the absence of divalent metal ions at high ionic strengths. This characteristic suggests that a functional group within the covalent structure of small ribozymes is facilitating catalysis. Structural and mechanistic analyses have demonstrated that the HDV ribozyme active site contains a cytosine with a perturbed pK(a) that serves as a general acid to protonate the leaving group. The reaction of the HDV ribozyme in monovalent cations alone never approaches the velocity of the Mg(2+)-dependent reaction, and there is significant biochemical evidence that a Mg(2+) ion participates directly in catalysis. A recent crystal structure of the HDV ribozyme revealed that there is a metal binding pocket in the HDV ribozyme active site. Modeling of the cleavage site into the structure suggested that this metal ion can interact directly with the scissile phosphate and the nucleophile. In this manner, the Mg(2+) ion can serve as a Lewis acid, facilitating deprotonation of the nucleophile and stabilizing the conformation of the cleavage site for in-line attack of the nucleophile at the scissile phosphate. This catalytic strategy had previously been observed only in much larger ribozymes. Thus, in contrast to most large and small ribozymes, the HDV ribozyme uses two distinct catalytic strategies in its cleavage reaction.

Structure-function analysis from the outside in: long-range tertiary contacts in RNA exhibit distinct catalytic roles

The conserved catalytic core of the Tetrahymena group I ribozyme is encircled by peripheral elements. We have conducted a detailed structure-function study of the five long-range tertiary contacts that fasten these distal elements together. Mutational ablation of each of the tertiary contacts destabilizes the folded ribozyme, indicating a role of the peripheral elements in overall stability. Once folded, three of the five tertiary contact mutants exhibit defects in overall catalysis that range from 20- to 100-fold. These and the subsequent results indicate that the structural ring of peripheral elements does not act as a unitary element; rather, individual connections have distinct roles as further revealed by kinetic and thermodynamic dissection of the individual reaction steps. Ablation of P14 or the metal ion core/metal ion core receptor (MC/MCR) destabilizes docking of the substrate-containing P1 helix into tertiary interactions with the ribozyme's conserved core. In contrast, ablation of the L9/P5 contact weakens binding of the guanosine nucleophile by slowing its association, without affecting P1 docking. The P13 and tetraloop/tetraloop receptor (TL/TLR) mutations had little functional effect and small, local structural changes, as revealed by hydroxyl radical footprinting, whereas the P14, MC/MCR, and L9/P5 mutants show structural changes distal from the mutation site. These changes extended into regions of the catalytic core involved in docking or guanosine binding. Thus, distinct allosteric pathways couple the long-range tertiary contacts to functional sites within the conserved core. This modular functional specialization may represent a fundamental strategy in RNA structure-function interrelationships.

Prevalence of syn nucleobases in the active sites of functional RNAs

Biological RNAs, like their DNA counterparts, contain helical stretches, which have standard Watson-Crick base pairs in the anti conformation. Most functional RNAs also adopt geometries with far greater complexity such as bulges, loops, and multihelical junctions. Occasionally, nucleobases in these regions populate the syn conformation wherein the base resides close to or over the ribose sugar, which leads to a more compact state. The importance of the syn conformation to RNA function is largely unknown. In this study, we analyze 51 RNAs with tertiary structure, including aptamers, riboswitches, ribozymes, and ribosomal RNAs, for number, location, and properties of syn nucleobases. These RNAs represent the set of nonoverlapping, moderate- to high-resolution structures available at present. We find that syn nucleobases are much more common among purines than pyrimidines, and that they favor C2'-endo-like conformations especially among those nucleobases in the intermediate syn conformation. Strikingly, most syn nucleobases participate in tertiary stacking and base-pairing interactions: Inspection of RNA structures revealed that the majority of the syn nucleobases are in regions assigned to function, with many syn nucleobases interacting directly with a ligand or ribozyme active site. These observations suggest that judicious placement of conformationally restricted nucleotides biased into the syn conformation could enhance RNA folding and catalysis. Such changes could also be useful for locking RNAs into functionally competent folds for use in X-ray crystallography and NMR.

RNA molecules with conserved catalytic cores but variable peripheries fold along unique energetically optimized pathways

Functional and kinetic constraints must be efficiently balanced during the folding process of all biopolymers. To understand how homologous RNA molecules with different global architectures fold into a common core structure we determined, under identical conditions, the folding mechanisms of three phylogenetically divergent group I intron ribozymes. These ribozymes share a conserved functional core defined by topologically equivalent tertiary motifs but differ in their primary sequence, size, and structural complexity. Time-resolved hydroxyl radical probing of the backbone solvent accessible surface and catalytic activity measurements integrated with structural-kinetic modeling reveal that each ribozyme adopts a unique strategy to attain the conserved functional fold. The folding rates are not dictated by the size or the overall structural complexity, but rather by the strength of the constituent tertiary motifs which, in turn, govern the structure, stability, and lifetime of the folding intermediates. A fundamental general principle of RNA folding emerges from this study: The dominant folding flux always proceeds through an optimally structured kinetic intermediate that has sufficient stability to act as a nucleating scaffold while retaining enough conformational freedom to avoid kinetic trapping. Our results also suggest a potential role of naturally selected peripheral A-minor interactions in balancing RNA structural stability with folding efficiency.