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.

Circular at the very beginning: on the initial genomes in the RNA world.

It is likely that an RNA world existed in early life, when RNA played both the roles of the genome and functional molecules, thereby undergoing Darwinian evolution. However, even with only one type of polymer, it seems quite necessary to introduce a labour division concerning these two roles because folding is required for functional molecules (ribozymes) but unfavourable for the genome (as a template in replication). Notably, while ribozymes tend to have adopted a linear form for folding without constraints, a circular form, which might have been topologically hindered in folding, seems more suitable for an RNA template. Another advantage of involving a circular genome could have been to resist RNA's end-degradation. Here, we explore the scenario of a circular RNA genome plus linear ribozyme(s) at the precellular stage of the RNA world through computer modelling. The results suggest that a one-gene scene could have been 'maintained', albeit with rather a low efficiency for the circular genome to produce the ribozyme, which required precise chain-break or chain-synthesis. This strict requirement may have been relieved by introducing a 'noncoding' sequence into the genome, which had the potential to derive a second gene through mutation. A two-gene scene may have 'run well' with the two corresponding ribozymes promoting the replication of the circular genome from different respects. Circular genomes with more genes might have arisen later in RNA-based protocells. Therefore, circular genomes, which are common in the modern living world, may have had their 'root' at the very beginning of life.

Gas-sensing riboceptors.

Understanding how cells sense gases or gaseous solutes is a fundamental question in biology and is pivotal for the evolution of molecular and organismal life. In numerous organisms, gases can diffuse into cells, be transported, generated, and sensed. Controlling gases in the cellular environment is essential to prevent cellular and molecular damage due to interactions with gas-dependent free radicals. Consequently, the mechanisms governing acute gas sensing are evolutionarily conserved and have been experimentally elucidated in various organisms. However, the scientific literature on direct gas sensing is largely based on hemoprotein-based gasoreceptors (or sensors). As RNA-based G-quadruplex (G4) structures can also bind to heme, I propose that some ribozymes can act as gas-sensing riboceptors (ribonucleic acid receptors). Additionally, I present a few other ideas for non-heme metal ion- or metal cluster-based gas-sensing riboceptors. Studying riboceptors can help understand the evolutionary origins of cellular and gasocrine signaling.

Fitness Landscapes and Evolution of Catalytic RNA.

The relationship between genotype and phenotype, or the fitness landscape, is the foundation of genetic engineering and evolution. However, mapping fitness landscapes poses a major technical challenge due to the amount of quantifiable data that is required. Catalytic RNA is a special topic in the study of fitness landscapes due to its relatively small sequence space combined with its importance in synthetic biology. The combination of in vitro selection and high-throughput sequencing has recently provided empirical maps of both complete and local RNA fitness landscapes, but the astronomical size of sequence space limits purely experimental investigations. Next steps are likely to involve data-driven interpolation and extrapolation over sequence space using various machine learning techniques. We discuss recent progress in understanding RNA fitness landscapes, particularly with respect to protocells and machine representations of RNA. The confluence of technical advances may significantly impact synthetic biology in the near future.

Computational Design of Allosteric Ribozymes via Genetic Algorithms.

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

In Vitro Amplification and Selection of Engineered RNase P Ribozyme for Gene Targeting Applications.

Ribozymes engineered from the RNase P catalytic RNA (M1 RNA) represent promising gene-targeting agents for clinical applications. We describe in this report an in vitro amplification and selection procedure for generating active RNase P ribozyme variants with improved catalytic efficiency. Using the amplification and selection procedure, we have previously generated ribozyme variants that were highly active in cleaving a herpes simplex virus 1-encoded mRNA in vitro and inhibiting its expression in virally infected human cells. In this chapter, we use an overlapping region of the mRNAs for the IE1 and IE2 proteins of human cytomegalovirus (HCMV) as a target substrate. We provide detailed protocols and include methods for establishing the procedure for the amplification and selection of active mRNA-cleaving RNase P ribozymes. The in vitro amplification and selection system represents an excellent approach for engineering highly active RNase P ribozymes that can be used in both basic research and clinical applications.

Cell-Free Biosensors Based on Modular Eukaryotic Riboswitches That Function in One Pot at Ambient Temperature.

Artificial riboswitches responsive to user-defined analytes can be constructed by successfully inserting in vitro selected aptamers, which bind to the analytes, into untranslated regions of mRNA. Among them, eukaryotic riboswitches are more promising as biosensors than bacterial ones because they function well at ambient temperature. In addition, cell-free expression systems allow the broader use of these riboswitches as cell-free biosensors in an environmentally friendly manner without cellular limitations. The current best cell-free eukaryotic riboswitch regulates eukaryotic canonical translation initiation through self-cleavage mediated by an implanted analyte-responsive ribozyme (i.e., an aptazyme, an aptamer-ribozyme fusion). However, it has critical flaws as a sensor: due to the less-active ribozyme used, self-cleavage and translation reactions must be conducted separately and sequentially, and a different aptazyme has to be selected to change the analyte specificity, even if an aptamer for the next analyte is available. We here stepwise engineered novel types of cell-free eukaryotic riboswitches that harness highly active self-cleavage and thus require no reaction partitioning. Despite the single-step and one-pot reaction, these riboswitches showed higher analyte dose dependency and sensitivities than the current best cell-free eukaryotic riboswitch requiring multistep reactions. In addition, the analyte specificity can be changed in an extremely facile way, simply by aptamer substitution (and the subsequent simple fine-tuning for giant aptamers). Given that cell-free systems can be lyophilized for storage and transport, the present one-pot and thus easy-to-handle cell-free biosensors utilizing eukaryotic riboswitches are expected to be widely used for on-the-spot sensing of analytes at ambient temperature.

A tribute to Sidney Altman, one of the architects of modern RNA biology.

This special issue of JBC pays tribute to Sidney Altman, whose discovery of a catalytic role for RNA, a breakthrough made independently by Thomas Cech, overturned the long-held dogma that only proteins can serve as catalysts in biological systems. The discovery of RNA catalysis galvanized biologists to think expansively in new directions and has given rise to a remarkable RNAissance in science and medicine. The collection of articles begins with the story of the discovery of RNase P and builds up to the emerging picture of an unexpectedly vast repertoire of RNase P variants in the three domains of life, including insights derived from recent high-resolution structures on how RNAs, ribonucleoproteins, or protein scaffolds can be used variably to generate an active site for catalyzing the same RNA processing reaction. The series of articles ends with a discussion of more recently discovered endonucleases (Argonautes, Cas), whose parallels with RNase P underscore recurring themes in diverse biological contexts.

Primordial aminoacyl-tRNA synthetases preferred minihelices to full-length tRNA.

Aminoacyl-tRNA synthetases (AARS) and tRNAs translate the genetic code in all living cells. Little is known about how their molecular ancestors began to enforce the coding rules for the expression of their own genes. Schimmel et al. proposed in 1993 that AARS catalytic domains began by reading an 'operational' code in the acceptor stems of tRNA minihelices. We show here that the enzymology of an AARS urzyme•TΨC-minihelix cognate pair is a rich in vitro realization of that idea. The TΨC-minihelixLeu is a very poor substrate for full-length Leucyl-tRNA synthetase. It is a superior RNA substrate for the corresponding urzyme, LeuAC. LeuAC active-site mutations shift the choice of both amino acid and RNA substrates. AARS urzyme•minihelix cognate pairs are thus small, pliant models for the ancestral decoding hardware. They are thus an ideal platform for detailed experimental study of the operational RNA code.

The RNA-DNA world and the emergence of DNA-encoded heritable traits.

The RNA world hypothesis confers a central role to RNA molecules in information encoding and catalysis. Even though evidence in support of this hypothesis has accumulated from both experiments and computational modelling, the transition from an RNA world to a world where heritable genetic information is encoded in DNA remains an open question. Recent experiments show that both RNA and DNA templates can extend complementary primers using free RNA/DNA nucleotides, either non-enzymatically or in the presence of a replicase ribozyme. Guided by these experiments, we analyse protocellular evolution with an expanded set of reaction pathways made possible through the presence of DNA nucleotides. By encapsulating these reactions inside three different types of protocellular compartments, each subject to distinct modes of selection, we show how protocells containing DNA-encoded replicases in low copy numbers and replicases in high copy numbers can dominate the population. This is facilitated by a reaction that leads to auto-catalytic synthesis of replicase ribozymes from DNA templates encoding the replicase after the chance emergence of a replicase through non-enzymatic reactions. Our work unveils a pathway for the transition from an RNA world to a mixed RNA-DNA world characterized by Darwinian evolution, where DNA sequences encode heritable phenotypes.

Emergent ribozyme behaviors in oxychlorine brines indicate a unique niche for molecular evolution on Mars.

Mars is a particularly attractive candidate among known astronomical objects to potentially host life. Results from space exploration missions have provided insights into Martian geochemistry that indicate oxychlorine species, particularly perchlorate, are ubiquitous features of the Martian geochemical landscape. Perchlorate presents potential obstacles for known forms of life due to its toxicity. However, it can also provide potential benefits, such as producing brines by deliquescence, like those thought to exist on present-day Mars. Here we show perchlorate brines support folding and catalysis of functional RNAs, while inactivating representative protein enzymes. Additionally, we show perchlorate and other oxychlorine species enable ribozyme functions, including homeostasis-like regulatory behavior and ribozyme-catalyzed chlorination of organic molecules. We suggest nucleic acids are uniquely well-suited to hypersaline Martian environments. Furthermore, Martian near- or subsurface oxychlorine brines, and brines found in potential lifeforms, could provide a unique niche for biomolecular evolution.

The discovery of a catalytic RNA within RNase P and its legacy.

Sidney Altman's discovery of the processing of one RNA by another RNA that acts like an enzyme was revolutionary in biology and the basis for his sharing the 1989 Nobel Prize in Chemistry with Thomas Cech. These breakthrough findings support the key role of RNA in molecular evolution, where replicating RNAs (and similar chemical derivatives) either with or without peptides functioned in protocells during the early stages of life on Earth, an era referred to as the RNA world. Here, we cover the historical background highlighting the work of Altman and his colleagues and the subsequent efforts of other researchers to understand the biological function of RNase P and its catalytic RNA subunit and to employ it as a tool to downregulate gene expression. We primarily discuss bacterial RNase P-related studies but acknowledge that many groups have significantly contributed to our understanding of archaeal and eukaryotic RNase P, as reviewed in this special issue and elsewhere.

From Enigma to Revelation: Unravelling Biological Functions of Ubiquitous Small Ribozymes.

RNA, widely recognized as an information-carrier molecule, is capable of catalyzing essential biological processes through ribozymes. Despite their ubiquity, specific functions in a biological context and phenotypes based on the ribozymes' activity are often unknown. Here, we present the discovery of a subgroup of minimal HDV-like ribozymes, which reside 3' to viral tRNAs and appear to cleave the 3'-trailers of viral premature tRNA transcripts. This proposed tRNA-processing function is unprecedented for any ribozymes, thus, we designate this subgroup as theta ribozymes. Most theta ribozymes were identified in Caudoviricetes bacteriophages, the main constituent (>90%) of the mammalian gut virome. Intriguingly, our findings further suggest the involvement of theta ribozymes in the transition of certain bacteriophages between distinct genetic codes, thus possibly contributing to the phage lysis trigger. Our discovery expands the limited repertoire of biological functions attributed to HDV-like ribozymes and provides insights into the fascinating world of RNA catalysis.

RNA targeting and cleavage by the type III-Dv CRISPR effector complex.

CRISPR-Cas are adaptive immune systems in bacteria and archaea that utilize CRISPR RNA-guided surveillance complexes to target complementary RNA or DNA for destruction1-5. Target RNA cleavage at regular intervals is characteristic of type III effector complexes6-8. Here, we determine the structures of the Synechocystis type III-Dv complex, an apparent evolutionary intermediate from multi-protein to single-protein type III effectors9,10, in pre- and post-cleavage states. The structures show how multi-subunit fusion proteins in the effector are tethered together in an unusual arrangement to assemble into an active and programmable RNA endonuclease and how the effector utilizes a distinct mechanism for target RNA seeding from other type III effectors. Using structural, biochemical, and quantum/classical molecular dynamics simulation, we study the structure and dynamics of the three catalytic sites, where a 2'-OH of the ribose on the target RNA acts as a nucleophile for in line self-cleavage of the upstream scissile phosphate. Strikingly, the arrangement at the catalytic residues of most type III complexes resembles the active site of ribozymes, including the hammerhead, pistol, and Varkud satellite ribozymes. Our work provides detailed molecular insight into the mechanisms of RNA targeting and cleavage by an important intermediate in the evolution of type III effector complexes.

Catalytic Metal Ion-Substrate Coordination during Nonenzymatic RNA Primer Extension.

Nonenzymatic template-directed RNA copying requires catalysis by divalent metal ions. The primer extension reaction involves the attack of the primer 3'-hydroxyl on the adjacent phosphate of a 5'-5'-imidazolium-bridged dinucleotide substrate. However, the nature of the interaction of the catalytic metal ion with the reaction center remains unclear. To explore the coordination of the catalytic metal ion with the imidazolium-bridged dinucleotide substrate, we examined catalysis by oxophilic and thiophilic metal ions with both diastereomers of phosphorothioate-modified substrates. We show that Mg2+ and Cd2+ exhibit opposite preferences for the two phosphorothioate substrate diastereomers, indicating a stereospecific interaction of the divalent cation with one of the nonbridging phosphorus substituents. High-resolution X-ray crystal structures of the products of primer extension with phosphorothioate substrates reveal the absolute stereochemistry of this interaction and indicate that catalysis by Mg2+ involves inner-sphere coordination with the nonbridging phosphate oxygen in the pro-SP position, while thiophilic cadmium ions interact with sulfur in the same position, as in one of the two phosphorothioate substrates. These results collectively suggest that during nonenzymatic RNA primer extension with a 5'-5'-imidazolium-bridged dinucleotide substrate the interaction of the catalytic Mg2+ ion with the pro-SP oxygen of the reactive phosphate plays a crucial role in the metal-catalyzed SN2(P) reaction.

Chemoenzymatic Installation of Site-Specific Chemical Groups on DNA Enhances the Catalytic Activity.

Functional DNAs are valuable molecular tools in chemical biology and analytical chemistry but suffer from low activities due to their limited chemical functionalities. Here, we present a chemoenzymatic method for site-specific installation of diverse functional groups on DNA, and showcase the application of this method to enhance the catalytic activity of a DNA catalyst. Through chemoenzymatic introduction of distinct chemical groups, such as hydroxyl, carboxyl, and benzyl, at specific positions, we achieve significant enhancements in the catalytic activity of the RNA-cleaving deoxyribozyme 10-23. A single carboxyl modification results in a 100-fold increase, while dual modifications (carboxyl and benzyl) yield an approximately 700-fold increase in activity when an RNA cleavage reaction is catalyzed on a DNA-RNA chimeric substrate. The resulting dually modified DNA catalyst, CaBn, exhibits a kobs of 3.76 min-1 in the presence of 1 mM Mg2+ and can be employed for fluorescent imaging of intracellular magnesium ions. Molecular dynamics simulations reveal the superior capability of CaBn to recruit magnesium ions to metal-ion-binding site 2 and adopt a catalytically competent conformation. Our work provides a broadly accessible strategy for DNA functionalization with diverse chemical modifications, and CaBn offers a highly active DNA catalyst with immense potential in chemistry and biotechnology.

RNA-catalyzed evolution of catalytic RNA.

An RNA polymerase ribozyme that was obtained by directed evolution can propagate a functional RNA through repeated rounds of replication and selection, thereby enabling Darwinian evolution. Earlier versions of the polymerase did not have sufficient copying fidelity to propagate functional information, but a new variant with improved fidelity can replicate the hammerhead ribozyme through reciprocal synthesis of both the hammerhead and its complement, with the products then being selected for RNA-cleavage activity. Two evolutionary lineages were carried out in parallel, using either the prior low-fidelity or the newer high-fidelity polymerase. The former lineage quickly lost hammerhead functionality as the population diverged toward random sequences, whereas the latter evolved new hammerhead variants with improved fitness compared to the starting RNA. The increase in fitness was attributable to specific mutations that improved the replicability of the hammerhead, counterbalanced by a small decrease in hammerhead activity. Deep sequencing analysis was used to follow the course of evolution, revealing the emergence of a succession of variants that progressively diverged from the starting hammerhead as fitness increased. This study demonstrates the critical importance of replication fidelity for maintaining heritable information in an RNA-based evolving system, such as is thought to have existed during the early history of life on Earth. Attempts to recreate RNA-based life in the laboratory must achieve further improvements in replication fidelity to enable the fully autonomous Darwinian evolution of RNA enzymes as complex as the polymerase itself.

High-Fidelity RNA Copying via 2',3'-Cyclic Phosphate Ligation.

Templated ligation offers an efficient approach to replicate long strands in an RNA world. The 2',3'-cyclic phosphate (>P) is a prebiotically available activation that also forms during RNA hydrolysis. Using gel electrophoresis and high-performance liquid chromatography, we found that the templated ligation of RNA with >P proceeds in simple low-salt aqueous solutions with 1 mM MgCl2 under alkaline pH ranging from 9 to 11 and temperatures from -20 to 25 °C. No additional catalysts were required. In contrast to previous reports, we found an increase in the number of canonical linkages to 50%. The reaction proceeds in a sequence-specific manner, with an experimentally determined ligation fidelity of 82% at the 3' end and 91% at the 5' end of the ligation site. With splinted oligomers, five ligations created a 96-mer strand, demonstrating a pathway for the ribozyme assembly. Due to the low salt requirements, the ligation conditions will be compatible with strand separation. Templated ligation mediated by 2',3'-cyclic phosphate in alkaline conditions therefore offers a performant replication and elongation reaction for RNA on early Earth.

Architecture of an RNA polymerase ribozyme illuminates the RNA World.

'Ribo-organisms' of the primordial RNA World would have needed ribozymes that catalyze RNA replication. McRae, Wan, Kristoffersen et al. recently revealed how these RNA replicases might have functioned by solving the structure of an artificial polymerase ribozyme. This work illustrates how complex RNA structures evolve, with implications for the origins of life.