Structure and function of the hairpin ribozyme

Simulation-Guided Conformational Space Exploration to Assess Reactive Conformations of a Ribozyme

Self-splicing ribozymes are small ribonucleic acid (RNA) enzymes that catalyze their own cleavage through a transphosphoesterification reaction. While this process is involved in some specific steps of viral RNA replication and splicing, it is also of importance in the context of the (putative) first autocatalytic RNA-based systems that could have preceded the emergence of modern life. The uncatalyzed phosphoester bond formation is thermodynamically very unfavorable, and many experimental studies have focused on understanding the molecular features of catalysis in these ribozymes. However, chemical reaction paths are short-lived and not easily characterized by experimental approaches, so molecular simulation approaches appear as an ideal tool to unveil the molecular details of the reaction. Here, we focus on the model hairpin ribozyme. We show that identifying a relevant initial conformation for reactivity studies, which is frequently overlooked in mixed quantum-classical studies that predominantly concentrate on the chemical reaction itself, can be highly challenging. These challenges stem from limitations in both available experimental structures (which are chemically altered to prevent self-cleavage) and the accuracy of force fields, together with the necessity for comprehensive sampling. We show that molecular dynamics simulations, combined with extensive conformational phase space exploration with Hamiltonian replica-exchange simulations, enable us to characterize the relevant conformational basins of the minimal hairpin ribozyme in the ligated state prior to self-cleavage. We find that what is usually considered a canonical reactive conformation with active site geometries and hydrogen-bond patterns that are optimal for the addition-elimination reaction with general acid/general base catalysis is metastable and only marginally populated. The thermodynamically stable conformation appears to be consistent with the expectations of a mechanism that does not require the direct participation of ribozyme residues in the reaction. While these observations may suffer from forcefield inaccuracies, all investigated forcefields lead to the same conclusions upon proper sampling, contrasting with previous investigations on shorter timescales suggesting that at least one reparametrization of the Amber99 forcefield allowed to stabilize aligned active site conformations. Our study demonstrates that identifying the most pertinent reactant state conformation holds equal importance alongside the accurate determination of the thermodynamics and kinetics of the chemical steps of the reaction.

RNA sequence to structure analysis from comprehensive pairwise mutagenesis of multiple self-cleaving ribozymes

Self-cleaving ribozymes are RNA molecules that catalyze the cleavage of their own phosphodiester backbones. These ribozymes are found in all domains of life and are also a tool for biotechnical and synthetic biology applications. Self-cleaving ribozymes are also an important model of sequence-to-function relationships for RNA because their small size simplifies synthesis of genetic variants and self-cleaving activity is an accessible readout of the functional consequence of the mutation. Here, we used a high-throughput experimental approach to determine the relative activity for every possible single and double mutant of five self-cleaving ribozymes. From this data, we comprehensively identified non-additive effects between pairs of mutations (epistasis) for all five ribozymes. We analyzed how changes in activity and trends in epistasis map to the ribozyme structures. The variety of structures studied provided opportunities to observe several examples of common structural elements, and the data was collected under identical experimental conditions to enable direct comparison. Heatmap-based visualization of the data revealed patterns indicating structural features of the ribozymes including paired regions, unpaired loops, non-canonical structures, and tertiary structural contacts. The data also revealed signatures of functionally critical nucleotides involved in catalysis. The results demonstrate that the data sets provide structural information similar to chemical or enzymatic probing experiments, but with additional quantitative functional information. The large-scale data sets can be used for models predicting structure and function and for efforts to engineer self-cleaving ribozymes.

A two-residue nascent-strand steric gate controls synthesis of 2'-O-methyl- and 2'-O-(2-methoxyethyl)-RNA

Steric exclusion is a key element of enzyme substrate specificity, including in polymerases. Such substrate specificity restricts the enzymatic synthesis of 2'-modified nucleic acids, which are of interest in nucleic-acid-based drug development. Here we describe the discovery of a two-residue, nascent-strand, steric control 'gate' in an archaeal DNA polymerase. We show that engineering of the gate to reduce steric bulk in the context of a previously described RNA polymerase activity unlocks the synthesis of 2'-modified RNA oligomers, specifically the efficient synthesis of both defined and random-sequence 2'-O-methyl-RNA (2'OMe-RNA) and 2'-O-(2-methoxyethyl)-RNA (MOE-RNA) oligomers up to 750 nt. This enabled the discovery of RNA endonuclease catalysts entirely composed of 2'OMe-RNA (2'OMezymes) for the allele-specific cleavage of oncogenic KRAS (G12D) and β-catenin CTNNB1 (S33Y) mRNAs, and the elaboration of mixed 2'OMe-/MOE-RNA aptamers with high affinity for vascular endothelial growth factor. Our results open up these 2'-modified RNAs-used in several approved nucleic acid therapeutics-for enzymatic synthesis and a wider exploration in directed evolution and nanotechnology.

A Nucleic Acid Sequence That is Catalytically Active in Both RNA and TNA Backbones

Threose nucleic acid (TNA) is considered a potential RNA progenitor due to its chemical simplicity, base pairing property, and capability of folding into a functional tertiary structure. However, it is unknown whether the functional property can be maintained during transition from TNA to RNA. Here, we use a toggle in vitro selection to identify nucleic acid catalyst sequences that are active in both TNA and RNA backbones. One such nucleic acid enzyme with exchangeable backbone (CAMELEON) catalyzes an RNA cleavage reaction when prepared as TNA (T) and RNA (R). Further biochemical characterization reveals that CAMELEON R and T exhibit different catalytic behaviors such as rate enhancement and magnesium dependence. Structural probing and mutagenesis experiments suggest that they likely fold into distinct tertiary structures. This work demonstrates that the catalytic activity can be preserved during backbone transition from TNA to RNA and provides further experimental support for TNA as an RNA precursor in evolution.

Cotranscriptionally encoded RNA strand displacement circuits

Engineered molecular circuits that process information in biological systems could address emerging human health and biomanufacturing needs. However, such circuits can be difficult to rationally design and scale. DNA-based strand displacement reactions have demonstrated the largest and most computationally powerful molecular circuits to date but are limited in biological systems due to the difficulty in genetically encoding components. Here, we develop scalable cotranscriptionally encoded RNA strand displacement (ctRSD) circuits that are rationally programmed via base pairing interactions. ctRSD circuits address the limitations of DNA-based strand displacement circuits by isothermally producing circuit components via transcription. We demonstrate circuit programmability in vitro by implementing logic and amplification elements, as well as multilayer cascades. Furthermore, we show that circuit kinetics are accurately predicted by a simple model of coupled transcription and strand displacement, enabling model-driven design. We envision ctRSD circuits will enable the rational design of powerful molecular circuits that operate in biological systems, including living cells.

Structural Insights Into Tautomeric Dynamics in Nucleic Acids and in Antiviral Nucleoside Analogs

DNA (2'-deoxyribonucleic acid) and RNA (ribonucleic acid) play diverse functional roles in biology and disease. Despite being comprised primarily of only four cognate nucleobases, nucleic acids can adopt complex three-dimensional structures, and RNA in particular, can catalyze biochemical reactions to regulate a wide variety of biological processes. Such chemical versatility is due in part to the phenomenon of nucleobase tautomerism, whereby the bases can adopt multiple, yet distinct isomeric forms, known as tautomers. For nucleobases, tautomers refer to structural isomers that differ from one another by the position of protons. By altering the position of protons on nucleobases, many of which play critical roles for hydrogen bonding and base pairing interactions, tautomerism has profound effects on the biochemical processes involving nucleic acids. For example, the transient formation of minor tautomers during replication could generate spontaneous mutations. These mutations could arise from the stabilization of mismatches, in the active site of polymerases, in conformations involving minor tautomers that are indistinguishable from canonical base pairs. In this review, we discuss the evidence for tautomerism in DNA, and its consequences to the fidelity of DNA replication. Also reviewed are RNA systems, such as the riboswitches and self-cleaving ribozymes, in which tautomerism plays a functional role in ligand recognition and catalysis, respectively. We also discuss tautomeric nucleoside analogs that are efficacious as antiviral drug candidates such as molnupiravir for coronaviruses and KP1212 for HIV. The antiviral efficacy of these analogs is due, in part, to their ability to exist in multiple tautomeric forms and induce mutations in the replicating viral genomes. From a technical standpoint, minor tautomers of nucleobases are challenging to identify directly because they are rare and interconvert on a fast, millisecond to nanosecond, time scale. Nevertheless, many approaches including biochemical, structural, computational and spectroscopic methods have been developed to study tautomeric dynamics in RNA and DNA systems, and in antiviral nucleoside analogs. An overview of these methods and their applications is included here.

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.

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.

The enzymes for genome size increase and maintenance of large (+)RNA viruses

With sizes <50 kb, viral RNA genomes are at the crossroads of genetic, biophysical, and biochemical stability in their host cell. Here, we analyze the enzymatic assets accompanying large RNA genome viruses, mostly based on recent scientific advances in Coronaviridae. We argue that, in addition to the presence of an RNA exonuclease (ExoN), two markers for the large size of viral RNA genomes are (i) the presence of one or more RNA methyltransferases (MTases) and (ii) a specific architecture of the RNA-dependent RNA polymerase active site. We propose that RNA genome expansion and maintenance are driven by an evolutionary ménage-à-trois made of fast and processive RNA polymerases, RNA repair ExoNs, and RNA MTases that relates to the transition between RNA- to DNA-based life.

Cloning and Detection of Aptamer-Ribozyme Conjugations

RNA aptamers can be used to target proteins or nucleic acids for therapeutic purposes and are candidates for RNA-mediated gene therapy. Like other small therapeutic RNAs, they can be expressed in cells from DNA templates that include a cellular promoter upstream of the RNA coding sequence. Secondary structures flanking aptamers can be used to enhance the activity or stability of these molecules. Notably, flanking self-cleaving ribozymes to remove extraneous nucleotides included during transcription as well as flanking hairpins to improve RNA stability have been used to increase the effect of therapeutic aptamers. Here we describe the cloning procedure of aptamers containing different flanking secondary structures and methods to compare their expression levels by a northern blot protocol optimized for the detection of small RNA molecules.

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.

Hepatitis delta virus-like circular RNAs from diverse metazoans encode conserved hammerhead ribozymes

Human hepatitis delta virus (HDV) is a unique infectious agent whose genome is composed of a small circular RNA. Recent data, however, have reported the existence of highly divergent HDV-like circRNAs in the transcriptomes of diverse vertebrate and invertebrate species. The HDV-like genomes described in amniotes such as birds and reptiles encode self-cleaving RNA motifs or ribozymes similar to the ones present in the human HDV, whereas no catalytic RNA domains have been reported for the HDV-like genomes detected in metagenomic data from some amphibians, fish, and invertebrates. Herein, we describe the self-cleaving motifs of the HDV-like genomes reported in newts and fish, which belong to the characteristic class of HDV ribozymes. Surprisingly, HDV-like genomes from a toad and a termite show conserved type III hammerhead ribozymes, which belong to an unrelated class of catalytic RNAs characteristic of plant genomes and plant subviral circRNAs, such as some viral satellites and viroids. Sequence analyses revealed the presence of similar HDV-like hammerhead ribozymes encoded in two termite genomes, but also in the genomes of several dipteran species. In vitro transcriptions confirmed the cleaving activity for these motifs, with moderate rates of self-cleavage. These data indicate that all described HDV-like agents contain self-cleaving motifs from either the HDV or the hammerhead class. Autocatalytic ribozymes in HDV-like genomes could be regarded as interchangeable domains and may have arisen from cellular transcriptomes, although we still cannot rule out some other evolutionary explanations.

Fundamental studies of functional nucleic acids: aptamers, riboswitches, ribozymes and DNAzymes

This review aims at juxtaposing common versus distinct structural and functional strategies that are applied by aptamers, riboswitches, and ribozymes/DNAzymes. Focusing on recently discovered systems, we begin our analysis with small-molecule binding aptamers, with emphasis on in vitro-selected fluorogenic RNA aptamers and their different modes of ligand binding and fluorescence activation. Fundamental insights are much needed to advance RNA imaging probes for detection of exo- and endogenous RNA and for RNA process tracking. Secondly, we discuss the latest gene expression-regulating mRNA riboswitches that respond to the alarmone ppGpp, to PRPP, to NAD⁺, to adenosine and cytidine diphosphates, and to precursors of thiamine biosynthesis (HMP-PP), and we outline new subclasses of SAM and tetrahydrofolate-binding RNA regulators. Many riboswitches bind protein enzyme cofactors that, in principle, can catalyse a chemical reaction. For RNA, however, only one system (glmS ribozyme) has been identified in Nature thus far that utilizes a small molecule - glucosamine-6-phosphate - to participate directly in reaction catalysis (phosphodiester cleavage). We wonder why that is the case and what is to be done to reveal such likely existing cellular activities that could be more diverse than currently imagined. Thirdly, this brings us to the four latest small nucleolytic ribozymes termed twister, twister-sister, pistol, and hatchet as well as to in vitro selected DNA and RNA enzymes that promote new chemistry, mainly by exploiting their ability for RNA labelling and nucleoside modification recognition. Enormous progress in understanding the strategies of nucleic acids catalysts has been made by providing thorough structural fundaments (e.g. first structure of a DNAzyme, structures of ribozyme transition state mimics) in combination with functional assays and atomic mutagenesis.

Biotechnological and Therapeutic Applications of Natural Nucleic Acid Structural Motifs

Genetic information and the blueprint of life are stored in the form of nucleic acids. The primary sequence of DNA, read from the canonical double helix, provides the code for RNA and protein synthesis. Yet these already-information-rich molecules have higher-order structures which play critical roles in transcription and translation. Uncovering the sequences, parameters, and conditions which govern the formation of these structural motifs has allowed researchers to study them and to utilize them in biotechnological and therapeutic applications in vitro and in vivo. This review covers both DNA and RNA structural motifs found naturally in biological systems including catalytic nucleic acids, non-coding RNA, aptamers, G-quadruplexes, i-motifs, and Holliday junctions. For each category, an overview of the structural characteristics, biological prevalence, and function will be discussed. The biotechnological and therapeutic applications of these structural motifs are highlighted. Future perspectives focus on the addition of proteins and unnatural modifications to enhance structural stability for greater applicability.

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.

Parsimonious Scenario for the Emergence of Viroid-Like Replicons De Novo

Viroids are small, non-coding, circular RNA molecules that infect plants. Different hypotheses for their evolutionary origin have been put forward, such as an early emergence in a precellular RNA World or several de novo independent evolutionary origins in plants. Here, we discuss the plausibility of de novo emergence of viroid-like replicons by giving theoretical support to the likelihood of different steps along a parsimonious evolutionary pathway. While Avsunviroidae-like structures are relatively easy to obtain through evolution of a population of random RNA sequences of fixed length, rod-like structures typical of Pospiviroidae are difficult to fix. Using different quantitative approaches, we evaluated the likelihood that RNA sequences fold into a rod-like structure and bear specific sequence motifs facilitating interactions with other molecules, e.g., RNA polymerases, RNases, and ligases. By means of numerical simulations, we show that circular RNA replicons analogous to Pospiviroidae emerge if evolution is seeded with minimal circular RNAs that grow through the gradual addition of nucleotides. Further, these rod-like replicons often maintain their structure if independent functional modules are acquired that impose selective constraints. The evolutionary scenario we propose here is consistent with the structural and biochemical properties of viroids described to date.

Intracellular Imaging with Genetically Encoded RNA-based Molecular Sensors

Genetically encodable sensors have been widely used in the detection of intracellular molecules ranging from metal ions and metabolites to nucleic acids and proteins. These biosensors are capable of monitoring in real-time the cellular levels, locations, and cell-to-cell variations of the target compounds in living systems. Traditionally, the majority of these sensors have been developed based on fluorescent proteins. As an exciting alternative, genetically encoded RNA-based molecular sensors (GERMS) have emerged over the past few years for the intracellular imaging and detection of various biological targets. In view of their ability for the general detection of a wide range of target analytes, and the modular and simple design principle, GERMS are becoming a popular choice for intracellular analysis. In this review, we summarize different design principles of GERMS based on various RNA recognition modules, transducer modules, and reporting systems. Some recent advances in the application of GERMS for intracellular imaging are also discussed. With further improvement in biostability, sensitivity, and robustness, GERMS can potentially be widely used in cell biology and biotechnology.

RNA Structural Dynamics As Captured by Molecular Simulations: A Comprehensive Overview

With both catalytic and genetic functions, ribonucleic acid (RNA) is perhaps the most pluripotent chemical species in molecular biology, and its functions are intimately linked to its structure and dynamics. Computer simulations, and in particular atomistic molecular dynamics (MD), allow structural dynamics of biomolecular systems to be investigated with unprecedented temporal and spatial resolution. We here provide a comprehensive overview of the fast-developing field of MD simulations of RNA molecules. We begin with an in-depth, evaluatory coverage of the most fundamental methodological challenges that set the basis for the future development of the field, in particular, the current developments and inherent physical limitations of the atomistic force fields and the recent advances in a broad spectrum of enhanced sampling methods. We also survey the closely related field of coarse-grained modeling of RNA systems. After dealing with the methodological aspects, we provide an exhaustive overview of the available RNA simulation literature, ranging from studies of the smallest RNA oligonucleotides to investigations of the entire ribosome. Our review encompasses tetranucleotides, tetraloops, a number of small RNA motifs, A-helix RNA, kissing-loop complexes, the TAR RNA element, the decoding center and other important regions of the ribosome, as well as assorted others systems. Extended sections are devoted to RNA-ion interactions, ribozymes, riboswitches, and protein/RNA complexes. Our overview is written for as broad of an audience as possible, aiming to provide a much-needed interdisciplinary bridge between computation and experiment, together with a perspective on the future of the field.

Structure-based mechanistic insights into catalysis by small self-cleaving ribozymes

Small self-cleaving ribozymes are widely distributed in nature and are essential for rolling-circle-based replication of satellite and pathogenic RNAs. Earlier structure-function studies on the hammerhead, hairpin, glmS, hepatitis delta virus and Varkud satellite ribozymes have provided insights into their overall architecture, their catalytic active site organization, and the role of nearby nucleobases and hydrated divalent cations in facilitating general acid-base and electrostatic-mediated catalysis. This review focuses on recent structure-function research on active site alignments and catalytic mechanisms of the Rzb hammerhead ribozyme, as well as newly-identified pistol, twister and twister-sister ribozymes. In contrast to an agreed upon mechanistic understanding of self-cleavage by Rzb hammerhead and pistol ribozymes, there exists a divergence of views as to the cleavage site alignments and catalytic mechanisms adopted by twister and twister-sister ribozymes. One approach to resolving this conundrum would be to extend the structural studies from currently available pre-catalytic conformations to their transition state mimic vanadate counterparts for both ribozymes.

Structural and Biochemical Properties of Novel Self-Cleaving Ribozymes

Fourteen well-defined ribozyme classes have been identified to date, among which nine are site-specific self-cleaving ribozymes. Very recently, small self-cleaving ribozymes have attracted renewed interest in their structure, biochemistry, and biological function since the discovery, during the last three years, of four novel ribozymes, termed twister, twister sister, pistol, and hatchet. In this review, we mainly address the structure, biochemistry, and catalytic mechanism of the novel ribozymes. They are characterized by distinct active site architectures and divergent, but similar, biochemical properties. The cleavage activities of the ribozymes are highly dependent upon divalent cations, pH, and base-specific mutations, which can cause changes in the nucleotide arrangement and/or electrostatic potential around the cleavage site. It is most likely that a guanine and adenine in close proximity of the cleavage site are involved in general acid-base catalysis. In addition, metal ions appear to play a structural rather than catalytic role although some of their crystal structures have shown a direct metal ion coordination to a non-bridging phosphate oxygen at the cleavage site. Collectively, the structural and biochemical data of the four newest ribozymes could contribute to advance our mechanistic understanding of how self-cleaving ribozymes accomplish their efficient site-specific RNA cleavages.

Unwinding the twister ribozyme: from structure to mechanism

The twister ribozyme motif has been identified by bioinformatic means very recently. Currently, four crystal structures with ordered active sites together with a series of chemical and biochemical data provide insights into how this RNA accomplishes its efficient self‐cleavage. Of particular interest for a mechanistic proposal are structural distinctions observed in the active sites that concern the conformation of the U‐A cleavage site dinucleotide (in‐line alignment of the attacking 2′‐O nucleophile to the to‐be‐cleaved P—O5′ bond versus suboptimal alignments) as well as the presence/absence of Mg²⁺ ions at the scissile phosphate. All structures support the notion that an active site guanine and the conserved adenine at the cleavage site are important contributors to cleavage chemistry, likely being involved in general acid base catalysis. Evidence for innersphere coordination of a Mg²⁺ ion to the pro‐S nonbridging oxygen of the scissile phosphate stems from two of the four crystal structures. Together with the finding of thio/rescue effects for phosphorothioate substrates, this suggests the participation of divalent ions in the overall catalytic strategy employed by twister ribozymes. In this context, it is notable that twister retains wild‐type activity when the phylogenetically conserved stem P1 is deleted, able to cleave a single nucleotide only. WIREs RNA 2017, 8:e1402. doi: 10.1002/wrna.1402

For further resources related to this article, please visit the WIREs website.

Intrinsic Base-Pair Rearrangement in the Hairpin Ribozyme Directs RNA Conformational Sampling and Tertiary Interface Formation

Dynamic fluctuations in RNA structure enable conformational changes that are required for catalysis and recognition. In the hairpin ribozyme, the catalytically active structure is formed as an intricate tertiary interface between two RNA internal loops. Substantial alterations in the structure of each loop are observed upon interface formation, or docking. The very slow on-rate for this relatively tight interaction has led us to hypothesize a double conformational capture mechanism for RNA-RNA recognition. We used extensive molecular dynamics simulations to assess conformational sampling in the undocked form of the loop domain containing the scissile phosphate (loop A). We observed several major accessible conformations with distinctive patterns of hydrogen bonding and base stacking interactions in the active-site internal loop. Several important conformational features characteristic of the docked state were observed in well-populated substates, consistent with the kinetic sampling of docking-competent states by isolated loop A. Our observations suggest a hybrid or multistage binding mechanism, in which initial conformational selection of a docking-competent state is followed by induced-fit adjustment to an in-line, chemically reactive state only after formation of the initial complex with loop B.

Pistol ribozyme adopts a pseudoknot fold facilitating site-specific in-line cleavage

The field of small self-cleaving nucleolytic ribozymes has been invigorated by the recent discovery of the twister, twister-sister, pistol and hatchet ribozymes. We report the crystal structure of a pistol ribozyme termed env25, which adopts a compact tertiary architecture stabilized by an embedded pseudoknot fold. The G-U cleavage site adopts a splayed-apart conformation with in-line alignment of the modeled 2'-O of G for attack on the adjacent to-be-cleaved P-O5' bond. Highly conserved residues G40 (N1 position) and A32 (N3 and 2'-OH positions) are aligned to act as a general base and a general acid, respectively, to accelerate cleavage chemistry, with their roles confirmed by cleavage assays on variants, and an increased pKa of 4.7 for A32. Our structure of the pistol ribozyme defined how the overall and local topologies dictate the in-line alignment at the G-U cleavage site, with cleavage assays on variants revealing key residues that participate in acid-base-catalyzed cleavage chemistry.

Mn(2+)-sensing mechanisms of yybP-ykoY orphan riboswitches

Gene regulation in cis by riboswitches is prevalent in bacteria. The yybP-ykoY riboswitch family is quite widespread, yet its ligand and function remained unknown. Here, we characterize the Lactococcus lactis yybP-ykoY orphan riboswitch as a Mn(2+)-dependent transcription-ON riboswitch, with a ∼30-40 μM affinity for Mn(2+). We further determined its crystal structure at 2.7 Å to elucidate the metal sensing mechanism. The riboswitch resembles a hairpin, with two coaxially stacked helices tethered by a four-way junction and a tertiary docking interface. The Mn(2+)-sensing region, strategically located at the highly conserved docking interface, has two metal binding sites. Whereas one site tolerates the binding of either Mg(2+) or Mn(2+), the other site strongly prefers Mn(2+) due to a direct contact from the N7 of an invariable adenosine. Mutagenesis and a Mn(2+)-free E. coli yybP-ykoY structure further reveal that Mn(2+) binding is coupled with stabilization of the Mn(2+)-sensing region and the aptamer domain.

A Mini-Twister Variant and Impact of Residues/Cations on the Phosphodiester Cleavage of this Ribozyme Class

Nucleolytic ribozymes catalyze site-specific cleavage of their phosphodiester backbones. A minimal version of the twister ribozyme is reported that lacks the phylogenetically conserved stem P1 while retaining wild-type activity. Atomic mutagenesis revealed that nitrogen atoms N1 and N3 of the adenine-6 at the cleavage site are indispensable for cleavage. By NMR spectroscopy, a pKa value of 5.1 was determined for a (13) C2-labeled adenine at this position in the twister ribozyme, which is significantly shifted compared to the pKa of the same adenine in the substrate alone. This finding pinpoints at a potential role for adenine-6 in the catalytic mechanism besides the previously identified invariant guanine-48 and a Mg(2+) ion, both of which are directly coordinated to the non-bridging oxygen atoms of the scissile phosphate; for the latter, additional evidence stems from the observation that Mn(2+) or Cd(2+) accelerated cleavage of phosphorothioate substrates. The relevance of this metal ion binding site is further emphasized by a new 2.6 Å X-ray structure of a 2'-OCH3 -U5 modified twister ribozyme.

Fitness Landscapes of Functional RNAs

The notion of fitness landscapes, a map between genotype and fitness, was proposed more than 80 years ago. For most of this time data was only available for a few alleles, and thus we had only a restricted view of the whole fitness landscape. Recently, advances in genetics and molecular biology allow a more detailed view of them. Here we review experimental and theoretical studies of fitness landscapes of functional RNAs, especially aptamers and ribozymes. We find that RNA structures can be divided into critical structures, connecting structures, neutral structures and forbidden structures. Such characterisation, coupled with theoretical sequence-to-structure predictions, allows us to construct the whole fitness landscape. Fitness landscapes then can be used to study evolution, and in our case the development of the RNA world.

Progress and Prospects of Anti-HBV Gene Therapy Development

Despite the availability of an effective vaccine against hepatitis B virus (HBV), chronic infection with the virus remains a major global health concern. Current drugs against HBV infection are limited by emergence of resistance and rarely achieve complete viral clearance. This has prompted vigorous research on developing better drugs against chronic HBV infection. Advances in understanding the life cycle of HBV and improvements in gene-disabling technologies have been impressive. This has led to development of better HBV infection models and discovery of new drug candidates. Ideally, a regimen against chronic HBV infection should completely eliminate all viral replicative intermediates, especially covalently closed circular DNA (cccDNA). For the past few decades, nucleic acid-based therapy has emerged as an attractive alternative that may result in complete clearance of HBV in infected patients. Several genetic anti-HBV strategies have been developed. The most studied approaches include the use of antisense oligonucleotides, ribozymes, RNA interference effectors and gene editing tools. This review will summarize recent developments and progress made in the use of gene therapy against HBV.

The dynamics of the RNA world: insights and challenges

The RNA world hypothesis of the origin of life, in which RNA emerged as both enzyme and information carrier, is receiving solid experimental support. The prebiotic synthesis of biomolecules, the catalytic aid offered by mineral surfaces, and the vast enzymatic repertoire of ribozymes are only pieces of the origin of life puzzle; the full picture can only emerge if the pieces fit together by either following from one another or coexisting with each other. Here, we review the theory of the origin, maintenance, and enhancement of the RNA world as an evolving population of dynamical systems. The dynamical view of the origin of life allows us to pinpoint the missing and the not fitting pieces: (1) How can the first self-replicating ribozyme emerge in the absence of template-directed information replication? (2) How can nucleotide replicators avoid competitive exclusion despite utilizing the very same resources (nucleobases)? (3) How can the information catastrophe be avoided? (4) How can enough genes integrate into a cohesive system in order to transition to a cellular stage? (5) How can the way information is stored and metabolic complexity coevolve to pave to road leading out of the RNA world to the present protein-DNA world?

(Dis)similar Analogues of Riboswitch Metabolites as Antibacterial Lead Compounds

The rise of antimicrobial resistance in human pathogenic bacteria has increased the necessity for the discovery of novel, yet unexplored antibacterial drug targets. Riboswitches, which are embedded in untranslated regions of bacterial messenger RNA (mRNA), represent such an interesting target structure. These RNA elements regulate gene expression upon binding to natural metabolites, second messengers, and inorganic ions, such as fluoride with high affinity and in a highly discriminative manner. Recently, efforts have been directed toward the identification of artificial riboswitch activators by establishing high-throughput screening assays, fragment-based screening, and structure-guided ligand design approaches. Emphasis in this review is placed on the special requirements and synthesis of new potential antibiotic drugs that target riboswitches in which dissimilarity is an important aspect in the design of potential lead compounds.

Role of tautomerism in RNA biochemistry

Heterocyclic nucleic acid bases and their analogs can adopt multiple tautomeric forms due to the presence of multiple solvent-exchangeable protons. In DNA, spontaneous formation of minor tautomers has been speculated to contribute to mutagenic mispairings during DNA replication, whereas in RNA, minor tautomeric forms have been proposed to enhance the structural and functional diversity of RNA enzymes and aptamers. This review summarizes the role of tautomerism in RNA biochemistry, specifically focusing on the role of tautomerism in catalysis of small self-cleaving ribozymes and recognition of ligand analogs by riboswitches. Considering that the presence of multiple tautomers of nucleic acid bases is a rare occurrence, and that tautomers typically interconvert on a fast time scale, methods for studying rapid tautomerism in the context of nucleic acids under biologically relevant aqueous conditions are also discussed.

In-line alignment and Mg²⁺ coordination at the cleavage site of the env22 twister ribozyme

Small self-cleaving nucleolytic ribozymes contain catalytic domains that accelerate site-specific cleavage/ligation of phosphodiester backbones. We report on the 2.9-Å crystal structure of the env22 twister ribozyme, which adopts a compact tertiary fold stabilized by co-helical stacking, double-pseudoknot formation and long-range pairing interactions. The U-A cleavage site adopts a splayed-apart conformation with the modelled 2'-O of U positioned for in-line attack on the adjacent to-be-cleaved P-O5' bond. Both an invariant guanosine and a Mg(2+) are directly coordinated to the non-bridging phosphate oxygens at the U-A cleavage step, with the former positioned to contribute to catalysis and the latter to structural integrity. The impact of key mutations on cleavage activity identified an invariant guanosine that contributes to catalysis. Our structure of the in-line aligned env22 twister ribozyme is compared with two recently reported twister ribozymes structures, which adopt similar global folds, but differ in conformational features around the cleavage site.

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.

Viroids: survivors from the RNA world?

Because RNA can be a carrier of genetic information and a biocatalyst, there is a consensus that it emerged before DNA and proteins, which eventually assumed these roles and relegated RNA to intermediate functions. If such a scenario--the so-called RNA world--existed, we might hope to find its relics in our present world. The properties of viroids that make them candidates for being survivors of the RNA world include those expected for primitive RNA replicons: (a) small size imposed by error-prone replication, (b) high G + C content to increase replication fidelity, (c) circular structure for assuring complete replication without genomic tags, (d) structural periodicity for modular assembly into enlarged genomes, (e) lack of protein-coding ability consistent with a ribosome-free habitat, and (f) replication mediated in some by ribozymes, the fingerprint of the RNA world. With the advent of DNA and proteins, those protoviroids lost some abilities and became the plant parasites we now know.

Comparison of ab Initio, DFT, and Semiempirical QM/MM Approaches for Description of Catalytic Mechanism of Hairpin Ribozyme

We have analyzed the capability of state-of-the-art multiscale computational approaches to provide atomic-resolution electronic structure insights into possible catalytic scenarios of the hairpin ribozyme by evaluating potential and free energy surfaces of the reactions by various hybrid QM/MM methods. The hairpin ribozyme is a unique catalytic RNA that achieves rate acceleration similar to other small self-cleaving ribozymes but without direct metal ion participation. Guanine 8 (G8) and adenine 38 (A38) have been identified as the catalytically essential nucleobases. However, their exact catalytic roles are still being investigated. In line with the available experimental data, we considered two reaction scenarios involving protonated A38H(+) as a general acid which is further assisted by either canonical G8 or deprotonated G8(-) forms. We used the spin-component scaled Møller-Plesset (SCS-MP2) method at the complete basis set limit as the reference method. The semiempirical AM1/d-PhoT and SCC-DFTBPR methods provided acceptable activation barriers with respect to the SCS-MP2 data but predicted significantly different reaction pathways. DFT functionals (BLYP and MPW1K) yielded the same reaction pathway as the SCS-MP2 method. The activation barriers were slightly underestimated by the GGA BLYP functional, although with accuracy comparable to the semiempirical methods. The SCS-MP2 method and hybrid MPW1K functional gave activation barriers that were closest to those derived from experimentally measured rate constants.

Cobalt(III)hexaammine-dependent photocrosslinks in the hairpin ribozyme

We have utilized the hairpin ribozyme, an RNA enzyme whose structure has been solved by high-resolution methods, to develop a new tool for mapping nucleobase-stacking interactions and potential metal-binding sites in RNA molecules. This tool involves the photoactivation of a specifically bound cobalt(III)hexaammine molecule at wavelengths corresponding to excitation of the metal ion complex only; no base excitation is involved. The photoexcitation initiates a process which strongly promotes the formation of a novel covalent bond or crosslink between one base (termed the "first base"), which is close in space to the excited cobalt(III)hexaammine complex, and another base upon which the first base is closely stacked. These crosslinked species can be isolated and sequenced; their activities can be analyzed to ensure that the crosslinked structures represent an active conformation of the molecule. We have shown that, as in electron transfer in DNA, several criteria must be met to result in the successful formation of these crosslinks. These include the appropriate oxidation potential of the first donor base, the stacking and close interaction of the two donor bases involved in the crosslink, and the binding of a specific cobalt(III)hexaammine molecule to the first donor base. Additionally, we have determined that this crosslinking is pH-sensitive, although the cause of this sensitivity remains unknown. This tool has proven useful in the past for the analysis of the hairpin ribozyme folded structure, and has been applied to identify potential metal-binding sites on the hairpin and extended hammerhead ribozymes.

RNAiFOLD: a constraint programming algorithm for RNA inverse folding and molecular design

Synthetic biology is a rapidly emerging discipline with long-term ramifications that range from single-molecule detection within cells to the creation of synthetic genomes and novel life forms. Truly phenomenal results have been obtained by pioneering groups--for instance, the combinatorial synthesis of genetic networks, genome synthesis using BioBricks, and hybridization chain reaction (HCR), in which stable DNA monomers assemble only upon exposure to a target DNA fragment, biomolecular self-assembly pathways, etc. Such work strongly suggests that nanotechnology and synthetic biology together seem poised to constitute the most transformative development of the 21st century. In this paper, we present a Constraint Programming (CP) approach to solve the RNA inverse folding problem. Given a target RNA secondary structure, we determine an RNA sequence which folds into the target structure; i.e. whose minimum free energy structure is the target structure. Our approach represents a step forward in RNA design--we produce the first complete RNA inverse folding approach which allows for the specification of a wide range of design constraints. We also introduce a Large Neighborhood Search approach which allows us to tackle larger instances at the cost of losing completeness, while retaining the advantages of meeting design constraints (motif, GC-content, etc.). Results demonstrate that our software, RNAiFold, performs as well or better than all state-of-the-art approaches; nevertheless, our approach is unique in terms of completeness, flexibility, and the support of various design constraints. The algorithms presented in this paper are publicly available via the interactive webserver http://bioinformatics.bc.edu/clotelab/RNAiFold; additionally, the source code can be downloaded from that site.

Recent advances in the rational design of silica-based nanoparticles for gene therapy

Gene therapy has attracted much attention in modern society and provides a promising approach for treating genetic disorders, diseases and cancers. Safe and effective vectors are vital tools to deliver genetic molecules to cells. This review summarizes recent advances in the rational design of silica-based nanoparticles and their applications in gene therapy. An overview of different types of genetic agents available for gene therapy is provided. The engineering of various silica nanoparticles is described, which can be used as versatile complexation tools for genetic agents and advanced gene therapy. Several challenges are raised and future research directions in the area of gene therapy using silica-based nanoparticles are proposed.

Click nucleic acid ligation: applications in biology and nanotechnology

Biochemical strategies that use a combination of synthetic oligonucleotides, thermostable DNA polymerases, and DNA ligases can produce large DNA constructs up to 1 megabase in length. Although these ambitious targets are feasible biochemically, comparable technologies for the chemical synthesis of long DNA strands lag far behind. The best available chemical approach is the solid-phase phosphoramidite method, which can be used to assemble DNA strands up to 150 bases in length. Beyond this point, deficiencies in the chemistry make it impossible to produce pure DNA. A possible alternative approach to the chemical synthesis of large DNA strands is to join together carefully purified synthetic oligonucleotides by chemical methods. Click ligation by the copper-catalyzed azide–alkyne (CuAAC) reaction could facilitate this process. In this Account, we describe the synthesis, characterization, and applications of oligonucleotides prepared by click ligation.

The alkyne and azide oligonucleotide strands can be prepared by standard protocols, and the ligation reaction is compatible with a wide range of chemical modifications to DNA and RNA. We have employed click ligation to synthesize DNA constructs up to 300 bases in length and much longer sequences are feasible. When the resulting triazole linkage is placed in a PCR template, various DNA polymerases correctly copy the entire base sequence. We have also successfully demonstrated both in vitro transcription and rolling circle amplification through the modified linkage. This linkage has shown in vivo biocompatibility: an antibiotic resistance gene containing triazole linkages functions in E. coli. Using click ligation, we have synthesized hairpin ribozymes up to 100 nucleotides in length and a hammerhead ribozyme with the triazole linkage located at the substrate cleavage site. At the opposite end of the length scale, click-ligated, cyclic mini-DNA duplexes have been used as models to study base pairing. Cyclic duplexes have potential therapeutic applications. They have extremely high thermodynamic stability, have increased resistance to enzymatic degradation, and have been investigated as decoys for regulatory proteins. For potential nanotechnology applications, we have synthesized double stranded DNA catenanes by click ligation. Other researchers have studied covalently fixed multistranded DNA constructs including triplexes and quadruplexes.

Functional dynamics of RNA ribozymes studied by NMR spectroscopy

Catalytic RNA motifs (ribozymes) are involved in various cellular processes. Although functional cleavage of the RNA phosphodiester backbone for self-cleaving ribozymes strongly differs with respect to sequence specificity, the structural context, and the underlying mechanism, these ribozyme motifs constitute evolved RNA molecules that carry out identical chemical functionality. Therefore, they represent ideal systems for detailed studies of the underlying structure-function relationship, illustrating the diversity of RNA's functional role in biology. Nuclear magnetic resonance (NMR) spectroscopic methods in solution allow investigation of structure and dynamics of functional RNA motifs at atomic resolution. In addition, characterization of RNA conformational transitions initiated either through addition of specific cofactors, as e.g. ions or small molecules, or by photo-chemical triggering of essential RNA functional groups provides insights into the reaction mechanism. Here, we discuss applications of static and time-resolved NMR spectroscopy connected with the design of suitable NMR probes that have been applied to characterize global and local RNA functional dynamics together with cleavage-induced conformational transitions of two RNA ribozyme motifs: a minimal hammerhead ribozyme and an adenine-dependent hairpin ribozyme.

QM/MM studies of hairpin ribozyme self-cleavage suggest the feasibility of multiple competing reaction mechanisms

The hairpin ribozyme is a prominent member of small ribozymes since it does not require metal ions to achieve catalysis. Guanine 8 (G8) and adenine 38 (A38) have been identified as key participants in self-cleavage and -ligation. We have carried out hybrid quantum-mechanical/molecular mechanical (QM/MM) calculations to evaluate the energy along several putative reaction pathways. The error of our DFT description of the QM region was tested and shown to be ~1 kcal/mol. We find that self-cleavage of the hairpin ribozyme may follow several competing microscopic reaction mechanisms, all with calculated activation barriers in good agreement with those from experiment (20-21 kcal/mol). The initial nucleophilic attack of the A-1(2'-OH) group on the scissile phosphate is predicted to be rate-limiting in all these mechanisms. An unprotonated G8(-) (together with A38H(+)) yields a feasible activation barrier (20.4 kcal/mol). Proton transfer to a nonbridging phosphate oxygen also leads to feasible reaction pathways. Finally, our calculations consider thio-substitutions of one or both nonbridging oxygens of the scissile phosphate and predict that they have only a negligible effect on the reaction barrier, as observed experimentally.

Thermodynamics and kinetics of the hairpin ribozyme from atomistic folding/unfolding simulations

We report a set of atomistic folding/unfolding simulations for the hairpin ribozyme using a Monte Carlo algorithm. The hairpin ribozyme folds in solution and catalyzes self-cleavage or ligation via a specific two-domain structure. The minimal active ribozyme has been studied extensively, showing stabilization of the active structure by cations and dynamic motion of the active structure. Here, we introduce a simple model of tertiary-structure formation that leads to a phase diagram for the RNA as a function of temperature and tertiary-structure strength. We then employ this model to capture many folding/unfolding events and to examine the transition-state ensemble (TSE) of the RNA during folding to its active "docked" conformation. The TSE is compact but with few tertiary interactions formed, in agreement with single-molecule dynamics experiments. To compare with experimental kinetic parameters, we introduce a novel method to benchmark Monte Carlo kinetic parameters to docking/undocking rates collected over many single molecular trajectories. We find that topology alone, as encoded in a biased potential that discriminates between secondary and tertiary interactions, is sufficient to predict the thermodynamic behavior and kinetic folding pathway of the hairpin ribozyme. This method should be useful in predicting folding transition states for many natural or man-made RNA tertiary structures.

Ribozymes: Flexible molecular devices at work

The discovery of ribozymes, RNAs with catalytic activity, revealed the extraordinary characteristic of this molecule, and corroborated the idea that RNA was the first informative polymer. The "RNA world" hypothesis asserts that the DNA/RNA/PROTEIN world arose from an earlier RNA world in which were present only RNA molecules able to perform both of the two functions performed separately by DNA and proteins in the present-day cells: the ability to transfer genetic information and to carry out catalytic activity. The catalytic properties of ribozymes are exclusively due to the capacity of RNA molecules to assume particular structures. Moreover, the structural versatility of RNA can allow to a single RNA sequence to fold in more than one structure, able to perform more than one function. In the first part of this work we will discuss the RNA plasticity, focusing on "bifunctional" ribozymes isolated by in vitro selection experiments, and on the consequences of this plasticity in the prospective of the emergence of new specific functions. The possibility that one sequence could have more than one structure/function, greatly increase the evolutionary potential of RNA, and the capacity of RNA to switch from a structure/function to another is probably one of the reasons of the evolutionary success also in modern-day cells. Naturally occurring ribozymes discovered in contemporary cells, demonstrate the crucial role that ribozymes still have in the modern protein world. In the second part of this paper we will discuss the capacity of natural ribozymes to modulate gene expression making use of their exclusive catalytic properties. Moreover, we will consider the possibility of their ancient origin.