Imidazole rescue of a cytosine mutation in a self-cleaving ribozyme

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.

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.

New Insights into the Dependence of CPEB3 Ribozyme Cleavage on Mn2+ and Mg2

CPEB3 ribozyme is a self-cleaving RNA that occurs naturally in mammals and requires divalent metal ions for efficient activity. Ribozymes exhibit preferences for specific metal ions, but the exact differences in the catalytic mechanisms of various metal ions on the CPEB3 ribozyme remain unclear. Our findings reveal that Mn2+ functions as a more effective cofactor for CPEB3 ribozyme catalysis compared to Mg2+, as confirmed by its stronger binding affinity to CPEB3 by EPR. Cleavage assays of CPEB3 mutants and molecular docking analyses further showed that excessive Mn2+ ions can bind to a second binding site near the catalytic site, hindering CPEB3 catalytic efficiency and contributing to the Mn2+ bell-shaped curve. These results implicate a pivotal role for the local nucleobase-Mn2+ interactions in facilitating RNA folding and modulating the directed attack of nucleophilic reagents. Our study provides new insights and experimental evidence for exploring the divalent cation dependent cleavage mechanism of the CPEB3 ribozyme.

Ceg1 depletion reveals mechanisms governing degradation of non-capped RNAs in Saccharomyces cerevisiae

Most functional eukaryotic mRNAs contain a 5' 7-methylguanosine (m7G) cap. Although capping is essential for many biological processes including mRNA processing, export and translation, the fate of uncapped transcripts has not been studied extensively. Here, we employed fast nuclear depletion of the capping enzymes in Saccharomyces cerevisiae to uncover the turnover of the transcripts that failed to be capped. We show that although the degradation of cap-deficient mRNA is dominant, the levels of hundreds of non-capped mRNAs increase upon depletion of the capping enzymes. Overall, the abundance of non-capped mRNAs is inversely correlated to the expression levels, altogether resembling the effects observed in cells lacking the cytoplasmic 5'-3' exonuclease Xrn1 and indicating differential degradation fates of non-capped mRNAs. The inactivation of the nuclear 5'-3' exonuclease Rat1 does not rescue the non-capped mRNA levels indicating that Rat1 is not involved in their degradation and consequently, the lack of the capping does not affect the distribution of RNA Polymerase II on the chromatin. Our data indicate that the cap presence is essential to initiate the Xrn1-dependent degradation of mRNAs underpinning the role of 5' cap in the Xrn1-dependent buffering of the cellular mRNA levels.

Structural basis for substrate binding and catalysis by a self-alkylating ribozyme

Ribozymes that react with small-molecule probes have important applications in transcriptomics and chemical biology, such as RNA labeling and imaging. Understanding the structural basis for these RNA-modifying reactions will enable the development of better tools for studying RNA. Nevertheless, high-resolution structures and underlying catalytic mechanisms for members of this ribozyme class remain elusive. Here, we focus on a self-alkylating ribozyme that catalyzes nitrogen-carbon bond formation between a specific guanine and a 2,3-disubstituted epoxide substrate and report the crystal structures of a self-alkylating ribozyme, including both alkylated and apo forms, at 1.71-Å and 2.49-Å resolution, respectively. The ribozyme assumes an elongated hairpin-like architecture preorganized to accommodate the epoxide substrate in a hook-shaped conformation. Observed reactivity of substrate analogs together with an inverse, log-linear pH dependence of the reaction rate suggests a requirement for epoxide protonation, possibly assisted by the ether oxygens within the substrate.

Asymmetric Roles of Two Histidine Residues in Streptococcus pyogenes Cas9 Catalytic Domains upon Chemical Rescue

CRISPR-Cas9 technology has been at the forefront of the field of biology. The Streptococcus pyogenes (SpyCas9) protein forms a complex with guide RNA and can recognize and cleave double-stranded DNA through hybridization based on 20 base pairings. SpyCas9 has two nuclease domains, HNH and RuvC, each of which cuts each DNA strand, and both contain critical histidine residues. Although previously reported crystal structures provide useful geometric information, the extent to which these residues functionally contribute to catalysis is unknown. Here, we mutated histidine residues on HNH and RuvC domains to alanine or glycine and attempted to rescue the enzymatic activity by adding the imidazole molecule, using an in vitro DNA cleavage assay. H840A and H840G exhibited rescued enzymatic activity on the HNH domain following imidazole addition, suggesting that H840 acts as a general base. We also tested various chemicals and found that the pK_a of imidazole derivatives, and not their molecular shape, correlated with the rescue effect. In contrast, both H983A and H983G on the RuvC domain did not exhibit a rescue effect following imidazole addition. Our chemical rescue approach will provide crucial insight into understanding Cas9 catalysis, complementing structural analyses.

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.

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.

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.

Evolutionary convergence in the biosyntheses of the imidazole moieties of histidine and purines

Background

The imidazole group is an ubiquitous chemical motif present in several key types of biomolecules. It is a structural moiety of purines, and plays a central role in biological catalysis as part of the side-chain of histidine, the amino acid most frequently found in the catalytic site of enzymes. Histidine biosynthesis starts with both ATP and the pentose phosphoribosyl pyrophosphate (PRPP), which is also the precursor for the de novo synthesis of purines. These two anabolic pathways are also connected by the imidazole intermediate 5-aminoimidazole-4-carboxamide ribotide (AICAR), which is synthesized in both routes but used only in purine biosynthesis. Rather surprisingly, the imidazole moieties of histidine and purines are synthesized by different, non-homologous enzymes. As discussed here, this phenomenon can be understood as a case of functional molecular convergence.

Results

In this work, we analyze these polyphyletic processes and argue that the independent origin of the corresponding enzymes is best explained by the differences in the function of each of the molecules to which the imidazole moiety is attached. Since the imidazole present in histidine is a catalytic moiety, its chemical arrangement allows it to act as an acid or a base. On the contrary, the de novo biosynthesis of purines starts with an activated ribose and all the successive intermediates are ribotides, with the key β-glycosidic bondage joining the ribose and the imidazole moiety. This prevents purine ribonucleotides to exhibit any imidazole-dependent catalytic activity, and may have been the critical trait for the evolution of two separate imidazole-synthesizing-enzymes. We also suggest that, in evolutionary terms, the biosynthesis of purines predated that of histidine.

Conclusions

As reviewed here, other biosynthetic routes for imidazole molecules are also found in extant metabolism, including the autocatalytic cyclization that occurs during the formation of creatinine from creatine phosphate, as well as the internal cyclization of the Ala-Ser-Gly motif of some members of the ammonia-lyase and aminomutase families, that lead to the MIO cofactor. The diversity of imidazole-synthesizing pathways highlights the biological significance of this key chemical group, whose biosyntheses evolved independently several times.

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.

Allosteric Modulation of the Faecalibacterium prausnitzii Hepatitis Delta Virus-like Ribozyme by Glucosamine 6-Phosphate: The Substrate of the Adjacent Gene Product

Self-cleaving ribozymes were discovered 30 years ago and have been found throughout nature, from bacteria to animals, but little is known about their biological functions and regulation, particularly how cofactors and metabolites alter their activity. A hepatitis delta virus-like self-cleaving ribozyme maps upstream of a phosphoglucosamine mutase (glmM) open reading frame in the genome of the human gut bacterium Faecalibacterium prausnitzii. The presence of a ribozyme in the untranslated region of glmM suggests a regulation mechanism of gene expression. In the bacterial hexosamine biosynthesis pathway, the enzyme glmM catalyzes the isomerization of glucosamine 6-phosphate into glucosamine 1-phosphate. In this study, we investigated the effect of these metabolites on the co-transcriptional self-cleavage rate of the ribozyme. Our results suggest that glucosamine 6-phosphate, but not glucosamine 1-phosphate, is an allosteric ligand that increases the self-cleavage rate of drz-Fpra-1, providing the first known example of allosteric modulation of a self-cleaving ribozyme by the substrate of the adjacent gene product. Given that the ribozyme is activated by the glmM substrate, but not the product, this allosteric modulation may represent a potential feed-forward mechanism of gene expression regulation in bacteria.

Can an Imidazole Be Formed from an Alanyl-Seryl-Glycine Tripeptide under Possible Prebiotic Conditions?

The five-membered heterocyclic imidazole group, which is an essential component of purines, histidine and many cofactors, has been abiotically synthesized in different model experiments that attempt to simulate the prebiotic environment. The evolutionary significance of imidazoles is highlighted not only by its presence in nucleic acid components and in histidine, but also by experimental reports of its ability to restore the catalytic activity of ribozymes. However, as of today there are no reports of histidine in carbonaceous chondrites, and although the abiotic synthesis of His reported by Shen et al. (1987, 1990a) proceeds via an Amadori rearrangement, like in the biosynthesis of histidine, neither the reactants nor the conditions are truly prebiotic. Based on the autocatalytic biosynthesis of 4-methylidene-imidazole-one (MIO), a cofactor of some members of the amino acid aromatic ammonia-lyases and aminomutases, which occur via the self-condensation of a simple Ala-Ser-Gly motif within the sequence of the enzymes, we propose a possible prebiotic synthesis of an imidazolide.

How RNA acts as a nuclease: some mechanistic comparisons in the nucleolytic ribozymes

Recent structural and mechanistic studies have shed considerable light on the catalytic mechanisms of nucleolytic ribozymes. The discovery of several new ribozymes in this class has now allowed comparisons to be made, and the beginnings of mechanistic groupings to emerge.

Inducible Bcl-2 gene RNA interference mediated by aptamer-integrated HDV ribozyme switch

The regulation of RNA interference (RNAi) could be a powerful method for the study of temporal and dose dependent effects of gene expression. In this study, we designed the hepatitis delta virus (HDV) ribozyme with an embedded theophylline aptamer as the sensor domain and the pri-miRNA of endogenous gene Bcl-2 as the effector domain to engineer an RNAi-regulatory device in MCF-7 cells. The system allowed us to control gene expression by adding theophylline into the culture media in a dose dependent fashion. This is the pioneering application of ribozyme switches to activate RNAi for modulating endogenous genes in mammalian cells. The platform sets the stage for investigations of other endogenous genes that regulate various biological functions such as differentiation, cell division or cell death, and provides a promising interface with other universal RNAi-based decision-making circuits that operate in mammalian cells. It can be used to study more genes associated with cancer and screen for potential drug targets for gene therapy.

DNA-Catalyzed Amide Hydrolysis

DNA catalysts (deoxyribozymes) for a variety of reactions have been identified by in vitro selection. However, for certain reactions this identification has not been achieved. One important example is DNA-catalyzed amide hydrolysis, for which a previous selection experiment instead led to DNA-catalyzed DNA phosphodiester hydrolysis. Subsequent efforts in which the selection strategy deliberately avoided phosphodiester hydrolysis led to DNA-catalyzed ester and aromatic amide hydrolysis, but aliphatic amide hydrolysis has been elusive. In the present study, we show that including modified nucleotides that bear protein-like functional groups (any one of primary amino, carboxyl, or primary hydroxyl) enables identification of amide-hydrolyzing deoxyribozymes. In one case, the same deoxyribozyme sequence without the modifications still retains substantial catalytic activity. Overall, these findings establish the utility of introducing protein-like functional groups into deoxyribozymes for identifying new catalytic function. The results also suggest the longer-term feasibility of deoxyribozymes as artificial proteases.

In vitro selection of deoxyribozymes active with Cd(2+) ions resulting in variants of DNAzyme 8-17

In vitro selection was performed to search for RNA-cleaving DNAzymes catalytically active with Cd(2+) ions from the oligonucleotide combinatorial library with a 23-nucleotide random region. All the selected, catalytically active variants turned out to belong to the 8-17 type DNAzyme. Three DNAzymes were prepared in shortened, cis-acting versions which were subjected to a detailed study of the kinetic properties and metal ion preferences. Although the selection protocol was designed for Cd(2+)-dependent DNAzymes, the variants showed broader metal ion specificity. They preferred Cd(2+) but were also active with Mn(2+) and Zn(2+), suggesting that binding of the catalytic ion does not require an extremely specific coordination pattern. The unexpected decrease of the catalytic activity of the variants along with the temperature increase suggested that some changes occurred in their structures or the rate-limiting step of the reaction was changed. Two elements of the catalytic core of DNAzyme 1/VIIWS, the nucleotide at position 12 and the three-base-pair hairpin motif, were mutated. The presence of a purine residue at position 12 was crucial for the catalytic activity but the changes at that position had a relatively small influence on the metal ion preferences of this variant. The middle base pair of the three-base-pair hairpin was changed from A-T to C-G interaction. The catalytic activity of the mutated variant was increased with Zn(2+), decreased with Mn(2+), and was not changed in the presence of Cd(2+) ions. Clearly, this base pair was important for defining the metal ion preferences of the DNAzyme 1/VIIWS.

Wobble pairs of the HDV ribozyme play specific roles in stabilization of active site dynamics

The hepatitis delta virus (HDV) is the only known human pathogen whose genome contains a catalytic RNA motif (ribozyme). The overall architecture of the HDV ribozyme is that of a double-nested pseudoknot, with two GU pairs flanking the active site. Although extensive studies have shown that mutation of either wobble results in decreased catalytic activity, little work has focused on linking these mutations to specific structural effects on catalytic fitness. Here we use molecular dynamics simulations based on an activated structure to probe the active site dynamics as a result of wobble pair mutations. In both wild-type and mutant ribozymes, the in-line fitness of the active site (as a measure of catalytic proficiency) strongly depends on the presence of a C75(N3H3+)N1(O5') hydrogen bond, which positions C75 as the general acid for the reaction. Our mutational analyses show that each GU wobble supports catalytically fit conformations in distinct ways; the reverse G25U20 wobble promotes high in-line fitness, high occupancy of the C75(N3H3+)G1(O5') general-acid hydrogen bond and stabilization of the G1U37 wobble, while the G1U37 wobble acts more locally by stabilizing high in-line fitness and the C75(N3H3+)G1(O5') hydrogen bond. We also find that stable type I A-minor and P1.1 hydrogen bonding above and below the active site, respectively, prevent local structural disorder from spreading and disrupting global conformation. Taken together, our results define specific, often redundant architectural roles for several structural motifs of the HDV ribozyme active site, expanding the known roles of these motifs within all HDV-like ribozymes and other structured RNAs.

Assessment of metal-assisted nucleophile activation in the hepatitis delta virus ribozyme from molecular simulation and 3D-RISM

The hepatitis delta virus ribozyme is an efficient catalyst of RNA 2'-O-transphosphorylation and has emerged as a key experimental system for identifying and characterizing fundamental features of RNA catalysis. Recent structural and biochemical data have led to a proposed mechanistic model whereby an active site Mg(2+) ion facilitates deprotonation of the O2' nucleophile, and a protonated cytosine residue (C75) acts as an acid to donate a proton to the O5' leaving group as noted in a previous study. This model assumes that the active site Mg(2+) ion forms an inner-sphere coordination with the O2' nucleophile and a nonbridging oxygen of the scissile phosphate. These contacts, however, are not fully resolved in the crystal structure, and biochemical data are not able to unambiguously exclude other mechanistic models. In order to explore the feasibility of this model, we exhaustively mapped the free energy surfaces with different active site ion occupancies via quantum mechanical/molecular mechanical (QM/MM) simulations. We further incorporate a three-dimensional reference interaction site model for the solvated ion atmosphere that allows these calculations to consider not only the rate associated with the chemical steps, but also the probability of observing the system in the presumed active state with the Mg(2+) ion bound. The QM/MM results predict that a pathway involving metal-assisted nucleophile activation is feasible based on the rate-controlling transition state barrier departing from the presumed metal-bound active state. However, QM/MM results for a similar pathway in the absence of Mg(2+) are not consistent with experimental data, suggesting that a structural model in which the crystallographically determined Mg(2+) is simply replaced with Na(+) is likely incorrect. It should be emphasized, however, that these results hinge upon the assumption of the validity of the presumed Mg(2+)-bound starting state, which has not yet been definitively verified experimentally, nor explored in depth computationally. Thus, further experimental and theoretical study is needed such that a consensus view of the catalytic mechanism emerges.

Facile and rapid deprotection conditions for the cleavage of synthetic oligonucleotides from 1,4-anhydroerythritol-based universal polymer support

In our previous report [Kumar, P.; Dhawan, G.; Chandra, R.; Gupta, K.C. Polyamine-assisted rapid and clean cleavage of oligonucleotides from cis-diol bearing universal support. Nucl. Acids Res. 2002, 30, e130 (1-8)], we demonstrated polyamine-mediated deprotection of oligonucleotides from cis-diol group bearing universal polymer support (I). However, vulnerability of the conventional dC(bz) to modifications under these conditions compelled us to employ dC(ac) during synthesis of oligonucleotide using conventional synthons. Here, a new set of simple and rapid deprotection conditions has been developed for the complete cleavage of oligonucleotides from the 1,4-anhydroerythritol-based universal polymer support employing conventional dC(bz) synthon. Using manganese-imidazole complex in aqueous ammonium hydroxide (∼ 30%), fully deprotected oligonucleotide sequences were obtained in 40 min, which were analyzed on reverse phase-HPLC and compared with the standard oligomers in terms of their retention time. Finally, their biological compatibility was established by analyzing PCR amplified products of npsA gene of N. meningitidis.

Transition State Features in the Hepatitis Delta Virus Ribozyme Reaction Revealed by Atomic Perturbations

Endonucleolytic ribozymes constitute a class of non-coding RNAs that catalyze single-strand RNA scission. With crystal structures available for all of the known ribozymes, a major challenge involves relating functional data to the physically observed RNA architecture. In the case of the hepatitis delta virus (HDV) ribozyme, there are three high-resolution crystal structures, the product state of the reaction and two precursor variants, with distinct mechanistic implications. Here, we develop new strategies to probe the structure and catalytic mechanism of a ribozyme. First, we use double-mutant cycles to distinguish differences in functional group proximity implicated by the crystal structures. Second, we use a corrected form of the Brønsted equation to assess the functional significance of general acid catalysis in the system. Our results delineate the functional relevance of atomic interactions inferred from structure, and suggest that the HDV ribozyme transition state resembles the cleavage product in the degree of proton transfer to the leaving group.

Mechanistic Studies Reveal Similar Catalytic Strategies for Phosphodiester Bond Hydrolysis by Protein-only and RNA-dependent Ribonuclease P

Ribonuclease P (RNase P) is an endonuclease that catalyzes the essential removal of the 5' end of tRNA precursors. Until recently, all identified RNase P enzymes were a ribonucleoprotein with a conserved catalytic RNA component. However, the discovery of protein-only RNase P (PRORP) shifted this paradigm, affording a unique opportunity to compare mechanistic strategies used by naturally evolved protein and RNA-based enzymes that catalyze the same reaction. Here we investigate the enzymatic mechanism of pre-tRNA hydrolysis catalyzed by the NYN (Nedd4-BP1, YacP nuclease) metallonuclease of Arabidopsis thaliana, PRORP1. Multiple and single turnover kinetic data support a mechanism where a step at or before chemistry is rate-limiting and provide a kinetic framework to interpret the results of metal alteration, mutations, and pH dependence. Catalytic activity has a cooperative dependence on the magnesium concentration (nH = 2) under kcat/Km conditions, suggesting that PRORP1 catalysis is optimal with at least two active site metal ions, consistent with the crystal structure. Metal rescue of Asp-to-Ala mutations identified two aspartates important for enhancing metal ion affinity. The single turnover pH dependence of pre-tRNA cleavage revealed a single ionization (pKa ∼ 8.7) important for catalysis, consistent with deprotonation of a metal-bound water nucleophile. The pH and metal dependence mirrors that observed for the RNA-based RNase P, suggesting similar catalytic mechanisms. Thus, despite different macromolecular composition, the RNA and protein-based RNase P act as dynamic scaffolds for the binding and positioning of magnesium ions to catalyze phosphodiester bond hydrolysis.

Inverse thio effects in the hepatitis delta virus ribozyme reveal that the reaction pathway is controlled by metal ion charge density

The hepatitis delta virus (HDV) ribozyme self-cleaves in the presence of a wide range of monovalent and divalent ions. Prior theoretical studies provided evidence that self-cleavage proceeds via a concerted or stepwise pathway, with the outcome dictated by the valency of the metal ion. In the present study, we measure stereospecific thio effects at the nonbridging oxygens of the scissile phosphate under a wide range of experimental conditions, including varying concentrations of diverse monovalent and divalent ions, and combine these with quantum mechanical/molecular mechanical (QM/MM) free energy simulations on the stereospecific thio substrates. The R_P substrate gives large normal thio effects in the presence of all monovalent ions. The S_P substrate also gives normal or no thio effects, but only for smaller monovalent and divalent cations, such as Li⁺, Mg²⁺, Ca²⁺, and Sr²⁺; in contrast, sizable inverse thio effects are found for larger monovalent and divalent cations, including Na⁺, K⁺, NH₄⁺, and Ba²⁺. Proton inventories are found to be unity in the presence of the larger monovalent and divalent ions, but two in the presence of Mg²⁺. Additionally, rate–pH profiles are inverted for the low charge density ions, and only imidazole plus ammonium ions rescue an inactive C75Δ variant in the absence of Mg²⁺. Results from the thio effect experiments, rate–pH profiles, proton inventories, and ammonium/imidazole rescue experiments, combined with QM/MM free energy simulations, support a change in the mechanism of HDV ribozyme self-cleavage from concerted and metal ion-stabilized to stepwise and proton transfer-stabilized as the charge density of the metal ion decreases.

The role of an active site Mg(2+) in HDV ribozyme self-cleavage: insights from QM/MM calculations

The hepatitis delta virus (HDV) ribozyme is a catalytic RNA motif embedded in the human pathogenic HDV RNA. It catalyzes self-cleavage of its sugar-phosphate backbone with direct participation of the active site cytosine C75. Biochemical and structural data support a general acid role of C75. Here, we used hybrid quantum mechanical/molecular mechanical (QM/MM) calculations to probe the reaction mechanism and changes in Gibbs energy along the ribozyme's reaction pathway with an N3-protonated C75H(+) in the active site, which acts as the general acid, and a partially hydrated Mg(2+) ion with one deprotonated, inner-shell coordinated water molecule that acts as the general base. We followed eight reaction paths with a distinct position and coordination of the catalytically important active site Mg(2+) ion. For six of them, we observed feasible activation barriers ranging from 14.2 to 21.9 kcal mol(-1), indicating that the specific position of the Mg(2+) ion in the active site is predicted to strongly affect the kinetics of self-cleavage. The deprotonation of the U-1(2'-OH) nucleophile and the nucleophilic attack of the resulting U-1(2'-O(-)) on the scissile phosphodiester are found to be separate steps, as deprotonation precedes the nucleophilic attack. This sequential mechanism of the HDV ribozyme differs from the concerted nucleophilic activation and attack suggested for the hairpin ribozyme. We estimate the pKa of the U-1(2'-OH) group to range from 8.8 to 11.2, suggesting that it is lowered by several units from that of a free ribose, comparable to and most likely smaller than the pKa of the solvated active site Mg(2+) ion. Our results thus support the notion that the structure of the HDV ribozyme, and particularly the positioning of the active site Mg(2+) ion, facilitate deprotonation and activation of the 2'-OH nucleophile.

New tools provide a second look at HDV ribozyme structure, dynamics and cleavage

The hepatitis delta virus (HDV) ribozyme is a self-cleaving RNA enzyme essential for processing viral transcripts during rolling circle viral replication. The first crystal structure of the cleaved ribozyme was solved in 1998, followed by structures of uncleaved, mutant-inhibited and ion-complexed forms. Recently, methods have been developed that make the task of modeling RNA structure and dynamics significantly easier and more reliable. We have used ERRASER and PHENIX to rebuild and re-refine the cleaved and cis-acting C75U-inhibited structures of the HDV ribozyme. The results correct local conformations and identify alternates for RNA residues, many in functionally important regions, leading to improved R values and model validation statistics for both structures. We compare the rebuilt structures to a higher resolution, trans-acting deoxy-inhibited structure of the ribozyme, and conclude that although both inhibited structures are consistent with the currently accepted hammerhead-like mechanism of cleavage, they do not add direct structural evidence to the biochemical and modeling data. However, the rebuilt structures (PDBs: 4PR6, 4PRF) provide a more robust starting point for research on the dynamics and catalytic mechanism of the HDV ribozyme and demonstrate the power of new techniques to make significant improvements in RNA structures that impact biologically relevant conclusions.

Disparate HDV ribozyme crystal structures represent intermediates on a rugged free-energy landscape

The hepatitis delta virus (HDV) ribozyme is a member of the class of small, self-cleaving catalytic RNAs found in a wide range of genomes from HDV to human. Both pre- and post-catalysis (precursor and product) crystal structures of the cis-acting genomic HDV ribozyme have been determined. These structures, together with extensive solution probing, have suggested that a significant conformational change accompanies catalysis. A recent crystal structure of a trans-acting precursor, obtained at low pH and by molecular replacement from the previous product conformation, conforms to the product, raising the possibility that it represents an activated conformer past the conformational change. Here, using fluorescence resonance energy transfer (FRET), we discovered that cleavage of this ribozyme at physiological pH is accompanied by a structural lengthening in magnitude comparable to previous trans-acting HDV ribozymes. Conformational heterogeneity observed by FRET in solution appears to have been removed upon crystallization. Analysis of a total of 1.8 µsec of molecular dynamics (MD) simulations showed that the crystallographically unresolved cleavage site conformation is likely correctly modeled after the hammerhead ribozyme, but that crystal contacts and the removal of several 2'-oxygens near the scissile phosphate compromise catalytic in-line fitness. A cis-acting version of the ribozyme exhibits a more dynamic active site, while a G-1 residue upstream of the scissile phosphate favors poor fitness, allowing us to rationalize corresponding changes in catalytic activity. Based on these data, we propose that the available crystal structures of the HDV ribozyme represent intermediates on an overall rugged RNA folding free-energy landscape.

Position-specific modification with imidazolyl group on10-23 DNAzyme realized catalytic activity enhancement

Nucleoside analogues with imidazolyl and histidinyl groups were synthesized for site-specific modification on the catalytic core of 10-23 DNAzyme. The distinct position-dependent effect of imidazolyl group was observed. Positive effect at A9 position was always observed. The pH- and Mg(2+)-dependence of the imidazolyl-modified DNAzymes suggested that imidazolyl group in 10-23 DNAzyme probably plays a dual role, its hydrogen bonding ability and spacial occupation play the favorable influence on the catalytic conformation of the modified DNAzymes. This research demonstrated that the catalytic performance of DNAzymes could be enhanced by incorporation of additional functional groups. Chemical modification is a feasible approach toward more efficient DNAzymes for therapeutic and biotechnological applications.

The noncoding RNA revolution-trashing old rules to forge new ones

Noncoding RNAs (ncRNAs) accomplish a remarkable variety of biological functions. They regulate gene expression at the levels of transcription, RNA processing, and translation. They protect genomes from foreign nucleic acids. They can guide DNA synthesis or genome rearrangement. For ribozymes and riboswitches, the RNA structure itself provides the biological function, but most ncRNAs operate as RNA-protein complexes, including ribosomes, snRNPs, snoRNPs, telomerase, microRNAs, and long ncRNAs. Many, though not all, ncRNAs exploit the power of base pairing to selectively bind and act on other nucleic acids. Here, we describe the pathway of ncRNA research, where every established "rule" seems destined to be overturned.

Replication and packaging of Turnip yellow mosaic virus RNA containing Flock house virus RNA1 sequence

Turnip yellow mosaic virus (TYMV) is a spherical plant virus that has a single 6.3 kb positive strand RNA as a genome. In this study, RNA1 sequence of Flock house virus (FHV) was inserted into the TYMV genome to test whether TYMV can accommodate and express another viral entity. In the resulting construct, designated TY-FHV, the FHV RNA1 sequence was expressed as a TYMV subgenomic RNA. Northern analysis of the Nicotiana benthamiana leaves agroinfiltrated with the TY-FHV showed that both genomic and subgenomic FHV RNAs were abundantly produced. This indicates that the FHV RNA1 sequence was correctly expressed and translated to produce a functional FHV replicase. Although these FHV RNAs were not encapsidated, the FHV RNA having a TYMV CP sequence at the 3'-end was efficiently encapsidated. When an eGFP gene was inserted into the B2 ORF of the FHV sequence, a fusion protein of B2-eGFP was produced as expected.

HDV family of self-cleaving ribozymes

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

Quantum mechanical/molecular mechanical free energy simulations of the self-cleavage reaction in the hepatitis delta virus ribozyme

The hepatitis delta virus (HDV) ribozyme catalyzes a self-cleavage reaction using a combination of nucleobase and metal ion catalysis. Both divalent and monovalent ions can catalyze this reaction, although the rate is slower with monovalent ions alone. Herein, we use quantum mechanical/molecular mechanical (QM/MM) free energy simulations to investigate the mechanism of this ribozyme and to elucidate the roles of the catalytic metal ion. With Mg²⁺ at the catalytic site, the self-cleavage mechanism is observed to be concerted with a phosphorane-like transition state and a free energy barrier of ∼13 kcal/mol, consistent with free energy barrier values extrapolated from experimental studies. With Na⁺ at the catalytic site, the mechanism is observed to be sequential, passing through a phosphorane intermediate, with free energy barriers of 2–4 kcal/mol for both steps; moreover, proton transfer from the exocyclic amine of protonated C75 to the nonbridging oxygen of the scissile phosphate occurs to stabilize the phosphorane intermediate in the sequential mechanism. To explain the slower rate observed experimentally with monovalent ions, we hypothesize that the activation of the O2′ nucleophile by deprotonation and orientation is less favorable with Na⁺ ions than with Mg²⁺ ions. To explore this hypothesis, we experimentally measure the pK_a of O2′ by kinetic and NMR methods and find it to be lower in the presence of divalent ions rather than only monovalent ions. The combined theoretical and experimental results indicate that the catalytic Mg²⁺ ion may play three key roles: assisting in the activation of the O2′ nucleophile, acidifying the general acid C75, and stabilizing the nonbridging oxygen to prevent proton transfer to it.

Controlling mammalian gene expression by allosteric hepatitis delta virus ribozymes

We engineered small molecule responsive allosteric ribozymes based on the genomic hepatitis delta virus (HDV) ribozyme by replacing the P4-L4 stem-loop with an RNA aptamer through a connector stem. When embedded in the 3′ untranslated region of a reporter gene mRNA, these RNA devices enabled regulation of cis-gene expression by theophylline and guanine by up to 29.5-fold in mammalian cell culture. Furthermore, a NOR logic gate device was constructed by placing two engineered ribozymes in tandem, demonstrating the modularity of the RNA devices. The significant improvement in the regulatory dynamic range (ON/OFF ratio) of the RNA devices based on the HDV ribozyme should provide new opportunities for practical applications.

Thio effects and an unconventional metal ion rescue in the genomic hepatitis delta virus ribozyme

Metal ion and nucleobase catalysis are important for ribozyme mechanism, but the extent to which they cooperate is unclear. A crystal structure of the hepatitis delta virus (HDV) ribozyme suggested that the pro-RP oxygen at the scissile phosphate directly coordinates a catalytic Mg(2+) ion and is within hydrogen bonding distance of the amine of the general acid C75. Prior studies of the genomic HDV ribozyme, however, showed neither a thio effect nor metal ion rescue using Mn(2+). Here, we combine experiment and theory to explore phosphorothioate substitutions at the scissile phosphate. We report significant thio effects at the scissile phosphate and metal ion rescue with Cd(2+). Reaction profiles with an SP-phosphorothioate substitution are indistinguishable from those of the unmodified substrate in the presence of Mg(2+) or Cd(2+), supporting the idea that the pro-SP oxygen does not coordinate metal ions. The RP-phosphorothioate substitution, however, exhibits biphasic kinetics, with the fast-reacting phase displaying a thio effect of up to 5-fold and the slow-reacting phase displaying a thio effect of ~1000-fold. Moreover, the fast- and slow-reacting phases give metal ion rescues in Cd(2+) of up to 10- and 330-fold, respectively. The metal ion rescues are unconventional in that they arise from Cd(2+) inhibiting the oxo substrate but not the RP substrate. This metal ion rescue suggests a direct interaction of the catalytic metal ion with the pro-RP oxygen, in line with experiments with the antigenomic HDV ribozyme. Experiments without divalent ions, with a double mutant that interferes with Mg(2+) binding, or with C75 deleted suggest that the pro-RP oxygen plays at most a redundant role in positioning C75. Quantum mechanical/molecular mechanical (QM/MM) studies indicate that the metal ion contributes to catalysis by interacting with both the pro-RP oxygen and the nucleophilic 2'-hydroxyl, supporting the experimental findings.

The DNA and RNA sugar-phosphate backbone emerges as the key player. An overview of quantum-chemical, structural biology and simulation studies

Knowledge of geometrical and physico-chemical properties of the sugar-phosphate backbone substantially contributes to the comprehension of the structural dynamics, function and evolution of nucleic acids. We provide a side by side overview of structural biology/bioinformatics, quantum chemical and molecular mechanical/simulation studies of the nucleic acids backbone. We highlight main features, advantages and limitations of these techniques, with a special emphasis given to their synergy. The present status of the research is then illustrated by selected examples which include classification of DNA and RNA backbone families, benchmark structure-energy quantum chemical calculations, parameterization of the dihedral space of simulation force fields, incorporation of arsenate into DNA, sugar-phosphate backbone self-cleavage in small RNA enzymes, and intricate geometries of the backbone in recurrent RNA building blocks. Although not apparent from the current literature showing limited overlaps between the QM, simulation and bioinformatics studies of the nucleic acids backbone, there in fact should be a major cooperative interaction between these three approaches in studies of the sugar-phosphate backbone.

Synthesis, properties, and applications of oligonucleotides containing an RNA dinucleotide phosphorothiolate linkage

RNA represents a prominent class of biomolecules. Present in all living systems, RNA plays many essential roles in gene expression, regulation, and development. Accordingly, many biological processes depend on the accurate enzymatic processing, modification, and cleavage of RNA. Understanding the catalytic mechanisms of these enzymes therefore represents an important goal in defining living systems at the molecular level. In this context, RNA molecules bearing 3'- or 5'-S-phosphorothiolate linkages comprise what are arguably among the most incisive mechanistic probes available. They have been instrumental in showing that RNA splicing systems are metalloenzymes and in mapping the ligands that reside within RNA active sites. The resulting models have in turn verified the functional relevance of crystal structures. In other cases, phosphorothiolates have offered an experimental strategy to circumvent the classic problem of kinetic ambiguity; mechanistic enzymologists have used this tool to assign precise roles to catalytic groups as general acids or bases. These insights into macromolecular function are enabled by the synthesis of nucleic acids bearing phosphorothiolate linkages and the unique chemical properties they impart. In this Account, we review the synthesis, properties, and applications of oligonucleotides and oligodeoxynucleotides containing an RNA dinucleotide phosphorothiolate linkage. Phosphorothioate linkages are structurally very similar to phosphorothiolate linkages, as reflected in the single letter of difference in nomenclature. Phosphorothioate substitutions, in which sulfur replaces one or both nonbridging oxygens within a phosphodiester linkage, are now widely available and are used routinely in numerous biochemical and medicinal applications. Indeed, synthetic phosphorothioate linkages can be introduced readily via a sulfurization step programmed into automated solid-phase oligonucleotide synthesizers. In contrast, phosphorothiolate oligonucleotides, in which sulfur replaces a specific 3'- or 5'-bridging oxygen, have presented a more difficult synthetic challenge, requiring chemical alterations to the attached sugar moiety. Here we begin by outlining the synthetic strategies used to access these phosphorothiolate RNA analogues. The Arbuzov reaction and phosphoramidite chemistry are often brought to bear in creating either 3'- or 5'-S-phosphorothiolate dinucleotides. We then summarize the responses of the phosphorothiolate derivatives to chemical and enzymatic cleavage agents, as well as mechanistic insights their use has engendered. They demonstrate particular utility as probes of metal-ion-dependent phosphotransesterification, general acid-base-catalyzed phosphotransesterification, and rate-limiting chemistry. The 3'- and 5'-S-phosphorothiolates have proven invaluable in elucidating the mechanisms of enzymatic and nonenzymatic phosphoryl transfer reactions. Considering that RNA cleavage represents a fundamental step in the maturation, degradation, and regulation of this important macromolecule, the significant synthetic challenges that remain offer rich research opportunities.

Two distinct catalytic strategies in the hepatitis δ virus ribozyme cleavage reaction

The hepatitis delta virus (HDV) ribozyme and related RNAs are widely dispersed in nature. This RNA is a small nucleolytic ribozyme that self-cleaves to generate products with a 2',3'-cyclic phosphate and a free 5'-hydroxyl. Although small ribozymes are dependent on divalent metal ions under biologically relevant buffer conditions, they function in the absence of divalent metal ions at high ionic strengths. This characteristic suggests that a functional group within the covalent structure of small ribozymes is facilitating catalysis. Structural and mechanistic analyses have demonstrated that the HDV ribozyme active site contains a cytosine with a perturbed pK(a) that serves as a general acid to protonate the leaving group. The reaction of the HDV ribozyme in monovalent cations alone never approaches the velocity of the Mg(2+)-dependent reaction, and there is significant biochemical evidence that a Mg(2+) ion participates directly in catalysis. A recent crystal structure of the HDV ribozyme revealed that there is a metal binding pocket in the HDV ribozyme active site. Modeling of the cleavage site into the structure suggested that this metal ion can interact directly with the scissile phosphate and the nucleophile. In this manner, the Mg(2+) ion can serve as a Lewis acid, facilitating deprotonation of the nucleophile and stabilizing the conformation of the cleavage site for in-line attack of the nucleophile at the scissile phosphate. This catalytic strategy had previously been observed only in much larger ribozymes. Thus, in contrast to most large and small ribozymes, the HDV ribozyme uses two distinct catalytic strategies in its cleavage reaction.

The glmS ribozyme tunes the catalytically critical pK(a) of its coenzyme glucosamine-6-phosphate

The glmS ribozyme riboswitch is the first known natural catalytic RNA that employs a small-molecule cofactor. Binding of glucosamine-6-phosphate (GlcN6P) uncovers the latent self-cleavage activity of the RNA, which adopts a catalytically competent conformation that is nonetheless inactive in the absence of GlcN6P. Structural and analogue studies suggest that the amine of GlcN6P functions as a general acid-base catalyst, while its phosphate is important for binding affinity. However, the solution pK(a) of the amine is 8.06 ± 0.05, which is not optimal for proton transfer. Here we used Raman crystallography directly to determine the pK(a)'s of GlcN6P bound to the glmS ribozyme. Binding to the RNA lowers the pK(a) of the amine of GlcN6P to 7.26 ± 0.09 and raises the pK(a) of its phosphate to 6.35 ± 0.09. Remarkably, the pK(a)'s of these two functional groups are unchanged from their values for free GlcN6P (8.06 ± 0.05 and 5.98 ± 0.05, respectively) when GlcN6P binds to the catalytically inactive but structurally unperturbed G40A mutant of the ribozyme, thus implicating the ribozyme active site guanine in pK(a) tuning. This is the first demonstration that a ribozyme can tune the pK(a) of a small-molecule ligand. Moreover, the anionic glmS ribozyme in effect stabilizes the neutral amine of GlcN6P by lowering its pK(a). This is unprecedented and illustrates the chemical sophistication of ribozyme active sites.

Base ionization and ligand binding: how small ribozymes and riboswitches gain a foothold in a protein world

Genome sequencing has produced thousands of nonprotein coding (nc)RNA sequences including new ribozymes and riboswitches. Such RNAs are notable for their extraordinary functionality, which entails exquisite folding that culminates in biocatalytic or ligand-binding capabilities. Here we discuss advances in relating ncRNA form to function with an emphasis on base pK(a) shifting by the hairpin and hepatitis delta virus ribozymes. We then describe ligand binding by the two smallest riboswitches, which target preQ(1) and S-adenosyl-(l)-homocysteine, followed by an analysis of a second-messenger riboswitch that binds cyclic-di-GMP. Each riboswitch is then compared to a protein that binds the same ligand to contrast binding properties. The results showcase the breadth of functionality attainable from ncRNAs, as well as molecular features notable for antibacterial design.

Metal binding motif in the active site of the HDV ribozyme binds divalent and monovalent ions

The hepatitis delta virus (HDV) ribozyme uses both metal ion and nucleobase catalysis in its cleavage mechanism. A reverse G·U wobble was observed in a recent crystal structure of the precleaved state. This unusual base pair positions a Mg(2+) ion to participate in catalysis. Herein, we used molecular dynamics (MD) and X-ray crystallography to characterize the conformation and metal binding characteristics of this base pair in product and precleaved forms. Beginning with a crystal structure of the product form, we observed formation of the reverse G·U wobble during MD trajectories. We also demonstrated that this base pair is compatible with the diffraction data for the product-bound state. During MD trajectories of the product form, Na(+) ions interacted with the reverse G·U wobble in the RNA active site, and a Mg(2+) ion, introduced in certain trajectories, remained bound at this site. Beginning with a crystal structure of the precleaved form, the reverse G·U wobble with bound Mg(2+) remained intact during MD simulations. When we removed Mg(2+) from the starting precleaved structure, Na(+) ions interacted with the reverse G·U wobble. In support of the computational results, we observed competition between Na(+) and Mg(2+) in the precleaved ribozyme crystallographically. Nonlinear Poisson-Boltzmann calculations revealed a negatively charged patch near the reverse G·U wobble. This anionic pocket likely serves to bind metal ions and to help shift the pK(a) of the catalytic nucleobase, C75. Thus, the reverse G·U wobble motif serves to organize two catalytic elements, a metal ion and catalytic nucleobase, within the active site of the HDV ribozyme.

Protonation of base pairs in RNA: context analysis and quantum chemical investigations of their geometries and stabilities

Base pairs involving protonated nucleobases play important roles in mediating global macromolecular conformational changes and in facilitation of catalysis in a variety of functional RNA molecules. Here we present our attempts at understanding the role of such base pairs by detecting possible protonated base pairs in the available RNA crystal structures using BPFind software, in their specific structural contexts, and by the characterization of their geometries, interaction energies, and stabilities using advanced quantum chemical computations. We report occurrences of 18 distinct protonated base pair combinations from a representative data set of RNA crystal structures and propose a theoretical model for one putative base pair combination. Optimization of base pair geometries was carried out at the B3LYP/cc-pVTZ level, and the BSSE corrected interaction energies were calculated at the MP2/aug-cc-pVDZ level of theory. The geometries for each of the base pairs were characterized in terms of H-bonding patterns observed, rmsd values observed on optimization, and base pair geometrical parameters. In addition, the intermolecular interaction in these complexes was also analyzed using Morokuma energy decomposition. The gas phase interaction energies of the base pairs range from -24 to -49 kcal/mol and reveal the dominance of Hartree-Fock component of interaction energy constituting 73% to 98% of the total interaction energy values. On the basis of our combined bioinformatics and quantum chemical analysis of different protonated base pairs, we suggest resolution of structural ambiguities and correlate their geometric and energetic features with their structural and functional roles. In addition, we also examine the suitability of specific base pairs as key elements in molecular switches and as nucleators for higher order structures such as base triplets and quartets.

DNA as a versatile chemical component for catalysis, encoding, and stereocontrol

DNA (deoxyribonucleic acid) is the genetic material common to all of Earth's organisms. Our biological understanding of DNA is extensive and well-exploited. In recent years, chemists have begun to develop DNA for nonbiological applications in catalysis, encoding, and stereochemical control. This Review summarizes key advances in these three exciting research areas, each of which takes advantage of a different subset of DNA's useful chemical properties.

Metal ions: supporting actors in the playbook of small ribozymes

Since the 1980s, several small RNA motifs capable of chemical catalysis have been discovered. These small ribozymes, composed of between approximately 40 and 200 nucleotides, have been found to play vital roles in the replication of subviral and viral pathogens, as well as in gene regulation in prokaryotes, and have recently been discovered in noncoding eukaryotic RNAs. All of the known natural small ribozymes - the hairpin, hammerhead, hepatitis delta virus, Varkud satellite, and glmS ribozymes--catalyze the same self-cleavage reaction as RNase A, resulting in two products, one bearing a 2'-3' cyclic phosphate and the other a 5'-hydroxyl group. Although originally thought to be obligate metalloenzymes like the group I and II self-splicing introns, the small ribozymes are now known to support catalysis in a wide variety of cations that appear to be only indirectly involved in catalysis. Nevertheless, under physiologic conditions, metal ions are essential for the proper folding and function of the small ribozymes, the most effective of these being magnesium. Metal ions contribute to catalysis in the small ribozymes primarily by stabilizing the catalytically active conformation, but in some cases also by activating RNA functional groups for catalysis, directly participating in catalytic acid-base chemistry, and perhaps by neutralizing the developing negative charge of the transition state. Although interactions between the small ribozymes and cations are relatively nonspecific, ribozyme activity is quite sensitive to the types and concentrations of metal ions present in solution, suggesting a close evolutionary relationship between cellular metal ion homeostasis and cation requirements of catalytic RNAs, and perhaps RNA in general.

The glmS ribozyme: use of a small molecule coenzyme by a gene-regulatory RNA

The glmS ribozyme is the first known example of a natural ribozyme that has evolved to require binding of an exogenous small molecule for activity. In Gram-positive bacteria, this RNA domain is part of the messenger RNA (mRNA) encoding the essential enzyme that synthesizes glucosamine-6-phosphate (GlcN6P). When present at physiologic concentration, this small molecule binds to the glmS ribozyme and uncovers a latent self-cleavage activity that ultimately leads to degradation of the mRNA. Biochemical and structural studies reveal that the RNA adopts a rigid fold stabilized by three pseudoknots and the packing of a peripheral domain against the ribozyme core. GlcN6P binding to this pre-organized RNA does not induce conformational changes; rather, the small molecule functions as a coenzyme, providing a catalytically essential amine group to the active site. The ribozyme is not a passive player, however. Active site functional groups are essential for catalysis, even in the presence of GlcN6P. In addition to being a superb experimental system with which to analyze how RNA catalysts can exploit small molecule coenzymes to broaden their chemical versatility, the presence of the glmS ribozyme in numerous pathogenic bacteria make this RNA an attractive target for the development of new antibiotics and antibacterial strategies.

Investigation of the catalytic mechanism of a synthetic DNAzyme with protein-like functionality: an RNaseA mimic?

The protein enzyme ribonuclease A (RNaseA) cleaves RNA with catalytic perfection, although with little sequence specificity, by a divalent metal ion (M(2+))-independent mechanism in which a pair of imidazoles provides general acid and base catalysis, while a cationic amine provides electrostatic stabilization of the transition state. Synthetic imitation of this remarkable organo-catalyst ("RNaseA mimicry") has been a longstanding goal in biomimetic chemistry. The 9(25)-11 DNAzyme contains synthetically modified nucleotides presenting both imidazole and cationic amine side chains, and catalyzes RNA cleavage with turnover in the absence of M(2+) similarly to RNaseA. Nevertheless, the catalytic roles, if any, of the "protein-like" functional groups have not been defined, and hence the question remains whether 9(25)-11 engages any of these functionalities to mimic aspects of the mechanism of RNaseA. To address this question, we report a mechanistic investigation of 9(25)-11 catalysis wherein we have employed a variety of experiments, such as DNAzyme functional group deletion, mechanism-based affinity labeling, and bridging and nonbridging phosphorothioate substitution of the scissile phosphate. Several striking parallels exist between the results presented here for 9(25)-11 and the results of analogous experiments applied previously to RNaseA. Specifically, our results implicate two particular imidazoles in general acid and base catalysis and suggest that a specific cationic amine stabilizes the transition state via diastereoselective interaction with the scissile phosphate. Overall, 9(25)-11 appears to meet the minimal criteria of an RNaseA mimic; this demonstrates how added synthetic functionality can expand the mechanistic repertoire available to a synthetic DNA-based catalyst.