A guanine nucleobase important for catalysis by the VS ribozyme

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.

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.

The Varkud Satellite Ribozyme: A Thirty-Year Journey through Biochemistry, Crystallography, and Computation

The discovery of catalytic RNAs or ribozymes introduced a new class of enzymes to biology. In addition to their increasingly important roles in modern life, ribozymes are key players in the RNA World hypothesis, which posits that life started or flourished with RNA supporting both genetic and enzymatic functions. Therefore, investigations into the mechanisms of ribozyme function provide an exciting opportunity to examine the foundational principles of biological catalysis. Ribozymes are also attractive model systems to investigate the relationship between structure and function in RNA. Endonucleolytic ribozymes represent the largest class of catalytic RNA, of which the Varkud satellite (VS) ribozyme is structurally the most complex. The last ribozyme to be discovered by accident, the VS ribozyme had eluded structural determination for over two decades. When we solved the first crystal structures of the VS ribozyme, an extensive body of biochemical and biophysical data had accumulated over the years with which we could evaluate the functional relevance of the structure. Conversely, the structures provided a new perspective from which to reexamine the functional data and test new hypotheses. The VS ribozyme is organized in a modular fashion where independently folding domains assemble into the active conformation of the ribozyme via three-way junctions. Structures of the VS ribozyme in complex with its substrate at different stages of activation enabled us to map the structural reorganization of the substrate that must precede catalysis. In addition to defining the global architecture of the RNA, the essential interactions between the substrate and catalytic domains, and the rearrangements in the substrate prior to catalysis, these structures provided detailed snapshots of the ribozyme active site, revealing potential catalytic interactions. High resolution structures of the active site bolstered the view that the catalytic mechanism involved nucleobase-mediated general acid-base catalysis and uncovered additional catalytic interactions between the cleavage site and catalytic residues. Informed by the crystal structures of the VS ribozyme, an integrated experimental and computational approach identified the key players and essential interactions that define the active site of the ribozyme. This confluence of biochemical, structural, and computational studies revealed the catalytic mechanism of the ribozyme at unprecedented detail. Additionally, comparative analyses of the active site structures of the VS ribozyme and other nucleic acid-based endoribonucleases revealed common architectural motifs and strikingly similar catalytic strategies. In this Account, we document the progress of VS ribozyme research starting from its discovery and extending to the elucidation of its detailed catalytic mechanism 30 years later.

The L-platform/L-scaffold framework: a blueprint for RNA-cleaving nucleic acid enzyme design

We develop an L-platform/L-scaffold framework we hypothesize may serve as a blueprint to facilitate site-specific RNA-cleaving nucleic acid enzyme design. Building on the L-platform motif originally described by Suslov and coworkers, we identify new critical scaffolding elements required to anchor a conserved general base guanine ("L-anchor") and bind functionally important metal ions at the active site ("L-pocket"). Molecular simulations, together with a broad range of experimental structural and functional data, connect the L-platform/L-scaffold elements to necessary and sufficient conditions for catalytic activity. We demonstrate that the L-platform/L-scaffold framework is common to five of the nine currently known naturally occurring ribozyme classes (Twr, HPr, VSr, HHr, Psr), and intriguingly from a design perspective, the framework also appears in an artificially engineered DNAzyme (8-17dz). The flexibility of the L-platform/L-scaffold framework is illustrated on these systems, highlighting modularity and trends in the variety of known general acid moieties that are supported. These trends give rise to two distinct catalytic paradigms, building on the classifications proposed by Wilson and coworkers and named for the implicated general base and acid. The "G + A" paradigm (Twr, HPr, VSr) exclusively utilizes nucleobase residues for chemistry, and the "G + M + " paradigm (HHr, 8-17dz, Psr) involves structuring of the "L-pocket" metal ion binding site for recruitment of a divalent metal ion that plays an active role in the chemical steps of the reaction. Finally, the modularity of the L-platform/L-scaffold framework is illustrated in the VS ribozyme where the "L-pocket" assumes the functional role of the "L-anchor" element, highlighting a distinct mechanism, but one that is functionally linked with the hammerhead ribozyme.

Confluence of theory and experiment reveals the catalytic mechanism of the Varkud satellite ribozyme

The Varkud satellite ribozyme catalyses site-specific RNA cleavage and ligation, and serves as an important model system to understand RNA catalysis. Here, we combine stereospecific phosphorothioate substitution, precision nucleobase mutation and linear free-energy relationship measurements with molecular dynamics, molecular solvation theory and ab initio quantum mechanical/molecular mechanical free-energy simulations to gain insight into the catalysis. Through this confluence of theory and experiment, we unify the existing body of structural and functional data to unveil the catalytic mechanism in unprecedented detail, including the degree of proton transfer in the transition state. Further, we provide evidence for a critical Mg²⁺ in the active site that interacts with the scissile phosphate and anchors the general base guanine in position for nucleophile activation. This novel role for Mg²⁺ adds to the diversity of known catalytic RNA strategies and unifies functional features observed in the Varkud satellite, hairpin and hammerhead ribozyme classes.

Nucleic Acid Catalysis under Potential Prebiotic Conditions

Catalysis by nucleic acids is indispensable for extant cellular life, and it is widely accepted that nucleic acid enzymes were crucial for the emergence of primitive life 3.5-4 billion years ago. However, geochemical conditions on early Earth must have differed greatly from the constant internal milieus of today's cells. In order to explore plausible scenarios for early molecular evolution, it is therefore essential to understand how different physicochemical parameters, such as temperature, pH, and ionic composition, influence nucleic acid catalysis and to explore to what extent nucleic acid enzymes can adapt to non-physiological conditions. In this article, we give an overview of the research on catalysis of nucleic acids, in particular catalytic RNAs (ribozymes) and DNAs (deoxyribozymes), under extreme and/or unusual conditions that may relate to prebiotic environments.

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.

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.

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.

Design and synthesis of a novel lanthanide fluorescent probe (TbIII-dtpa-bis(2,6-diaminopurine)) and its application to the detection of uric acid in urine sample

In this study, a novel fluorescent probe, Tb^III-dtpa-bis(2,6-diaminopurine) (Tb-dtpa-bdap), is designed based on the principle of complementary base pairing and synthesized for uric acid detection. The synthesized fluorescent probe is characterized by ¹H NMR, ¹³C NMR, infra-red (IR) spectrum and ultraviolet-visible (UV-vis) spectra. It is found that the fluorescence of Tb-dtpa-bdap solution can be quenched obviously in the presence of uric acid. The affecting factors, including solution acidity, uric acid concentration and interfering substances, on the detection of uric acid using this probe are examined. Under optimized conditions, the fluorescence intensities of Tb-dtpa-bdap solution towards different uric acid concentrations show a linear response in the range from 1.00 × 10^-5 mol·L^-1 to 5.00 × 10^-5 mol·L^-1 with a linear correlation coefficient (R²) of 0.9877. And the obtained limit of detection (LOD) is about 5.80 × 10^-6 mol·L^-1, which is lower than the level of uric acid in actual urine. The mechanism on the detection of uric acid by using Tb-dtpa-bdap is inferred from the experimental results. The facts demonstrate that the proposed fluorescent probe can be successfully applied for the determination of uric acid in human urine samples.

Complexity in pH-Dependent Ribozyme Kinetics: Dark pKa Shifts and Wavy Rate-pH Profiles

Charged bases occur in RNA enzymes, or ribozymes, where they play key roles in catalysis. Cationic bases donate protons and perform electrostatic catalysis, while anionic bases accept protons. We previously published simulations of rate-pH profiles for ribozymes in terms of species plots for the general acid and general base that have been useful for understanding how ribozymes respond to pH. In that study, we did not consider interaction between the general acid and general base or interaction with other species on the RNA. Since that report, diverse small ribozyme classes have been discovered, many of which have charged nucleobases or metal ions in the active site that can either directly interact and participate in catalysis or indirectly interact as "influencers". Herein, we simulate experimental rate-pH profiles in terms of species plots in which reverse protonated charged nucleobases interact. These analyses uncover two surprising features of pH-dependent enzyme kinetics. (1) Cooperativity between the general acid and general base enhances population of the functional forms of a ribozyme and manifests itself as hidden or "dark" pK_a shifts, real pK_a shifts that accelerate the reaction but are not readily observed by standard experimental approaches, and (2) influencers favorably shift the pK_as of proton-transferring nucleobases and manifest themselves as "wavy" rate-pH profiles. We identify parallels with the protein enzyme literature, including reverse protonation and wavelike behavior, while pointing out that RNA is more prone to reverse protonation. The complexities uncovered, which arise from simple pairwise interactions, should aid deconvolution of complex rate-pH profiles for RNA and protein enzymes and suggest veiled catalytic devices for promoting catalysis that can be tested by experiment and calculation.

Structural Basis for Substrate Helix Remodeling and Cleavage Loop Activation in the Varkud Satellite Ribozyme

The Varkud satellite (VS) ribozyme catalyzes site-specific RNA cleavage and ligation reactions. Recognition of the substrate involves a kissing loop interaction between the substrate and the catalytic domain of the ribozyme, resulting in a rearrangement of the substrate helix register into a so-called "shifted" conformation that is critical for substrate binding and activation. We report a 3.3 Å crystal structure of the complete ribozyme that reveals the active, shifted conformation of the substrate, docked into the catalytic domain of the ribozyme. Comparison to previous NMR structures of isolated, inactive substrates provides a physical description of substrate remodeling, and implicates roles for tertiary interactions in catalytic activation of the cleavage loop. Similarities to the hairpin ribozyme cleavage loop activation suggest general strategies to enhance fidelity in RNA folding and ribozyme cleavage.

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.

Cooperative Interactions in the Hammerhead Ribozyme Drive pKa Shifting of G12 and Its Stacked Base C17

General acid-base catalysis is a key mechanistic strategy in protein and RNA enzymes. Ribozymes use hydrated metal ions, nucleobases, and organic cofactors to carry this out. In most small ribozymes, a guanosine is positioned to participate in proton transfer with the nucleophilic 2'-OH. The unshifted pK_a values for nucleobases and solvated metal ions are far from neutrality, however, and thus nonideal for general acid-base catalysis. Herein, evidence is provided for cooperative interaction in the hammerhead ribozyme among the guanine that interacts with the nucleophilic 2'-OH, G12, the -1 nucleobase C17, and Mg²⁺ ions. We introduce global fitting for analyzing ribozyme rate-pH data parametric in Mg²⁺ concentration and benchmark this method on data from the hepatitis delta virus ribozyme. We then apply global fitting to new rate-pH data for the hammerhead ribozyme using a minimal three-dimensional, four-channel cooperative model. The value for the pK_a of G12 that we obtain is channel-dependent and varies from 8.1 to 9.9, shifting closest toward neutrality in the presence of two cationic species: C17H⁺ and a Mg²⁺ ion. The value for the pK_a of the -1 nucleotide, C17, is increased a remarkable 3.5-5 pK_a units toward neutrality. Shifting of the pK_a of C17 appears to be driven by an electrostatic sandwich of C17 between carbonyl groups of the 5'-neighboring U and of G12 and involves cation-π interactions. Rate-pH profiles reveal that the major reactive channel under biological Mg²⁺ and pH involves a cationic C17 rather than a second metal ion. Substitution of a cationic base for a metal underscores the versatility of RNA.

The structure of a nucleolytic ribozyme that employs a catalytic metal ion

The TS ribozyme (originally called "twister sister") is a catalytic RNA. We present a crystal structure of the ribozyme in a pre-reactive conformation. Two co-axial helical stacks are organized by a three-way junction and two tertiary contacts. Five divalent metal ions are directly coordinated to RNA ligands, making important contributions to the RNA architecture. The scissile phosphate lies in a quasihelical loop region that is organized by a network of hydrogen bonding. A divalent metal ion is directly bound to the nucleobase 5' to the scissile phosphate, with an inner-sphere water molecule positioned to interact with the O2' nucleophile. The rate of ribozyme cleavage correlated in a log-linear manner with divalent metal ion pKa, consistent with proton transfer in the transition state, and we propose that the bound metal ion is a likely general base for the cleavage reaction. Our data indicate that the TS ribozyme functions predominantly as a metalloenzyme.

Structural and Biochemical Properties of Novel Self-Cleaving Ribozymes

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

Insights into RNA structure and dynamics from recent NMR and X-ray studies of the Neurospora Varkud satellite ribozyme

Despite the large number of noncoding RNAs and their importance in several biological processes, our understanding of RNA structure and dynamics at atomic resolution is still very limited. Like many other RNAs, the Neurospora Varkud satellite (VS) ribozyme performs its functions through dynamic exchange of multiple conformational states. More specifically, the VS ribozyme recognizes and cleaves its stem-loop substrate via a mechanism that involves several structural transitions within its stem-loop substrate. The recent publications of high-resolution structures of the VS ribozyme, obtained by NMR spectroscopy and X-ray crystallography, offer an opportunity to integrate the data and closely examine the structural and dynamic properties of this model RNA system. Notably, these investigations provide a valuable example of the divide-and-conquer strategy for structural and dynamic characterization of a large RNA, based on NMR structures of several individual subdomains. The success of this divide-and-conquer approach reflects the modularity of RNA architecture and the great care taken in identifying the independently-folding modules. Together with previous biochemical and biophysical characterizations, the recent NMR and X-ray studies provide a coherent picture into how the VS ribozyme recognizes its stem-loop substrate. Such in-depth characterization of this RNA enzyme will serve as a model for future structural and engineering studies of dynamic RNAs and may be particularly useful in planning divide-and-conquer investigations. WIREs RNA 2017, 8:e1421. doi: 10.1002/wrna.1421 For further resources related to this article, please visit the WIREs website.

Molecular Dynamics Study of Twister Ribozyme: Role of Mg(2+) Ions and the Hydrogen-Bonding Network in the Active Site

The recently discovered twister ribozyme is thought to utilize general acid-base catalysis in its self-cleavage mechanism, but the roles of nucleobases and metal ions in the mechanism are unclear. Herein, molecular dynamics simulations of the env22 twister ribozyme are performed to elucidate the structural and equilibrium dynamical properties, as well as to examine the role of Mg(2+) ions and possible candidates for the general base and acid in the self-cleavage mechanism. The active site region and the ends of the pseudoknots were found to be less mobile than other regions of the ribozyme, most likely providing structural stability and possibly facilitating catalysis. A purported catalytic Mg(2+) ion and the closest neighboring Mg(2+) ion remained chelated and relatively immobile throughout the microsecond trajectories, although removal of these Mg(2+) ions did not lead to any significant changes in the structure or equilibrium motions of the ribozyme on the microsecond time scale. In addition, a third metal ion, a Na(+) ion remained close to A1(O5'), the leaving group atom, during the majority of the microsecond trajectories, suggesting that it might stabilize the negative charge on A1(O5') during self-cleavage. The locations of these cations and their interactions with key nucleotides in the active site suggest that they may be catalytically relevant. The P1 stem is partially melted at its top and bottom in the crystal structure and further unwinds in the trajectories. The simulations also revealed an interconnected network comprised of hydrogen-bonding and π-stacking interactions that create a relatively rigid network around the self-cleavage site. The nucleotides involved in this network are among the highly conserved nucleotides in twister ribozymes, suggesting that this interaction network may be important to structure and function.

The Novel Chemical Mechanism of the Twister Ribozyme

We describe the multifactorial origins of catalysis by the twister ribozyme. We provide evidence that the adenine immediately 3' to the scissile phosphate (A1) acts as a general acid. Substitution of ring nitrogen atoms indicates that very unusually the N3 of A1 is the proton donor to the oxyanion leaving group. A1 is accommodated in a specific binding pocket that raises its pKa toward neutrality, juxtaposes its N3 with the O5' to be protonated, and helps create the in-line trajectory required for nucleophilic attack. A1 performs general acid catalysis while G33 acts as a general base. A 100-fold stereospecific phosphorothioate effect at the scissile phosphate is consistent with a significant stabilization of the transition state by the ribozyme, and functional group substitution at G33 indicates that its exocyclic N2 interacts directly with the scissile phosphate. A model of the ribozyme active site is proposed that accommodates these catalytic strategies.

Rational engineering of the Neurospora VS ribozyme to allow substrate recognition via different kissing-loop interactions

The Neurospora VS ribozyme is a catalytic RNA that has the unique ability to specifically recognize and cleave a stem-loop substrate through formation of a highly stable kissing-loop interaction (KLI). In order to explore the engineering potential of the VS ribozyme to cleave alternate substrates, we substituted the wild-type KLI by other known KLIs using an innovative engineering method that combines rational and combinatorial approaches. A bioinformatic search of the protein data bank was initially performed to identify KLIs that are structurally similar to the one found in the VS ribozyme. Next, substrate/ribozyme (S/R) pairs that incorporate these alternative KLIs were kinetically and structurally characterized. Interestingly, several of the resulting S/R pairs allowed substrate cleavage with substantial catalytic efficiency, although with reduced activity compared to the reference S/R pair. Overall, this study describes an innovative approach for RNA engineering and establishes that the KLI of the trans VS ribozyme can be adapted to cleave other folded RNA substrates.

In Vitro Selection of a DNAzyme Cooperatively Binding Two Lanthanide Ions for RNA Cleavage

Trivalent lanthanide ions (Ln(3+)) were recently employed to select RNA-cleaving DNAzymes, and three new DNAzymes have been reported so far. In this work, dysprosium (Dy(3+)) was used with a library containing 50 random nucleotides. After six rounds of in vitro selection, a new DNAzyme named Dy10a was obtained and characterized. Dy10a has a bulged hairpin structure cleaving a RNA/DNA chimeric substrate. Dy10a is highly active in the presence of the five Ln(3+) ions in the middle of the lanthanide series (Sm(3+), Eu(3+), Gd(3+), Tb(3+), and Dy(3+)), while its activity descends on the two sides. The cleavage rate reaches 0.6 min(-1) at pH 6 with just 200 nM Sm(3+), which is the fastest among all known Ln(3+)-dependent enzymes. Dy10a binds two Ln(3+) ions cooperatively. When a phosphorothioate (PS) modification is introduced at the cleavage junction, the activity decreases by >2500-fold for both the Rp and Sp diastereomers, and thiophilic Cd(2+) cannot rescue the activity. The pH-rate profile has a slope of 0.37 between pH 4.2 and 5.2, and the slope was even lower at higher pH. On the basis of these data, a model of metal binding is proposed. Finally, a catalytic beacon sensor that can detect Ho(3+) down to 1.7 nM is constructed.

Chemistry and Biology of Self-Cleaving Ribozymes

Self-cleaving ribozymes were discovered 30 years ago, but their biological distribution and catalytic mechanisms are only beginning to be defined. Each ribozyme family is defined by a distinct structure, with unique active sites accelerating the same transesterification reaction across the families. Biochemical studies show that general acid-base catalysis is the most common mechanism of self-cleavage, but metal ions and metabolites can be used as cofactors. Ribozymes have been discovered in highly diverse genomic contexts throughout nature, from viroids to vertebrates. Their biological roles include self-scission during rolling-circle replication of RNA genomes, co-transcriptional processing of retrotransposons, and metabolite-dependent gene expression regulation in bacteria. Other examples, including highly conserved mammalian ribozymes, suggest that many new biological roles are yet to be discovered.

Crystal structure of the Varkud satellite ribozyme

Varkud Satellite (VS) ribozyme mediates rolling circle replication of a plasmid found in the Neurospora mitochondria. We report crystal structures of this ribozyme at 3.1Å resolution, revealing an intertwined dimer formed by an exchange of substrate helices. Within each protomer, an arrangement of three-way helical junctions organizes seven helices into a global fold that creates a docking site for the substrate helix of the other protomer, resulting in the formation of two active sites in trans. This mode of RNA-RNA association resembles the process of domain swapping in proteins and has implications for RNA regulation and evolution. Within each active site, adenine and guanine nucleobases abut the scissile phosphate, poised to serve direct roles in catalysis. Similarities to the active sites of the hairpin and hammerhead ribozymes highlight the functional significance of active site features, underscore the ability of RNA to access functional architectures from distant regions of sequence space, and suggest convergent evolution.

The NMR structure of the II-III-VI three-way junction from the Neurospora VS ribozyme reveals a critical tertiary interaction and provides new insights into the global ribozyme structure

As part of an effort to structurally characterize the complete Neurospora VS ribozyme, NMR solution structures of several subdomains have been previously determined, including the internal loops of domains I and VI, the I/V kissing-loop interaction and the III-IV-V junction. Here, we expand this work by determining the NMR structure of a 62-nucleotide RNA (J236) that encompasses the VS ribozyme II-III-VI three-way junction and its adjoining stems. In addition, we localize Mg(2+)-binding sites within this structure using Mn(2+)-induced paramagnetic relaxation enhancement. The NMR structure of the J236 RNA displays a family C topology with a compact core stabilized by continuous stacking of stems II and III, a cis WC/WC G•A base pair, two base triples and two Mg(2+) ions. Moreover, it reveals a remote tertiary interaction between the adenine bulges of stems II and VI. Additional NMR studies demonstrate that both this bulge-bulge interaction and Mg(2+) ions are critical for the stable folding of the II-III-VI junction. The NMR structure of the J236 RNA is consistent with biochemical studies on the complete VS ribozyme, but not with biophysical studies performed with a minimal II-III-VI junction that does not contain the II-VI bulge-bulge interaction. Together with previous NMR studies, our findings provide important new insights into the three-dimensional architecture of this unique ribozyme.

Fitness Landscapes of Functional RNAs

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

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.

Nuclear magnetic resonance structure of the III-IV-V three-way junction from the Varkud satellite ribozyme and identification of magnesium-binding sites using paramagnetic relaxation enhancement

The VS ribozyme is a catalytic RNA found within some natural isolates of Neurospora that is being used as a model system to improve our understanding of RNA structure, catalysis, and engineering. The catalytic domain contains five helical domains (SLII-SLVI) that are organized by two three-way junctions. The III-IV-V junction is required for high-affinity binding of the substrate domain (SLI) through formation of a kissing loop interaction with SLV. Here, we determine the high-resolution nuclear magnetic resonance (NMR) structure of a 47-nucleotide RNA containing the III-IV-V junction (J345). The J345 RNA adopts a Y-shaped fold typical of the family C three-way junctions, with coaxial stacking between stems III and IV and an acute angle between stems III and V. The NMR structure reveals that the core of the III-IV-V junction contains four stacked base triples, a U-turn motif, a cross-strand stacking interaction, an A-minor interaction, and a ribose zipper. In addition, the NMR structure shows that the cCUUGg tetraloop used to stabilize stem IV adopts a novel RNA tetraloop fold, different from the known gCUUGc tetraloop structure. Using Mn(2+)-induced paramagnetic relaxation enhancement, we identify six Mg(2+)-binding sites within J345, including one associated with the cCUUGg tetraloop and two with the junction core. The NMR structure of J345 likely represents the conformation of the III-IV-V junction in the context of the active VS ribozyme and suggests that this junction functions as a dynamic hinge that contributes to substrate recognition and catalysis. Moreover, this study highlights a new role for family C three-way junctions in long-range tertiary interactions.

A remarkably stable kissing-loop interaction defines substrate recognition by the Neurospora Varkud Satellite ribozyme

Kissing loops are tertiary structure elements that often play key roles in functional RNAs. In the Neurospora VS ribozyme, a kissing-loop interaction between the stem-loop I (SLI) substrate and stem-loop V (SLV) of the catalytic domain is known to play an important role in substrate recognition. In addition, this I/V kissing-loop interaction is associated with a helix shift in SLI that activates the substrate for catalysis. To better understand the role of this kissing-loop interaction in substrate recognition and activation by the VS ribozyme, we performed a thermodynamic characterization by isothermal titration calorimetry using isolated SLI and SLV stem-loops. We demonstrate that preshifted SLI variants have higher affinity for SLV than shiftable SLI variants, with an energetic cost of 1.8-3 kcal/mol for the helix shift in SLI. The affinity of the preshifted SLI for SLV is remarkably high, the interaction being more stable by 7-8 kcal/mol than predicted for a comparable duplex containing three Watson-Crick base pairs. The structural basis of this remarkable stability is discussed in light of previous NMR studies. Comparative thermodynamic studies reveal that kissing-loop complexes containing 6-7 Watson-Crick base pairs are as stable as predicted from comparable RNA duplexes; however, those with 2-3 Watson-Crick base pairs are more stable than predicted. Interestingly, the stability of SLI/ribozyme complexes is similar to that of SLI/SLV complexes. Thus, the I/V kissing loop interaction represents the predominant energetic contribution to substrate recognition by the trans-cleaving VS ribozyme.

NMR localization of divalent cations at the active site of the Neurospora VS ribozyme provides insights into RNA-metal-ion interactions

Metal cations represent key elements of RNA structure and function. In the Neurospora VS ribozyme, metal cations play diverse roles; they are important for substrate recognition, formation of the active site, and shifting the pKa's of two key nucleobases that contribute to the general acid-base mechanism. Recently, we determined the NMR structure of the A730 loop of the VS ribozyme active site (SLVI) that contributes the general acid (A756) in the enzymatic mechanism of the cleavage reaction. Our studies showed that magnesium (Mg(2+)) ions are essential to stabilize the formation of the S-turn motif within the A730 loop that exposes the A756 nucleobase for catalysis. In this article, we extend these NMR investigations by precisely mapping the Mg(2+)-ion binding sites using manganese-induced paramagnetic relaxation enhancement and cadmium-induced chemical-shift perturbation of phosphorothioate RNAs. These experiments identify five Mg(2+)-ion binding sites within SLVI. Four Mg(2+) ions in SLVI are associated with known RNA structural motifs, including the G-U wobble pair and the GNRA tetraloop, and our studies reveal novel insights about Mg(2+) ion binding to these RNA motifs. Interestingly, one Mg(2+) ion is specifically associated with the S-turn motif, confirming its structural role in the folding of the A730 loop. This Mg(2+) ion is likely important for formation of the active site and may play an indirect role in catalysis.

General acid-base catalysis mediated by nucleobases in the hairpin ribozyme

The catalytic mechanism by which the hairpin ribozyme accelerates cleavage or ligation of the phosphodiester backbone of RNA has been incompletely understood. There is experimental evidence for an important role for an adenine (A38) and a guanine (G8), and it has been proposed that these act in general acid-base catalysis. In this work we show that a large reduction in cleavage rate on substitution of A38 by purine (A38P) can be reversed by replacement of the 5'-oxygen atom at the scissile phosphate by sulfur (5'-PS), which is a much better leaving group. This is consistent with A38 acting as the general acid in the unmodified ribozyme. The rate of cleavage of the 5'-PS substrate by the A38P ribozyme increases with pH log-linearly, indicative of a requirement for a deprotonated base with a relatively high pK(a). On substitution of G8 by diaminopurine, the 5'-PS substrate cleavage rate at first increases with pH and then remains at a plateau, exhibiting an apparent pK(a) consistent with this nucleotide acting in general base catalysis. Alternative explanations for the pH dependence of hairpin ribozyme reactivity are discussed, from which we conclude that general acid-base catalysis by A38 and G8 is the simplest and most probable explanation consistent with all the experimental data.

The chemical versatility of RNA

The ability of RNA to both store genetic information and catalyse chemical reactions has led to the hypothesis that it predates DNA and proteins. While there is no doubt that RNA is capable of storing the genetic information of a primitive organism, only two classes of reactions-phosphoryl transfer and peptide bond formation-have been observed to be catalysed by RNA in nature. However, these naturally occurring ribozymes use a wide range of catalytic strategies that could be applied to other reactions. Furthermore, RNA can bind several cofactors that are used by protein enzymes to facilitate a wide variety of chemical processes. Despite its limited functional groups, these observations indicate RNA is a versatile molecule that could, in principle, catalyse the myriad reactions necessary to sustain life.

Mechanisms of RNA catalysis

Ribozymes are RNA molecules that act as chemical catalysts. In contemporary cells, most known ribozymes carry out phosphoryl transfer reactions. The nucleolytic ribozymes comprise a class of five structurally-distinct species that bring about site-specific cleavage by nucleophilic attack of the 2'-O on the adjacent 3'-P to form a cyclic 2',3'-phosphate. In general, they will also catalyse the reverse reaction. As a class, all these ribozymes appear to use general acid-base catalysis to accelerate these reactions by about a million-fold. In the Varkud satellite ribozyme, we have shown that the cleavage reaction is catalysed by guanine and adenine nucleobases acting as general base and acid, respectively. The hairpin ribozyme most probably uses a closely similar mechanism. Guanine nucleobases appear to be a common choice of general base, but the general acid is more variable. By contrast, the larger ribozymes such as the self-splicing introns and RNase P act as metalloenzymes.

Crystal structure of an RNA polymerase ribozyme in complex with an antibody fragment

All models of the RNA world era invoke the presence of ribozymes that can catalyse RNA polymerization. The class I ligase ribozyme selected in vitro 15 years ago from a pool of random RNA sequences catalyses formation of a 3',5'-phosphodiester linkage analogous to a single step of RNA polymerization. Recently, the three-dimensional structure of the ligase was solved in complex with U1A RNA-binding protein and independently in complex with an antibody fragment. The RNA adopts a tripod arrangement and appears to use a two-metal ion mechanism similar to protein polymerases. Here, we discuss structural implications for engineering a true polymerase ribozyme and describe the use of the antibody framework both as a portable chaperone for crystallization of other RNAs and as a platform for exploring steps in evolution from the RNA world to the RNA-protein world.

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.

Helix-length compensation studies reveal the adaptability of the VS ribozyme architecture

Compensatory mutations in RNA are generally regarded as those that maintain base pairing, and their identification forms the basis of phylogenetic predictions of RNA secondary structure. However, other types of compensatory mutations can provide higher-order structural and evolutionary information. Here, we present a helix-length compensation study for investigating structure–function relationships in RNA. The approach is demonstrated for stem-loop I and stem-loop V of the Neurospora VS ribozyme, which form a kissing–loop interaction important for substrate recognition. To rapidly characterize the substrate specificity (k_cat/K_M) of several substrate/ribozyme pairs, a procedure was established for simultaneous kinetic characterization of multiple substrates. Several active substrate/ribozyme pairs were identified, indicating the presence of limited substrate promiscuity for stem Ib variants and helix-length compensation between stems Ib and V. 3D models of the I/V interaction were generated that are compatible with the kinetic data. These models further illustrate the adaptability of the VS ribozyme architecture for substrate cleavage and provide global structural information on the I/V kissing–loop interaction. By exploring higher-order compensatory mutations in RNA our approach brings a deeper understanding of the adaptability of RNA structure, while opening new avenues for RNA research.

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.

Arginine cofactors on the polymerase ribozyme

The RNA world hypothesis states that the early evolution of life went through a stage in which RNA served both as genome and as catalyst. The central catalyst in an RNA world organism would have been a ribozyme that catalyzed RNA polymerization to facilitate self-replication. An RNA polymerase ribozyme was developed previously in the lab but it is not efficient enough for self-replication. The factor that limits its polymerization efficiency is its weak sequence-independent binding of the primer/template substrate. Here we tested whether RNA polymerization could be improved by a cationic arginine cofactor, to improve the interaction with the substrate. In an RNA world, amino acid-nucleic acid conjugates could have facilitated the emergence of the translation apparatus and the transition to an RNP world. We chose the amino acid arginine for our study because this is the amino acid most adept to interact with RNA. An arginine cofactor was positioned at ten different sites on the ribozyme, using conjugates of arginine with short DNA or RNA oligonucleotides. However, polymerization efficiency was not increased in any of the ten positions. In five of the ten positions the arginine reduced or modulated polymerization efficiency, which gives insight into the substrate-binding site on the ribozyme. These results suggest that the existing polymerase ribozyme is not well suited to using an arginine cofactor.

Modeling the RNA 2'OH activation: possible roles of metal ion and nucleobase as catalysts in self-cleaving ribozymes

The RNA 2'OH activation as taking place in the first chemical step of self-cleaving ribozymes is studied theoretically by DFT and MP2 methods using a continuum solvation model (CPCM). The reaction of proton transfer is studied in the presence of two kinds of catalysts: a fully hydrated metal ion (Mg(2+)) or partially hydrated nucleobase (guanine), taken separately or together leading to three different modes of activation. The metal ion is either directly bound (inner-sphere) or indirectly bound (outer-sphere) to the 2'OH group and a hydroxide ion acts as a general or specific base; the nucleobase is taken in anionic or in neutral enol-tautomeric forms playing itself the role of general base. The presence of a close metal ion (outer-sphere) lowers the pK(a) value of the 2'OH group by several log units in both metal-ion and nuleobase catalysis. The direct metal coordination to the 2'OH group (inner-sphere) further stabilizes the developing negative charge on the nucleophile. The switching from the inner-sphere to the outer-sphere coordination appears to be driven by the energy cost for reorganizing the first coordination shell rather than by the electrostatic repulsion between the ligands. The metal-ion catalysis is more effective with a specific base in the dianionic mechanism. On the other hand, the nucleobase catalysis is more effective in the monoanionic mechanism and in the presence of a metal ion acting as a cofactor through nonspecific electrostatic interactions. The results establish a baseline to study the possible roles of metal and nucleobase catalysts and their environment in more realistic models for self-cleaving ribozymes.

glmS Riboswitch binding to the glucosamine-6-phosphate α-anomer shifts the pKa toward neutrality

The glmS riboswitch regulates gene expression through a self-cleavage activity. The reaction is catalyzed with the assistance of the metabolite cofactor glucosamine-6-phosphate (GlcN6P), whose amino group is proposed to serve as the general acid during the reaction. This reaction is pH-dependent with a pK(a) that is lower than the observed pK(a) for the amine of GlcN6P in solution. GlcN6P, like other pyranose sugars, undergoes spontaneous and rapid interconversion between the α and β anomers at the C1 position. Here we demonstrate by NMR that the Bacillus anthracis glmS riboswitch selectively binds the α-anomer of GlcN6P with a maximum binding affinity of 0.36 mM and that binding is pH-dependent. We also report that the anomeric ratio between α and β is pH-dependent and the pK(a)s of the two amines differ by 0.5 pH units, α being the higher of the two (pK(a)=8.3). The pH dependence of binding reveals a pK(a) of 6.7, suggesting that the glmS RNA reduces the pK(a) of the GlcN6P amine by 1.6 units in the ground state. We reevaluated previously obtained kinetic data and found the reaction pK(a) is 6.9, within error of the binding data. The data support a model where the reaction pK(a) corresponds to that of the GlcN6P amine. This observation has broader relevance for considering how the microenvironment of an RNA, despite its anionic character, can reduce the pK(a)s of functional groups for use in catalysis.