Solvent structure and hammerhead ribozyme catalysis

Although the hammerhead ribozyme is regarded as a prototype for understanding RNA catalysis, the mechanistic roles of associated metal ions and water molecules in the cleavage reaction remain controversial. We have investigated the catalytic potential of observed divalent metal ions and water molecules bound to a 2 A structure of the full-length hammerhead ribozyme by using X-ray crystallography in combination with molecular dynamics simulations. A single Mn(2+) is observed to bind directly to the A9 phosphate in the active site, accompanying a hydrogen-bond network involving a well-ordered water molecule spanning N1 of G12 (the general base) and 2'-O of G8 (previously implicated in general acid catalysis) that we propose, based on molecular dynamics calculations, facilitates proton transfer in the cleavage reaction. Phosphate-bridging metal interactions and other mechanistic hypotheses are also tested with this approach.     

Probing general base catalysis in the hammerhead ribozyme

Recent structural and computational studies have shed new light on the catalytic mechanism and active site structure of the RNA cleaving hammerhead ribozyme. Consequently, specific ribozyme functional groups have been hypothesized to be directly involved in general/acid base catalysis. In order to test this hypothesis, we have developed an affinity label to identify the functional general base in the S. mansoni hammerhead ribozyme. The ribozyme was reacted with a substrate analogue bearing a 2'-bromoacetamide group in place of the nucleophilic 2'-hydroxyl group which would normally be deprotonated by a general base. The electrophilic 2'-bromoacetamide group is poised to alkylate the general base, which is subsequently identified by footprinting analysis. Herein, we demonstrate alkylation of N1 of G12 in the hammerhead ribozyme in a pH and [Mg(2+)] dependent manner that is consistent with the native cleavage reaction. These results provide substantial evidence that deprotonated N1 of G12 functions directly as a general base in the hammerhead ribozyme; moreover, our experiments provide evidence that the pKa of G12 is perturbed downward in the context of the active site structure. We also observed other pH-independent alkylations, which do not appear to reflect the catalytic mechanism, but offer further insight into ribozyme conformation and structure.     

Diels-Alder ribozyme catalysis: a computational approach

A computational comparison of the Diels-Alder reaction of a maleimide and an anthracene in water and the active site of the ribozyme Diels-Alderase is reported. During the course of the catalyzed reaction, the maleimide is held in the hydrophobic pocket while the anthracene approaches to the maleimide through the back passage of the active site. The active site is so narrow that the anthracene has to adopt a tilted approach angle toward maleimide. The conformation of the active site changes marginally at different states of the reaction. Active site dynamics contribution to catalysis has been ruled out. The active site stabilizes the product more than the transition state (TS). The reaction coordinates of the ribozyme reaction in TS, RC1-CD1 and RC4-CD2, are 2.35 and 2.33 A, respectively, compared to 2.37 and 2.36 A in water. The approach angle of anthracene toward maleimide is twisted by 18 degrees in the TS structure of ribozyme reaction while no twisted angle is found in TS of the reaction in water. The free energy barriers for reactions in both ribozyme and water were obtained by umbrella sampling combined with SCCDFTB/MM. The calculated free energy barriers for the ribozyme and water reactions are in good agreement with the experimental values. As expected, Mulliken charges of the atoms involved in the ribozyme reaction change in a similar manner as that of the reaction in water. The proficiency of the Diels-Alder ribozyme reaction originates from the active site holding the two reactants in reactive conformations, in which the reacting atoms are brought together in van der Waals distances and reactants approach to each other at an appropriate angle.     

Allosterically activated Diels-Alder catalysis by a ribozyme

We describe the allosteric control of Diels-Alder reactions by a small organic effector, theophylline. This is achieved by converting a Diels-Alder ribozyme into an allosterically regulated system. In contrast to other published systems, we have a bond-forming reaction with two small-molecule substrates and multiple turnover. This system could be very attractive for the development of assays for a variety of analytes and can be regarded as a prototype of fully synthetic signaling cascades.     

Promiscuous catalysis by the tetrahymena group I ribozyme

Catalytic promiscuity, the ability of an enzyme to catalyze alternative reactions, has been suggested to have played an important role in the evolution of new catalytic activities in protein enzymes. Similarly, promiscuous activities may have been advantageous in an earlier RNA world. The Tetrahymena Group I ribozyme naturally catalyzes the site-specific guanosine attack on an anionic phosphate diester and has been shown to also catalyze aminoacyl transfer to water, albeit with a small rate acceleration (<10-fold). This inefficient catalysis could be due to the differences in charge and/or geometry requirements for the two reactions. Herein, we describe a new promiscuous activity of this ribozyme, the site-specific guanosine attack on a neutral phosphonate diester. This alternative substrate lacks the negative charge at the reaction center but, in contrast to the aminoacyl substrate, can undergo nucleophilic attack with the same geometry as the natural substrate. Our results show that the neutral phosphonate reaction is catalyzed about 1 x 106-fold, substantially better than the acyl transfer but far below the normal anionic substrate. We conclude that both charge and geometry are important factors for catalysis of the normal reaction and that promiscuous catalytic activities of ribozymes could have been created or enhanced by reorienting and swapping RNA domains.     

Essential role of an active-site guanine in glmS ribozyme catalysis

The glmS ribozyme is a catalytic riboswitch that is activated for endonucleolytic cleavage by the coenzyme glucosamine-6-phosphate. Using kinetic assays and X-ray crystallography, we identify an active-site mutation of a conserved guanine that abolishes catalysis without perturbing coenzyme binding. Our results provide evidence that coenzyme function requires a specific nucleobase to interact with the nucleophile of the cleavage reaction.     

An active-site guanine participates in glmS ribozyme catalysis in its protonated state

Active-site guanines that occupy similar positions have been proposed to serve as general base catalysts in hammerhead, hairpin, and glmS ribozymes, but no specific roles for these guanines have been demonstrated conclusively. Structural studies place G33(N1) of the glmS ribozyme of Bacillus anthracis within hydrogen-bonding distance of the 2'-OH nucleophile. Apparent pK(a) values determined from the pH dependence of cleavage kinetics for wild-type and mutant glmS ribozymes do not support a role for G33, or any other active-site guanine, in general base catalysis. Furthermore, discrepancies between apparent pK(a) values obtained from functional assays and microscopic pK(a) values obtained from pH-fluorescence profiles with ribozymes containing a fluorescent guanosine analogue, 8-azaguanosine, at position 33 suggest that the pH-dependent step in catalysis does not involve G33 deprotonation. These results point to an alternative model in which G33(N1) in its neutral, protonated form donates a hydrogen bond to stabilize the transition state.     

An unusual pH-independent and metal-ion-independent mechanism for hairpin ribozyme catalysis

Background: Hairpin ribozymes (RNA enzymes) catalyze the same chemical reaction as ribonuclease A and yet RNAs do not usually have functional groups analogous to the catalytically essential histidine and lysine sidechains of protein ribonucleases. Some RNA enzymes appear to recruit metal ions to act as Lewis acids in charge stabilization and metal-bound hydroxide for general base catalysis, but it has been reported that the hairpin ribozyme functions in the presence of metal ion chelators. This led us to investigate whether the hairpin ribozyme exploits a metal-ion-independent catalytic strategy.Results: Substitution of sulfur for nonbridging oxygens of the reactive phosphate of the hairpin ribozyme has small, stereospecific and metal-ionindependent effects on cleavage and ligation mediated by this ribozyme. Cobalt hexammine, an exchange-inert metal complex, supports full hairpin ribozyme activity, and the ribozyme's catalytic rate constants display only a shallow dependence on pH.Conclusions: Direct metal ion coordination to phosphate oxygens is not essential for hairpin ribozyme catalysis and metal-bound hydroxide does not serve as the general base in this catalysis. Several models might account for the unusual pH and metal ion independence: hairpin cleavage and ligation might be limited by a slow conformational change; a pH-independent or metalcation-independent chemical step, such as breaking the 5′ oxygen-phosphorus bond, might be rate determining; or finally, functional groups within the ribozyme might participate directly in catalytic chemistry. Whichever the case, the hairpin ribozyme appears to employ a unique strategy for RNA catalysis.     

Role of Mg2+ in hammerhead ribozyme catalysis from molecular simulation

Molecular dynamics simulations have been performed to investigate the role of Mg2+ in the full-length hammerhead ribozyme cleavage reaction. In particular, the aim of this work is to characterize the binding mode and conformational events that give rise to catalytically active conformations and stabilization of the transition state. Toward this end, a series of eight 12 ns molecular dynamics simulations have been performed with different divalent metal binding occupations for the reactant, early and late transition state using recently developed force field parameters for metal ions and reactive intermediates in RNA catalysis. In addition, hybrid QM/MM calculations of the early and late transition state were performed to study the proton-transfer step in general acid catalysis that is facilitated by the catalytic Mg2+ ion. The simulations suggest that Mg2+ is profoundly involved in the hammerhead ribozyme mechanism both at structural and catalytic levels. Binding of Mg2+ in the active site plays a key structural role in the stabilization of stem I and II and to facilitate formation of near attack conformations and interactions between the nucleophile and G12, the implicated general base catalyst. In the transition state, Mg2+ binds in a bridging position where it stabilizes the accumulated charge of the leaving group while interacting with the 2'OH of G8, the implicated general acid catalyst. The QM/MM simulations provide support that, in the late transition state, the 2'OH of G8 can transfer a proton to the leaving group while directly coordinating the bridging Mg2+ ion. The present study provides evidence for the role of Mg2+ in hammerhead ribozyme catalysis. The proposed simulation model reconciles the interpretation of available experimental structural and biochemical data, and provides a starting point for more detailed investigation of the chemical reaction path with combined QM/MM methods.     

Separate metal requirements for loop interactions and catalysis in the extended hammerhead ribozyme

An extended hammerhead ribozyme derived from Schistosoma mansoni, including conserved loops in stems I and II, has been examined to directly monitor the relationship between docking of loops and its activity using site-directed spin labeling (SDSL) and EPR spectroscopy. Dynamics with EPR spectroscopy and fast-quench kinetics measurements have shown that the docking of stems I and II occurs at low Mg2+ concentrations ([Mg2+]1/2,dock = 0.7 mM, 0.1 M NaCl), but a much weaker Mg2+ interaction ([Mg2+]1/2,cat approximately 90 mM) increases activity to very high maximum rates of approximately 1 s-1 at 0.1 M Na+ and pH 7.0.     

Phase-transfer catalysis in reactions with the participation of polyvinyl chloride

When polyvinyl chloride is dehydrochlorinated under conditions of phase transfer catalysis, the nature of the interface determines the supramolecular structure of the synthesized polyacetylene, the packing density of crystalline regions in the polyacetylene, and the chemical properties of the polymer. The ionic mechanism of catalysis of these reactions has been established.This review is based on materials reported at the Conference Phase-Transfer Catalysis: New Ideas and Methods(March, 1994).Translated fromIzvestiya Akademii Nauk. Seriya Khimicheskaya, No. 10, pp. 1886–1894, October, 1995.     

Pseudorotation barriers of biological oxyphosphoranes: a challenge for simulations of ribozyme catalysis

Pseudorotation reactions of biologically relevant oxyphosphoranes were studied by using density functional and continuum solvation methods. A series of 16 pseudorotation reactions involving acyclic and cyclic oxyphosphoranes in neutral and monoanionic (singly deprotonated) forms were studied, in addition to pseudorotation of PF5. The effect of solvent was treated by using three different solvation models for comparison. The barriers to pseudorotation ranged from 1.5 to 8.1 kcal mol(-1) and were influenced systematically by charge state, apicophilicity of ligands, intramolecular hydrogen bonding, cyclic structure and solvation. Barriers to pseudorotation for monoanionic phosphoranes occur with the anionic oxo ligand as the pivotal atom, and are generally lower than for neutral phosphoranes. The OCH3 groups were observed to be more apicophilic than OH groups, and hence pseudorotations that involve axial OCH3/equatorial OH exchange had higher reaction and activation free energy values. Solvent generally lowered barriers relative to the gas-phase reactions. These results, together with isotope 18O exchange experiments, support the assertion that dianionic phosphoranes are not sufficiently long-lived to undergo pseudorotation. Comparison of the density functional results with those from several semiempirical quantum models highlight a challenge for new-generation hybrid quantum mechanical/molecular mechanical potentials for non-enzymatic and enzymatic phosphoryl transfer reactions: the reliable modeling of pseudorotation processes.     

Quantum mechanical/molecular mechanical simulation study of the mechanism of hairpin ribozyme catalysis

The molecular mechanism of hairpin ribozyme catalysis is studied with molecular dynamics simulations using a combined quantum mechanical and molecular mechanical (QM/MM) potential with a recently developed semiempirical AM1/d-PhoT model for phosphoryl transfer reactions. Simulations are used to derive one- and two-dimensional potentials of mean force to examine specific reaction paths and assess the feasibility of proposed general acid and base mechanisms. Density-functional calculations of truncated active site models provide complementary insight to the simulation results. Key factors utilized by the hairpin ribozyme to enhance the rate of transphosphorylation are presented, and the roles of A38 and G8 as general acid and base catalysts are discussed. The computational results are consistent with available experimental data, provide support for a general acid/base mechanism played by functional groups on the nucleobases, and offer important insight into the ability of RNA to act as a catalyst without explicit participation by divalent metal ions.     

A unique mechanism for RNA catalysis: the role of metal cofactors in hairpin ribozyme cleavage

Background: Ribozymes are biological catalysts that promote the hydrolysis and transesterification of phosphate diesters of RNA. They typically require divalent magnesium ions for activation, although it has proven difficult to differentiate structural from catalytic roles for the magnesium ions and to identify the molecular mechanism of catalysis. Direct inner-sphere coordination is usually invoked in the catalytic step, although there is no evidence to support the generality of such a pathway for all ribozymes.Results: We studied the catalytic pathway for the hairpin class of ribozyme. The substitutionally inert transition metal complex cobalt hexaammine [Co(NH₃)₆³⁺) was shown to be as active as Mg²⁺(aq) in promoting hairpin ribozyme activity, demonstrating that inner-sphere pathways are not used by this class of ribozyme. These results were confirmed by studies with R_P- and S_P-phosphorothioate substrate analogs which show a similar reactivity to that of the native substrate towards the magnesium-activated ribozyme. Monovalent cations enhance the activity of Co(NH₃)₆³⁺-promoted reactions, but inhibit Mg²⁺-activated catalysis, demonstrating a requirement for hydrated cations at several key sites in the ribozyme.Conclusions: These results provide clear support for a model of RNA catalysis that does not involve direct coordination of magnesium to the phosphate ester, nor activation of a bound water molecule. A mechanism in which catalysis is carried out by functional groups on the RNA ribozyme itself is possible; such functional groups are likely to have pK_a values that are appropriate for carrying out this catalysis. The metal cofactor would then serve to define the architecture of the catalytic pocket and contribute to the stabilization of transient species, as has been described earlier. Hydrolytic pathways in nucleic acid reactions are apparently more diverse than was previously thought, and the hairpin ribozyme falls into a mechanistically distinct class from the Tetrahymena and the hammerhead ribozymes.     

Synthesis of novel C4-linked imidazole ribonucleoside phosphoramidites for probing general acid and base catalysis in ribozyme

We describe the synthesis of novel C4-linked imidazole ribonucleoside phosphoramidite (PA) 1a by which the imidazole moiety is incorporated into VS ribozyme to study its role in general acid and base catalysis. Investigation of protecting groups for the imidazole-N first indicated that pivaloyloxymethyl (POM) was adequate as an N-protecting group for the imidazole nucleoside, which could be readily removed under mild basic conditions. Further, the synthetic method was extended to synthesis of 2′-deoxy- and 2′-O-allyl nucleoside PAs 1b and 1c.     

Functional Systems Derived from Nucleobase Self-assembly

Dynamic and reversible non-covalent interactions endow synthetic systems and materials with smart adaptive functions that allow them to response to diverse stimuli, interact with external agents, or repair structural defects. Inspired by the outstanding performance and selectivity of DNA in living systems, scientists are increasingly employing Watson−Crick nucleobase pairing to control the structure and properties of self-assembled materials. Two sets of complementary purine-pyrimidine pairs (guanine:cytosine and adenine:thymine(uracil)) are available that provide selective and directional H-bonding interactions, present multiple metal-coordination sites, and exhibit rich redox chemistry. In this review, we highlight several recent examples that profit from these features and employ nucleobase interactions in functional systems and materials, covering the fields of energy/electron transfer, charge transport, adaptive nanoparticles, porous materials, macromolecule self-assembly, or polymeric materials with adhesive or self-healing ability.     

Rescue of an abasic hairpin ribozyme by cationic nucleobases: evidence for a novel mechanism of RNA catalysis

The hairpin ribozyme catalyzes a reversible phosphodiester cleavage reaction. We examined the roles of conserved nucleobases in catalysis using an abasic ribozyme rescue strategy. Loss of the active site G8 nucleobase reduced the cleavage rate constant by 350-fold while loss of A9 and A10 nucleobases reduced activity less than 10-fold. Certain heterocyclic amines restored partial activity when provided in solution to the variant lacking G8. Heterocyclic amines that were capable of rescue shared the exocyclic amine and cyclic amide in common with the Watson-Crick hydrogen bonding face of guanine. In contrast to the shallow pH dependence of unmodified ribozyme activity, rescue activity increased sharply with decreasing pH. These results support a novel model for RNA catalysis in which a cationic nucleobase contributes electrostatic stabilization to negative charge developing in the transition state.     

Origin of mutational effects at the C3 and G8 positions on hammerhead ribozyme catalysis from molecular dynamics simulations

A series of ten 60 ns molecular dynamics (MD) simulations of the native and mutated full length hammerhead ribozymes in the reactant state and in an activated precursor state (G8:2'OH deprotonated) are reported. Mutant simulations include the C3U, G8A, and G8I single mutants and a C3U/G8A double mutant that exhibits an experimental rescue effect. The results provide critical details into the origin of the observed mutation effects and support a mechanism where the 2'OH of G8 acts as a general acid catalyst that is held in position through Watson-Crick hydrogen bonding between G8 and C3.