Structural diversity of self-cleaving ribozymes

Evolution Driven Microscale Combinatorial Chemistry in Intracellular Mimicking Droplets to Engineer Thermostable RNA for Cellular Imaging

Fluorescent light-up aptamer/fluorogen pairs are powerful tools for tracking RNA in the cell, however limitations in thermostability and fluorescence intensity exist. Current in vitro selection techniques struggle to mimic complex intracellular environments, limiting in vivo biomolecule functionality. Taking inspiration from microenvironment-dependent RNA folding observed in cells and organelle-mimicking droplets, an efficient system is created that uses microscale heated water droplets to simulate intracellular conditions, effectively replicating the intracellular RNA folding landscape. This system is integrated with microfluidic droplet sorting to evolve RNA aptamers. Through this approach, an RNA aptamer is engineered with improved fluorescence activity by exploring the chemical fitness landscape under biomimetic conditions. The enhanced RNA aptamer named eBroccoli has increased fluorescence intensity and thermal stability, both in vitro and in vivo in bacterial and mammalian cells. In mammalian cell culture conditions, a fluorescence improvement of 3.9-times is observed and biological thermal stability up to 45 °C is observed in bacterial systems. eBroccoli enable real-time visualization of nanoscale stress granule formation in mammalian cells during heat shock at 42 °C. By introducing the concept of "biomimetic equivalence" based on RNA folding, the platform offers a simple yet effective strategy to mimic intracellular complexity in evolution-based engineering.

Computational Design of Allosteric Ribozymes via Genetic Algorithms

In vitro selection of allosteric ribozymes has many challenges, such as complex and time-consuming experimental procedures, uncertain results, and the unwanted functionality of the enriched sequences. The precise computational design of allosteric ribozymes is achievable using RNA secondary structure folding principles. The computational design of allosteric ribozymes is based on experimentally validated EAs, random search algorithms, and a partition function for RNA folding. The in silico design achieves an accuracy exceeding 90%. Various algorithms with different logic gates have been automated via computer programs that can quickly create many allosteric sequences. This can eliminate the need for in vitro selection of allosteric ribozymes, thus vastly reducing the time and cost required.

A comparative survey of the influence of small self-cleaving ribozymes on gene expression in human cell culture

Self-cleaving ribozymes are versatile tools for synthetic biologists when it comes to controlling gene expression. Up to date, 12 different classes are known, and over the past decades more and more details about their structure, cleavage mechanisms and natural environments have been uncovered. However, when these motifs are applied to mammalian gene expression constructs, the outcome can often be unexpected. A variety of factors, such as surrounding sequences and positioning of the ribozyme influences the activity and hence performance of catalytic RNAs. While some information about the efficiency of individual ribozymes (each tested in specific contexts) is known, general trends obtained from standardized, comparable experiments are lacking, complicating decisions such as which ribozyme to choose and where to insert it into the target mRNA. In many cases, application-specific optimization is required, which can be very laborious. Here, we systematically compared different classes of ribozymes within the 3'-UTR of a given reporter gene. We then examined position-dependent effects of the best-performing ribozymes. Moreover, we tested additional variants of already widely used hammerhead ribozymes originating from various organisms. We were able to identify functional structures suited for aptazyme design and generated highly efficient hammerhead ribozyme variants originating from the human genome. The present dataset will aide decisions about how to apply ribozymes for affecting gene expression as well as for developing ribozyme-based switches for controlling gene expression in human cells.

Fast Transport and Transformation of Biomacromolecular Substances via Thermo-Stimulated Active "Inhalation-Exhalation" Cycles of Hierarchically Structured Smart pNIPAM-DNA Hydrogels

Although smart hydrogels hold great promise in biosensing and biomedical applications, their response to external stimuli is governed by the passive diffusion-dependent substance transport between hydrogels and environments and within the 3D hydrogel matrices, resulting in slow response to biomacromolecules and limiting their extensive applications. Herein, inspired by the respiration systems of organisms, an active strategy to achieve highly efficient biomolecular substance transport through the thermo-stimulated "inhalation-exhalation" cycles of hydrogel matrices is demonstrated. The cryo-structured poly(N-isopropylacrylamide) (pNIPAM)-DNA hydrogels, composed of functional DNA-tethered pNIPAM networks and free-water-containing macroporous channels, exhibit thermally triggered fast and reversible shrinking/swelling cycles with high-volume changes, which drive the formation of dynamic water stream to accelerate the intake of external substances and expelling of endogenous substances, thus promoting the functional properties of hydrogel systems. Demonstrated by catalytic DNAzyme and CRISPR-Cas12a-incorporating hydrogels, significantly enhanced catalytic efficiency with up to 280% and 390% is achieved, upon the introduction of active "inhalation-exhalation" cycles, respectively. Moreover, remotely near-infrared (NIR)-triggering of "inhalation-exhalation" cycles is achieved after the introduction of NIR-responsive MXene nanosheets into the hydrogel matrix. These hydrogel systems with enhanced substance transport and transformation properties hold promise in the development of more effective biosensing and therapeutic systems.

Inspecting the Ribozyme Region of Hepatitis Delta Virus Genotype 1: Conservation and Variability

The hepatitis delta virus (HDV) genome has an autocatalytic region called the ribozyme, which is essential for viral replication. The aim of this study was to use next-generation sequencing (NGS) to analyze the ribozyme quasispecies (QS) in order to study its evolution and identify highly conserved regions potentially suitable for a gene-silencing strategy. HDV RNA was extracted from 2 longitudinal samples of chronic HDV patients and the ribozyme (nucleotide, nt 688-771) was analyzed using NGS. QS conservation, variability and genetic distance were analyzed. Mutations were identified by aligning sequences with their specific genotype consensus. The main relevant mutations were tested in vitro. The ribozyme was conserved overall, with a hyper-conserved region between nt 715-745. No difference in QS was observed over time. The most variable region was between nt 739-769. Thirteen mutations were observed, with three showing a higher frequency: T23C, T69C and C64 deletion. This last strongly reduced HDV replication by more than 1 log in vitro. HDV Ribozyme QS was generally highly conserved and was maintained during follow-up. The most conserved portion may be a valuable target for a gene-silencing strategy. The presence of the C64 deletion may strongly impair viral replication, as it is a potential mechanism of viral persistence.

A mechanistic view of long noncoding RNAs in cancer

Long noncoding RNAs (lncRNAs) have emerged as important modulators of a wide range of biological processes in normal and disease states. In particular, lncRNAs have garnered significant interest as novel players in the molecular pathology of cancer, spurring efforts to define the functions, and elucidate the mechanisms through which cancer‐associated lncRNAs operate. In this review, we discuss the prevalent mechanisms employed by lncRNAs, with a critical assessment of the methodologies used to determine each molecular function. We survey the abilities of cancer‐associated lncRNAs to enact diverse trans functions throughout the nucleus and in the cytoplasm and examine the local roles of cis‐acting lncRNAs in modulating the expression of neighboring genes. In linking lncRNA functions and mechanisms to their roles in cancer biology, we contend that a detailed molecular understanding of lncRNA functionality is key to elucidating their contributions to tumorigenesis and to unlocking their therapeutic potential.

This article is categorized under:

Regulatory RNAs/RNAi/Riboswitches > Regulatory RNAs

RNA in Disease and Development > RNA in Disease

The hammerhead self-cleaving motif as a precursor to complex endonucleolytic ribozymes

Connections between distinct catalytic RNA motifs through networks of mutations that retain catalytic function (neutral networks) were likely central to the evolution of biocatalysis. Despite suggestions that functional RNAs collectively form an interconnected web of neutral networks, little evidence has emerged to demonstrate the existence of such intersecting networks in naturally occurring RNAs. Here we show that neutral networks of two naturally occurring, seemingly unrelated endonucleolytic ribozymes, the hammerhead (HH) and hairpin (HP), intersect. Sequences at the intersection of these networks exhibit catalytic functions corresponding to both ribozymes by potentially populating both catalytic folds and enable a smooth crossover between the two. Small and structurally simple endonucleolytic motifs like the HH ribozyme could, through mutational walks along their neutral networks, encounter novel catalytic phenotypes, and structurally flexible, bifunctional sequences at the intersection of these networks could have acted as nodes for evolutionary diversification in an RNA world. Considering the simplicity and small size of the HH ribozyme, we propose that this self-cleaving motif could have been a precursor to other more complex endonucleolytic ribozymes. More generally, our results suggest that RNAs that possess distinct sequences, structures, and catalytic functions, can potentially share evolutionary history through mutational connections in sequence space.

Non-associative phase separation in an evaporating droplet as a model for prebiotic compartmentalization

The synthetic pathways of life’s building blocks are envisaged to be through a series of complex prebiotic reactions and processes. However, the strategy to compartmentalize and concentrate biopolymers under prebiotic conditions remains elusive. Liquid-liquid phase separation is a mechanism by which membraneless organelles form inside cells, and has been hypothesized as a potential mechanism for prebiotic compartmentalization. Associative phase separation of oppositely charged species has been shown to partition RNA, but the strongly negative charge exhibited by RNA suggests that RNA-polycation interactions could inhibit RNA folding and its functioning inside the coacervates. Here, we present a prebiotically plausible pathway for non-associative phase separation within an evaporating all-aqueous sessile droplet. We quantitatively investigate the kinetic pathway of phase separation triggered by the non-uniform evaporation rate, together with the Marangoni flow-driven hydrodynamics inside the sessile droplet. With the ability to undergo liquid-liquid phase separation, the drying droplets provide a robust mechanism for formation of prebiotic membraneless compartments, as demonstrated by localization and storage of nucleic acids, in vitro transcription, as well as a three-fold enhancement of ribozyme activity. The compartmentalization mechanism illustrated in this model system is feasible on wet organophilic silica-rich surfaces during early molecular evolution.

Prebiotic compartmentalization could prove essential for the evolution of life. Guo et al. show that liquid-liquid separation in an aqueous two-phase system driven by evaporation may already suffice to facilitate chemical processes required for the RNA world hypothesis.

The case of the missing allosteric ribozymes

The RNA World theory encompasses the hypothesis that sophisticated ribozymes and riboswitches were the primary drivers of metabolic processes in ancient organisms. Several types of catalytic RNAs and many classes of ligand-sensing RNA switches still exist in modern cells. Curiously, allosteric ribozymes formed by the merger of RNA enzyme and RNA switch components are largely absent in today's biological systems. This is true despite the striking abundances of various classes of both self-cleaving ribozymes and riboswitch aptamers. Here we present the known types of ligand-controlled ribozymes and riboswitches and discuss the possible reasons why fused ribozyme-aptamer constructs have been disfavored through evolution.

B2 and ALU retrotransposons are self-cleaving ribozymes whose activity is enhanced by EZH2

Transposable elements make up half of the mammalian genome. One of the most abundant is the short interspersed nuclear element (SINE). Among their million copies, B2 accounts for ∼350,000 in the mouse genome and has garnered special interest because of emerging roles in epigenetic regulation. Our recent work demonstrated that B2 RNA binds stress genes to retard transcription elongation. Although epigenetically silenced, B2s become massively up-regulated during thermal and other types of stress. Specifically, an interaction between B2 RNA and the Polycomb protein, EZH2, results in cleavage of B2 RNA, release of B2 RNA from chromatin, and activation of thermal stress genes. Although an established RNA-binding protein and histone methyltransferase, EZH2 is not known to be a nuclease. Here, we provide evidence for the surprising conclusion that B2 is a self-cleaving ribozyme. Ribozyme activity depends on Mg⁺² and monovalent cations but is resistant to protease treatment. However, contact with EZH2 accelerates cleavage rate by >100-fold, suggesting that EZH2 promotes a cleavage-competent RNA conformation. B2 modification-interference analysis demonstrates that phosphorothioate changes at A and C nucleotides can substitute for EZH2. B2 nucleotides 45 to 55 and 100 to 101 are essential for activity. Finally, another family of SINEs, the human ALU element, also produces a self-cleaving RNA and is cleaved during T-cell activation as well as thermal and endoplasmic reticulum (ER) stress. Thus, B2/ALU SINEs may be classified as "epigenetic ribozymes" that function as transcriptional switches during stress. Given their high copy numbers, B2 and ALU may represent the predominant ribozyme activity in mammalian cells.

Novel ribozymes: discovery, catalytic mechanisms, and the quest to understand biological function

Small endonucleolytic ribozymes promote the self-cleavage of their own phosphodiester backbone at a specific linkage. The structures of and the reactions catalysed by members of individual families have been studied in great detail in the past decades. In recent years, bioinformatics studies have uncovered a considerable number of new examples of known catalytic RNA motifs. Importantly, entirely novel ribozyme classes were also discovered, for most of which both structural and biochemical information became rapidly available. However, for the majority of the new ribozymes, which are found in the genomes of a variety of species, a biological function remains elusive. Here, we concentrate on the different approaches to find catalytic RNA motifs in sequence databases. We summarize the emerging principles of RNA catalysis as observed for small endonucleolytic ribozymes. Finally, we address the biological functions of those ribozymes, where relevant information is available and common themes on their cellular activities are emerging. We conclude by speculating on the possibility that the identification and characterization of proteins that we hypothesize to be endogenously associated with catalytic RNA might help in answering the ever-present question of the biological function of the growing number of genomically encoded, small endonucleolytic ribozymes.

Evolution of ribozymes in the presence of a mineral surface

Mineral surfaces are often proposed as the sites of critical processes in the emergence of life. Clay minerals in particular are thought to play significant roles in the origin of life including polymerizing, concentrating, organizing, and protecting biopolymers. In these scenarios, the impact of minerals on biopolymer folding is expected to influence evolutionary processes. These processes include both the initial emergence of functional structures in the presence of the mineral and the subsequent transition away from the mineral-associated niche. The initial evolution of function depends upon the number and distribution of sequences capable of functioning in the presence of the mineral, and the transition to new environments depends upon the overlap between sequences that evolve on the mineral surface and sequences that can perform the same functions in the mineral's absence. To examine these processes, we evolved self-cleaving ribozymes in vitro in the presence or absence of Na-saturated montmorillonite clay mineral particles. Starting from a shared population of random sequences, RNA populations were evolved in parallel, along separate evolutionary trajectories. Comparative sequence analysis and activity assays show that the impact of this clay mineral on functional structure selection was minimal; it neither prevented common structures from emerging, nor did it promote the emergence of new structures. This suggests that montmorillonite does not improve RNA's ability to evolve functional structures; however, it also suggests that RNAs that do evolve in contact with montmorillonite retain the same structures in mineral-free environments, potentially facilitating an evolutionary transition away from a mineral-associated niche.

Allosteric DNA nanoswitches for controlled release of a molecular cargo triggered by biological inputs

A rationally designed new class of DNA-based nanoswitches allosterically regulated by specific biological targets, antibodies and transcription factors, can load and release a molecular cargo in a controlled fashion.

Here we demonstrate the rational design of a new class of DNA-based nanoswitches which are allosterically regulated by specific biological targets, antibodies and transcription factors, and are able to load and release a molecular cargo (i.e. doxorubicin) in a controlled fashion. In our first model system we rationally designed a stem-loop DNA-nanoswitch that adopts two mutually exclusive conformations: a “Load” conformation containing a doxorubicin-intercalating domain and a “Release” conformation containing a duplex portion recognized by a specific transcription-factor (here Tata Binding Protein). The binding of the transcription factor pushes this conformational equilibrium towards the “Release” state thus leading to doxorubicin release from the nanoswitch. In our second model system we designed a similar stem-loop DNA-nanoswitch for which conformational change and subsequent doxorubicin release can be triggered by a specific antibody. Our approach augments the current tool kit of smart drug release mechanisms regulated by different biological inputs.

Cellular Delivery of RNA Nanoparticles

RNA nanostructures can be programmed to exhibit defined sizes, shapes and stoichiometries from naturally occurring or de novo designed RNA motifs. These constructs can be used as scaffolds to attach functional moieties, such as ligand binding motifs or gene expression regulators, for nanobiology applications. This review is focused on four areas of importance to RNA nanotechnology: the types of RNAs of particular interest for nanobiology, the assembly of RNA nanoconstructs, the challenges of cellular delivery of RNAs in vivo, and the delivery carriers that aid in the matter. The available strategies for the design of nucleic acid nanostructures, as well as for formulation of their carriers, make RNA nanotechnology an important tool in both basic research and applied biomedical science.

Chemistry and Biology of Self-Cleaving Ribozymes

Self-cleaving ribozymes were discovered 30 years ago, but their biological distribution and catalytic mechanisms are only beginning to be defined. Each ribozyme family is defined by a distinct structure, with unique active sites accelerating the same transesterification reaction across the families. Biochemical studies show that general acid-base catalysis is the most common mechanism of self-cleavage, but metal ions and metabolites can be used as cofactors. Ribozymes have been discovered in highly diverse genomic contexts throughout nature, from viroids to vertebrates. Their biological roles include self-scission during rolling-circle replication of RNA genomes, co-transcriptional processing of retrotransposons, and metabolite-dependent gene expression regulation in bacteria. Other examples, including highly conserved mammalian ribozymes, suggest that many new biological roles are yet to be discovered.

Single-Molecule Visualization of the Activity of a Zn(2+)-Dependent DNAzyme

We demonstrate the single-molecule imaging of the catalytic reaction of a Zn(2+)-dependent DNAzyme in a DNA origami nanostructure. The single-molecule catalytic activity of the DNAzyme was examined in the designed nanostructure, a DNA frame. The DNAzyme and a substrate strand attached to two supported dsDNA molecules were assembled in the DNA frame in two different configurations. The reaction was monitored by observing the configurational changes of the incorporated DNA strands in the DNA frame. This configurational changes were clearly observed in accordance with the progress of the reaction. The separation processes of the dsDNA molecules, as induced by the cleavage by the DNAzyme, were directly visualized by high-speed atomic force microscopy (AFM). This nanostructure-based AFM imaging technique is suitable for the monitoring of various chemical and biochemical catalytic reactions at the single-molecule level.

A general approach to the design of allosteric, transcription factor-regulated DNAzymes

Here we explore a general strategy for the rational design of nucleic acid catalysts that can be allosterically activated by specific nucleic-acid binding proteins. To demonstrate this we have combined a catalytic DNAzyme sequence and the consensus sequence recognized by specific transcription factors to create a construct exhibiting two low-energy conformations: a more stable conformation lacking catalytic activity and lacking the transcription factor binding site, and a less stable conformation that is both catalytically active and competent to bind the transcription factor. The presence of the target transcription factor pushes the equilibrium between these states towards the latter conformation, concomitantly activating catalysis. To demonstrate this we have designed and characterized two peroxidase-like DNAzymes whose activities are triggered upon binding either TATA binding protein or the microphthalmia-associated transcription factor. Our approach augments the current tool kit for the allosteric control of DNAzymes and ribozymes and, because transcription factors control many key biological functions, could have important clinical and diagnostic applications.

New classes of self-cleaving ribozymes revealed by comparative genomics analysis

Enzymes made of RNA catalyze reactions that are essential for protein synthesis and RNA processing. However, such natural ribozymes are exceedingly rare, as evident by the fact that the discovery rate for new classes has dropped to one per decade from about one per year during the 1980s. Indeed, only 11 distinct ribozyme classes have been experimentally validated to date. Recently, we recognized that self-cleaving ribozymes frequently associate with certain types of genes from bacteria. Herein this synteny was exploited to identify divergent architectures for two previously known ribozyme classes and to discover additional noncoding RNA motifs that are self-cleaving RNA candidates. Three new self-cleaving classes, named twister sister, pistol and hatchet, have been identified from this collection, suggesting that even more ribozymes remain hidden in modern cells.

In vitro evolution of distinct self-cleaving ribozymes in diverse environments

In vitro evolution experiments have long been used to evaluate the roles of RNA in both modern and ancient biology, and as a tool for biotechnology applications. The conditions under which these experiments have been conducted, however, do not reflect the range of cellular environments in modern biology or our understanding of chemical environments on the early earth, when the atmosphere and oceans were largely anoxic and soluble Fe²⁺ was abundant. To test the impact of environmental factors relevant to RNA's potential role in the earliest forms of life, we evolved populations of self-cleaving ribozymes in an anoxic atmosphere with varying pH in the presence of either Fe²⁺ or Mg²⁺. Populations evolved under these different conditions are dominated by different sequences and secondary structures, demonstrating global differences in the underlying fitness landscapes. Comparisons between evolutionary outcomes and catalytic activities also indicate that Mg²⁺ can readily take the place of Fe²⁺ in supporting the catalysis of RNA cleavage at neutral pH, but not at lower pH. These results highlight the importance of considering the specific environments in which functional biopolymers evolve when evaluating their potential roles in the origin of life, extant biology, or biotechnology.

Using droplet-based microfluidics to improve the catalytic properties of RNA under multiple-turnover conditions

In vitro evolution methodologies are powerful approaches to identify RNA with new functionalities. While Systematic Evolution of Ligands by Exponential enrichment (SELEX) is an efficient approach to generate new RNA aptamers, it is less suited for the isolation of efficient ribozymes as it does not select directly for the catalysis. In vitro compartmentalization (IVC) in aqueous droplets in emulsions allows catalytic RNAs to be selected under multiple-turnover conditions but suffers severe limitations that can be overcome using the droplet-based microfluidics workflow described in this paper. Using microfluidics, millions of genes in a library can be individually compartmentalized in highly monodisperse aqueous droplets and serial operations performed on them. This allows the different steps of the evolution process (gene amplification, transcription, and phenotypic assay) to be uncoupled, making the method highly flexible, applicable to the selection and evolution of a variety of RNAs, and easily adaptable for evolution of DNA or proteins. To demonstrate the method, we performed cycles of random mutagenesis and selection to evolve the X-motif, a ribozyme which, like many ribozymes selected using SELEX, has limited multiple-turnover activity. This led to the selection of variants, likely to be the optimal ribozymes that can be generated using point mutagenesis alone, with a turnover number under multiple-turnover conditions, k(ss) cat, ∼ 28-fold higher than the original X-motif, primarily due to an increase in the rate of product release, the rate-limiting step in the multiple-turnover reaction.

Evolution of the division of labor between genes and enzymes in the RNA world

The RNA world is a very likely interim stage of the evolution after the first replicators and before the advent of the genetic code and translated proteins. Ribozymes are known to be able to catalyze many reaction types, including cofactor-aided metabolic transformations. In a metabolically complex RNA world, early division of labor between genes and enzymes could have evolved, where the ribozymes would have been transcribed from the genes more often than the other way round, benefiting the encapsulating cells through this dosage effect. Here we show, by computer simulations of protocells harboring unlinked RNA replicators, that the origin of replicational asymmetry producing more ribozymes from a gene template than gene strands from a ribozyme template is feasible and robust. Enzymatic activities of the two modeled ribozymes are in trade-off with their replication rates, and the relative replication rates compared to those of complementary strands are evolvable traits of the ribozymes. The degree of trade-off is shown to have the strongest effect in favor of the division of labor. Although some asymmetry between gene and enzymatic strands could have evolved even in earlier, surface-bound systems, the shown mechanism in protocells seems inevitable and under strong positive selection. This could have preadapted the genetic system for transcription after the subsequent origin of chromosomes and DNA.

Local neutral networks help maintain inaccurately replicating ribozymes

The error threshold of replication limits the selectively maintainable genome size against recurrent deleterious mutations for most fitness landscapes. In the context of RNA replication a distinction between the genotypic and the phenotypic error threshold has been made; where the latter concerns the maintenance of secondary structure rather than sequence. RNA secondary structure is treated as a proxy for function. The phenotypic error threshold allows higher per digit mutation rates than its genotypic counterpart, and is known to increase with the frequency of neutral mutations in sequence space. Here we show that the degree of neutrality, i.e. the frequency of nearest-neighbour (one-step) neutral mutants is a remarkably accurate proxy for the overall frequency of such mutants in an experimentally verifiable formula for the phenotypic error threshold; this we achieve by the full numerical solution for the concentration of all sequences in mutation-selection balance up to length 16. We reinforce our previous result that currently known ribozymes could be selectively maintained by the accuracy known from the best available polymerase ribozymes. Furthermore, we show that in silico stabilizing selection can increase the mutational robustness of ribozymes due to the fact that they were produced by artificial directional selection in the first place. Our finding offers a better understanding of the error threshold and provides further insight into the plausibility of an ancient RNA world.

HDV family of self-cleaving ribozymes

The hepatitis delta virus (HDV) ribozymes are catalytic RNAs capable of cleaving their own sugar-phosphate backbone. The HDV virus possesses the ribozymes in both sense and antisense genomic transcripts, where they are essential for processing during replication. These ribozymes have been the subject of intense biochemical scrutiny and have yielded a wealth of mechanistic insights. In recent years, many HDV-like ribozymes have been identified in nearly all branches of life. The ribozymes are implicated in a variety of biological events, including episodic memory in mammals and retrotransposition in many eukaryotes. Detailed analysis of additional HDV-like ribozyme isolates will likely reveal many more biological functions and provide information about the evolution of this unique RNA.

A widespread self-cleaving ribozyme class is revealed by bioinformatics

Ribozymes are noncoding RNAs that promote chemical transformations with rate enhancements approaching those of protein enzymes. Although ribozymes are likely to have been abundant during the RNA world era, only ten classes are known to exist among contemporary organisms. We report the discovery and analysis of an additional self-cleaving ribozyme class, called twister, which is present in many species of bacteria and eukarya. Nearly 2,700 twister ribozymes were identified that conform to a secondary structure consensus that is small yet complex, with three stems conjoined by internal and terminal loops. Two pseudoknots provide tertiary structure contacts that are critical for catalytic activity. The twister ribozyme motif provides another example of a natural RNA catalyst and calls attention to the potentially varied biological roles of this and other classes of widely distributed self-cleaving RNAs.

Genome-wide probing of RNA structure reveals active unfolding of mRNA structures in vivo

RNA plays a dual role as an informational molecule and a direct effector of biological tasks. The latter function is enabled by RNA’s ability to adopt complex secondary and tertiary folds and thus has motivated extensive computational^1–2 and experimental^3–8 efforts for determining RNA structures. Existing approaches for evaluating RNA structure have been largely limited to in vitro systems, yet the thermodynamic forces which drive RNA folding in vitro may not be sufficient to predict stable RNA structures in vivo⁵. Indeed, the presence of RNA binding proteins and ATP-dependent helicases can influence which structures are present inside cells. Here we present an approach for globally monitoring RNA structure in native conditions in vivo with single nucleotide precision. This method is based on in vivo modification with dimethyl sulfate (DMS), which reacts with unpaired adenine and cytosine residues⁹, followed by deep sequencing to monitor modifications. Our data from yeast and mammalian cells are in excellent agreement with known mRNA structures and with the high-resolution crystal structure of the Saccharomyces cerevisiae ribosome¹⁰. Comparison between in vivo and in vitro data reveals that in rapidly dividing cells there are vastly fewer structured mRNA regions in vivo than in vitro. Even thermostable RNA structures are often denatured in cells, highlighting the importance of cellular processes in regulating RNA structure. Indeed, analysis of mRNA structure under ATP-depleted conditions in yeast reveals that energy-dependent processes strongly contribute to the predominantly unfolded state of mRNAs inside cells. Our studies broadly enable the functional analysis of physiological RNA structures and reveal that, in contrast to the Anfinsen view of protein folding, thermodynamics play an incomplete role in determining mRNA structure in vivo.

An in vitro evolved glmS ribozyme has the wild-type fold but loses coenzyme dependence

Uniquely among known ribozymes, the glmS ribozyme-riboswitch requires a small-molecule coenzyme, glucosamine-6-phosphate (GlcN6P). Although consistent with its gene-regulatory function, the use of GlcN6P is unexpected because all of the other characterized self-cleaving ribozymes use RNA functional groups or divalent cations for catalysis. To determine what active site features make this ribozyme reliant on GlcN6P and to evaluate whether it might have evolved from a coenzyme-independent ancestor, we isolated a GlcN6P-independent variant through in vitro selection. Three active site mutations suffice to generate a highly reactive RNA that adopts the wild-type fold but uses divalent cations for catalysis and is insensitive to GlcN6P. Biochemical and crystallographic comparisons of wild-type and mutant ribozymes show that a handful of functional groups fine-tune the RNA to be either coenzyme or cation dependent. These results indicate that a few mutations can confer new biochemical activities on structured RNAs. Thus, families of structurally related ribozymes with divergent function may exist.

Evolution of the R2 retrotransposon ribozyme and its self-cleavage site

R2 is a non-long terminal repeat retrotransposon that inserts site-specifically in the tandem 28S rRNA genes of many animals. Previously, R2 RNA from various species of Drosophila was shown to self-cleave from the 28S rRNA/R2 co-transcript by a hepatitis D virus (HDV)-like ribozyme encoded at its 5' end. RNA cleavage was at the precise 5' junction of the element with the 28S gene. Here we report that RNAs encompassing the 5' ends of R2 elements from throughout its species range fold into HDV-like ribozymes. In vitro assays of RNA self-cleavage conducted in many R2 lineages confirmed activity. For many R2s, RNA self-cleavage was not at the 5' end of the element but at 28S rRNA sequences up to 36 nucleotides upstream of the junction. The location of cleavage correlated well with the types of endogenous R2 5' junctions from different species. R2 5' junctions were uniform for most R2s in which RNA cleavage was upstream in the rRNA sequences. The 28S sequences remaining on the first DNA strand synthesized during retrotransposition are postulated to anneal to the target site and uniformly prime second strand DNA synthesis. In species where RNA cleavage occurred at the R2 5' end, the 5' junctions were variable. This junction variation is postulated to result from the priming of second strand DNA synthesis by chance microhomologies between the target site and the first DNA strand. Finally, features of R2 ribozyme evolution, especially changes in cleavage site and convergence on the same active site sequences, are discussed.

Small, highly active DNAs that hydrolyze DNA

DNA phosphoester bonds are exceedingly resistant to hydrolysis in the absence of chemical or enzymatic catalysts. This property is particularly important for organisms with large genomes, as resistance to hydrolytic degradation permits the long-term storage of genetic information. Here we report the creation and analysis of two classes of engineered deoxyribozymes that selectively and rapidly hydrolyze DNA. Members of class I deoxyribozymes carry a catalytic core composed of only 15 conserved nucleotides and attain an observed rate constant (k(obs)) of ~1 min(-1) when incubated near neutral pH in the presence of Zn(2+). Natural DNA sequences conforming to the class I consensus sequence and structure were found that undergo hydrolysis under selection conditions (2 mM Zn(2+), pH 7), which demonstrates that the inherent structure of certain DNA regions might promote catalytic reactions, leading to genomic instability.

In-silico design of computational nucleic acids for molecular information processing

Within recent years nucleic acids have become a focus of interest for prototype implementations of molecular computing concepts. During the same period the importance of ribonucleic acids as components of the regulatory networks within living cells has increasingly been revealed. Molecular computers are attractive due to their ability to function within a biological system; an application area extraneous to the present information technology paradigm. The existence of natural information processing architectures (predominately exemplified by protein) demonstrates that computing based on physical substrates that are radically different from silicon is feasible. Two key principles underlie molecular level information processing in organisms: conformational dynamics of macromolecules and self-assembly of macromolecules. Nucleic acids support both principles, and moreover computational design of these molecules is practicable. This study demonstrates the simplicity with which one can construct a set of nucleic acid computing units using a new computational protocol. With the new protocol, diverse classes of nucleic acids imitating the complete set of boolean logical operators were constructed. These nucleic acid classes display favourable thermodynamic properties and are significantly similar to the approximation of successful candidates implemented in the laboratory. This new protocol would enable the construction of a network of interconnecting nucleic acids (as a circuit) for molecular information processing.

Computational tests of a thermal cycling strategy to isolate more complex functional nucleic acid motifs from random sequence pools by in vitro selection

The dual information-function nature of nucleic acids has been exploited in the laboratory to isolate novel receptors and catalysts from random DNA and RNA sequences by cycles of in vitro selection and amplification. This strategy is particularly effective because, unlike polypeptides with random amino acid sequences, nucleic acids with random base sequences are often capable of stably folding into defined three-dimensional structures. However, the pervasive base-pairing potential of nucleic acids is also known to lead to kinetic traps in their folding landscapes. That is, the same DNA or RNA sequence can often adopt alternative base-paired structures that are local energy minima, and these folds may interconvert very slowly. We have used simulations with nucleic acid folding algorithms to evaluate the effect of misfolding on in vitro selection experiments. We demonstrate that kinetic traps can prevent the recovery of novel families of complex functional motifs by two mechanisms. First, misfolding can lead to the stochastic loss of unique sequences in the first round of selection. Second, frequent misfolding can reduce the average activity of multiple copies of a sequence to such an extent that it will be outcompeted after multiple rounds of selection. In these simulations, adding thermal cycling to sample multiple folds of one sequence during a selection for a self-modifying catalytic activity can improve the recovery of rare examples of more complex structures. Although newly isolated sequences may fold poorly, they can represent footholds in sequence space that can be improved to reliably fold after a few mutations. Thus, it is plausible that thermal cycling by day-night cycles or other mechanisms on the primordial earth may have been important for the evolution of the first RNA catalysts, and a fold sampling strategy might be used to search for more effective nucleic acid catalysts in the laboratory today.

Processing and translation initiation of non-long terminal repeat retrotransposons by hepatitis delta virus (HDV)-like self-cleaving ribozymes

Many non-long terminal repeat (non-LTR) retrotransposons lack internal promoters and are co-transcribed with their host genes. These transcripts need to be liberated before inserting into new loci. Using structure-based bioinformatics, we show that several classes of retrotransposons in phyla-spanning arthropods, nematodes, and chordates utilize self-cleaving ribozymes of the hepatitis delta virus (HDV) family for processing their 5' termini. Ribozyme-terminated retrotransposons include rDNA-specific R2, R4, and R6, telomere-specific SART, and Baggins and RTE. The self-scission of the R2 ribozyme is strongly modulated by the insertion site sequence in the rDNA, with the most common insertion sequences promoting faster processing. The ribozymes also promote translation initiation of downstream open reading frames in vitro and in vivo. In some organisms HDV-like and hammerhead ribozymes appear to be dedicated to processing long and short interspersed elements, respectively. HDV-like ribozymes serve several distinct functions in non-LTR retrotransposition, including 5' processing, translation initiation, and potentially trans-templating.

Identification of hammerhead ribozymes in all domains of life reveals novel structural variations

Hammerhead ribozymes are small self-cleaving RNAs that promote strand scission by internal phosphoester transfer. Comparative sequence analysis was used to identify numerous additional representatives of this ribozyme class than were previously known, including the first representatives in fungi and archaea. Moreover, we have uncovered the first natural examples of "type II" hammerheads, and our findings reveal that this permuted form occurs in bacteria as frequently as type I and III architectures. We also identified a commonly occurring pseudoknot that forms a tertiary interaction critical for high-speed ribozyme activity. Genomic contexts of many hammerhead ribozymes indicate that they perform biological functions different from their known role in generating unit-length RNA transcripts of multimeric viroid and satellite virus genomes. In rare instances, nucleotide variation occurs at positions within the catalytic core that are otherwise strictly conserved, suggesting that core mutations are occasionally tolerated or preferred.

Structure-based search reveals hammerhead ribozymes in the human microbiome

Deep sequencing of viral or bacterial nucleic acids monitors the presence and diversity of microbes in select populations and locations. Metagenomic study of mammalian viromes can help trace paths of viral transmissions within or between species. High throughput sequencing of patient and untreated sewage microbiomes showed many sequences with no similarity to genomic sequences of known function or origin. To estimate the distribution of functional RNAs in these microbiomes, we used the hammerhead ribozyme (HHR) motif to search for sequences capable of assuming its three-way junction fold. Although only two of the three possible natural HHR topologies had been known, our analysis revealed highly active ribozymes that terminated in any of the three stems. The most abundant of these are type II HHRs, one of which is the fastest natural cis-acting HHR yet discovered. Altogether, 13 ribozymes were confirmed in vitro, but only one showed sequence similarity to previously described HHRs. Sequences surrounding the ribozymes do not generally show similarity to known genes, except in one case, where a ribozyme is immediately preceded by a bacterial RadC gene. We demonstrate that a structure-based search for a known functional RNA is a powerful tool for analysis of metagenomic datasets, complementing sequence alignments.

Computational approaches for RNA energy parameter estimation

Methods for efficient and accurate prediction of RNA structure are increasingly valuable, given the current rapid advances in understanding the diverse functions of RNA molecules in the cell. To enhance the accuracy of secondary structure predictions, we developed and refined optimization techniques for the estimation of energy parameters. We build on two previous approaches to RNA free-energy parameter estimation: (1) the Constraint Generation (CG) method, which iteratively generates constraints that enforce known structures to have energies lower than other structures for the same molecule; and (2) the Boltzmann Likelihood (BL) method, which infers a set of RNA free-energy parameters that maximize the conditional likelihood of a set of reference RNA structures. Here, we extend these approaches in two main ways: We propose (1) a max-margin extension of CG, and (2) a novel linear Gaussian Bayesian network that models feature relationships, which effectively makes use of sparse data by sharing statistical strength between parameters. We obtain significant improvements in the accuracy of RNA minimum free-energy pseudoknot-free secondary structure prediction when measured on a comprehensive set of 2518 RNA molecules with reference structures. Our parameters can be used in conjunction with software that predicts RNA secondary structures, RNA hybridization, or ensembles of structures. Our data, software, results, and parameter sets in various formats are freely available at http://www.cs.ubc.ca/labs/beta/Projects/RNA-Params.

The RNA structure alignment ontology

Multiple sequence alignments are powerful tools for understanding the structures, functions, and evolutionary histories of linear biological macromolecules (DNA, RNA, and proteins), and for finding homologs in sequence databases. We address several ontological issues related to RNA sequence alignments that are informed by structure. Multiple sequence alignments are usually shown as two-dimensional (2D) matrices, with rows representing individual sequences, and columns identifying nucleotides from different sequences that correspond structurally, functionally, and/or evolutionarily. However, the requirement that sequences and structures correspond nucleotide-by-nucleotide is unrealistic and hinders representation of important biological relationships. High-throughput sequencing efforts are also rapidly making 2D alignments unmanageable because of vertical and horizontal expansion as more sequences are added. Solving the shortcomings of traditional RNA sequence alignments requires explicit annotation of the meaning of each relationship within the alignment. We introduce the notion of "correspondence," which is an equivalence relation between RNA elements in sets of sequences as the basis of an RNA alignment ontology. The purpose of this ontology is twofold: first, to enable the development of new representations of RNA data and of software tools that resolve the expansion problems with current RNA sequence alignments, and second, to facilitate the integration of sequence data with secondary and three-dimensional structural information, as well as other experimental information, to create simultaneously more accurate and more exploitable RNA alignments.

Evolutionary origins and directed evolution of RNA

In vitro selection experiments show first and foremost that it is possible that functional nucleic acids can arise from random sequence libraries. Indeed, even simple sequence and structural motifs can prove to be robust binding species and catalysts, indicating that it may have been possible to transition from even the earliest self-replicators to a nascent, RNA-catalyzed metabolism. Because of the diversity of aptamers and ribozymes that can be selected, it is possible to construct a 'fossil record' of the evolution of the RNA world, with in vitro selected catalysts filling in as doppelgangers for molecules long gone. In this way a plausible pathway from simple oligonucleotide replicators to genomic polymerases can be imagined, as can a pathway from basal ribozyme activities to the ribosome. Most importantly, though, in vitro selection experiments can give a true and quantitative idea of the likelihood that these scenarios could have played out in the RNA world. Simple binding species and catalysts could have evolved into other structures and functions. As replicating sequences grew longer, new, more complex functions or faster catalytic activities could have been accessed. Some activities may have been isolated in sequence space, but others could have been approached along large, interconnected neutral networks. As the number, type, and length of ribozymes increased, RNA genomes would have evolved and eventually there would have been no area in a fitness landscape that would have been inaccessible. Self-replication would have inexorably led to life.

Searching genomes for ribozymes and riboswitches

A discussion of experimental approaches and theoretical difficulties in the identification of ribozymes with novel catalytic functions.

New regulatory RNAs with complex structures have recently been discovered, among them the first catalytic riboswitch, a gene-regulatory RNA sequence with catalytic activity. Here we discuss some of the experimental approaches and theoretical difficulties attached to the identification of new ribozymes in genomes.

Multiple pattern matching: a Markov chain approach

RNA motifs typically consist of short, modular patterns that include base pairs formed within and between modules. Estimating the abundance of these patterns is of fundamental importance for assessing the statistical significance of matches in genomewide searches, and for predicting whether a given function has evolved many times in different species or arose from a single common ancestor. In this manuscript, we review in an integrated and self-contained manner some basic concepts of automata theory, generating functions and transfer matrix methods that are relevant to pattern analysis in biological sequences. We formalize, in a general framework, the concept of Markov chain embedding to analyze patterns in random strings produced by a memoryless source. This conceptualization, together with the capability of automata to recognize complicated patterns, allows a systematic analysis of problems related to the occurrence and frequency of patterns in random strings. The applications we present focus on the concept of synchronization of automata, as well as automata used to search for a finite number of keywords (including sets of patterns generated according to base pairing rules) in a general text.

Evolution of High-Branching Deoxyribozymes from a Catalytic DNA with a Three-Way Junction

Here, we report the evolution of two star-shaped (five-way junction) deoxyribozymes from a catalytic DNA containing a three-way junction scaffold. The transition was shown to be a switch rather than a gradual progression. The star-shaped motifs, surprisingly, only took five selection cycles to be detected, and another four to dominate the evolving population. Chemical probing experiments indicated that the two deoxyribozymes belong to the same family despite noticeable variations in both the primary sequence and the secondary structure. Our findings not only describe the evolution of high-branching nucleic acid structures from a low-branching catalytic module, but they also illustrate the idea of deriving a rare structural motif by sampling the sequence variants of a given functional nucleic acid.

RNA synthetic biology

RNA molecules play important and diverse regulatory roles in the cell by virtue of their interaction with other nucleic acids, proteins and small molecules. Inspired by this natural versatility, researchers have engineered RNA molecules with new biological functions. In the last two years efforts in synthetic biology have produced novel, synthetic RNA components capable of regulating gene expression in vivo largely in bacteria and yeast, setting the stage for scalable and programmable cellular behavior. Immediate challenges for this emerging field include determining how computational and directed-evolution techniques can be implemented to increase the complexity of engineered RNA systems, as well as determining how such systems can be broadly extended to mammalian systems. Further challenges include designing RNA molecules to be sensors of intracellular and environmental stimuli, probes to explore the behavior of biological networks and components of engineered cellular control systems.

Characterization of long RNA-cleaving deoxyribozymes with short catalytic cores: the effect of excess sequence elements on the outcome of in vitro selection

We previously conducted an in vitro selection experiment for RNA-cleaving deoxyribozymes, using a combinatorial DNA library containing 80 random nucleotides. Ultimately, 110 different sequence classes were isolated, but the vast majority contained a short14-15 nt catalytic DNA motif commonly known as 8-17. Herein, we report extensive truncation experiments conducted on multiple sequence classes to confirm the suspected catalytic role played by 8-17 and to determine the effect of excess sequence elements on the activity of this motif and the outcome of selection. Although we observed beneficial, detrimental and neutral consequences for activity, the magnitude of the effect rarely exceeded 2-fold. These deoxyribozymes appear to have survived increasing selection pressure despite the presence of additional sequence elements, rather than because of them. A new deoxyribozyme with comparable activity, called G15-30, was approximately 2.5-fold larger and experienced a approximately 4-fold greater inhibitory effect from excess sequence elements than the average 8-17 motif. Our results suggest that 8-17 may be less susceptible to the potential inhibitory effects of excess arbitrary sequence than larger motifs, which represents a previously unappreciated selective advantage that may contribute to its widespread recurrence.