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

Catalytic Metal Ion-Substrate Coordination during Nonenzymatic RNA Primer Extension

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

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.

Dissociative Transition State in Hepatitis Delta Virus Ribozyme Catalysis

Ribonucleases and small nucleolytic ribozymes are both able to catalyze RNA strand cleavage through 2'-O-transphosphorylation, provoking the question of whether protein and RNA enzymes facilitate mechanisms that pass through the same or distinct transition states. Here, we report the primary and secondary 18O kinetic isotope effects for hepatitis delta virus ribozyme catalysis that reveal a dissociative, metaphosphate-like transition state in stark contrast to the late, associative transition states observed for reactions catalyzed by specific base, Zn2+ ions, or ribonuclease A. This new information provides evidence for a discrete ribozyme active site design that modulates the RNA cleavage pathway to pass through an altered transition state.

Phosphorothioate Substitutions in RNA Structure Studied by Molecular Dynamics Simulations, QM/MM Calculations, and NMR Experiments

Phosphorothioates (PTs) are important chemical modifications of the RNA backbone where a single nonbridging oxygen of the phosphate is replaced with a sulfur atom. PT can stabilize RNAs by protecting them from hydrolysis and is commonly used as a tool to explore their function. It is, however, unclear what basic physical effects PT has on RNA stability and electronic structure. Here, we present molecular dynamics (MD) simulations, quantum mechanical (QM) calculations, and NMR spectroscopy measurements, exploring the effects of PT modifications in the structural context of the neomycin-sensing riboswitch (NSR). The NSR is the smallest biologically functional riboswitch with a well-defined structure stabilized by a U-turn motif. Three of the signature interactions of the U-turn: an H-bond, an anion-π interaction, and a potassium binding site; are formed by RNA phosphates, making the NSR an ideal model for studying how PT affects RNA structure and dynamics. By comparing with high-level QM calculations, we reveal the distinct physical properties of the individual interactions facilitated by the PT. The sulfur substitution, besides weakening the direct H-bond interaction, reduces the directionality of H-bonding while increasing its dispersion and induction components. It also reduces the induction and increases the dispersion component of the anion-π stacking. The sulfur force-field parameters commonly employed in the literature do not reflect these distinctions, leading to the unsatisfactory description of PT in simulations of the NSR. We show that it is not possible to accurately describe the PT interactions using one universal set of van der Waals sulfur parameters and provide suggestions for improving the force-field performance.

In vitro Selection of Chemically Modified DNAzymes

DNAzymes are in vitro selected DNA oligonucleotides with catalytic activities. RNA cleavage is one of the most extensively studied DNAzyme reactions. To expand the chemical functionality of DNA, various chemical modifications have been made during and after selection. In this review, we summarize examples of RNA-cleaving DNAzymes and focus on those modifications introduced during in vitro selection. By incorporating various modified nucleotides via polymerase chain reaction (PCR) or primer extension, a few DNAzymes were obtained that can be specifically activated by metal ions such as Zn2+ and Hg2+. In addition, some modifications were introduced to mimic RNase A that can cleave RNA substrates in the absence of divalent metal ions. In addition, single modifications at the fixed regions of DNA libraries, especially at the cleavage junctions, have been tested, and examples of DNAzymes with phosphorothioate and histidine-glycine modified tertiary amine were successfully obtained specific for Cu2+, Cd2+, Zn2+, and Ni2+. Labeling fluorophore/quencher pair right next to the cleavage junction was also used to obtain signaling DNAzymes for detecting various metal ions and cells. Furthermore, we reviewed work on the cleavage of 2'-5' linked RNA and L-RNA substrates. Finally, applications of these modified DNAzymes as biosensors, RNases, and biochemical probes are briefly described with a few future research opportunities outlined at the end.

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

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

Crucial Roles of Two Hydrated Mg2+ Ions in Reaction Catalysis of the Pistol Ribozyme

Pistol ribozymes constitute a new class of small self‐cleaving RNAs. Crystal structures have been solved, providing three‐dimensional snapshots along the reaction coordinate of pistol phosphodiester cleavage, corresponding to the pre‐catalytic state, a vanadate mimic of the transition state, and the product. The results led to the proposed underlying chemical mechanism. Importantly, a hydrated Mg²⁺ ion remains innersphere‐coordinated to N7 of G33 in all three states, and is consistent with its likely role as acid in general acid base catalysis (δ and β catalysis). Strikingly, the new structures shed light on a second hydrated Mg²⁺ ion that approaches the scissile phosphate from its binding site in the pre‐cleavage state to reach out for water‐mediated hydrogen bonding in the cyclophosphate product. The major role of the second Mg²⁺ ion appears to be the stabilization of product conformation. This study delivers a mechanistic understanding of ribozyme‐catalyzed backbone cleavage.

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.

From general base to general acid catalysis in a sodium-specific DNAzyme by a guanine-to-adenine mutation

Recently, a few Na+-specific RNA-cleaving DNAzymes were reported, where nucleobases are likely to play critical roles in catalysis. The NaA43 and NaH1 DNAzymes share the same 16-nt Na+-binding motif, but differ in one or two nucleotides in a small catalytic loop. Nevertheless, they display an opposite pH-dependency, implicating distinct catalytic mechanisms. In this work, rational mutation studies locate a catalytic adenine residue, A22, in NaH1, while previous studies found a guanine (G23) to be important for the catalysis of NaA43. Mutation with pKa-perturbed analogs, such as 2-aminopurine (∼3.8), 2,6-diaminopurine (∼5.6) and hypoxanthine (∼8.7) affected the overall reaction rate. Therefore, we propose that the N1 position of G23 (pKa ∼6.6) in NaA43 functions as a general base, while that of A22 (pKa ∼6.3) in NaH1 as a general acid. Further experiments with base analogs and a phosphorothioate-modified substrate suggest that the exocyclic amine in A22 and both of the non-bridging oxygens at the scissile phosphate are important for catalysis for NaH1. This is an interesting example where single point mutations can change the mechanism of cleavage from general base to general acid, and it can also explain this Na+-dependent DNAzyme scaffold being sensitive to a broad range of metal ions and molecules.

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.

Unveiling the Sulfur-Sulfur Bridge: Accurate Structural and Energetic Characterization of a Homochalcogen Intermolecular Bond

By combining rotational spectroscopy in supersonic expansion with the capability of state-of-the-art quantum-chemical computations in accurately determining structural and energetic properties, the genuine nature of a sulfur-sulfur chalcogen bond between dimethyl sulfide and sulfur dioxide has been unveiled in a gas-jet environment free from collision, solvent and matrix perturbations. A SAPT analysis pointed out that electrostatic S⋅⋅⋅S interactions play the dominant role in determining the stability of the complex, largely overcoming dispersion and C-H⋅⋅⋅O hydrogen-bond contributions. Indeed, in agreement with the analysis of the quadrupole-coupling constants and of the methyl internal rotation barrier, the NBO and NOCV/CD approaches show a marked charge transfer between the sulfur atoms. Based on the assignment of the rotational spectra for 7 isotopologues, an accurate semi-experimental equilibrium structure for the heavy-atom backbone of the molecular complex has been determined, which is characterized by a S⋅⋅⋅S distance (2.947(3) Å) well below the sum of van der Waals radii.

Instantaneous Iodine-Assisted DNAzyme Cleavage of Phosphorothioate RNA

Metal ions play a critical role in the RNA-cleavage reaction by interacting with the scissile phosphate and stabilizing the highly negatively charged transition state. Many metal-dependent DNAzymes have been selected for RNA cleavage. Herein, we report that the Ce13d DNAzyme can use nonmetallic iodine (I₂) to cleave a phosphorothioate (PS)-modified substrate. The cleavage yield exceeded 60% for both the R_p and S_p stereoisomers in 10 s, while the yield without the enzyme strand was only ∼10%. The Ce13d cleavage with I₂ also required Na⁺, consistent with the property of Ce13d and confirming the similar role of I₂ as a metal ion. Ce13d had the highest yield among eight tested DNAzymes, with the second highest DNAzyme showing only 20% cleavage. The incomplete cleavage was due to competition from desulfurization and isomerization reactions. This DNAzyme was engineered for fluorescence-based I₂ detection. With EDTA for masking metal ions, I₂ was selectively detected down to 4.7 nM. Oxidation of I^- with Fe³⁺ produced I₂ in situ, allowing detection of Fe³⁺ down to 78 nM. By harnessing nonelectrostatic interactions, such as the I₂/sulfur interaction observed here, more nonmetal species might be discovered to assist DNAzyme-based RNA cleavage.

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

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

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

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

Multi-metal-dependent nucleic acid enzymes

Nucleic acid enzymes require metal ions for activity, and many recently discovered enzymes can use multiple metals, either binding to the scissile phosphate or also playing an allosteric role.

Structure-based insights into self-cleavage by a four-way junctional twister-sister ribozyme

Here we report on the crystal structure and cleavage assays of a four-way junctional twister-sister self-cleaving ribozyme. Notably, 11 conserved spatially separated loop nucleotides are brought into close proximity at the ribozyme core through long-range interactions mediated by hydrated Mg²⁺ cations. The C62–A63 step at the cleavage site adopts a splayed-apart orientation, with flexible C62 directed outwards, whereas A63 is directed inwards and anchored by stacking and hydrogen-bonding interactions. Structure-guided studies of key base, sugar, and phosphate mutations in the twister-sister ribozyme, suggest contributions to the cleavage chemistry from interactions between a guanine at the active site and the non-bridging oxygen of the scissile phosphate, a feature found previously also for the related twister ribozyme. Our four-way junctional pre-catalytic structure differs significantly in the alignment at the cleavage step (splayed-apart vs. base-stacked) and surrounding residues and hydrated Mg²⁺ ions relative to a reported three-way junctional pre-catalytic structure of the twister-sister ribozyme.

Twister-sister is a self-cleaving ribozyme. Here, the authors report the 2.0 Å crystal structure of the four-way junctional twister-sister ribozyme in the pre-catalytic state and discuss mechanistic implications based on their mutagenesis experiments and comparisons with other ribozyme structures.

Two Completely Different Mechanisms for Highly Specific Na+ Recognition by DNAzymes

Our view of the interaction between Na⁺ and nucleic acids was changed by a few recently discovered Na⁺ -specific RNA-cleaving DNAzymes. In addition to nonspecific electrostatic interactions, highly specific recognition is also possible. Herein, two such DNAzymes, named EtNa and Ce13d, are compared to elucidate their mechanisms of Na⁺ binding. Mutation studies indicate that they have different sequence requirements. Phosphorothioate (PS) substitution at the scissile phosphate drops the activity of EtNa 140-fold, and it cannot be rescued by thiophilic Cd²⁺ or Mn²⁺ , whereas the activity of PS-modified Ce13d can be rescued. Na⁺ -dependent activity assays indicate that two Na⁺ ions bind cooperatively in EtNa, and each Na⁺ likely interacts with a nonbridging oxygen atom in the scissile phosphate, whereas Ce13d binds only one Na⁺ ion in a well-defined Na⁺ aptamer, and this Na⁺ ion does not directly interact with the scissile phosphate. Both DNAzymes display a normal pH-rate profile, with a single deprotonation reaction required for catalysis. For EtNa, Na⁺ fails to protect the conserved nucleotides from dimethyl sulfate attack, and no specific Na⁺ binding is detected by 2-aminopurine fluorescence, both of which are different from those observed for Ce13d. This work suggests that EtNa binds Na⁺ mainly through its scissile phosphate without significant involvement of the nucleotides in the enzyme strand, whereas Ce13d has a well-defined aptamer for Na⁺ binding. Therefore, DNA has at least two distinct ways to achieve highly selective Na⁺ binding.

Folding of the silver aptamer in a DNAzyme probed by 2-aminopurine fluorescence

The RNA-cleaving Ag10c DNAzyme was recently isolated via in vitro selection and it can bind two Ag⁺ ions for activity. The Ag10c contains a well-defined Ag⁺ binding aptamer as indicated by DMS footprinting. Since aptamer binding is often accompanied with conformational changes, we herein used 2-aminopurine (2AP) to probe its folding in the presence of Ag⁺. The Ag10c was respectively labeled with 2AP at three different positions, both in the substrate strand and in the enzyme strand, one at a time. Ag⁺-induced folding was observed at the substrate cleavage junction and the A9 position of the enzyme strand, consistent with aptamer binding. The measured K_d at the A9 position was 18 μM Ag⁺ with a Hill coefficient of 2.17, similar to those obtained from the previous cleavage activity based assays. However, labeling a 2AP at the A2 position inhibited the activity and folding. Compared to other metal ions, Ag⁺ has a unique sigmoidal folding profile indicative of multiple silver binding cooperatively. This suggests that Ag⁺ can induce a local folding in the enzyme loop and this folding is important for activity. This study provides important biophysical insights into this new DNAzyme, suggesting the possibility of designing folding-based biosensors for Ag⁺.

A Silver-Specific DNAzyme with a New Silver Aptamer and Salt-Promoted Activity

Most RNA-cleaving DNAzymes require a metal ion to interact with the scissile phosphate for activity. Therefore, few unmodified DNAzymes work with thiophilic metals because of their low affinity for phosphate. Recently, an Ag⁺-specific Ag10c DNAzyme was reported via in vitro selection. Herein, Ag10c is characterized to rationalize the role of the strongly thiophilic Ag⁺. Systematic mutation studies indicate that Ag10c is a highly conserved DNAzyme and its Ag⁺ binding is unrelated to C-Ag⁺-C interaction. Its activity is enhanced by increasing Na⁺ concentrations in buffer. At the same metal concentration, activity decreases in the following order: Li⁺ > Na⁺ > K⁺. Ag10c binds one Na⁺ ion and two Ag⁺ ions for catalysis. The pH-rate profile has a slope of ∼1, indicating a single deprotonation step. Phosphorothioate substitution at the scissile phosphate suggests that Na⁺ interacts with the pro-R_p oxygen of the phosphate, and dimethyl sulfate footprinting indicates that the DNAzyme loop is a silver aptamer binding two Ag⁺ ions. Therefore, Ag⁺ exerts its function allosterically, while the scissile phosphate interacts with Na⁺, Li⁺, Na⁺, or Mg²⁺. This work suggests the possibility of isolating thiophilic metal aptamers based on DNAzyme selection, and it also demonstrates a new Ag⁺ aptamer.

Probing fast ribozyme reactions under biological conditions with rapid quench-flow kinetics

Reaction kinetics on the millisecond timescale pervade the protein and RNA fields. To study such reactions, investigators often perturb the system with abiological solution conditions or substrates in order to slow the rate to timescales accessible by hand mixing; however, such perturbations can change the rate-limiting step and obscure key folding and chemical steps that are found under biological conditions. Mechanical methods for collecting data on the millisecond timescale, which allow these perturbations to be avoided, have been developed over the last few decades. These methods are relatively simple and can be conducted on affordable and commercially available instruments. Here, we focus on using the rapid quench-flow technique to study the fast reaction kinetics of RNA enzymes, or ribozymes, which often react on the millisecond timescale under biological conditions. Rapid quench of ribozymes is completely parallel to the familiar hand-mixing approach, including the use of radiolabeled RNAs and fractionation of reactions on polyacrylamide gels. We provide tips on addressing and preventing common problems that can arise with the rapid-quench technique. Guidance is also offered on ensuring the ribozyme is properly folded and fast-reacting. We hope that this article will facilitate the broader use of rapid-quench instrumentation to study fast-reacting ribozymes under biological reaction conditions.

Topological constraints of structural elements in regulation of catalytic activity in HDV-like self-cleaving ribozymes

Self-cleaving ribozymes fold into intricate structures, which orient active site groups into catalytically competent conformations. Most ribozyme families have distinct catalytic cores stabilized by tertiary interactions between domains peripheral to those cores. We show that large hepatitis delta virus (HDV)-like ribozymes are activated by peripheral domains that bring two helical segments, P1 and P2, into proximity – a “pinch” that results in rate acceleration by almost three orders of magnitude. Kinetic analysis of ribozymes with systematically altered length and stability of the peripheral domain revealed that about one third of its free energy of formation is used to lower an activation energy barrier, likely related to a rate-limiting conformational change leading to the pre-catalytic state. These findings provide a quantitative view of enzyme regulation by peripheral domains and may shed light on the energetics of allosteric regulation.

Secondary structure confirmation and localization of Mg2+ ions in the mammalian CPEB3 ribozyme

Most of today's knowledge of the CPEB3 ribozyme, one of the few small self-cleaving ribozymes known to occur in humans, is based on comparative studies with the hepatitis delta virus (HDV) ribozyme, which is highly similar in cleavage mechanism and probably also in structure. Here we present detailed NMR studies of the CPEB3 ribozyme in order to verify the formation of the predicted nested double pseudoknot in solution. In particular, the influence of Mg(2+), the ribozyme's crucial cofactor, on the CPEB3 structure is investigated. NMR titrations, Tb(3+)-induced cleavage, as well as stoichiometry determination by hydroxyquinoline sulfonic acid fluorescence and equilibrium dialysis, are used to evaluate the number, location, and binding mode of Mg(2+)ions. Up to eight Mg(2+)ions interact site-specifically with the ribozyme, four of which are bound with high affinity. The global fold of the CPEB3 ribozyme, encompassing 80%-90% of the predicted base pairs, is formed in the presence of monovalent ions alone. Low millimolar concentrations of Mg(2+)promote a more compact fold and lead to the formation of additional structures in the core of the ribozyme, which contains the inner small pseudoknot and the active site. Several Mg(2+)binding sites, which are important for the functional fold, appear to be located in corresponding locations in the HDV and CPEB3 ribozyme, demonstrating the particular relevance of Mg(2+)for the nested double pseudoknot structure.

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.

An Ultrasensitive Light-up Cu(2+) Biosensor Using a New DNAzyme Cleaving a Phosphorothioate-Modified Substrate

Cu(2+) is a very important metal ion in biology, environmental science, and industry. Developing biosensors for Cu(2+) is a key topic in analytical chemistry. DNAzyme-based sensors are highly attractive for their excellent sensitivity, stability, and programmability. In the past decade, a few Cu(2+) biosensors were reported using DNAzymes with DNA cleavage or DNA ligation activity. However, they require unstable ascorbate or imidazole activation. So far, no RNA-cleaving DNAzymes specific for Cu(2+) are known. In this work, a phosphorothioate (PS) RNA-containing library was used for in vitro selection, and a few new Cu(2+)-specific RNA-cleaving DNAzymes were isolated. Among them, a DNAzyme named PSCu10 was studied further. It has only eight nucleotides in the enzyme loop with a cleavage rate of 0.1 min(-1) in the presence of 1 μM Cu(2+) at pH 6.0 (its optimal pH). Between the two diastereomers of the PS RNA chiral center, the R(p) isomer is 37 times more active than the S(p) one. Among the other divalent metal ions, only Hg(2+) can cleave the substrate due to its extremely high thiophilicity. A catalytic beacon sensor was designed with a detection limit of 1.6 nM Cu(2+) and extremely high selectivity. PSCu10 is specific for Cu(2+), and it has no cleavage in the presence of ascorbate, which reduces Cu(2+) to Cu(+).

A Two-Metal-Ion-Mediated Conformational Switching Pathway for HDV Ribozyme Activation

RNA enzymes serve as a potentially powerful platform from which to design catalysts and engineer new biotechnology. A fundamental understanding of these systems provides insight to guide design. The hepatitis delta virus ribozyme (HDVr) is a small, self-cleaving RNA motif widely distributed in nature, that has served as a paradigm for understanding basic principles of RNA catalysis. Nevertheless, questions remain regarding the precise roles of divalent metal ions and key nucleotides in catalysis. In an effort to establish a reaction mechanism model consistent with available experimental data, we utilize molecular dynamics simulations to explore different conformations and metal ion binding modes along the HDVr reaction path. Building upon recent crystallographic data, our results provide a dynamic model of the HDVr reaction mechanism involving a conformational switch between multiple non-canonical G25:U20 base pair conformations in the active site. These local nucleobase dynamics play an important role in catalysis by modulating the metal binding environments of two Mg²⁺ ions that support catalysis at different steps of the reaction pathway. The first ion plays a structural role by inducing a base pair flip necessary to obtain the catalytic fold in which C75 moves towards to the scissile phosphate in the active site. Ejection of this ion then permits a second ion to bind elsewhere in the active site and facilitate nucleophile activation. The simulations collectively describe a mechanistic scenario that is consistent with currently available experimental data from crystallography, phosphorothioate substitutions, and chemical probing studies. Avenues for further experimental verification are suggested.

Force Field for Mg(2+), Mn(2+), Zn(2+), and Cd(2+) Ions That Have Balanced Interactions with Nucleic Acids

Divalent metal ions are of fundamental importance to the function and folding of nucleic acids. Divalent metal ion-nucleic acid interactions are complex in nature and include both territorial and site specific binding. Commonly employed nonbonded divalent ion models, however, are often parametrized against bulk ion properties and are subsequently utilized in biomolecular simulations without considering any data related to interactions at specific nucleic acid sites. Previously, we assessed the ability of 17 different nonbonded Mg(2+) ion models to reproduce different properties of Mg(2+) in aqueous solution including radial distribution functions, solvation free energies, water exchange rates, and translational diffusion coefficients. In the present work, we depart from the recently developed 12-6-4 potential models for divalent metal ions developed by Li and Merz and tune the pairwise parameters for Mg(2+), Mn(2+), Zn(2+), and Cd(2+) binding dimethyl phosphate, adenosine, and guanosine in order to reproduce experimental site specific binding free energies derived from potentiometric pH titration data. We further apply these parameters to investigate a metal ion migration previously proposed to occur during the catalytic reaction of the hammerhead ribozyme. The new parameters are shown to be accurate and balanced for nucleic acid binding in comparison with available experimental data and provide an important tool for molecular dynamics and free energy simulations of nucleic acids where these ions may exhibit different binding modes.

Biochemical analysis of hatchet self-cleaving ribozymes

Hatchet RNAs are members of a novel self-cleaving ribozyme class that was recently discovered by using a bioinformatics search strategy. The consensus sequence and secondary structure of this class includes 13 highly conserved and numerous other modestly conserved nucleotides interspersed among bulges linking four base-paired substructures. A representative hatchet ribozyme from a metagenomic source requires divalent ions such as Mg(2+) to promote RNA strand scission with a maximum rate constant of ∼4 min(-1). As with all other small self-cleaving ribozymes discovered to date, hatchet ribozymes employ a general mechanism for catalysis involving the nucleophilic attack of a ribose 2'-oxygen atom on an adjacent phosphorus center. Kinetic characteristics of the reaction demonstrate that members of this ribozyme class have an essential requirement for divalent metal ions and that they might have a complex active site that employs multiple catalytic strategies to accelerate RNA cleavage by internal phosphoester transfer.

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

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

Desulfurization Activated Phosphorothioate DNAzyme for the Detection of Thallium

Thallium (Tl) is a highly toxic heavy metal situated between mercury and lead in the periodic table. While its neighbors have been thoroughly studied for DNA-based sensing, little is known about thallium detection. In this work, in vitro selection of RNA-cleaving DNAzymes is carried out using Tl(3+) as the target metal cofactor. Both normal DNA and phosphorothioate (PS)-modified DNA are tested for this purpose. While no Tl(3+)-dependent DNAzymes are obtained, a DNA oligonucleotide containing a single PS-modified RNA nucleotide is found to cleave by ∼7% by Tl(3+) at the RNA position. The remaining 93% are desulfurized. By hybridization of this PS-modified oligonucleotide with the Tm7 DNAzyme, the cleavage yield increases to ∼40% in the presence of Tl(3+) and Er(3+). Tm7 is an Er(3+)-dependent RNA-cleaving DNAzyme. It cleaves only the normal substrate but is completely inactive using the PS-modified substrate. Tl(3+) desulfurizes the PS substrate to the normal substrate to be cleaved by Tm7 and Er(3+). This system is engineered into a catalytic beacon for Tl(3+) with a detection limit of 1.5 nM, which is below its maximal contamination limit defined by the U.S. Environmental Protection Agency (10 nM).

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.

Understanding thio-effects in simple phosphoryl systems: role of solvent effects and nucleophile charge

Recent experimental work (J. Org. Chem., 2012, 77, 5829) demonstrated pronounced differences in measured thio-effects for the hydrolysis of (thio)phosphodichloridates by water and hydroxide nucleophiles. In the present work, we have performed detailed quantum chemical calculations of these reactions, with the aim of rationalizing the molecular bases for this discrimination. The calculations highlight the interplay between nucleophile charge and transition state solvation in SN2(P) mechanisms as the basis of these differences, rather than a change in mechanism.

Rational evolution of Cd2+-specific DNAzymes with phosphorothioate modified cleavage junction and Cd2+ sensing

In vitro selection of RNA-cleaving DNAzymes is a powerful method for isolating metal-specific DNA. A few successful examples are known, but it is still difficult to target some thiophilic metals such as Cd²⁺ due to limited functional groups in DNA. While using modified bases expands the chemical functionality of DNA, a single phosphorothioate modification might boost its affinity for thiophilic metals without complicating the selection process or using bases that are not commercially available. In this work, the first such in vitro selection for Cd²⁺ is reported. After using a blocking DNA and negative selections to rationally direct the library outcome, a highly specific DNAzyme with only 12 nucleotides in the catalytic loop is isolated. This DNAzyme has a cleavage rate of 0.12 min⁻¹ with 10 μM Cd²⁺ at pH 6.0. The R_p form of the substrate is cleaved ∼100-fold faster than the S_p form. The DNAzyme is most active with Cd²⁺ and its selectivity against Zn²⁺ is over 100 000-fold. Its application in detecting Cd²⁺ is also demonstrated. The idea of introducing single modifications in the fixed region expands the scope of DNA/metal interactions with minimal perturbation of DNA structure and property.

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.

Lanthanide-dependent RNA-cleaving DNAzymes as metal biosensors

Lanthanides represent a group of very important but challenging analytes for biosensor development. These 15 elements are very similar in their chemical properties. So far, limited success has been realized using the rational ligand design approach. My laboratory has successfully accomplished the task of carrying out combinatorial selection to isolate lanthanide-dependent RNA-cleaving DNAzymes. We report two new DNAzymes, each discovered in a different selection condition and both are highly specific to lanthanides. When both DNAzymes are used together, it is possible to identify the last few heavy lanthanides. Upon introducing a phosphorothioate modification, one of the abovementioned DNAzymes becomes highly active with many toxic heavy metals. With the selection of more DNAzymes with different activity patterns cross the lanthanide series, a sensor array might be produced for identifying each ion. This article is a minireview of the current developments on this topic and some of the historical aspects. It reflects the main content of the Fred Beamish Award presentation delivered at the 2014 Canadian Society for Chemistry Conference in Vancouver. Future directions in this area are also discussed.

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.

A new heavy lanthanide-dependent DNAzyme displaying strong metal cooperativity and unrescuable phosphorothioate effect

In vitro selection of RNA-cleaving DNAzymes was performed using three heavy lanthanide ions (Ln³⁺): Ho³⁺, Er³⁺ and Tm³⁺. The resulting sequences were aligned together and about half of the library contained a new family of DNAzyme. These DNAzymes have a simple loop structure, and they are active only with the seven heavy Ln³⁺. Among the tested non-lanthanide ions, only Y³⁺ induced cleavage and even Pb²⁺ failed to cleave, suggesting a very high specificity. A representative DNAzyme, Tm7, has a sigmoidal metal binding curve with a Hill coefficient of 3, indicating that three metal ions are involved in the catalytic step. Its pH-rate profile has a slope of 1, suggesting a single deprotonation step is involved in the rate-limiting step. Tm7 has a cleavage rate of 1.6 min⁻¹ at pH 7.8 with 10 μM Er³⁺. Phosphorothioate substitution at the cleavage junction completely inhibits the activity, which cannot be rescued by Cd²⁺ alone, or by a mixture of Er³⁺ and Cd²⁺, suggesting that two interacting metal ions are involved in direct bonding to both non-bridging oxygen atoms. A new model involving three lanthanide ions is proposed based on this study. A biosensor is engineered using Tm7 to detect Dy³⁺ down to 14 nM.

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.

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

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

Identification of minimal HDV-like ribozymes with unique divalent metal ion dependence in the human microbiome

HDV-like self-cleaving ribozymes have been found in a wide variety of organisms, implicated in diverse biological processes, and their activity typically shows a strong divalent metal dependence, but little metal specificity. Recent studies suggested that very short variants of these ribozymes exist in nature, but their distribution and biochemical properties have not been established. To map out the distribution of small HDV-like ribozymes, the drz-Spur-3 sequence was minimized to yield a core construct for structure-based bioinformatic searches. These searches revealed several microbial ribozymes, particularly in the human microbiome. Kinetic profile of the smallest ribozyme revealed two distinct metal binding sites, only one of which promotes fast catalysis. Furthermore, this ribozyme showed markedly reduced activity in Ca(2+), even in the presence of physiological Mg(2+) concentrations. Our study substantially expands the number of microbial HDV-like ribozymes and provides an example of cleavage regulation by divalent metals.

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.

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.

A catalytic role for methionine revealed by a combination of computation and experiments on phosphite dehydrogenase

Novel role for methionine in enzyme catalysis.