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

Identification of Hammerhead-variant ribozyme sequences in SARS-CoV-2

The SARS-CoV-2 RNA virus and variants, responsible for the COVID-19 pandemic has become endemic, raised a need for further understanding of the viral genome and biology. Despite vast research on SARS-CoV-2, no ribozymes have been found in the virus genome. Here we report the identification of 39 Hammerhead-variant ribozyme sequences (CoV-HHRz) in SARS-CoV-2. These sequences are highly conserved within SARS-CoV-2 variants but show large diversity among other coronaviruses. In vitro CoV-HHRz sequences possess the characteristics of typical ribozymes; cleavage is pH and ion dependent, although their activity is relatively low and Mn2+ is required for cleavage. The cleavage sites of four CoV-HHRz coincide with the breakpoint of expressed subgenomic RNA (sgRNAs) in SARS-CoV-2 transcriptome data suggesting in vivo activity. The CoV-HHRz are involved in processing sgRNAs for ORF7b, ORF 10 and ORF1ab nsp13 which are essential for viral packaging and life cycle.

Identification of HDV-like theta ribozymes involved in tRNA-based recoding of gut bacteriophages

Trillions of microorganisms, collectively known as the microbiome, inhabit our bodies with the gut microbiome being of particular interest in biomedical research. Bacteriophages, the dominant virome constituents, can utilize suppressor tRNAs to switch to alternative genetic codes (e.g., the UAG stop-codon is reassigned to glutamine) while infecting hosts with the standard bacterial code. However, what triggers this switch and how the bacteriophage manipulates its host is poorly understood. Here, we report the discovery of a subgroup of minimal hepatitis delta virus (HDV)-like ribozymes - theta ribozymes - potentially involved in the code switch leading to the expression of recoded lysis and structural phage genes. We demonstrate their HDV-like self-scission behavior in vitro and find them in an unreported context often located with their cleavage site adjacent to tRNAs, indicating a role in viral tRNA maturation and/or regulation. Every fifth associated tRNA is a suppressor tRNA, further strengthening our hypothesis. The vast abundance of tRNA-associated theta ribozymes - we provide 1753 unique examples - highlights the importance of small ribozymes as an alternative to large enzymes that usually process tRNA 3'-ends. Our discovery expands the short list of biological functions of small HDV-like ribozymes and introduces a previously unknown player likely involved in the code switch of certain recoded gut bacteriophages.

Flanking Sequence Cotranscriptionally Regulates Twister Ribozyme Activity

Small nucleolytic ribozymes are RNAs that cleave their own phosphodiester backbone. While proteinaceous enzymes are regulated by a variety of known mechanisms, methods of regulation for ribozymes remain unclear. Twister is one ribozyme class for which many structural and catalytic properties have been elucidated. However, few studies have analyzed the activity of twister ribozymes in the context of a native flanking sequence, even though ribozymes as transcribed in nature do not exist in isolation. Interactions between the ribozyme and its neighboring sequences can induce conformational changes that inhibit self-cleavage, providing a regulatory mechanism that could naturally determine ribozyme activity in vivo and in synthetic applications. To date, eight twister ribozymes have been identified within the staple crop rice (Oryza sativa). Herein, we select several twister ribozymes from rice and show that they are differentially regulated by their flanking sequence using published RNA-seq data sets, structure probing, and cotranscriptional cleavage assays. We found that the Osa 1-2 ribozyme does not interact with its flanking sequences. However, sequences flanking the Osa 1-3 and Osa 1-8 ribozymes form inactive conformations, referred to here as "ribozymogens", that attenuate ribozyme self-cleavage activity. For the Osa 1-3 ribozyme, we show that activity can be rescued upon addition of a complementary antisense oligonucleotide, suggesting ribozymogens can be controlled via external signals. In all, our data provide a plausible mechanism wherein flanking sequence differentially regulates ribozyme activity in vivo. More broadly, the ability to regulate ribozyme behavior locally has potential applications in control of gene expression and synthetic biology.

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.

Evolution of the RNA world: From signals to codes

The accumulated material in evolutionary biology, greatly enhanced by the achievements of modern synthetic biology, allows us to envision certain key hypothetical stages of prebiotic (chemical) evolution. This is often understood as the further evolution in the RNA World towards the RNA-protein World. It is a path towards the emergence of translation and the genetic code (I), signaling pathways with signaling molecules (II), and the appearance of RNA-based components of future gene regulatory networks (III). We believe that these evolutionary paths can be constructively viewed from the perspective of the concept of biological codes (Barbieri, 2003). Crucial evolutionary events in these directions would involve the emergence of RNA-based adaptors. Such adaptors connect two families of functionally and chemically distinct molecules into one functional entity. The emergence of primitive translation processes is undoubtedly the major milestone in the evolutionary path towards modern life. The key aspect here is the appearance of adaptors between amino acids and their cognate triplet codons. The initial steps are believed to involve the emergence of proto-transfer RNAs capable of self-aminoacylation. The second significant evolutionary breakthrough is the development of biochemical regulatory networks based on signaling molecules of the RNA World (ribonucleotides and their derivatives), as well as receptors and effectors (riboswitches) for these messengers. Some authors refer to this as the "lost language of the RNA World." The third evolutionary step is the emergence of signal sequences for ribozymes on the molecules of their RNA targets. This level of regulation in the RNA World is comparable to the gene regulatory networks of modern organisms. We believe that the signal sequences on target molecules have been rediscovered and developed by evolution into the gene regulatory networks of modern cells. In conclusion, the immense diversity of modern biological codes, in some of its key characteristics, can be traced back to the achievements of prebiotic evolution.

Astrovirology: how viruses enhance our understanding of life in the Universe

Viruses are the most numerically abundant biological entities on Earth. As ubiquitous replicators of molecular information and agents of community change, viruses have potent effects on the life on Earth, and may play a critical role in human spaceflight, for life-detection missions to other planetary bodies and planetary protection. However, major knowledge gaps constrain our understanding of the Earth's virosphere: (1) the role viruses play in biogeochemical cycles, (2) the origin(s) of viruses and (3) the involvement of viruses in the evolution, distribution and persistence of life. As viruses are the only replicators that span all known types of nucleic acids, an expanded experimental and theoretical toolbox built for Earth's viruses will be pivotal for detecting and understanding life on Earth and beyond. Only by filling in these knowledge and technical gaps we will obtain an inclusive assessment of how to distinguish and detect life on other planetary surfaces. Meanwhile, space exploration requires life-support systems for the needs of humans, plants and their microbial inhabitants. Viral effects on microbes and plants are essential for Earth's biosphere and human health, but virus-host interactions in spaceflight are poorly understood. Viral relationships with their hosts respond to environmental changes in complex ways which are difficult to predict by extrapolating from Earth-based proxies. These relationships should be studied in space to fully understand how spaceflight will modulate viral impacts on human health and life-support systems, including microbiomes. In this review, we address key questions that must be examined to incorporate viruses into Earth system models, life-support systems and life detection. Tackling these questions will benefit our efforts to develop planetary protection protocols and further our understanding of viruses in astrobiology.

Using 5' 32 P-labeled Primer and Reverse Transcription to Probe RNA Structure

RNA molecules perform numerous cellular functions necessary for cell viability, some of which can depend on the RNA's structure. Therefore, it is important to study RNA structure and RNA folding to better understand the molecular basis of these functions. RNA chemical mapping strategies to elucidate RNA structural changes involve using chemical reagents that form adducts or cleave RNA. Selective 2'-hydroxyl acylation analyzed by primer extension (SHAPE) measures RNA flexibility by modification of the 2' hydroxyl groups of flexible nucleotides. These RNA adducts can be detected using ³² P-labeled primers and reverse transcription (RT) followed by PAGE analysis. This strategy can reveal the base-paired regions of the RNA and provide insight into tertiary structure and solvent accessibility. This protocol provides a method to interrogate RNA structure using furoyl acylimidazole (FAI). © 2023 Wiley Periodicals LLC. Basic Protocol 1: Reverse transcription (RT) primer labeling with ³² P radionuclide Basic Protocol 2: Characterization of RNA structure with radiolabeled primer and reverse transcription (RT).

Rapid Kinetics of Pistol Ribozyme: Insights into Limits to RNA Catalysis

Pistol ribozyme (Psr) is a distinct class of small endonucleolytic ribozymes, which are important experimental systems for defining fundamental principles of RNA catalysis and designing valuable tools in biotechnology. High-resolution structures of Psr, extensive structure-function studies, and computation support a mechanism involving one or more catalytic guanosine nucleobases acting as a general base and divalent metal ion-bound water acting as an acid to catalyze RNA 2'-O-transphosphorylation. Yet, for a wide range of pH and metal ion concentrations, the rate of Psr catalysis is too fast to measure manually and the reaction steps that limit catalysis are not well understood. Here, we use stopped-flow fluorescence spectroscopy to evaluate Psr temperature dependence, solvent H/D isotope effects, and divalent metal ion affinity and specificity unconstrained by limitations due to fast kinetics. The results show that Psr catalysis is characterized by small apparent activation enthalpy and entropy changes and minimal transition state H/D fractionation, suggesting that one or more pre-equilibrium steps rather than chemistry is rate limiting. Quantitative analyses of divalent ion dependence confirm that metal aquo ion pK_a correlates with higher rates of catalysis independent of differences in ion binding affinity. However, ambiguity regarding the rate-limiting step and similar correlation with related attributes such as ionic radius and hydration free energy complicate a definitive mechanistic interpretation. These new data provide a framework for further interrogation of Psr transition state stabilization and show how thermal instability, metal ion insolubility at optimal pH, and pre-equilibrium steps such as ion binding and folding limit the catalytic power of Psr suggesting potential strategies for further optimization.

Catalytic mechanism and pH dependence of a methyltransferase ribozyme (MTR1) from computational enzymology

A methyltransferase ribozyme (MTR1) was selected in vitro to catalyze alkyl transfer from exogenous O6-methylguanine (O6mG) to a target adenine N1, and recently, high-resolution crystal structures have become available. We use a combination of classical molecular dynamics, ab initio quantum mechanical/molecular mechanical (QM/MM) and alchemical free energy (AFE) simulations to elucidate the atomic-level solution mechanism of MTR1. Simulations identify an active reactant state involving protonation of C10 that hydrogen bonds with O6mG:N1. The deduced mechanism involves a stepwise mechanism with two transition states corresponding to proton transfer from C10:N3 to O6mG:N1 and rate-controlling methyl transfer (19.4  kcal·mol-1 barrier). AFE simulations predict the pKa for C10 to be 6.3, close to the experimental apparent pKa of 6.2, further implicating it as a critical general acid. The intrinsic rate derived from QM/MM simulations, together with pKa calculations, enables us to predict an activity-pH profile that agrees well with experiment. The insights gained provide further support for a putative RNA world and establish new design principles for RNA-based biochemical tools.

Signal amplification strategies in biosensing of extracellular vesicles (EVs)

Extracellular vesicles (EVs) are membrane-enclosed vesicles secreted from mammalian cells. EVs act as multicomponent delivery vehicles to carry a wide variety of biological molecular information and participate in intercellular communications. Since elevated levels of EVs are associated with some pathological states such as inflammatory diseases and cancers, probing circulating EVs holds a great potential for early diagnostics. To this end, several detection methods have been developed in which biosensors have attracted great attentions in identification of EVs due to their simple instrumentation, versatile design and portability for point-of-care applications. The concentrations of EVs in bodily fluids are extremely low (i.e. 1-100 per μl) at early stages of a disease, which necessitates the use of signal amplification strategies for EVs detection. In this way, this review presents and discusses various amplification strategies for EVs biosensors based on detection modalities including surface plasmon resonance (SPR), calorimetry, fluorescence, electrochemical and electrochemiluminescence (ECL). In addition, microfluidic systems employed for signal amplification are reviewed and discussed in terms of their design and integration with the detection methods.

Tapping the potential of synthetic riboswitches: reviewing the versatility of the tetracycline aptamer

Synthetic riboswitches are a versatile class of regulatory elements that are becoming increasingly established in synthetic biology applications. They are characterized by their compact size and independence from auxiliary protein factors. While naturally occurring riboswitches were mostly discovered in bacteria, synthetic riboswitches have been designed for all domains of life. Published design strategies far exceed the number of riboswitches found in nature. A core element of any riboswitch is a binding domain, called an aptamer, which is characterized by high specificity and affinity for its ligand. Aptamers can be selected de novo, allowing the design of synthetic riboswitches against a broad spectrum of targets. The tetracycline aptamer has proven to be well suited for riboswitch engineering. Since its selection, it has been used in a variety of applications and is considered to be well established and characterized. Using the tetracycline aptamer as an example, we aim to discuss a large variety of design approaches for synthetic riboswitch engineering and their application. We aim to demonstrate the versatility of riboswitches in general and the high potential of synthetic RNA devices for creating new solutions in both the scientific and medical fields.

Novel Viroid-Like RNAs Naturally Infect a Filamentous Fungus

To date, viroids have been found to naturally infect only plants, resulting in substantial losses for some crops. Whether viroids or viroid-like RNAs naturally infect non-plant hosts remains unknown. Here the existence of a set of exogenous, single-stranded circular RNAs, ranging in size from 157 to 450 nucleotides, isolated from the fungus Botryosphaeria dothidea and nominated B. dothidea RNAs (BdcRNAs) is reported. BdcRNAs replicate autonomously in the nucleus via a rolling-circle mechanism following a symmetric pathway. BdcRNA infection induces symptoms, because BdcRNAs can apparently modulate, to different degrees, specific biological traits (e.g., alter morphology, decrease growth rate, attenuate virulence, and increase or decrease tolerance to osmotic and oxidative stress) of the host fungus. Overall, BdcRNAs have genome characteristics similar to those of viroids and exhibit pathogenic effects on fungal hosts. It is proposed that these novel fungus infecting RNAs should be termed mycoviroids. BdcRNA(s) may be considered additional inhabitants at the frontier of life in terms of genomic complexity, and represent a new class of acellular entities endowed with regulatory functions, and novel epigenomic carriers of biological information.

The Impact of DFT Functional, Cluster Model Size, and Implicit Solvation on the Structural Description of Single-Metal-Mediated DNA Phosphodiester Bond Cleavage: The Case Study of APE1

Phosphodiester bond hydrolysis in nucleic acids is a ubiquitous reaction that can be facilitated by enzymes called nucleases, which often use metal ions to achieve catalytic function. While a two-metal-mediated pathway has been well established for many enzymes, there is growing support that some enzymes require only one metal for the catalytic step. Using human apurinic/apyrimidinic endonuclease (APE1) as a prototypical example and cluster models, this study clarifies the impact of DFT functional, cluster model size, and implicit solvation on single-metal-mediated phosphodiester bond cleavage and provides insight into how to efficiently model this chemistry. Initially, a model containing 69 atoms built from a high-resolution X-ray crystal structure is used to explore the reaction pathway mapped by a range of DFT functionals and basis sets, which provides support for the use of standard functionals (M06-2X and B3LYP-D3) to study this reaction. Subsequently, systematically increasing the model size to 185 atoms by including additional amino acids and altering residue truncation points highlights that small models containing only a few amino acids or β carbon truncation points introduce model strains and lead to incorrect metal coordination. Indeed, a model that contains all key residues (general base and acid, residues that stabilize the substrate, and amino acids that maintain the metal coordination) is required for an accurate structural depiction of the one-metal-mediated phosphodiester bond hydrolysis by APE1, which results in 185 atoms. The additional inclusion of the broader enzyme environment through continuum solvation models has negligible effects. The insights gained in the present work can be used to direct future computational studies of other one-metal-dependent nucleases to provide a greater understanding of how nature achieves this difficult chemistry.

RNA World Modeling: A Comparison of Two Complementary Approaches

Simple Summary

Despite years of dedicated research, scientists are still not sure what the first ”living” cell would have looked like. One of the most well-known hypotheses is the RNA world hypothesis, which assumes that, in the beginning, life relied on RNA molecules instead of DNA as information carriers and primitive enzymes. The population of such RNAs is made up of self-replicating molecules (replicases) that could make copies of themselves and parasite molecules that could only be copied by replicases. In this study, we further investigated the interplay between these hypothetical prebiotic RNA species, since it plays a crucial role in generating diversity and complexity in prebiotic molecular evolution. We compared two approaches that are commonly used to investigate such simple prebiotic systems, representing different modeling and observation scales—namely, microscopic and macroscopic. In both cases, we were able to obtain consistent results.

Abstract

The origin of life remains one of the major scientific questions in modern biology. Among many hypotheses aiming to explain how life on Earth started, RNA world is probably the most extensively studied. It assumes that, in the very beginning, RNA molecules served as both enzymes and as genetic information carriers. However, even if this is true, there are many questions that still need to be answered—for example, whether the population of such molecules could achieve stability and retain genetic information for many generations, which is necessary in order for evolution to start. In this paper, we try to answer this question based on the parasite–replicase model (RP model), which divides RNA molecules into enzymes (RNA replicases) capable of catalyzing replication and parasites that do not possess replicase activity but can be replicated by RNA replicases. We describe the aforementioned system using partial differential equations and, based on the analysis of the simulation, surmise general rules governing its evolution. We also compare this approach with one where the RP system is modeled and implemented using a multi-agent modeling technique. We show that approaching the description and analysis of the RP system from different perspectives (microscopic represented by MAS and macroscopic depicted by PDE) provides consistent results. Therefore, applying MAS does not lead to erroneous results and allows us to study more complex situations where many cases are concerned, which would not be possible through the PDE model.

Overview of Methods for Large-Scale RNA Synthesis

In recent years, it has become clear that RNA molecules are involved in almost all vital cellular processes and pathogenesis of human disorders. The functional diversity of RNA comes from its structural richness. Although composed of only four nucleotides, RNA molecules present a plethora of secondary and tertiary structures critical for intra and intermolecular contacts with other RNAs and ligands (proteins, small metabolites, etc.). In order to fully understand RNA function it is necessary to define its spatial structure. Crystallography, nuclear magnetic resonance and cryogenic electron microscopy have demonstrated considerable success in determining the structures of biologically important RNA molecules. However, these powerful methods require large amounts of sample. Despite their limitations, chemical synthesis and in vitro transcription are usually employed to obtain milligram quantities of RNA for structural studies, delivering simple and effective methods for large-scale production of homogenous samples. The aim of this paper is to provide an overview of methods for large-scale RNA synthesis with emphasis on chemical synthesis and in vitro transcription. We also present our own results of testing the efficiency of these approaches in order to adapt the material acquisition strategy depending on the desired RNA construct.

Application of RtcB ligase to monitor self-cleaving ribozyme activity by RNA-seq

Self-cleaving ribozymes are catalytic RNAs and can be found in all domains of life. They catalyze a site-specific cleavage that results in a 5' fragment with a 2',3' cyclic phosphate (2',3' cP) and a 3' fragment with a 5' hydroxyl (5' OH) end. Recently, several strategies to enrich self-cleaving ribozymes by targeted biochemical methods have been introduced by us and others. Here, we develop an alternative strategy in which 5' OH RNAs are specifically ligated by RtcB ligase, which first guanylates the 3' phosphate of the adapter and then ligates it directly to RNAs with 5' OH ends. Our results demonstrate that adapter ligation to highly structured ribozyme fragments is much more efficient using the thermostable RtcB ligase from Pyrococcus horikoshii than the broadly applied Escherichia coli enzyme. Moreover, we investigated DNA, RNA and modified RNA adapters for their suitability in RtcB ligation reactions. We used the optimized RtcB-mediated ligation to produce RNA-seq libraries and captured a spiked 3' twister ribozyme fragment from E. coli total RNA. This RNA-seq-based method is applicable to detect ribozyme fragments as well as other cellular RNAs with 5' OH termini from total RNA.

cyPhyRNA-seq: a genome-scale RNA-seq method to detect active self-cleaving ribozymes by capturing RNAs with 2',3' cyclic phosphates and 5' hydroxyl ends

Self-cleaving ribozymes are catalytically active RNAs that cleave themselves into a 5′-fragment with a 2′,3′-cyclic phosphate and a 3′-fragment with a 5′-hydroxyl. They are widely applied for the construction of synthetic RNA devices and RNA-based therapeutics. However, the targeted discovery of self-cleaving ribozymes remains a major challenge. We developed a transcriptome-wide method, called cyPhyRNA-seq, to screen for ribozyme cleavage fragments in total RNA extract. This approach employs the specific ligation-based capture of ribozyme 5′-fragments using a variant of the Arabidopsis thaliana tRNA ligase we engineered. To capture ribozyme 3′-fragments, they are enriched from total RNA by enzymatic treatments. We optimized and enhanced the individual steps of cyPhyRNA-seq in vitro and in spike-in experiments. Then, we applied cyPhyRNA-seq to total RNA isolated from the bacterium Desulfovibrio vulgaris and detected self-cleavage of the three predicted type II hammerhead ribozymes, whose activity had not been examined to date. cyPhyRNA-seq can be used for the global analysis of active self-cleaving ribozymes with the advantage to capture both ribozyme cleavage fragments from total RNA. Especially in organisms harbouring many self-cleaving RNAs, cyPhyRNA-seq facilitates the investigation of cleavage activity. Moreover, this method has the potential to be used to discover novel self-cleaving ribozymes in different organisms.

Enzyme Nomenclature and Classification: the State of the Art

The IUBMB Enzyme classification system, available at the IUBMB ExplorEnz website, uses a four component number (the EC number) that identifies an enzyme in terms of reaction catalysed. There were originally six recognised groups of enzymes: Oxidoreductases (EC 1), Transferases (EC 2), Hydrolases, Lyases (EC 4), Isomerases (EC 5) and Ligases (EC 6). Of these the lyases, which are defined as "enzymes that cleave C-C, C-O, C-N and other bonds by means other than by hydrolysis or oxidation" present particular recognition and classification problems. Recently a new class, the Translocases (EC 7) has been added, which incorporates enzymes that catalyse the movement of ions or molecules across membranes or their separation within membranes. A new subclass of the isomerases has also been included for those enzymes that alter the conformations of proteins and nucleic acids. Newly reported enzymes are being regularly added to the list after validation and where new information affects the classification of an existing entry, a new EC number is created, but the old one is not re-used.

Streamlined CRISPR genome engineering in wild-type bacteria using SIBR-Cas

CRISPR-Cas is a powerful tool for genome editing in bacteria. However, its efficacy is dependent on host factors (such as DNA repair pathways) and/or exogenous expression of recombinases. In this study, we mitigated these constraints by developing a simple and widely applicable genome engineering tool for bacteria which we termed SIBR-Cas (Self-splicing Intron-Based Riboswitch-Cas). SIBR-Cas was generated from a mutant library of the theophylline-dependent self-splicing T4 td intron that allows for tight and inducible control over CRISPR-Cas counter-selection. This control delays CRISPR-Cas counter-selection, granting more time for the editing event (e.g. by homologous recombination) to occur. Without the use of exogenous recombinases, SIBR-Cas was successfully applied to knock-out several genes in three wild-type bacteria species (Escherichia coli MG1655, Pseudomonas putida KT2440 and Flavobacterium IR1) with poor homologous recombination systems. Compared to other genome engineering tools, SIBR-Cas is simple, tightly regulated and widely applicable for most (non-model) bacteria. Furthermore, we propose that SIBR can have a wider application as a simple gene expression and gene regulation control mechanism for any gene or RNA of interest in bacteria.

The function of twister ribozyme variants in non-LTR retrotransposition in Schistosoma mansoni

The twister ribozyme is widely distributed over numerous organisms and is especially abundant in Schistosoma mansoni, but has no confirmed biological function. Of the 17 non-LTR retrotransposons known in S. mansoni, none have thus far been associated with ribozymes. Here we report the identification of novel twister variant (T-variant) ribozymes and their function in S. mansoni non-LTR retrotransposition. We show that T-variant ribozymes are located at the 5′ end of Perere-3 non-LTR retrotransposons in the S. mansoni genome. T-variant ribozymes were demonstrated to be catalytically active in vitro. In reporter constructs, T-variants were shown to cleave in vivo, and cleavage of T-variants was sufficient for the translation of downstream reporter genes. Our analysis shows that the T-variants and Perere-3 are transcribed together. Target site duplications (TSDs); markers of target-primed reverse transcription (TPRT) and footmarks of retrotransposition, are located adjacent to the T-variant cleavage site and suggest that T-variant cleavage has taken place inS. mansoni. Sequence heterogeneity in the TSDs indicates that Perere-3 retrotransposition is not site-specific. The TSD sequences contribute to the 5′ end of the terminal ribozyme helix (P1 stem). Based on these results we conclude that T-variants have a functional role in Perere-3 retrotransposition.

Self-cleaving ribozymes: substrate specificity and synthetic biology applications

Various self-cleaving ribozymes appearing in nature catalyze the sequence-specific intramolecular cleavage of RNA and can be engineered to catalyze cleavage of appropriate substrates in an intermolecular fashion, thus acting as true catalysts. The mechanisms of the small, self-cleaving ribozymes have been extensively studied and reviewed previously. Self-cleaving ribozymes can possess high catalytic activity and high substrate specificity; however, substrate specificity is also engineerable within the constraints of the ribozyme structure. While these ribozymes share a common fundamental catalytic mechanism, each ribozyme family has a unique overall architecture and active site organization, indicating that several distinct structures yield this chemical activity. The multitude of catalytic structures, combined with some flexibility in substrate specificity within each family, suggests that such catalytic RNAs, taken together, could access a wide variety of substrates. Here, we give an overview of 10 classes of self-cleaving ribozymes and capture what is understood about their substrate specificity and synthetic applications. Evolution of these ribozymes in an RNA world might be characterized by the emergence of a new ribozyme family followed by rapid adaptation or diversification for specific substrates.

Self-cleaving ribozymes have become important tools of synthetic biology. Here we summarize the substrate specificity and applications of the main classes of these ribozymes.

Satellite-Like W-Elements: Repetitive, Transcribed, and Putative Mobile Genetic Factors with Potential Roles for Biology and Evolution of Schistosoma mansoni

A large portion of animal and plant genomes consists of noncoding DNA. This part includes tandemly repeated sequences and gained attention because it offers exciting insights into genome biology. We investigated satellite-DNA elements of the platyhelminth Schistosoma mansoni, a parasite with remarkable biological features. Schistosoma mansoni lives in the vasculature of humans causing schistosomiasis, a disease of worldwide importance. Schistosomes are the only trematodes that have evolved separate sexes, and the sexual maturation of the female depends on constant pairing with the male. The schistosome karyotype comprises eight chromosome pairs, males are homogametic (ZZ) and females are heterogametic (ZW). Part of the repetitive DNA of S. mansoni are W-elements (WEs), originally discovered as female-specific satellite DNAs in the heterochromatic block of the W-chromosome. Based on new genome and transcriptome data, we performed a reanalysis of the W-element families (WEFs). Besides a new classification of 19 WEFs, we provide first evidence for stage-, sex-, pairing-, gonad-, and strain-specific/preferential transcription of WEs as well as their mobile nature, deduced from autosomal copies of full-length and partial WEs. Structural analyses suggested roles as sources of noncoding RNA-like hammerhead ribozymes, for which we obtained functional evidence. Finally, the variable WEF occurrence in different schistosome species revealed remarkable divergence. From these results, we propose that WEs potentially exert enduring influence on the biology of S. mansoni. Their variable occurrence in different strains, isolates, and species suggests that schistosome WEs may represent genetic factors taking effect on variability and evolution of the family Schistosomatidae.

Mechanisms of catalytic RNA molecules

Ribozymes are folded catalytic RNA molecules that perform important biological functions. Since the discovery of the first RNA with catalytic activity in 1982, a large number of ribozymes have been reported. While most catalytic RNA molecules act alone, some RNA-based catalysts, such as RNase P, the ribosome, and the spliceosome, need protein components to perform their functions in the cell. In the last decades, the structure and mechanism of several ribozymes have been studied in detail. Aside from the ribosome, which catalyzes peptide bond formation during protein synthesis, the majority of known ribozymes carry out mostly phosphoryl transfer reactions, notably trans-esterification or hydrolysis reactions. In this review, we describe the main features of the mechanisms of various types of ribozymes that can function with or without the help of proteins to perform their biological functions.

Identification of over 200-fold more hairpin ribozymes than previously known in diverse circular RNAs

Self-cleaving ribozymes are catalytic RNAs that cut themselves at a specific inter-nucleotide linkage. They serve as a model of RNA catalysis, and as an important tool in biotechnology. For most of the nine known structural classes of self-cleaving ribozymes, at least hundreds of examples are known, and some are present in multiple domains of life. By contrast, only four unique examples of the hairpin ribozyme class are known, despite its discovery in 1986. We bioinformatically predicted 941 unique hairpin ribozymes of a different permuted form from the four previously known hairpin ribozymes, and experimentally confirmed several diverse predictions. These results profoundly expand the number of natural hairpin ribozymes, enabling biochemical analysis based on natural sequences, and suggest that a distinct permuted form is more biologically relevant. Moreover, all novel hairpins were discovered in metatranscriptomes. They apparently reside in RNA molecules that vary both in size—from 381 to 5170 nucleotides—and in protein content. The RNA molecules likely replicate as circular single-stranded RNAs, and potentially provide a dramatic increase in diversity of such RNAs. Moreover, these organisms have eluded previous attempts to isolate RNA viruses from metatranscriptomes—suggesting a significant untapped universe of viruses or other organisms hidden within metatranscriptome sequences.

siRNA: Mechanism of action, challenges, and therapeutic approaches

Owing to specific and compelling gene silencing, RNA interference (RNAi) is expected to become an essential approach in treating a variety of infectious, hemato-oncological, cardiovascular, and neurodegenerative conditions. The mechanism of action of small interfering RNA (siRNA) is based on post-transcriptional gene silencing. siRNA molecules are usually specific and efficient in the knockdown of disease-related genes. However, they are characterized by low cellular uptake and are susceptible to nuclease-mediated degradation. Therefore, siRNAs require a carrier for their protection and efficient delivery into target cells. The current review highlights the siRNA-based mechanism of action, challanges, and recent advances in clinical applications.

Hovlinc is a recently evolved class of ribozyme found in human lncRNA

Although naturally occurring catalytic RNA molecules-ribozymes-have attracted a great deal of research interest, very few have been identified in humans. Here, we developed a genome-wide approach to discovering self-cleaving ribozymes and identified a naturally occurring ribozyme in humans. The secondary structure and biochemical properties of this ribozyme indicate that it belongs to an unidentified class of small, self-cleaving ribozymes. The sequence of the ribozyme exhibits a clear evolutionary path, from its appearance between ~130 and ~65 million years ago (Ma), to acquiring self-cleavage activity very recently, ~13-10 Ma, in the common ancestors of humans, chimpanzees and gorillas. The ribozyme appears to be functional in vivo and is embedded within a long noncoding RNA belonging to a class of very long intergenic noncoding RNAs. The presence of a catalytic RNA enzyme in lncRNA creates the possibility that these transcripts could function by carrying catalytic RNA domains.

A cell-based ribozyme reporter system employing a chromosomally-integrated 5' exonuclease gene

Background

Bioinformatic genome surveys indicate that self-cleaving ribonucleic acids (ribozymes) appear to be widespread among all domains of life, although the functions of only a small number have been validated by biochemical methods. Alternatively, cell-based reporter gene assays can be used to validate ribozyme function. However, reporter activity can be confounded by phenomena unrelated to ribozyme-mediated cleavage of RNA.

Results

We established a ribozyme reporter system in Escherichia coli in which a significant reduction of reporter activity is manifest when an active ribozyme sequence is fused to the reporter gene and the expression of a foreign Bacillus subtilis RNaseJ1 5′ exonuclease is induced from a chromosomally-integrated gene in the same cell.

Conclusions

The reporter system could be useful for validating ribozyme function in candidate sequences identified from bioinformatics.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12860-021-00357-7.

Experimental Resurrection of Ancestral Mammalian CPEB3 Ribozymes Reveals Deep Functional Conservation

Self-cleaving ribozymes are genetic elements found in all domains of life, but their evolution remains poorly understood. A ribozyme located in the second intron of the cytoplasmic polyadenylation binding protein 3 gene (CPEB3) shows high sequence conservation in mammals, but little is known about the functional conservation of self-cleaving ribozyme activity across the mammalian tree of life or during the course of mammalian evolution. Here, we use a phylogenetic approach to design a mutational library and a deep sequencing assay to evaluate the in vitro self-cleavage activity of numerous extant and resurrected CPEB3 ribozymes that span over 100 My of mammalian evolution. We found that the predicted sequence at the divergence of placentals and marsupials is highly active, and this activity has been conserved in most lineages. A reduction in ribozyme activity appears to have occurred multiple different times throughout the mammalian tree of life. The in vitro activity data allow an evaluation of the predicted mutational pathways leading to extant ribozyme as well as the mutational landscape surrounding these ribozymes. The results demonstrate that in addition to sequence conservation, the self-cleavage activity of the CPEB3 ribozyme has persisted over millions of years of mammalian evolution.

An RNA Triangle with Six Ribozyme Units Can Promote a Trans-Splicing Reaction through Trimerization of Unit Ribozyme Dimers

Ribozymes are catalytic RNAs that are attractive platforms for the construction of nanoscale objects with biological functions. We designed a dimeric form of the Tetrahymena group I ribozyme as a unit structure in which two ribozymes were connected in a tail-to-tail manner with a linker element. We introduced a kink-turn motif as a bent linker element of the ribozyme dimer to design a closed trimer with a triangular shape. The oligomeric states of the resulting ribozyme dimers (kUrds) were analyzed biochemically and observed directly by atomic force microscopy (AFM). Formation of kUrd oligomers also triggered trans-splicing reactions, which could be monitored with a reporter system to yield a fluorescent RNA aptamer as the trans-splicing product.

Virus-associated ribozymes and nano carriers against COVID-19

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is a zoo tonic, highly pathogenic virus. The new type of coronavirus with contagious nature spread from Wuhan (China) to the whole world in a very short time and caused the new coronavirus disease (COVID-19). COVID-19 has turned into a global public health crisis due to spreading by close person-to-person contact with high transmission capacity. Thus, research about the treatment of the damages caused by the virus or prevention from infection increases everyday. Besides, there is still no approved and definitive, standardized treatment for COVID-19. However, this disaster experienced by human beings has made us realize the significance of having a system ready for use to prevent humanity from viral attacks without wasting time. As is known, nanocarriers can be targeted to the desired cells in vitro and in vivo. The nano-carrier system targeting a specific protein, containing the enzyme inhibiting the action of the virus can be developed. The system can be used by simple modifications when we encounter another virus epidemic in the future. In this review, we present a potential treatment method consisting of a nanoparticle-ribozyme conjugate, targeting ACE-2 receptors by reviewing the virus-associated ribozymes, their structures, types and working mechanisms.

Biochemical analysis of cleavage and ligation activities of the pistol ribozyme from Paenibacillus polymyxa

Nine distinct classes of self-cleaving ribozymes are known to date, of which the pistol ribozyme class was discovered only 5 years ago. Self-cleaving ribozymes are able to cleave their own phosphodiester backbone at a specific site with rates much higher than those of spontaneous RNA degradation. Our study focuses on a bioinformatically predicted pistol ribozyme from the bacterium Paenibacillus polymyxa. We provide a biochemical characterization of this ribozyme, which includes an investigation of the effect of various metal ions on ribozyme cleavage and a kinetic analysis of ribozyme activity under increasing Mg²⁺ concentrations and pH. Based on the obtained results, we discuss a possible catalytic role of divalent metal ions. Moreover, we investigated the ligation activity of the P. polymyxa pistol ribozyme - an aspect that has not been previously analysed for this ribozyme class. We determined that the P. polymyxa pistol ribozyme is almost fully cleaved at equilibrium with the ligation rate constant being nearly 30-fold lower than the cleavage rate constant. In summary, we have characterized an additional representative of this recently discovered ribozyme class isolated from P. polymyxa. We expect that our biochemical characterization of a pistol representative in a cultivatable, genetically tractable organism will support our future investigation of the biological roles of this ribozyme class in bacteria.

Oligonucleotide-Based Approaches to Inhibit Dengue Virus Replication

Dengue fever is one of the most common viral infections affecting humans. It is an expanding public health problem, particularly in tropical and subtropical regions. No effective vaccine or antiviral therapies against Dengue virus (DENV) infection are available. Therefore, there is a strong need to develop safe and effective therapeutic strategies that can reduce the burden and duration of hospitalizations due to this life-threatening disease. Oligonucleotide-based strategies are considered as an attractive means of inhibiting viral replication since oligonucleotides can be designed to interact with any viral RNA, provided its sequence is known. The resultant targeted destruction of viral RNA interferes with viral replication without inducing any adverse effects on cellular processes. In this review, we elaborate the ribozymes, RNA interference, CRISPR, aptamer and morpholino strategies for the inhibition of DENV replication and discuss the challenges involved in utilizing such approaches.

Porphyrin metalation catalyzed by DNAzymes and nanozymes

In this review, DNA and nanomaterial based catalysts for porphyrin metalation reactions are summarized, including the selection of DNAzymes, choice of nanomaterials, their catalytic mechanisms, and applications of the reactions.

Recent progress in non-native nucleic acid modifications

While Nature harnesses RNA and DNA to store, read and write genetic information, the inherent programmability, synthetic accessibility and wide functionality of these nucleic acids make them attractive tools for use in a vast array of applications.

Small circRNAs with self-cleaving ribozymes are highly expressed in diverse metazoan transcriptomes

Ribozymes are catalytic RNAs present in modern genomes but regarded as remnants of a prebiotic RNA world. The paradigmatic hammerhead ribozyme (HHR) is a small self-cleaving motif widespread from bacterial to human genomes. Here, we report that most of the classical type I HHRs frequently found in the genomes of animals are contained within a novel family of non-autonomous non-LTR retrotransposons of the retrozyme class. These retroelements are expressed as abundant linear and circular RNAs of ∼170-400 nt in different animal tissues. Bioinformatic and in vitro analyses indicate an efficient self-cleavage of the HHRs harboured in most invertebrate retrozymes, whereas HHRs in retrozymes of vertebrates, such as the axolotl and other amphibians, require to act as dimeric motifs to reach higher self-cleavage rates. Ligation assays of retrozyme RNAs with a protein ligase versus HHR self-ligation indicate that, most likely, tRNA ligases and not the ribozymes are involved in the step of RNA circularization. Altogether, these results confirm the existence of a new and conserved pathway in animals and, likely, eukaryotes in general, for the efficient biosynthesis of RNA circles through small ribozymes, which opens the door for the development of new tools in the emerging field of study of circRNAs.

Imaginary Ribozymes

Biocatalysis is dominated by protein enzymes, and only a few classes of ribozymes are known to contribute to the task of promoting biochemical transformations. The RNA World theory encompasses the notion that earlier forms of life made use of a much greater diversity of ribozymes and other functional RNAs to guide complex metabolic states long before proteins had emerged in evolution. In recent years, the discoveries of various classes of ribozymes, riboswitches, and other noncoding RNAs in bacteria have provided additional support for the hypothesis that RNA molecules indeed have the catalytic competence to promote diverse chemical reactions without the aid of protein enzymes. Herein, some of the most striking observations made from examinations of natural riboswitches that bind small ligands are highlighted and used as a basis to imagine the characteristics and functions of long-extinct ribozymes from the RNA World.

Molecular dynamics analysis of Mg2+ -dependent cleavage of a pistol ribozyme reveals a fail-safe secondary ion for catalysis

Pistol ribozymes comprise a class of small, self-cleaving RNAs discovered via comparative genomic analysis. Prior work in the field has probed the kinetics of the cleavage reaction, as well as the influence of various metal ion cofactors that accelerate the process. In the current study, we performed unbiased and unconstrained molecular dynamics simulations from two current high-resolution pistol crystal structures, and we analyzed trajectory data within the context of the currently accepted ribozyme mechanistic framework. Root-mean-squared deviations, radial distribution functions, and distributions of nucleophilic angle-of-attack reveal insights into the potential roles of three magnesium ions with respect to catalysis and overall conformational stability of the molecule. A series of simulation trajectories containing in silico mutations reveal the relatively flexible and partially interchangeable roles of two particular magnesium ions within solvated hydrogen-bonding distances from the catalytic center.

Light-controlled twister ribozyme with single-molecule detection resolves RNA function in time and space

Small ribozymes such as Oryza sativa twister spontaneously cleave their own RNA when the ribozyme folds into its active conformation. The coupling between twister folding and self-cleavage has been difficult to study, however, because the active ribozyme rapidly converts to product. Here, we describe the synthesis of a photocaged nucleotide that releases guanosine within microseconds upon photosolvolysis with blue light. Application of this tool to O. sativa twister achieved the spatial (75 µm) and temporal (≤30 ms) control required to resolve folding and self-cleavage events when combined with single-molecule fluorescence detection of the ribozyme folding pathway. Real-time observation of single ribozymes after photo-deprotection showed that the precleaved folded state is unstable and quickly unfolds if the RNA does not react. Kinetic analysis showed that Mg²⁺ and Mn²⁺ ions increase ribozyme efficiency by making transitions to the high energy active conformation more probable, rather than by stabilizing the folded ground state or the cleaved product. This tool for light-controlled single RNA folding should offer precise and rapid control of other nucleic acid systems.

Single Mutation in Hammerhead Ribozyme Favors Cleavage Activity with Manganese over Magnesium

Hammerhead ribozymes are one of the most studied classes of ribozymes so far, from both the structural and biochemical point of views. The activity of most hammerhead ribozymes is cation-dependent. Mg2+ is one of the most abundant divalent cations in the cell and therefore plays a major role in cleavage activity for most hammerhead ribozymes. Besides Mg2+, cleavage can also occur in the presence of other cations such as Mn2+. The catalytic core of hammerhead ribozymes is highly conserved, which could contribute to a preference of hammerhead ribozymes toward certain cations. Here, we show a naturally occurring variation in the catalytic core of hammerhead ribozymes, A6C, that can favor one metallic ion, Mn2+, over several other cations.

Catalytic Nucleic Acids: Biochemistry, Chemical Biology, Biosensors, and Nanotechnology

Since the initial discovery of ribozymes in the early 1980s, catalytic nucleic acids have been used in different areas. Compared with protein enzymes, catalytic nucleic acids are programmable in structure, easy to modify, and more stable especially for DNA. We take a historic view to summarize a few main interdisciplinary areas of research on nucleic acid enzymes that may have broader impacts. Early efforts on ribozymes in the 1980s have broken the notion that all enzymes are proteins, supplying new evidence for the RNA world hypothesis. In 1994, the first catalytic DNA (DNAzyme) was reported. Since 2000, the biosensor applications of DNAzymes have emerged and DNAzymes are particularly useful for detecting metal ions, a challenging task for enzymes and antibodies. Combined with nanotechnology, DNAzymes are key building elements for switches allowing dynamic control of materials assembly. The search for new DNAzymes and ribozymes is facilitated by developments in DNA sequencing and computational algorithms, further broadening our fundamental understanding of their biochemistry.

Molecular crowding and RNA catalysis

Molecular crowding promotes RNA folding and catalysis and could have played vital roles in the evolution of primordial ribozymes and protocells.