Role of an active site adenine in hairpin ribozyme catalysis

The hairpin ribozyme is a small catalytic RNA that accelerates reversible cleavage of a phosphodiester bond. Structural and mechanistic studies suggest that divalent metals stabilize the functional structure but do not participate directly in catalysis. Instead, two active site nucleobases, G8 and A38, appear to participate in catalytic chemistry. The features of A38 that are important for active site structure and chemistry were investigated by comparing cleavage and ligation reactions of ribozyme variants with A38 modifications. An abasic substitution of A38 reduced cleavage and ligation activity by 14,000-fold and 370,000-fold, respectively, highlighting the critical role of this nucleobase in ribozyme function. Cleavage and ligation activity of unmodified ribozymes increased with increasing pH, evidence that deprotonation of some functional group with an apparent pK(a) value near 6 is important for activity. The pH-dependent transition in activity shifted by several pH units in the basic direction when A38 was substituted with an abasic residue, or with nucleobase analogs with very high or low pK(a) values that are expected to retain the same protonation state throughout the experimental pH range. Certain exogenous nucleobases that share the amidine group of adenine restored activity to abasic ribozyme variants that lack A38. The pH dependence of chemical rescue reactions also changed according to the intrinsic basicity of the rescuing nucleobase, providing further evidence that the protonation state of the N1 position of purine analogs is important for rescue activity. These results are consistent with models of the hairpin ribozyme catalytic mechanism in which interactions with A38 provide electrostatic stabilization to the transition state.     

The hairpin ribozyme

The hairpin ribozyme is a naturally occurring RNA that catalyzes sequence-specific cleavage and ligation of RNA. It has been the subject of extensive biochemical and structural studies, perhaps the most detailed for any catalytic RNA to date. Comparison of the structures of its constituent domains free and fully assembled demonstrates that the RNA undergoes extensive structural rearrangement. This rearrangement results in a distortion of the substrate RNA that primes it for cleavage. This ribozyme is known to achieve catalysis employing exclusively RNA functional groups. Metal ions or other catalytic cofactors are not used. Current experimental evidence points to a combination of at least four mechanistic strategies by this RNA: (1) precise substrate orientation, (2) preferential transition state binding, (3) electrostatic catalysis, and (4) general acid base catalysis.     

Electrostatic interactions in the hairpin ribozyme account for the majority of the rate acceleration without chemical participation by nucleobases

Molecular dynamics simulations using a combined quantum mechanical/molecular mechanical potential are used to determine the two-dimensional free energy profiles for the mechanism of RNA transphosphorylation in solution and catalyzed by the hairpin ribozyme. A mechanism is explored whereby the reaction proceeds without explicit chemical participation by conserved nucleobases in the active site. The ribozyme lowers the overall free energy barrier by up to 16 kcal/mol, accounting for the majority of the observed rate enhancement. The barrier reduction in this mechanism is achieved mainly by the electrostatic environment provided by the ribozyme without recruitment of active site nucleobases as acid or base catalysts. The results establish a baseline mechanism that invokes only the solvation and specific hydrogen-bonding interactions present in the ribozyme active site and provide a departure point for the exploration of alternate mechanisms where nucleobases play an active chemical role.     

Structural effects of nucleobase variations at key active site residue Ade38 in the hairpin ribozyme

The hairpin ribozyme requires functional groups from Ade38 to achieve efficient bond cleavage or ligation. To identify molecular features that contribute to catalysis, structures of position 38 base variants 2,6-diaminopurine (DAP), 2-aminopurine (AP), cytosine (Cyt), and guanine (Gua) were determined between 2.2 and 2.8 A resolution. For each variant, two substrate modifications were compared: (1) a 2'-O-methyl-substituent at Ade-1 was used in lieu of the nucleophile to mimic the precatalytic state, and (2) a 3'-deoxy-2',5'-phosphodiester linkage between Ade-1 and Gua+1 was used to mimic a reaction-intermediate conformation. While the global fold of each variant remained intact, the results revealed the importance of Ade38 N1 and N6 groups. Absence of N6 resulting from AP38 coincided with failure to localize the precatalytic scissile phosphate. Cyt38 severely impaired catalysis in a prior study, and its structures here indicated an anti base conformation that sequesters the imino moiety from the scissile bond. Gua38 was shown to be even more deleterious to activity. Although the precatalytic structure was nominally affected, the reaction-intermediate conformation indicated a severe electrostatic clash between the Gua38 keto oxygen and the pro-Rp oxygen of the scissile bond. Overall, position 38 modifications solved in the presence of 2'-OMe Ade-1 deviated from in-line geometry, whereas variants with a 2',5' linkage exhibited S-turn destabilization, as well as base conformational changes from syn to anti. These findings demonstrate the importance of the Ade38 Watson-Crick face in attaining a reaction-intermediate state and the sensitivity of the RNA fold to restructuring when electrostatic and shape features fail to complement.     

An important role of G638 in the cis-cleavage reaction of the Neurospora VS ribozyme revealed by a novel nucleotide analog incorporation method

We describe a chemical coupling procedure that allows joining of two RNAs, one of which contains a site-specific base analog substitution, in the absence of divalent ions. This method allows incorporation of nucleotide analogs at specific positions even into large, cis-cleaving ribozymes. Using this method we have studied the effects of substitution of G638 in the cleavage site loop of the VS ribozyme with a variety of purine analogs having different functional groups and pK(a) values. Cleavage rate versus pH profiles combined with kinetic solvent isotope experiments indicate an important role for G638 in proton transfer during the rate-limiting step of the cis-cleavage reaction.     

The pH dependence of hairpin ribozyme catalysis reflects ionization of an active site adenine

Understanding how self-cleaving ribozymes mediate catalysis is crucial in light of compelling evidence that human and bacterial gene expression can be regulated through RNA self-cleavage. The hairpin ribozyme catalyzes reversible phosphodiester bond cleavage through a mechanism that does not require divalent metal cations. Previous structural and biochemical evidence implicated the amidine group of an active site adenosine, A38, in a pH-dependent step in catalysis. We developed a way to determine microscopic pK(a) values in active ribozymes based on the pH-dependent fluorescence of 8-azaadenosine (8azaA). We compared the microscopic pK(a) for ionization of 8azaA at position 38 with the apparent pK(a) for the self-cleavage reaction in a fully functional hairpin ribozyme with a unique 8azaA at position 38. Microscopic and apparent pK(a) values were virtually the same, evidence that A38 protonation accounts for the decrease in catalytic activity with decreasing pH. These results implicate the neutral unprotonated form of A38 in a transition state that involves formation of the 5'-oxygen-phosphorus bond.     

The ionic environment determines ribozyme cleavage rate by modulation of nucleobase pK a

Several small ribozymes employ general acid–base catalysis as a mechanism to enhance site-specific RNA cleavage, even though the functional groups on the ribonucleoside building blocks of RNA have pK_a values far removed from physiological pH. The rate of the cleavage reaction is strongly affected by the identity of the metal cation present in the reaction solution; however, the mechanism(s) by which different cations contribute to rate enhancement has not been determined. Using the Neurospora VS ribozyme, we provide evidence that different cations confer particular shifts in the apparent pK_a values of the catalytic nucleobases, which in turn determines the fraction of RNA in the protonation state competent for general acid–base catalysis at a given pH, which determines the observed rate of the cleavage reaction. Despite large differences in observed rates of cleavage in different cations, mathematical models of general acid–base catalysis indicate that k₁, the intrinsic rate of the bond-breaking step, is essentially constant irrespective of the identity of the cation(s) in the reaction solution. Thus, in contrast to models that invoke unique roles for metal ions in ribozyme chemical mechanisms, we find that most, and possibly all, of the ion-specific rate enhancement in the VS ribozyme can be explained solely by the effect of the ions on nucleobase pK_a. The inference that k₁ is essentially constant suggests a resolution of the problem of kinetic ambiguity in favor of a model in which the lower pK_a is that of the general acid and the higher pK_a is that of the general base.     

Formation of an active site in trans by interaction of two complete Varkud Satellite ribozymes

The complete VS ribozyme comprises seven helical segments, connected by three three-way RNA junctions. In the presence of Mg²⁺ ions, cleavage occurs within the internal loop of helix I. This requires the participation of a guanine (G638) within the helix I loop, and a remote adenine (A756) within an internal loop of helix VI. Previous structural studies have suggested that helix I docks into the fold of the remaining part of the ribozyme, bringing A756 and G638 close to the scissile phosphate to allow the cleavage reaction to proceed. We show here that while either A756C or G638A individually exhibit very low cleavage activity, a mixture of the two variants leads to cleavage of the A756C RNA, but not the G638A RNA. The rate of cleavage depends on the concentration of the VS G638A RNA, as expected for a bimolecular interaction. This regaining of cleavage activity by complementation indicates that helix I of one VS RNA can interact with another VS RNA molecule to generate a functional active site in trans.     

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.     

A pH-dependent conformational change, rather than the chemical step, appears to be rate-limiting in the hammerhead ribozyme cleavage reaction.

We have investigated the chemical basis for a previously observed 7.8 A conformational change in the hammerhead ribozyme that positions the substrate for in-line attack. We have found that the conformational change can only be observed at or above pH 8.5 (in the presence of Co(2+)) and requires the presence of an ionizable 2'-OH at the cleavage site, and note that this observed apparent pK(a) of 8.5 for the conformational change is within experimental error (+/-0.5) of the previously reported apparent kinetic pK(a) of 8.5 for the hammerhead ribozyme in the presence of Co(2+). We have solved two crystal structures of hammerhead ribozymes having 2'-OCH(3) or 2'-F substitutions at the cleavage site and have found that these will not undergo a conformational change equivalent to that observed for the hammerhead ribozyme having an unmodified attacking nucleophile under otherwise identical conditions. We have also characterized the kinetics of cleavage in the crystal. In addition to verifying that the particular sequence of RNA that we crystallized cleaves faster in the crystal than in solution, we also find that the extent of cleavage in the crystal is complete, unlike in solution where this and most other hammerhead ribozyme substrates are cleaved only to about 70 % completion. The initial cleavage rate in the crystal obeys the expected log-linear relation between cleavage-rate and pH with a slope of 0.7, as has been observed for other hammerhead ribozyme sequences in solution, indicating that in both the crystal and in solution the pH-dependent step is rate-limiting. However, the cleavage rate in the crystal is biphasic, with the most dramatic distinction between initial (slower) and final (faster) phases appearing at pH 6.0. The initial phase corresponds to the pH-dependent cleavage rate observed in solution, but the second, faster phase is roughly pH-independent and closely parallels the cleavage rate observed at pH 8 (0.4/minute). This result is particularly remarkable because it entails that the rapidly cleaving phase at pH 6 is comparable to the cleavage rate for the fastest cleaving hammerhead ribozymes at pH 6. Based upon these observations, we conclude that the pH-dependent conformational change is the rate-determining step under standard conditions for the hammerhead ribozyme self-cleavage reaction, and that an ionizable 2'-proton at cleavage site is required for this conformational change. We further hypothesize that deprotonation of the cleavage-site 2'-oxygen drives this conformational change.     

Functional group requirements in the probable active site of the VS ribozyme

The VS ribozyme catalyses the site-specific cleavage of a phosphodiester linkage by a transesterification reaction that entails the attack of the neighbouring 2'-oxygen with departure of the 5'-oxygen. We have previously suggested that the A730 loop is an important component of the active site of the ribozyme, and that A756 is especially important in the cleavage reaction. Functional group modification experiments reported here indicate that the base of A756 is more important than its ribose for catalysis. A number of changes to the base, including complete ablation, lead to cleavage rates that are reduced 1000-fold, while removal of the 2'-hydroxyl group from the ribose results in tenfold slower cleavage. 2-Aminopurine fluorescence experiments indicate that this 2'-hydroxyl group is important for the structure of the A730 loop. Catalytic activity is especially sensitive to changes involving the exocyclic amine of A756; by contrast, the cleavage activity is only weakly sensitive to modification at the 7-position of the purine nucleus. These results suggest that the Watson-Crick edge of the adenine base is important in ribozyme function. We sought to test the possibility of a direct role of the nucleobase in the chemistry of the cleavage reaction. Addition of imidazole base in the medium failed to restore the activity of a ribozyme from which the nucleobase of A756 was removed. However, no restoration was obtained with exogenous adenine base either, indicating that the cavity that might result from ablation of the base was closed.     

Catalytic roles for proton transfer and protonation in ribozymes

Utilization of proton transfer in catalysis, which is well known in the mechanisms of protein enzymes, has been described only relatively recently for RNA enzymes. In this article, we present a current understanding of proton transfer by nucleic acids. Rate enhancement and specificity conferred by general acid-base catalysis are discussed. We also present possibilities for electrostatic catalysis from general acids and bases as well as cationic base pairs. The microenvironments of a large RNA provide the possibility of histidine-like pK(a)s for proton transfer, as well as lysine- and arginine-like pK(a)s for electrostatic catalysis. Discussion on proton transfer focuses on the hepatitis delta virus (HDV) and hairpin ribozymes, with select examples drawn from the protein literature. Discussion on electrostatic catalysis also draws on these two ribozymes, and a postulate for electrostatic catalysis by a cationic base pair in the mechanism of peptidyl transfer in the ribosome is presented. We also provide a perspective on possibilities for phosphoryl transfer mechanisms involving phosphorane intermediates and unusual tautomeric forms of the bases. Lastly, a distinction is made between ground state and "transition state" pK(a)s. We favor a model in which changes in pH lead to changes in the distribution of reactive and nonreactive ionizations of the ribozyme molecules in the ground state, and therefore suggest that "pK(a) changes in the transition state" do not provide an acceptable explanation for observed pH-rate profiles.     

Role of SLV in SLI substrate recognition by the Neurospora VS ribozyme

Substrate recognition by the VS ribozyme involves a magnesium-dependent loop/loop interaction between the SLI substrate and the SLV hairpin from the catalytic domain. Recent NMR studies of SLV demonstrated that magnesium ions stabilize a U-turn loop structure and trigger a conformational change for the extruded loop residue U700, suggesting a role for U700 in SLI recognition. Here, we kinetically characterized VS ribozyme mutants to evaluate the contribution of U700 and other SLV loop residues to SLI recognition. To help interpret the kinetic data, we structurally characterized the SLV mutants by NMR spectroscopy and generated a three-dimensional model of the SLI/SLV complex by homology modeling with MC-Sym. We demonstrated that the mutation of U700 by A, C, or G does not significantly affect ribozyme activity, whereas deletion of U700 dramatically impairs this activity. The U700 backbone is likely important for SLI recognition, but does not appear to be required for either the structural integrity of the SLV loop or for direct interactions with SLI. Thus, deletion of U700 may affect other aspects of SLI recognition, such as magnesium ion binding and SLV loop dynamics. As part of our NMR studies, we developed a convenient assay based on detection of unusual (31)P and (15)N N7 chemical shifts to probe the formation of U-turn structures in RNAs. Our model of the SLI/SLV complex, which is compatible with biochemical data, leads us to propose novel interactions at the loop I/loop V interface.     

The active site histidines of creatine kinase. A critical role of His 61 situated on a flexible loop.

A histidine residue with a pKa of 7 has been inferred to act as a general acid-base catalyst for the reaction of creatine kinase (CK), catalyzing the reversible phosphorylation of creatine by ATP. The chicken sarcomeric muscle mitochondrial isoenzyme Mib-CK contains several histidine residues that are conserved throughout the family of creatine kinases. By X-ray crystal structure analysis, three of them (His 61, His 92, and His 186) were recently shown to be located close to the active site of the enzyme. These residues were exchanged against alanine or aspartate by in vitro mutagenesis, and the six mutant proteins were expressed in E. coli and purified. Structural integrity of the mutant proteins was checked by small-angle X-ray scattering. Kinetic analysis showed the mutant His 61 Asp to be completely inactive in the direction of ATP consumption while exhibiting a residual activity of 1.7% of the wild-type (wt) activity in the reverse direction. The respective His to Ala mutant of residue 61 showed approximately 1% wt activity in the forward and 10% wt activity in the reverse reaction. All other mutants showed near wt activities. Changes in the kinetic parameters K(m) or Vmax, as well as a significant loss of synergism in substrate binding, could be observed with all active mutants. These effects were most pronounced for the binding of creatine and phosphocreatine, whereas ATP or ADP binding were less severely affected. Based on our results, we assume that His 92 and His 186 are involved in the binding of creatine and ATP in the active site, whereas His 61 is of importance for the catalytic reaction but does not serve as an acid-base catalyst in the transphosphorylation of creatine and ATP. In addition, our data support the idea that the flexible loop bearing His 61 is able to move towards the active site and to participate in catalysis.     

Screening of the active site from water by the incoming ligand triggers catalysis and inhibition in serine proteases

The pKa of the catalytic His57 N(epsilon)H in the tetrahedral complex (TC) of chymotrypsin with trifluoromethyl ketone inhibitors is 4-5 units higher relative to the free enzyme (FE). Such stable TC's, formed with transition state (TS) analog inhibitors, are topologically similar to the catalytic TS. Thus, analysis of this pKa shift may shed light on the role of water solvation in the general base catalysis by histidine. We applied our QM/SCRF(VS) approach to study this shift. The method enables explicit quantum mechanical DFT calculations of large molecular clusters that simulate chemical reactions at the active site (AS) of water solvated enzymes. We derived an analytical expression for the pKa dependence on the degree of water exposure of the ionizable group, and on the total charge in the enzyme AS, Q(A) and Q(B), when the target ionizable functional group (His57 in this study) is in the acidic (A) and basic (B) forms, respectively. Q2(B) > Q2(A) both in the FE and in the TC of chymotrypsin. Therefore, water solvation decreases the relative stability of the protonated histidine in both. Ligand binding reduces the degree of water solvation of the imidazole ring, and consequently elevates the histidine pKa. Thus, the binding of the ligand plays a triggering role that switches on the cascade of catalytic reactions in serine proteases.     

The catalytic mechanism of Drosophila alcohol dehydrogenase: evidence for a proton relay modulated by the coupled ionization of the active site Lysine/Tyrosine pair and a NAD+ ribose OH switch

The ionization properties of the active site residues in Drosophila lebanonensis alcohol dehydrogenase (DADH) were investigated theoretically by using an approach developed to account for multiple locations of the hydrogen atoms of the titratable and polar groups. The electrostatic calculations show that (a) the protonation/deprotonation transition of the binary complex of DADH is related to the coupled ionization of Tyr151 and Lys155 in the active site and (b) the pH dependence of the proton abstraction is correlated with a reorganization of the hydrogen bond network in the active site. On this basis, a proton relay mechanism for substrate dehydrogenation is proposed in which the O2' ribose hydroxyl group from the coenzyme has a key role and acts as a switch. The proton relay chain includes the active site catalytic residues, as well as a chain of eight water molecules that connects the active site with the bulk solvent.     

Probing the location and function of the conserved histidine residue of phosphoglucose isomerase by using an active site directed inhibitor N-bromoacetylethanolamine phosphate

Phosphoglucose isomerase (EC 5.3.1.9) catalyzes the interconversion of D-glucopyranose-6-phosphate and D-fructofuranose-6-phosphate by promoting an intrahydrogen transfer between C1 and C2. A conserved histidine exists throughout all phosphoglucose isomerases and was hypothesized to be the base catalyzing the isomerization reaction. In the present study, this conserved histidine, His311, of the enzyme from Bacillus stearothermophilus was subjected to mutational analysis, and the mutational effect on the inactivation kinetics by N-bromoacetylethanolamine phosphate was investigated. The substitution of His311 with alanine, asparagine, or glutamine resulted in the decrease of activity, in k(cat)/K(M), by a factor of 10(3), indicating the importance of this residue. N-bromoacetylethanolamine phosphate inactivated irreversibly the activity of wild-type phosphoglucose isomerase; however, His311 --> Ala became resistant to this inhibitor, indicating that His311 is located in the active site and is responsible for the inactivation of the enzyme by this active site-directed inhibitor. The pKa of His311 was estimated to be 6.31 according to the pH dependence of the inactivation. The proximity of this value with the pKa value of 6.35, determined from the pH dependence of k(cat)/K(M), supports a role of His311 as a general base in the catalysis.     

The crystal structure of a (-) gamma-lactamase from an Aureobacterium species reveals a tetrahedral intermediate in the active site

The structure of the recombinant (-) gamma-lactamase from an Aureobacterium species has been solved at 1.73A resolution in the cubic space group F23 with unit cell parameters a=b=c=240.6A. The trimeric enzyme has an alpha/beta hydrolase fold and closely resembles the cofactor free haloperoxidases. The structure has been solved in complex with a covalently bound ligand originating from the host cell and also in the unligated form. The associated density in the former structure has been interpreted as the two-ring ligand (3aR,7aS)-3a,4,7,7a-tetrahydro-benzo [1,3] dioxol-2-one which forms a tetrahedral complex with OG of the catalytic Ser98. Soaks of these crystals with the industrial substrate gamma-lactam or its structural analogue, norcamphor, result in the displacement of the ligand from the enzyme active site, thereby allowing determination of the unligated structure. The presence of the ligand in the active site protects the enzyme from serine hydrolase inhibitors. Cyclic ethylene carbonate, the first ring of the ligand, was shown to be a substrate of the enzyme.     

The role of intramolecular hydrogen bonding on nucleobase acidification following metal coordination: possible implications of an “indirect” role of metals in acid–base catalysis of nucleic acids

The acidifying effect of Pt^II on nucleobase –NH and –NH₂ groups depends both on the site of metal coordination and on the efficiency of stabilization of the deprotonated nucleobase via intracomplex hydrogen bonding. Weakly acidic nucleobase protons with pK _a values between 9 and 17 can be acidified by a single Pt^II to have pK _a values which are well within the physiological pH range. This could open the possibility of an acid–base catalysis occurring at pH 7, with the metal–nucleobase entity functioning either as an acid or a base. Examples of Pt^II complexes studied here include, among others, mixed nucleobase systems of 1-methylcytosine and 1,9-dimethyladenine as well as a complex of the rare iminooxo tautomer of 1-methylcytosine having the metal bonded at N4.