General acid-base catalysis in the mechanism of a hepatitis delta virus ribozyme

RNA catalysis moving towards metabolic reactions: progress with ribozyme catalyzed alkyl transfer

The RNA world hypothesis proposes that the early stages of the emergence of life on Earth comprised primitive cells in which RNA acted both to store genetic information and catalyze chemical reactions as RNA enzymes (ribozymes). Most contemporary ribozymes catalyze phosphoryl transfer reactions, but early ribozymes would have been required to catalyze a broader range of metabolic interconversions. None has been found in modern cells, yet ribozymes have been generated by in vitro evolution to catalyze several different chemical reactions, providing proof of principle of RNA-catalyzed metabolism. Recently, several different ribozymes that accelerate methyl or alkyl transfer have been isolated. As we discuss here, one of these, MTR1, uses a remarkably sophisticated catalytic mechanism involving nucleobase-mediated general acid catalysis.

Direct testing of natural twister ribozymes from over a thousand organisms reveals a broad tolerance for structural imperfections

Twister ribozymes are an extensively studied class of nucleolytic RNAs. Thousands of natural twisters have been proposed using sequence homology and structural descriptors. Yet, most of these candidates have not been validated experimentally. To address this gap, we developed Cleavage High-Throughput Assay (CHiTA), a high-throughput pipeline utilizing massively parallel oligonucleotide synthesis and next-generation sequencing to test putative ribozymes en masse in a scarless fashion. As proof of principle, we applied CHiTA to a small set of known active and mutant ribozymes. We then used CHiTA to test two large sets of naturally occurring twister ribozymes: over 1600 previously reported putative twisters and ∼1000 new candidate twisters. The new candidates were identified computationally in ∼1000 organisms, representing a massive increase in the number of ribozyme-harboring organisms. Approximately 94% of the twisters we tested were active and cleaved site-specifically. Analysis of their structural features revealed that many substitutions and helical imperfections can be tolerated. We repeated our computational search with structural descriptors updated from this analysis, whereupon we identified and confirmed the first intrinsically active twister ribozyme in mammals. CHiTA broadly expands the number of active twister ribozymes found in nature and provides a powerful method for functional analyses of other RNAs.

The Role of General Acid Catalysis in the Mechanism of an Alkyl Transferase Ribozyme

MTR1 is an in vitro-selected alkyl transferase ribozyme that transfers an alkyl group from O 6-alkylguanine to N1 of the target adenine in the RNA substrate (A63). The structure of the ribozyme suggested a mechanism in which a cytosine (C10) acts as a general acid to protonate O 6-alkylguanine N1. Here, we have analyzed the role of the C10 general acid and the A63 nucleophile by atomic mutagenesis and computation. C10 was substituted by n1c and n1c, c5n variants. The n1c variant has an elevated pK a (11.4 as the free nucleotide) and leads to a 104-fold lower activity that is pH-independent. Addition of the second c5n substitution with a lower pK a restored both the rate and pH dependence of alkyl transfer. Quantum mechanical calculations indicate that protonation of O 6-alkylguanine lowers the barrier to alkyl transfer and that there is a significantly elevated barrier to proton transfer for the n1c single substitution. The calculated pK a values are in good agreement with the apparent values from measured rates. Increasing the pK a of the nucleophile by A63 n7c substitution led to a 6-fold higher rate. The increased reactivity of the nucleophile corresponds to a βnuc of ∼0.5, indicating significant C-N bond formation in the transition state. Taken together, these results are consistent with a two-step mechanism comprising protonation of the O 6-alkylguanine followed by alkyl transfer.

Exploring the catalytic mechanism of the 10-23 DNAzyme: insights from pH-rate profiles

The 10-23 DNAzyme, a catalytic DNA molecule with RNA-cleaving activity, has garnered significant interest for its potential therapeutic applications as a gene-silencing agent. However, the lack of a detailed understanding about its mechanism has hampered progress. A recent structural analysis has revealed a highly organized conformation thanks to the stabilization of specific interactions within the catalytic core of the 10-23 DNAzyme, which facilitate the cleavage of RNA. In this configuration, it has been shown that G14 is in good proximity to the cleavage site which suggests its role as a general base, by activating the 2'-OH nucleophile, in the catalysis of the 10-23 DNAzyme. Also, the possibility of a hydrated metal acting as a general acid has been proposed. In this study, through activity assays, we offer evidence of the involvement of general acid-base catalysis in the mechanism of the 10-23 DNAzyme by analyzing its pH-rate profiles and the role of G14, and metal cofactors like Mg2+ and Pb2+. By substituting G14 with its analogue 2-aminopurine and examining the resultant pH-rate profiles, we propose the participation of G14 in a catalytically relevant proton transfer event, acting as a general base. Further analysis, using Pb2+ as a cofactor, suggests the capability of the hydrated metal ion to act as a general acid. These functional results provide critical insights into the catalytic strategies of RNA-cleaving DNAzymes, revealing common mechanisms among nucleic acid enzymes that cleave RNA.

Structure and sequence at an RNA template 5' end influence insertion of transgenes by an R2 retrotransposon protein

R2 non-long terminal repeat retrotransposons insert site-specifically into ribosomal RNA genes (rDNA) in a broad range of multicellular eukaryotes. R2-encoded proteins can be leveraged to mediate transgene insertion at 28S rDNA loci in cultured human cells. This strategy, precise RNA-mediated insertion of transgenes (PRINT), relies on the codelivery of an mRNA encoding R2 protein and an RNA template encoding a transgene cassette of choice. Here, we demonstrate that the PRINT RNA template 5' module, which as a complementary DNA 3' end will generate the transgene 5' junction with rDNA, influences the efficiency and mechanism of gene insertion. Iterative design and testing identified optimal 5' modules consisting of a hepatitis delta virus-like ribozyme fold with high thermodynamic stability, suggesting that RNA template degradation from its 5' end may limit transgene insertion efficiency. We also demonstrate that transgene 5' junction formation can be either precise, formed by annealing the 3' end of first-strand complementary DNA with the upstream target site, or imprecise, by end-joining, but this difference in junction formation mechanism is not a major determinant of insertion efficiency. Sequence characterization of imprecise end-joining events indicates surprisingly minimal reliance on microhomology. Our findings expand the current understanding of the role of R2 retrotransposon transcript sequence and structure, and especially the 5' ribozyme fold, for retrotransposon mobility and RNA-templated gene synthesis in cells.

Direct testing of natural twister ribozymes from over a thousand organisms reveals a broad tolerance for structural imperfections

Twister ribozymes are an extensively studied class of nucleolytic RNAs. Thousands of natural twisters have been proposed using sequence homology and structural descriptors. Yet, most of these candidates have not been validated experimentally. To address this gap, we developed CHiTA (Cleavage High-Throughput Assay), a high-throughput pipeline utilizing massively parallel oligonucleotide synthesis and next-generation sequencing to test putative ribozymes en masse in a scarless fashion. As proof of principle, we applied CHiTA to a small set of known active and mutant ribozymes. We then used CHiTA to test two large sets of naturally occurring twister ribozymes: over 1, 600 previously reported putative twisters and ∼1, 000 new candidate twisters. The new candidates were identified computationally in ∼1, 000 organisms, representing a massive increase in the number of ribozyme-harboring organisms. Approximately 94% of the twisters we tested were active and cleaved site-specifically. Analysis of their structural features revealed that many substitutions and helical imperfections can be tolerated. We repeated our computational search with structural descriptors updated from this analysis, whereupon we identified and confirmed the first intrinsically active twister ribozyme in mammals. CHiTA broadly expands the number of active twister ribozymes found in nature and provides a powerful method for functional analyses of other RNAs.

GRAPHICAL ABSTRACT

Assessment of Nucleobase Protomeric and Tautomeric States in Nucleic Acid Structures for Interaction Analysis and Structure-Based Ligand Design

With increasing interest in RNA as a therapeutic and a potential target, the role of RNA structures has become more important. Even slight changes in nucleobases, such as modifications or protomeric and tautomeric states, can have a large impact on RNA structure and function, while local environments in turn affect protonation and tautomerization. In this work, the application of empirical tools for pKa and tautomer prediction for RNA modifications was elucidated and compared with ab initio quantum mechanics (QM) methods and expanded toward macromolecular RNA structures, where QM is no longer feasible. In this regard, the Protonate3D functionality within the molecular operating environment (MOE) was expanded for nucleobase protomer and tautomer predictions and applied to reported examples of altered protonation states depending on the local environment. Overall, observations of nonstandard protomers and tautomers were well reproduced, including structural C+G:C(A) and A+GG motifs, several mismatches, and protonation of adenosine or cytidine as the general acid in nucleolytic ribozymes. Special cases, such as cobalt hexamine-soaked complexes or the deprotonation of guanosine as the general base in nucleolytic ribozymes, proved to be challenging. The collected set of examples shall serve as a starting point for the development of further RNA protonation prediction tools, while the presented Protonate3D implementation already delivers reasonable protonation predictions for RNA and DNA macromolecules. For cases where higher accuracy is needed, like following catalytic pathways of ribozymes, incorporation of QM-based methods can build upon the Protonate3D-generated starting structures. Likewise, this protonation prediction can be used for structure-based RNA-ligand design approaches.

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.

Characterizing Protonation-Coupled Conformational Ensembles in RNA via pH-Differential Mutational Profiling with DMS Probing

RNA molecules undergo conformational transitions in response to cellular and environmental stimuli. Site-specific protonation, a fundamental chemical property, can alter the conformational landscape of RNA to regulate their functions. However, characterizing protonation-coupled RNA conformational ensembles on a large scale remains challenging. Here, we present pH-differential mutational profiling (PD-MaP) with dimethyl sulfate probing for high-throughput detection of protonation-coupled conformational ensembles in RNA. We demonstrated this approach on microRNA-21 precursor (pre-miR-21) and recapitulated a previously discovered A+-G-coupled conformational ensemble. Additionally, we identified a secondary protonation event involving an A+-C mismatch. We validated the occurrence of both protonation-coupled ensembles in pre-miR-21 using NMR relaxation dispersion spectroscopy. Furthermore, the application of PD-MaP on a library of well-annotated human primary microRNAs uncovered widespread protonation-coupled conformational ensembles, suggesting their potentially broad functions in biology.

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.

Prebiotic Foam Environments to Oligomerize and Accumulate RNA

When water interacts with porous rocks, its wetting and surface tension properties create air bubbles in large number. To probe their relevance as a setting for the emergence of life, we microfluidically created foams that were stabilized with lipids. A persistent non-equilibrium setting was provided by a thermal gradient. The foam's large surface area triggers capillary flows and wet-dry reactions that accumulate, aggregate and oligomerize RNA, offering a compelling habitat for RNA-based early life as it offers both wet and dry conditions in direct neighborhood. Lipids were screened to stabilize the foams. The prebiotically more probable myristic acid stabilized foams over many hours. The capillary flow created by the evaporation at the water-air interface provided an attractive force for molecule localization and selection for molecule size. For example, self-binding oligonucleotide sequences accumulated and formed micrometer-sized aggregates which were shuttled between gas bubbles. The wet-dry cycles at the foam bubble interfaces triggered a non-enzymatic RNA oligomerization from 2',3'-cyclic CMP and GMP which despite the small dry reaction volume was superior to the corresponding dry reaction. The found characteristics make heated foams an interesting, localized setting for early molecular evolution.

In vitro evolution of ribonucleases from expanded genetic alphabets

The ability of nucleic acids to catalyze reactions (as well as store and transmit information) is important for both basic and applied science, the first in the context of molecular evolution and the origin of life and the second for biomedical applications. However, the catalytic power of standard nucleic acids (NAs) assembled from just four nucleotide building blocks is limited when compared with that of proteins. Here, we assess the evolutionary potential of libraries of nucleic acids with six nucleotide building blocks as reservoirs for catalysis. We compare the outcomes of in vitro selection experiments toward RNA-cleavage activity of two nucleic acid libraries: one built from the standard four independently replicable nucleotides and the other from six, with the two added nucleotides coming from an artificially expanded genetic information system (AEGIS). Results from comparative experiments suggest that DNA libraries with increased chemical diversity, higher information density, and larger searchable sequence spaces are one order of magnitude richer reservoirs of molecules that catalyze the cleavage of a phosphodiester bond in RNA than DNA libraries built from a standard four-nucleotide alphabet. Evolved AEGISzymes with nitro-carrying nucleobase Z appear to exploit a general acid–base catalytic mechanism to cleave that bond, analogous to the mechanism of the ribonuclease A family of protein enzymes and heavily modified DNAzymes. The AEGISzyme described here represents a new type of catalysts evolved from libraries built from expanded genetic alphabets.

A Cofactor-Based Mechanism for the Origin of the Genetic Code

The origin of the genetic code is probably the central problem of the studies on the origin of life. The key question to answer is the molecular mechanism that allows the association of the amino acids with their triplet codons. We proposed that the codon-anticodon duplex located in the acceptor stem of primitive tRNAs would facilitate the chemical reactions required to synthesize cognate amino acids from simple amino acids (glycine, valine, and aspartic acid) linked to the 3' acceptor end. In our view, various nucleotide-A-derived cofactors (with reactive chemical groups) may be attached to the codon-anticodon duplex, which allows group-transferring reactions from cofactors to simple amino acids, thereby producing the final amino acid. The nucleotide-A-derived cofactors could be incorporated into the RNA duplex (helix) by docking Adenosine (cofactor) into the minor groove via an interaction similar to the A-minor motif, forming a base triple between Adenosine and one complementary base pair of the duplex. Furthermore, we propose that this codon-anticodon duplex could initially catalyze a self-aminoacylation reaction with a simple amino acid. Therefore, the sequence of bases in the codon-anticodon duplex would determine the reactions that occurred during the formation of new amino acids for selective binding of nucleotide-A-derived cofactors.

The role of Na+ in catalysis by the 8-17 DNAzyme

The 8-17 DNAzyme is the most studied deoxyribozyme in terms of its molecular mechanism; hence it has become a model system to understand the basis behind DNA catalysis. New functional studies and the recent attainment of high-resolution X-ray structures, in addition to theoretical calculations have offered a great opportunity to gain a broader comprehension of its mechanism; however many aspects are unclear yet, especially regarding the precise role of metal ions in catalysis. Recently, molecular dynamics simulations have suggested for the first time a specific and dynamical participation of Na⁺ in the mechanism through the reaction pathway, besides the roles proposed for divalent metal cofactors. Herein, we present experimental evidence of a cooperative role of the monovalent cation Na⁺ in catalysis that is in line with these theoretical suggestions. Our findings show a clear influence of the concentration of Na⁺ on the activity of the 8-17 DNAzyme when Pb²⁺ is used as the cofactor. Interestingly, this effect is not noticed with Mg²⁺, indicating a particular contribution of the monovalent ion to catalysis that would operate preferentially with Pb²⁺. We have also found that Na⁺ affects the pK_a of the general base and the general acid, indicating its influence on general acid-base catalysis, already identified as part of the mechanism of the 8-17 DNAzyme. Finally, our results emphasize the need to consider Na⁺ carefully in the design and analysis of functional studies of catalytic DNAs and its possible specific role in their mechanisms.

Calculations of pKa Values for a Series of Naturally Occurring Modified Nucleobases

Modified nucleobases are found in functionally important regions of RNA and are often responsible for essential structural roles. Many of these nucleobase modifications are dynamically regulated in nature, with each modification having a different biological role in RNA. Despite the high abundance of modifications, many of their characteristics are still poorly understood. One important property of a nucleobase is its pK_a value, which has been widely studied for unmodified nucleobases, but not for the modified versions. In this study, the pK_a values of modified nucleobases were determined by performing ab initio quantum mechanical calculations using a B3LYP density functional with the 6-31+G(d,p) basis set and a combination of implicit-explicit solvation systems. This method, which was previously employed to determine the pK_a values of unmodified nucleobases, is applicable to a variety of modified nucleobases. Comparisons of the pK_a values of modified nucleobases give insight into their structural and energetic impacts within nucleic acids.

The development and improvement of ribonucleic acid therapy strategies

The biological understanding of RNA has evolved since the discovery of catalytic RNAs in the early 1980s and the establishment of RNA interference (RNAi) in the 1990s. RNA is no longer seen as the simple mid-product between transcription and translation but as potential molecules to be developed as RNA therapeutic drugs. RNA-based therapeutic drugs have gained recognition because of their ability to regulate gene expression and perform cellular functions. Various nucleobase, backbone, and sugar-modified oligonucleotides have been synthesized, as natural oligonucleotides have some limitations such as poor low nuclease resistance, binding affinity, poor cellular uptake, and toxicity, which affect their use as RNA therapeutic drugs. In this review, we briefly discuss different RNA therapeutic drugs and their internal connections, including antisense oligonucleotides, small interfering RNAs (siRNAs) and microRNAs (miRNAs), aptamers, small activating RNAs (saRNAs), and RNA vaccines. We also discuss the important roles of RNA vaccines and their use in the fight against COVID-19. In addition, various chemical modifications and delivery systems used to improve the performance of RNA therapeutic drugs and overcome their limitations are discussed.

Investigation of the pKa of the Nucleophilic O2' of the Hairpin Ribozyme

Small ribozymes cleave their RNA phosphodiester backbone by catalyzing a transphosphorylation reaction wherein a specific O2' functions as the nucleophile. While deprotonation of this alcohol through its acidification would increase its nucleophilicity, little is known about the pK_a of this O2' in small ribozymes, in part because high pK_a's are not readily accessible experimentally. Herein, we turn to molecular dynamics to calculate the pK_a of the nucleophilic O2' in the hairpin ribozyme and to study interactions within the active site that may impact its value. We estimate the pK_a of the nucleophilic O2' in the wild-type hairpin ribozyme to be 18.5 ± 0.8, which is higher than the reference compound, and identify a correlation between proper positioning of the O2' for nucleophilic attack and elevation of its pK_a. We find that monovalent ions may play a role in depression of the O2' pK_a, while the exocyclic amine appears to be important for organizing the ribozyme active site. Overall, this study suggests that the pK_a of the O2' is raised in the ground state and lowers during the course of the reaction owing to positioning and metal ion interactions.

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.

Beyond the Plateau: pL Dependence of Proton Inventories as a Tool for Studying Ribozyme and Ribonuclease Catalysis

Acid/base catalysis is an important catalytic strategy used by ribonucleases and ribozymes; however, understanding the number and identity of functional groups involved in proton transfer remains challenging. The proton inventory (PI) technique analyzes the dependence of the enzyme reaction rate on the ratio of D₂O to H₂O and can provide information about the number of exchangeable sites that produce isotope effects and their magnitude. The Gross-Butler (GB) equation is used to evaluate H/D fractionation factors from PI data typically collected under conditions (i.e., a "plateau" in the pH-rate profile) assuming minimal change in active site residue ionization. However, restricting PI analysis to these conditions is problematic for many ribonucleases, ribozymes, and their variants due to ambiguity in the roles of active site residues, the lack of a plateau within the accessible pL range, or cooperative interactions between active site functional groups undergoing ionization. Here, we extend the integration of species distributions for alternative enzyme states in noncooperative models of acid/base catalysis into the GB equation, first used by Bevilacqua and colleagues for the HDV ribozyme, to develop a general population-weighted GB equation that allows simulation and global fitting of the three-dimensional relationship of the D₂O ratio (n) versus pL versus k_n/k₀. Simulations using the GPW-GB equation of PI results for RNase A, HDVrz, and VSrz illustrate that data obtained at multiple selected pL values across the pL-rate profile can assist in the planning and interpreting of solvent isotope effect experiments to distinguish alternative mechanistic models.

Theoretical Studies of the Self Cleavage Pistol Ribozyme Mechanism

Ribozymes are huge complex biological catalysts composed of a combination of RNA and proteins. Nevertheless, there is a reduced number of small ribozymes, the self-cleavage ribozymes, that are formed just by RNA and, apparently, they existed in cells of primitive biological systems. Unveiling the details of these “fossils” enzymes can contribute not only to the understanding of the origins of life but also to the development of new simplified artificial enzymes. A computational study of the reactivity of the pistol ribozyme carried out by means of classical MD simulations and QM/MM hybrid calculations is herein presented to clarify its catalytic mechanism. Analysis of the geometries along independent MD simulations with different protonation states of the active site basic species reveals that only the canonical system, with no additional protonation changes, renders reactive conformations. A change in the coordination sphere of the Mg2+ ion has been observed during the simulations, which allows proposing a mechanism to explain the unique mode of action of the pistol ribozyme by comparison with other ribozymes. The present results are at the center of the debate originated from recent experimental and theoretical studies on pistol ribozyme.

Selection of a self-cleaving ribozyme activated in a chemically and thermally denaturing environment

A new self-cleaving ribozyme was obtained from in vitro selection, displaying site-specific cleavage activity under denaturing conditions, such as high temperatures up to 95 °C, and denaturing solvents (20 M formamide). Adding salt such as Mg2+ also inhibited its activity. The conserved ribozyme sequence was found in the genome of several extremophiles, suggesting its potential biological relevance. This study provides an example of a ribozyme working under exotic conditions, which may expand the application of ribozymes in non-biological environments.

Functional Roles of Chelated Magnesium Ions in RNA Folding and Function

RNA regulates myriad cellular events such as transcription, translation, and splicing. To perform these essential functions, RNA often folds into complex tertiary structures in which its negatively charged ribose-phosphate backbone interacts with metal ions. Magnesium, the most abundant divalent metal ion in cells, neutralizes the backbone, thereby playing essential roles in RNA folding and function. This has been known for more than 50 years, and there are now thousands of in vitro studies, most of which have used ≥10 mM free Mg²⁺ ions to achieve optimal RNA folding and function. In the cell, however, concentrations of free Mg²⁺ ions are much lower, with most Mg²⁺ ions chelated by metabolites. In this Perspective, we curate data from a number of sources to provide extensive summaries of cellular concentrations of metabolites that bind Mg²⁺ and to estimate cellular concentrations of metabolite-chelated Mg²⁺ species, in the representative prokaryotic and eukaryotic systems Escherichia coli, Saccharomyces cerevisiae, and iBMK cells. Recent research from our lab and others has uncovered the fact that such weakly chelated Mg²⁺ ions can enhance RNA function, including its thermodynamic stability, chemical stability, and catalysis. We also discuss how metabolite-chelated Mg²⁺ complexes may have played roles in the origins of life. It is clear from this analysis that bound Mg²⁺ should not be simply considered non-RNA-interacting and that future RNA research, as well as protein research, could benefit from considering chelated magnesium.

Hydrated metal ion as a general acid in the catalytic mechanism of the 8-17 DNAzyme

The RNA-cleaving 8-17 DNAzyme, which is a metalloenzyme that depends on divalent metal ions for its function, is the most studied catalytic DNA in terms of its mechanism. By the end of 2017, a report of the crystal structure of the enzyme-substrate complex in the presence of Pb²⁺ probed some of the previous findings and opened new questions, especially around the participation of the metal ion in the catalytic mechanism and the promiscuity exhibited by the enzyme in terms of the metal cofactor required for catalysis. In this article we explore the role of the divalent metal ion in the mechanism of the 8-17 DNAzyme as a general acid, by measuring the influence of pH over the activity of a slower variant of the enzyme in the presence of Pb²⁺. We replaced G14, which has been identified as a general base in the mechanism of the enzyme, by the unnatural analog 2-aminopurine, with a lower pK_a value of the N1 group. With this approach, we obtained a bell-shaped pH-rate profile with experimental pK_a values of 5.4 and 7.0. Comparing these results with previous pH-rate profiles in the presence of Mg²⁺, our findings suggest the stabilization of the 5'-O leaving group by the hydrated metal ion acting as a general acid, in addition to the activation of the 2'-OH nucleophile by the general base G14.

Extended magnesium and calcium force field parameters for accurate ion-nucleic acid interactions in biomolecular simulations

Magnesium and calcium play an essential role in the folding and function of nucleic acids. To correctly describe their interactions with DNA and RNA in biomolecular simulations, an accurate parameterization is crucial. In most cases, the ion parameters are optimized based on a set of experimental solution properties such as solvation free energies, radial distribution functions, water exchange rates, and activity coefficient derivatives. However, the transferability of such bulk-optimized ion parameters to quantitatively describe biomolecular systems is limited. Here, we extend the applicability of our previous bulk-optimized parameters by including experimental binding affinities toward the phosphate oxygen on nucleic acids. In particular, we systematically adjust the combination rules that are an integral part of the pairwise interaction potentials of classical force fields. This allows us to quantitatively describe specific ion binding to nucleic acids without changing the solution properties in the most simple and efficient way. We show the advancement of the optimized Lorentz combination rule for two representative nucleic acid systems. For double-stranded DNA, the optimized combination rule for Ca²⁺ significantly improves the agreement with experiments, while the standard combination rule leads to unrealistically distorted DNA structures. For the add A-riboswitch, the optimized combination rule for Mg²⁺ improves the structure of two specifically bound Mg²⁺ ions as judged by the experimental distance to the binding site. Including experimental binding affinities toward specific ion binding sites on biomolecules, therefore, provides a promising perspective to develop a more accurate description of metal cations for biomolecular simulations.

Engineered extracellular vesicles as versatile ribonucleoprotein delivery vehicles for efficient and safe CRISPR genome editing

Transient delivery of CRISPR-based genome editing effectors is important to reduce off-target effects and immune responses. Recently extracellular vesicles (EVs) have been explored for Cas9 ribonucleoprotein (RNP) delivery. However, lack of mechanisms to enrich RNPs into EVs limited the efficiency of EVs as a RNP delivery vehicle. Here we describe a mechanism to actively enrich RNPs into EVs. We used the specific interaction between RNA aptamer and aptamer-binding protein (ABP) to enrich RNPs into EVs. We inserted RNA aptamer com into single guide RNA (sgRNA), and fused com-binding ABP Com to both termini of tetraspan protein CD63 that is abundant in exosomes. We found that the Com/com interaction enriched Cas9 and adenine base editor (ABE) RNPs into EVs, via forming a three-component complex including CD63-Com fusion protein, com-modified sgRNA and Cas9 or ABE. The RNP enriched EVs are efficient in genome editing and transiently expressed. The system is capable of delivering RNPs targeting multiple loci for multiplex genome editing. In addition, Cas9 from different species can be used together. The EV-delivered RNPs are active in vivo. The data show that the aptamer and ABP interactions can be utilized to actively enrich RNPs into EVs for improved genome editing efficiency and safety.

Photon Dissipation as the Origin of Information Encoding in RNA and DNA

Ultraviolet light incident on organic material can initiate its spontaneous dissipative structuring into chromophores which can catalyze their own replication. This may have been the case for one of the most ancient of all chromophores dissipating the Archean UVC photon flux, the nucleic acids. Oligos of nucleic acids with affinity to particular amino acids which foment UVC photon dissipation would most efficiently catalyze their own reproduction and thus would have been selected through non-equilibrium thermodynamic imperatives which favor dissipation. Indeed, we show here that those amino acids with characteristics most relevant to fomenting UVC photon dissipation are precisely those with greatest chemical affinity to their codons or anticodons. This could provide a thermodynamic basis for the specificity in the amino acid-nucleic acid interaction and an explanation for the accumulation of information in nucleic acids since this information is relevant to the optimization of dissipation of the externally imposed thermodynamic potentials. The accumulation of information in this manner provides a link between evolution and entropy production.

A Computational and Modeling Study of the Reaction Mechanism of Staphylococcus aureus Monoglycosyltransferase Reveals New Insights on the GT51 Family of Enzymes

Bacterial glycosyltransferases of the GT51 family are key enzymes in bacterial cell wall synthesis. Inhibiting cell wall synthesis is a very effective approach for development of antibiotics, as this can lead to either bacteriostatic or bactericidal effects. Even though the existence of this family has been known for over 50 years, only one potent inhibitor exists, which is an analog of the lipid IV product and derived from a natural product. Drug development focused on bacterial transglycosylase has been hampered due to little being know about its structure and reaction mechanism. In this study, Staphylococcus aureus monoglycosyltransferase was investigated at an atomistic level using computational methods. Classical molecular dynamics simulations were used to reveal information about the large-scale dynamics of the enzyme-substrate complex and the importance of magnesium in structure and function of the protein, while mixed mode quantum mechanics/molecular mechanics calculations unveiled a novel hypothesis for the reaction mechanism. From these results, we present a new model for the binding mode of lipid II and the reaction mechanism of the GT51 glycosyltransferases. A metal-bound hydroxide catalyzed reaction mechanism yields an estimated free energy barrier of 16.1 ± 1.0 kcal/mol, which is in line with experimental values. The importance of divalent cations is also further discussed. These findings could significantly aid targeted drug design, particularly the efficient development of transition state analogues as potential inhibitors for the GT51 glycosyltransferases.

CRISPR-Cas9 System for Plant Genome Editing: Current Approaches and Emerging Developments

Targeted genome editing using CRISPR-Cas9 has been widely adopted as a genetic engineering tool in various biological systems. This editing technology has been in the limelight due to its simplicity and versatility compared to other previously known genome editing platforms. Several modifications of this editing system have been established for adoption in a variety of plants, as well as for its improved efficiency and portability, bringing new opportunities for the development of transgene-free improved varieties of economically important crops. This review presents an overview of CRISPR-Cas9 and its application in plant genome editing. A catalog of the current and emerging approaches for the implementation of the system in plants is also presented with details on the existing gaps and limitations. Strategies for the establishment of the CRISPR-Cas9 molecular construct such as the selection of sgRNAs, PAM compatibility, choice of promoters, vector architecture, and multiplexing approaches are emphasized. Progress in the delivery and transgene detection methods, together with optimization approaches for improved on-target efficiency are also detailed in this review. The information laid out here will provide options useful for the effective and efficient exploitation of the system for plant genome editing and will serve as a baseline for further developments of the system. Future combinations and fine-tuning of the known parameters or factors that contribute to the editing efficiency, fidelity, and portability of CRISPR-Cas9 will indeed open avenues for new technological advancements of the system for targeted gene editing in plants.

Visualizing group II intron dynamics between the first and second steps of splicing

Group II introns are ubiquitous self-splicing ribozymes and retrotransposable elements evolutionarily and chemically related to the eukaryotic spliceosome, with potential applications as gene-editing tools. Recent biochemical and structural data have captured the intron in multiple conformations at different stages of catalysis. Here, we employ enzymatic assays, X-ray crystallography, and molecular simulations to resolve the spatiotemporal location and function of conformational changes occurring between the first and the second step of splicing. We show that the first residue of the highly-conserved catalytic triad is protonated upon 5’-splice-site scission, promoting a reversible structural rearrangement of the active site (toggling). Protonation and active site dynamics induced by the first step of splicing facilitate the progression to the second step. Our insights into the mechanism of group II intron splicing parallels functional data on the spliceosome, thus reinforcing the notion that these evolutionarily-related molecular machines share the same enzymatic strategy.

Group II introns are self-splicing ribozymes. Here, the authors employ enzymatic assay, X-ray crystallography and molecular dynamics simulations to show that protonation of the group II intron catalytic triad plays an important role for the transition from the first to the second step of splicing.

Insights into DNA catalysis from structural and functional studies of the 8-17 DNAzyme

The review examines functional knowledge gathered over two decades of research on the 8-17 DNAzyme, focusing on three aspects: the structural requirements for catalysis, the role of metal ions and the participation of general acid-base catalysis.

Evaluating the Catalytic Potential of a General RNA-Cleaving FANA Enzyme

The discovery of synthetic genetic polymers (XNAs) with catalytic activity demonstrates that natural genetic polymers are not unique in their ability to function as enzymes. However, all known examples of in vitro selected XNA enzymes function with lower activity than their natural counterparts, suggesting that XNAs might be limited in their ability to fold into structures with high catalytic activity. To explore this problem, we evaluated the catalytic potential of FANAzyme 12-7, an RNA-cleaving catalyst composed entirely of 2'-fluoroarabino nucleic acid (FANA) that was evolved to cleave RNA at a specific phosphodiester bond located between an unpaired guanine and a paired uracil in the substrate recognition arm. Here, we show that this activity extends to chimeric DNA substrates that contain a central riboguanosine (riboG) residue at the cleavage site. Surprisingly, FANAzyme 12-7 rivals known DNAzymes that were previously evolved to cleave chimeric DNA substrates under physiological conditions. These data provide convincing evidence that FANAzyme 12-7 maintains the catalytic potential of equivalent DNAzymes, which has important implications for the evolution of XNA catalysts and their contributions to future applications in synthetic biology.

Small Molecule Rescue and Glycosidic Conformational Analysis of the Twister Ribozyme

The number of self-cleaving ribozymes has increased sharply in recent years, giving rise to elaborations of the four known ribozyme catalytic strategies, α, β, γ, and δ. One such extension is utilized by the twister ribozyme, which is hypothesized to conduct δ, or general acid catalysis, via N3 of the syn adenine +1 nucleobase indirectly via buffer catalysis at biological pH and directly at lower pH. Herein, we test the δ catalysis role of A1 via chemical rescue and the catalytic relevance of the syn orientation of the nucleobase by conformational analysis. Using inhibited twister ribozyme variants with A1(N3) deaza or A1 abasic modifications, we observe >100-fold chemical rescue effects in the presence of protonatable biological small molecules such as imidazole and histidine, similar to observed rescue values previously reported for C75U/C76Δ in the HDV ribozyme. Brønsted plots for the twister variants support a model in which small molecules rescue catalytic activity via a proton transfer mechanism, suggesting that A1 in the wild type is involved in proton transfer, most likely general acid catalysis. Additionally, through glycosidic conformational analysis in an appropriate background that accommodates the bromine atom, we observe that an 8BrA1-modified twister ribozyme is up to 10-fold faster than a nonmodified A1 ribozyme, supporting crystallographic data that show that A1 is syn when conducting proton transfer. Overall, this study provides functional evidence that the nucleotide immediately downstream of the cleavage site participates directly or indirectly in general acid-base catalysis in the twister ribozyme while occupying the syn conformation.

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.

Cellular Concentrations of Nucleotide Diphosphate-Chelated Magnesium Ions Accelerate Catalysis by RNA and DNA Enzymes

RNAs are involved in myriad cellular events. In general, RNA function is affected by cellular conditions. For instance, molecular crowding promotes RNA folding through compaction of the RNA. Metabolites generally destabilize RNA secondary structure, which improves RNA folding cooperativity and increases the fraction of functional RNA. Our recent studies demonstrate that cellular concentrations of amino acid-chelated magnesium (aaCM) stimulate RNA folding and catalysis while protecting RNAs from magnesium ion-induced degradation. However, effects of other cellular magnesium ion chelators on RNA function have not been tested. Herein, we report that nucleotide diphosphate-chelated magnesium, which is of intermediate strength, promotes RNA catalysis much like aaCM. Nucleotides are some of the major metabolites in cells and have one to three phosphates, which have increasingly tight binding of magnesium. On the basis of binding calculations, ∼85% ATP, ∼40% ADP, and only 5% AMP are estimated to possess a magnesium ion under cellular conditions of 0.50 mM Mg²⁺_free. We tested the self-cleaving activity of the hammerhead ribozyme in the presence of these chelated magnesium species. Our results indicate that NTP-chelated magnesium and NMP-chelated magnesium do not appreciably stimulate RNA catalysis, whereas NDP-chelated magnesium promotes RNA catalysis up to 6.5-fold. Inspired by NDP, we observed similar stimulatory effects for several other Mg²⁺ diphosphate-containing metabolites, including UDP-GlcNAC and UDP-Glc; in addition, we found similar effects for a DNAzyme. Thus, rate stimulatory effects are general with respect to the diphosphate and nucleic acid enzyme. These results implicate magnesium-chelated diphosphate metabolites as general facilitators of RNA function inside cells.

The effect of adenine protonation on RNA phosphodiester backbone bond cleavage elucidated by deaza-nucleobase modifications and mass spectrometry

The catalytic strategies of small self-cleaving ribozymes often involve interactions between nucleobases and the ribonucleic acid (RNA) backbone. Here we show that multiply protonated, gaseous RNA has an intrinsic preference for the formation of ionic hydrogen bonds between adenine protonated at N3 and the phosphodiester backbone moiety on its 5'-side that facilitates preferential phosphodiester backbone bond cleavage upon vibrational excitation by low-energy collisionally activated dissociation. Removal of the basic N3 site by deaza-modification of adenine was found to abrogate preferential phosphodiester backbone bond cleavage. No such effects were observed for N1 or N7 of adenine. Importantly, we found that the pH of the solution used for generation of the multiply protonated, gaseous RNA ions by electrospray ionization affects phosphodiester backbone bond cleavage next to adenine, which implies that the protonation patterns in solution are at least in part preserved during and after transfer into the gas phase. Our study suggests that interactions between protonated adenine and phosphodiester moieties of RNA may play a more important mechanistic role in biological processes than considered until now.

Design of highly active double-pseudoknotted ribozymes: a combined computational and experimental study

Design of RNA sequences that adopt functional folds establishes principles of RNA folding and applications in biotechnology. Inverse folding for RNAs, which allows computational design of sequences that adopt specific structures, can be utilized for unveiling RNA functions and developing genetic tools in synthetic biology. Although many algorithms for inverse RNA folding have been developed, the pseudoknot, which plays a key role in folding of ribozymes and riboswitches, is not addressed in most algorithms. For the few algorithms that attempt to predict pseudoknot-containing ribozymes, self-cleavage activity has not been tested. Herein, we design double-pseudoknot HDV ribozymes using an inverse RNA folding algorithm and test their kinetic mechanisms experimentally. More than 90% of the positively designed ribozymes possess self-cleaving activity, whereas more than 70% of negative control ribozymes, which are predicted to fold to the necessary structure but with low fidelity, do not possess it. Kinetic and mutation analyses reveal that these RNAs cleave site-specifically and with the same mechanism as the WT ribozyme. Most ribozymes react just 50- to 80-fold slower than the WT ribozyme, and this rate can be improved to near WT by modification of a junction. Thus, fast-cleaving functional ribozymes with multiple pseudoknots can be designed computationally.

The Role of Proton Transfer on Mutations

Hydrogen bonds play a critical role in nucleobase studies as they encode genes, map protein structures, provide stability to the base pairs, and are involved in spontaneous and induced mutations. Proton transfer mechanism is a critical phenomenon that is related to the acid-base characteristics of the nucleobases in Watson-Crick base pairs. The energetic and dynamical behavior of the proton can be depicted from these characteristics and their adjustment to the water molecules or the surrounding ions. Further, new pathways open up in which protonated nucleobases are generated by proton transfer from the ionized water molecules and elimination of a hydroxyl radical in this review, the analysis will be focused on understanding the mechanism of untargeted mutations in canonical, wobble, Hoogsteen pairs, and mutagenic tautomers through the non-covalent interactions. Further, rare tautomer formation through the single proton transfer (SPT) and the double proton transfer (DPT), quantum tunneling in nucleobases, radiation-induced bystander effects, role of water in proton transfer (PT) reactions, PT in anticancer drugs-DNA interaction, displacement and oriental polarization, possible models for mutations in DNA, genome instability, and role of proton transfer using kinetic parameters for RNA will be discussed.

Classification of the nucleolytic ribozymes based upon catalytic mechanism

The nucleolytic ribozymes carry out site-specific RNA cleavage reactions by nucleophilic attack of the 2'-oxygen atom on the adjacent phosphorus with an acceleration of a million-fold or greater. A major part of this arises from concerted general acid-base catalysis. Recent identification of new ribozymes has expanded the group to a total of nine and this provides a new opportunity to identify sub-groupings according to the nature of the general base and acid. These include nucleobases, hydrated metal ions, and 2'-hydroxyl groups. Evolution has selected a number of different combinations of these elements that lead to efficient catalysis. These differences provide a new mechanistic basis for classifying these ribozymes.

Comparison of the Structures and Mechanisms of the Pistol and Hammerhead Ribozymes

Comparison of the secondary and three-dimensional structures of the hammerhead and pistol ribozymes reveals many close similarities, so in this work we have asked if they are mechanistically identical. We have determined a new crystal structure of the pistol ribozyme and have shown that G40 acts as general base in the cleavage reaction. The conformation in the active site ensures an in-line attack of the O2' nucleophile, and the conformation at the scissile phosphate and the position of the general base are closely similar to those in the hammerhead ribozyme. However, the two ribozymes differ in the nature of the general acid. 2'-Amino substitution experiments indicate that the general acid of the hammerhead ribozyme is the O2' of G8, while that of the pistol ribozyme is a hydrated metal ion. The two ribozymes are related but mechanistically distinct.

Evidence for Hidden Involvement of N3-Protonated Guanine in RNA Structure and Function

Charged nucleobases have been found to occur in several known RNA molecules and are considered essential for their structure and function. The mechanism of their involvement is however not yet fully understood. Revelation of the role of N7-protonated guanine, in modulating the geometry and stability of noncanonical base pairs formed through its unprotonated edges [Watson-Crick (WC) and sugar], has triggered the need to evaluate the feasibility of similar roles of other protonated nucleobases [Halder et al., Phys Chem Chem Phys, 2015, 17, 26249]. In this context, N3 protonation of guanine makes an interesting case as its influence on the charge distribution of the WC edge is similar to that of N7 protonation, though its thermodynamic cost of protonation is significantly higher. In this work, we have carried out structural bioinformatics analyses and quantum mechanics-based calculations to show that N3 protonation of guanine may take place in a cellular environment, at least in the G:C W:W Trans and G:G W:H Cis base pairs. Our results provide a reasonable starting point for future investigations in order to address the larger mechanistic question.

Template-directed RNA polymerization and enhanced ribozyme catalysis inside membraneless compartments formed by coacervates

Membraneless compartments, such as complex coacervates, have been hypothesized as plausible prebiotic micro-compartments due to their ability to sequester RNA; however, their compatibility with essential RNA World chemistries is unclear. We show that such compartments can enhance key prebiotically-relevant RNA chemistries. We demonstrate that template-directed RNA polymerization is sensitive to polycation identity, with polydiallyldimethylammonium chloride (PDAC) outperforming poly(allylamine), poly(lysine), and poly(arginine) in polycation/RNA coacervates. Differences in RNA diffusion rates between PDAC/RNA and oligoarginine/RNA coacervates imply distinct biophysical environments. Template-directed RNA polymerization is relatively insensitive to Mg²⁺ concentration when performed in PDAC/RNA coacervates as compared to buffer, even enabling partial rescue of the reaction in the absence of magnesium. Finally, we show enhanced activities of multiple nucleic acid enzymes including two ribozymes and a deoxyribozyme, underscoring the generality of this approach, in which functional nucleic acids like aptamers and ribozymes, and in some cases key cosolutes localize within the coacervate microenvironments.

Membraneless compartments have been theorized to be prebiotic micro-compartments as they spontaneously encapsulate RNA and proteins. Here, the authors report membraneless compartments can enhance RNA chemistries, affecting template directed RNA polymerization and stimulating nucleic acid enzymes.

Protonation of Nucleobases in Single- and Double-Stranded DNA

Single-stranded model oligodeoxyribonucleotides, each containing a single protonatable base-cytosine, adenine, guanine, or 5-methylcytosine-centrally located in a background of non-protonatable thymine residues, were acid-titrated in aqueous solution, with UV monitoring. The basicity of the central base was shown to depend on the type of the central base and its nearest neighbours and to rise with increasing oligonucleotide length and decreasing ionic strength of the solution. More complex model oligonucleotides, each containing a centrally located 5-methylcytosine base, were comparatively evaluated in single-stranded and double-stranded form, by UV spectroscopy and high-field NMR. The N³ protonation of the 5-methylcytosine moiety in the double-stranded case occurred at much lower pH, at which the duplex was already experiencing general dissociation, than in the single-stranded case. The central guanine:5-methylcytosine base pair remained intact up to this point, possibly due to an unusual alternative protonation on O² of the 5-methylcytosine moiety, already taking place at neutral or weakly basic pH, as indicated by UV spectroscopy, thus suggesting that 5-methylcytosine sites in double-stranded DNA might be protonated to a significant extent under physiological conditions.

Cellular Small Molecules Contribute to Twister Ribozyme Catalysis

The number of self-cleaving small ribozymes has increased sharply in recent years. Advances have been made in describing these ribozymes in terms of four catalytic strategies: α describes in-line attack, β describes neutralization of the nonbridging oxygens, γ describes activation of the nucleophile, and δ describes stabilization of the leaving group. Current literature presents the rapid self-cleavage of the twister ribozyme in terms of all four of these classic catalytic strategies. Herein, we describe the nonspecific contribution of small molecules to ribozyme catalysis. At biological pH, the rate of the wild-type twister ribozyme is enhanced up to 5-fold in the presence of moderate buffer concentrations, similar to the 3-5-fold effects reported previously for buffer catalysis for protein enzymes. We observe this catalytic enhancement not only with standard laboratory buffers, but also with diverse biological small molecules, including imidazole, amino acids, and amino sugars. Brønsted plots suggest that small molecules assist in proton transfer, most likely with δ catalysis. Cellular small molecules provide a simple way to overcome the limited functional diversity of RNA and have the potential to participate in the catalytic mechanisms of many ribozymes in vivo.

Evidence That Nucleophile Deprotonation Exceeds Bond Formation in the HDV Ribozyme Transition State

Steric constraints imposed by the active sites of protein and RNA enzymes pose major challenges to the investigation of structure-function relationships within these systems. As a strategy to circumvent such constraints in the HDV ribozyme, we have synthesized phosphoramidites from propanediol derivatives and incorporated them at the 5'-termini of RNA and DNA oligonucleotides to generate a series of novel substrates with nucleophiles perturbed electronically through geminal fluorination. In nonenzymatic, hydroxide-catalyzed intramolecular transphosphorylation of the DNA substrates, pH-rate profiles revealed that fluorine substitution reduces the maximal rate and the kinetic p K_a, consistent with the expected electron-withdrawing effect. In HDV ribozyme reactions, we observed that the RNA substrates undergo transphosphorylation relatively efficiently, suggesting that the conformational constraints imposed by a ribofuranose ring are not strictly required for ribozyme catalysis. In contrast to the nonenzymatic reactions, however, substrate fluorination modestly increases the ribozyme reaction rate, consistent with a mechanism in which (1) the 2'-hydroxyl nucleophile exists predominantly in its neutral, protonated form in the ground state and (2) the 2'-hydroxyl bears some negative charge in the rate-determining step, consistent with a transition state in which the extent of 2'-OH deprotonation exceeds the extent of P-O bond formation.

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.

Nucleic acids: function and potential for abiogenesis

The emergence of functional cooperation between the three main classes of biomolecules - nucleic acids, peptides and lipids - defines life at the molecular level. However, how such mutually interdependent molecular systems emerged from prebiotic chemistry remains a mystery. A key hypothesis, formulated by Crick, Orgel and Woese over 40 year ago, posits that early life must have been simpler. Specifically, it proposed that an early primordial biology lacked proteins and DNA but instead relied on RNA as the key biopolymer responsible not just for genetic information storage and propagation, but also for catalysis, i.e. metabolism. Indeed, there is compelling evidence for such an 'RNA world', notably in the structure of the ribosome as a likely molecular fossil from that time. Nevertheless, one might justifiably ask whether RNA alone would be up to the task. From a purely chemical perspective, RNA is a molecule of rather uniform composition with all four bases comprising organic heterocycles of similar size and comparable polarity and pK a values. Thus, RNA molecules cover a much narrower range of steric, electronic and physicochemical properties than, e.g. the 20 amino acid side-chains of proteins. Herein we will examine the functional potential of RNA (and other nucleic acids) with respect to self-replication, catalysis and assembly into simple protocellular entities.

Evidence of a General Acid-Base Catalysis Mechanism in the 8-17 DNAzyme

DNAzymes are catalytic DNA molecules that can perform a variety of reactions. Although advances have been made in obtaining DNAzymes via in vitro selection and many of them have been developed into sensors and imaging agents for metal ions, bacteria, and other molecules, the structural features responsible for these enzymatic reactions are still not well understood. Previous studies of the 8-17 DNAzyme have suggested conserved guanines close to the phosphodiester transfer site may play a role in the catalytic reaction. To identify the specific guanine and functional group of the guanine responsible for the reaction, we herein report the effects of replacing G1.1 and G14 (G; p Ka,N1 = 9.4) with analogues with a different p Ka at the N1 position, such as inosine (G14I; p Ka,N1 = 8.7), 2,6-diaminopurine (G14diAP; p Ka,N1 = 5.6), and 2-aminopurine (G14AP; p Ka,N1 = 3.8) on pH-dependent reaction rates. A comparison of the pH dependence of the reaction rates of these DNAzymes demonstrated that G14 in the bulge loop next to the cleavage site, is involved in proton transfer at the catalytic site. In contrast, we did not find any evidence of G1.1 being involved in acid-base catalysis. These results support general acid-base catalysis as a feasible strategy used in DNA catalysis, as in RNA and protein enzymes.