Theoretical studies on the hydrolysis of phosphate diesters in the gas phase,solution, and RNase A

Density functional theory, polarizable continuum models and semiempirical hybrid quantum mechanical/molecular mechanical (QM/MM) calculations were applied to the hydrolysis of phosphate diesters in the gas phase, in solution, and in the enzyme RNase A. Neutralization of the negative charge of the pentacovalent phosphorane intermediates provides a substantial stabilization of the transition‐state structures in the gas phase. Inclusion of solvent effects on the phosphate/phosphorane species was critical to reproducing the trends in reactivity observed experimentally. Finally, the catalytic mechanism for the hydrolysis of uridine 2′,3′‐cyclic phosphate by RNase A was studied by QM/MM calculations. Our results suggest that the rate‐limiting transition state of the reaction corresponds to the approach of a water molecule to the phosphate and its activation by His119. Thus, His119 acts as a generalized base for the reaction. The water attack leads to a pentacovalent phosphorane transition state of formal charge ?2; this excess of negative charge in the transition state is stabilized by a number of positively charged residues including His12 and Lys41. In the second stage of the reaction, the phosphorane is converted into products. This part of the reaction proceeds without a detectable barrier, and it is facilitated by a proton transfer from Lys41 to the departing O_2′. © 2001 John Wiley & Sons, Inc. Int J Quantum Chem, 2001     

The active site of cellobiohydrolase Cel6A from Trichoderma reesei: the roles of aspartic acids D221 and D175

Trichoderma reesei cellobiohydrolase Cel6A is an inverting glycosidase. Structural studies have established that the tunnel-shaped active site of Cel6A contains two aspartic acids, D221 and D175, that are close to the glycosidic oxygen of the scissile bond and at hydrogen-bonding distance from each other. Here, site-directed mutagenesis, X-ray crystallography, and enzyme kinetic studies have been used to confirm the role of residue D221 as the catalytic acid. D175 is shown to affect protonation of D221 and to contribute to the electrostatic stabilization of the partial positive charge in the transition state. Structural and modeling studies suggest that the single-displacement mechanism of Cel6A may not directly involve a catalytic base. The value of (D2O)(V) of 1.16 +/- 0.14 for hydrolysis of cellotriose suggests that the large direct effect expected for proton transfer from the nucleophilic water through a water chain (Grotthus mechanism) is offset by an inverse effect arising from reversibly breaking the short, tight hydrogen bond between D221 and D175 before catalysis.     

Transition state differences in hydrolysis reactions of alkyl versus aryl phosphate monoester monoanions

Although aryl phosphates have been the subject of numerous experimental studies, far less data bearing on the mechanism and transition states for alkyl phosphate reactions have been presented. Except for esters with very good leaving groups such as 2,4-dinitrophenol, the monoanion of phosphate esters is more reactive than the dianion. Several mechanisms have been proposed for the hydrolysis of the monoanion species. (18)O kinetic isotope effects in the nonbridging oxygen atoms and in the P-O(R) ester bond, and solvent deuterium isotope effects, have been measured for the hydrolysis of m-nitrobenzyl phosphate. The results rule out a proposed mechanism in which the phosphoryl group deprotonates water and then undergoes attack by hydroxide. The results are most consistent with a preequilibrium proton transfer from the phosphoryl group to the ester oxygen atom, followed by rate-limiting P-O bond fission, as originally proposed by Kirby and co-workers in 1967. The transition state for m-nitrobenzyl phosphate (leaving group pK(a) 14.9) exhibits much less P-O bond fission than the reaction of the more labile p-nitrophenyl phosphate (leaving group pK(a) = 7.14). This seemingly anti-Hammond behavior results from weakening of the P-O(R) ester bond resulting from protonation, an effect which calculations have shown is much more pronounced for aryl phosphates than for alkyl ones.     

Evidence for direct attack by hydroxide in phosphodiester hydrolysis

Phosphodiester hydrolysis has been the subject of intense study due to its importance in biology. Despite the numerous significant analyses of phosphodiester cleavage mechansim, comparatively little is known about the nucleophiles in these reactions. To determine whether hydroxide acts as a nucleophile or a general base in the hydrolysis of thymidine-5'-p-nitrophenyl phosphate,we determined solvent deuterium isotope effects (D2Ok), ionic strength effects, and 18O isotope effects on the solvent nucleophile (18knuc). The D2Ok for hydroxide-catalyzed phosphodiester hydrolysis is slightly inverse (0.9 +/- 0.1), suggesting that a proton transfer does not occur in the transition state. A significant alpha effect is observed with hydroperoxide, demonstrating that oxyanions can act as nucleophiles in the reaction. Additionally, the ionic strength dependencies of hydroxide and hydroperoxide catalysis are indistinguishable, suggesting that they act by the same mechanism. Finally, the 18knuc for the hydroxide-catalyzed reaction is 1.068 +/- 0.007, well in excess of the equilibrium 18O isotope effect between water and hydroxide (1.040 +/- 0.003). Together, the data are most consistent with direct nucleophilic attack by hydroxide. From the observed 18knuc and the known equilibrium component, the kinetic component of the isotope effect was calculated to be 1.027 +/- 0.010. This large kinetic component suggests that little bond order to the nucleophile occurs in the transition state.     

Mechanism of sulfide oxidations by peroxymonocarbonate

A detailed mechanism for the oxidation of aryl sulfides by peroxymonocarbonate ion in cosolvent/water media is described. Kinetic studies were performed to characterize the transition state, including a Hammett correlation and variation of solvent composition. The results are consistent with a charge-separated transition state relative to the reactants, with an increase of positive charge on the sulfur following nucleophilic attack of the sulfide at the electrophilic oxygen of peroxymonocarbonate. In addition, an average solvent isotope effect of 1.5 +/- 0.2 for most aryl sulfide oxidations is consistent with proton transfer in the transition state of the rate-determining step. Activation parameters for oxidation of ethyl phenyl sulfide in tert-butyl alcohol/water are reported. From the pH dependence of oxidation rates and (13)C NMR equilibrium experiments, the estimated pK(a) of peroxymonocarbonate was found to be approximately 10.6.     

Characterization of the transition state of functional enzyme dynamics

Through characterization of the solvent isotope effect on protein dynamics, we have examined determinants of the rate limitation to enzyme catalysis. A global conformational change in Ribonuclease A limits the overall rate of catalytic turnover. Here we show that this motion is sensitive to solvent deuterium content; the isotope effect is 2.2, a value equivalent to the isotope effect on the catalytic rate constant. We further demonstrate that the protein motion possesses a linear proton inventory plot, indicating that a single proton is transferred in the transition state. These results provide compelling evidence for close coupling between enzyme dynamics and function and demonstrate that characterization of the transition state for protein motion in atomic detail is experimentally accessible.     

Transition-state analysis of the DNA repair enzyme MutY

The transition state (TS) structure of MutY-catalyzed DNA hydrolysis was solved using multiple kinetic isotope effect (KIE) measurements. MutY is a base excision repair enzyme which cleaves adenine from 8-oxo-G:A mismatches in vivo, and also from G:A mismatches in vitro. TS analysis of G:A-DNA hydrolysis revealed a stepwise S(N)1 (D(N)*A(N)(double dagger)) mechanism proceeding through a highly reactive oxacarbenium ion intermediate which would have a lifetime in solution of <10(-10) s. C-N bond cleavage is reversible; the N-glycoside bond breaks and reforms repeatedly before irreversible water attack on the oxacarbenium ion. KIEs demonstrated that MutY uses general acid catalysis by protonating N7. It enforces a 3'-exo sugar ring conformation and other sugar ring distortions to stabilize the oxacarbenium ion. Combining the experimental TS structure with the previously reported crystal structure of an abortive Michaelis complex elucidates the step-by-step catalytic sequence.     

A comprehensive theoretical study on the hydrolysis of carbonyl sulfide in the neutral water

The detailed hydration mechanism of carbonyl sulfide (COS) in the presence of up to five water molecules has been investigated at the level of HF and MP2 with the basis set of 6-311++G(d, p). The nucleophilic addition of water molecule occurs in a concerted way across the C==S bond of COS rather than across the C==O bond. This preferential reaction mechanism could be rationalized in terms of Fukui functions for the both nucleophilic and electrophilic attacks. The activation barriers, DeltaH( not equal) (298), for the rate-determining steps of one up to five-water hydrolyses of COS across the C==S bond are 199.4, 144.4, 123.0, 115.5, and 107.9 kJ/mol in the gas phase, respectively. The most favorable hydrolysis path of COS involves a sort of eight-membered ring transition structure and other two water molecules near to the nonreactive oxygen atom but not involved in the proton transfer, suggesting that the hydrolysis of COS can be significantly mediated by the water molecule(s) and the cooperative effects of the water molecule(s) in the nonreactive region. The catalytic effect of water molecule(s) due to the alleviation of ring strain in the proton transfer process may result from the synergistic effects of rehybridization and charge reorganization from the precoordination complex to the rate-determining transition state structure induced by water molecule. The studies on the effect of temperature on the hydrolysis of COS show that the higher temperature is unfavorable for the hydrolysis of COS. PCM solvation models almost do not modify the calculated energy barriers in a significant way.     

New insights into the mechanism of the Schiff base hydrolysis catalyzed by type I dehydroquinate dehydratase from S. enterica: a theoretical study

The reaction pathway of Schiff base hydrolysis catalyzed by type I dehydroquinate dehydratase (DHQD) from S. enterica has been studied by performing molecular dynamics (MD) simulations and density functional theory (DFT) calculations and the corresponding potential energy profile has also been identified. On the basis of the results, the catalytic hydrolysis process for the wild-type enzyme consists of three major reaction steps, including nucleophilic attack on the carbon atom involved in the carbon-nitrogen double bond of the Schiff base intermediate by a water molecule, deprotonation of the His143 residue, and dissociation between the product and the Lys170 residue of the enzyme. The remarkable difference between this and the previously proposed reaction mechanism is that the second step here, absent in the previously proposed reaction mechanism, plays an important role in facilitating the reaction through a key proton transfer by the His143 residue, resulting in a lower energy barrier. Comparison with our recently reported results on the Schiff base formation and dehydration processes clearly shows that the Schiff base hydrolysis is rate-determining in the overall reaction catalyzed by type I DHQD, consistent with the experimental prediction, and the calculated energy barrier of ～16.0 kcal mol(-1) is in good agreement with the experimentally derived activation free energy of ～14.3 kcal mol(-1). When the imidazole group of His143 residue is missing, the Schiff base hydrolysis is initiated by a hydroxide ion in the solution, rather than a water molecule, and both the reaction mechanism and the kinetics of Schiff base hydrolysis have been remarkably changed, clearly elucidating the catalytic role of the His143 residue in the reaction. The new mechanistic insights obtained here will be valuable for the rational design of high-activity inhibitors of type I DHQD as non-toxic antimicrobials, anti-fungals, and herbicides.     

Semiempirical Molecular Orbital Studies of the Acylation Step in the Lipase‐Catalyzed Ester Hydrolysis

In this study, we present the results from the semiempirical molecular orbital calculations for the acylation step in the lipase‐catalyzed ester hydrolysis. The results reveal that the lowest energy path for the formation of the tetrahedral intermediate is for the serine residue of the catalytic triad to attack the substrate, followed by coupling heavy atom movement and proton transfer. The calculations of four active site models show that the cooperation of the aspartate group and the oxyanion hole is capable of lowering the activation energy by about 16 kcalmol^?1. Our results further suggest that the lipase‐catalyzed ester hydrolysis adopts the single proton transfer mechanism.     

Analysis of classical and quantum paths for deprotonation of methylamine by methylamine dehydrogenase.

The hydrogen-transfer reaction catalysed by methylamine dehydrogenase (MADH) with methylamine (MA) as substrate is a good model system for studies of proton tunnelling in enzyme reactions--an area of great current interest--for which atomistic simulations will be vital. Here, we present a detailed analysis of the key deprotonation step of the MADH/MA reaction and compare the results with experimental observations. Moreover, we compare this reaction with the related aromatic amine dehydrogenase (AADH) reaction with tryptamine, recently studied by us, and identify possible causes for the differences observed in the measured kinetic isotope effects (KIEs) of the two systems. We have used combined quantum mechanics/molecular mechanics (QM/MM) techniques in molecular dynamics simulations and variational transition state theory with multidimensional tunnelling calculations averaged over an ensemble of paths. The results reveal important mechanistic complexity. We calculate activation barriers and KIEs for the two possible proton transfers identified-to either of the carboxylate oxygen atoms of the catalytic base (Asp428beta)-and analyse the contributions of quantum effects. The activation barriers and tunnelling contributions for the two possible proton transfers are similar and lead to a phenomenological activation free energy of 16.5+/-0.9 kcal mol(-1) for transfer to either oxygen (PM3-CHARMM calculations applying PM3-SRP specific reaction parameters), in good agreement with the experimental value of 14.4 kcal mol(-1). In contrast, for the AADH system, transfer to the equivalent OD1 was found to be preferred. The structures of the enzyme complexes during reaction are analysed in detail. The hydrogen bond of Thr474beta(MADH)/Thr172beta(AADH) to the catalytic carboxylate group and the nonconserved active site residue Tyr471beta(MADH)/Phe169beta(AADH) are identified as important factors in determining the preferred oxygen acceptor. The protein environment has a significant effect on the reaction energetics and hence on tunnelling contributions and KIEs. These environmental effects, and the related clearly different preferences for the two carboxylate oxygen atoms (with different KIEs) in MADH/MA and AADH/tryptamine, are possible causes of the differences observed in the KIEs between these two important enzyme reactions.     

A theoretical study on the catalytic mechanism of Mus musculus adenosine deaminase

The catalytic mechanism of Mus musculus adenosine deaminase (ADA) has been studied by quantum mechanics and two‐layered ONIOM calculations. Our calculations show that the previously proposed mechanism, involving His238 as the general base to activate the Zn‐bound water, has a high activation barrier of about 28 kcal/mol at the proposed rate‐determining nucleophilic addition step, and the corresponding calculated kinetic isotope effects are significantly different from the recent experimental observations. We propose a revised mechanism based on calculations, in which Glu217 serves as the general base to abstract the proton of the Zn‐bound water, and the protonated Glu217 then activates the substrate for the subsequent nucleophilic addition. The rate‐determining step is the proton transfer from Zn‐OH to 6‐NH₂ of the tetrahedral intermediate, in which His238 serves as a proton shuttle for the proton transfer. The calculated kinetic isotope effects agree well with the experimental data, and calculated activation energy is also consistent with the experimental reaction rate. © 2010 Wiley Periodicals, Inc. J Comput Chem, 2010     

Neutral hydrolyses of carbon disulfide: An ab initio study of water catalysis

The water-catalyzed hydrolysis reaction of carbon disulfide (CS(2)) has been investigated at the levels of HF and MP2 with the basis set of 6-311++G(d,p) using the combined supramolecular/continuum models, in which up to six water molecules are involved in the hydrolysis and the effect of water bulk solvent is taken into account according to the polarizable continuum model (PCM). The activation Gibbs free energies in water solution, DeltaG(sol) (not equal) (298 K), for the rate-determining steps of one up to six water hydrolyses are 247.9, 184.2, 152.3, 141.8, 134.4, and 118.9 kJ/mol, respectively. The most favorable hydrolysis path of CS(2) involves a sort of eight-membered ring transition structure formed by six water molecules, among which three water molecules are not involved in the proton transfer, two near to the nonreactive sulfur atom, and one below the parent carbon disulfide. This suggests that the hydrolysis of CS(2) can be mediated with the water molecule(s) and be significantly facilitated by the cooperative effects of the water molecule(s) in the nonreactive region. The catalytic effects of water molecule(s) due to the alleviation of ring strain in the proton transfer process may result from the synergistic effects of rehybridization and charge reorganization from the prereaction complex to the rate-determining transition state structure induced by water molecule(s). PCM solvation models could significantly lower the rate-determining activation Gibbs free energies by 20-38 kJ/mol when two up to six explicit water molecules involved in the neutral hydrolysis of CS(2).     

The role of the putative catalytic base in the phosphoryl transfer reaction in a protein kinase: first-principles calculations

Protein kinases are important enzymes controlling the majority of cellular signaling events via a transfer of the gamma-phosphate of ATP to a target protein. Even after many years of study, the mechanism of this reaction is still poorly understood. Among many factors that may be responsible for the 1011-fold rate enhancement due to this enzyme, the role of the conserved aspartate (Asp166) has been given special consideration. While the essential presence of Asp166 has been established by mutational studies, its function is still debated. The general base catalyst role assigned to Asp166 on the basis of its position in the active site has been brought into question by the pH dependence of the reaction rate, isotope measurements, and pre-steady-state kinetics. Recent semiempirical calculations have added to the controversy surrounding the role of Asp166 in the catalytic mechanism. No major role for Asp166 has been found in these calculations, which have predicted the reaction process consisting of an early transfer of a substrate proton onto the phosphate group. These conclusions were inconsistent with experimental observations. To address these differences between experimental results and theory with a more reliable computational approach and to provide a theoretical platform for understanding catalysis in this important enzyme family, we have carried out first-principles structural and dynamical calculations of the reaction process in cAPK kinase. To preserve the essential features of the reaction, representations of all of the key conserved residues (82 atoms) were included in the calculation. The structural calculations were performed using the local basis density functional (DFT) approach with both hybrid B3LYP and PBE96 generalized gradient approximations. This kind of calculation has been shown to yield highly accurate structural information for a large number of systems. The optimized reactant state structure is in good agreement with X-ray data. In contrast to semiempirical methods, the lowest energy product state places the substrate proton on Asp166. First-principles molecular dynamics simulations provide additional support for the stability of this product state. The latter also demonstrate that the proton transfer to Asp166 occurs at a point in the reaction where bond cleavage at the PO bridging position is already advanced. This mechanism is further supported by the calculated structure of the transition state in which the substrate hydroxyl group is largely intact. A metaphoshate-like structure is present in the transition state, which is consistent with the X-ray structures of transition state mimics. On the basis of the calculated structure of the transition state, it is estimated to be 85% dissociative. Our analysis also indicates an increase in the hydrogen bond strength between Asp166 and substrate hydroxyl and a small decrease in the bond strength of the latter in the transition state. In summary, our calculations demonstrate the importance of Asp166 in the enzymatic mechanism as a proton acceptor. However, the proton abstraction from the substrate occurs late in the reaction process. Thus, in the catalytic mechanism of cAPK protein kinase, Asp166 plays a role of a "proton trap" that locks the transferred phosphoryl group to the substrate. These results resolve prior inconsistencies between theory and experiment and bring new understanding of the role of Asp166 in the protein kinase catalytic mechanism.     

Theoretical study of general base-catalyzed hydrolysis of aryl esters and implications for enzymatic reactions

In this work, the mechanism of general base-catalyzed hydrolysis of aryl esters is investigated in vacuo with density functional theory and in solutions with a polarized continuum model. The hydrolysis is found to proceed via a concerted mechanism featuring simultaneous addition and elimination steps accompanied by proton transfers, consistent with experimental evidence. Reasonable agreement with measured kinetic isotope effects provides additional validation. It is found that solvation substantially lowers the transition state energy, but has a small effect on the reaction exothermicity. An enzyme oxyanion hole, modeled by an ammonia molecule hydrogen bonded to the acyl carbonyl oxygen, is found to stabilize the near-tetrahedral transition state. Implications of these findings for the hydrolysis step of the dehalogenation reaction catalyzed by 4-chlorobenzoyl-CoA dehalogenase are discussed.     

Probing the origin of the compromised catalysis of E. coli alkaline phosphatase in its promiscuous sulfatase reaction

The catalytic promiscuity of E. coli alkaline phosphatase (AP) and many other enzymes provides a unique opportunity to dissect the origin of enzymatic rate enhancements via a comparative approach. Here, we use kinetic isotope effects (KIEs) to explore the origin of the 109-fold greater catalytic proficiency by AP for phosphate monoester hydrolysis relative to sulfate monoester hydrolysis. The primary 18O KIEs for the leaving group oxygen atoms in the AP-catalyzed hydrolysis of p-nitrophenyl phosphate (pNPP) and p-nitrophenylsulfate (pNPS) decrease relative to the values observed for nonenzymatic hydrolysis reactions. Prior linear free energy relationship results suggest that the transition states for AP-catalyzed reactions of phosphate and sulfate esters are "loose" and indistinguishable from that in solution, suggesting that the decreased primary KIEs do not reflect a change in the nature of the transition state but rather a strong interaction of the leaving group oxygen atom with an active site Zn2+ ion. Furthermore, the primary KIEs for the two reactions are identical within error, suggesting that the differential catalysis of these reactions cannot be attributed to differential stabilization of the leaving group. In contrast, AP perturbs the KIE for the nonbridging oxygen atoms in the reaction of pNPP but not pNPS, suggesting a differential interaction with the transferred group in the transition state. These and prior results are consistent with a strong electrostatic interaction between the active site bimetallo Zn2+ cluster and one of the nonbridging oxygen atoms on the transferred group. We suggest that the lower charge density of this oxygen atom on a transferred sulfuryl group accounts for a large fraction of the decreased stabilization of the transition state for its reaction relative to phosphoryl transfer.     

Reaction of imidazole with toluene-4-sulfonate salts of substituted phenyl N-methylpyridinium-4-carboxylate esters: special base catalysis by imidazole

The reaction of imidazole in aqueous solution with toluene-4-sulfonate salts of substituted phenyl N-methylpyridinium-4-carboxylate esters obeys the rate law: k(obs) - k(background) = k2[Im] + k3[Im]2 where [Im] is the imidazole concentration present as free base. The parameters k2 and k3 fit Br?nsted type free energy correlations against the pKa of the leaving phenol with betaLg values of -0.65 and -0.42 respectively. The imidazolysis is insensitive to catalysis by general bases and yet k3 for the 3-cyanophenyl ester possesses a deuterium oxide solvent isotope effect of 4.43 consistent with rate limiting proton transfer. A special catalytic function is proposed for decomposition of the tetrahedral addition intermediate (T+/-) via k3 whereby the catalytic imidazole interacts electrophilically with the leaving phenolate ion and removes a proton from the nitrogen in the rate limiting step with subsequent non-rate limiting ArO-C bond fission. This is consistent with the change in effective charge on the leaving oxygen in the transition structure of k3 which is more positive (-0.42) than that expected (-0.60) for the equilibrium formation of the zwitterion intermediate. The catalytic function at the leaving oxygen is likely to be an electrophilic role of the NH as a hydrogen bond donor. In the k2 step the deuterium oxide solvent isotope effect of 1.51 for the 3-cyanophenyl ester and the betaLg of -0.65 are consistent with rate limiting expulsion of the phenolate ion from the T+/- intermediate. The absence of general base catalysis of imidazolysis rules out the established mechanism for aminolysis of esters where T+/- is stabilised by a standard rate limiting proton transfer. The kinetically equivalent term for k3 where T- reacts with the imidazolium ion as an acid catalyst would require this step to be rate limiting and involve proton transfer not consistent with departure of the good aryl oxide leaving group.     

Tunneling and classical paths for proton transfer in an enzyme reaction dominated by tunneling: oxidation of tryptamine by aromatic amine dehydrogenase

Proton tunneling dominates the oxidative deamination of tryptamine catalyzed by the enzyme aromatic amine dehydrogenase. For reaction with the fast substrate tryptamine, a H/D kinetic isotope effect (KIE) of 55 +/- 6 has been reported-one of the largest observed in an enzyme reaction. We present here a computational analysis of this proton-transfer reaction, applying combined quantum mechanics/molecular mechanics (QM/MM) methods (PM3-SRP//PM3/CHARMM22). In particular, we extend our previous computational study (Masgrau et al. Science 2006, 312, 237) by using improved energy corrections, high-level QM/MM methods, and an ensemble of paths to estimate the tunneling contributions. We have carried out QM/MM molecular dynamics simulations and variational transition state theory calculations with small-curvature tunneling corrections. The results provide detailed insight into the processes involved in the reaction. Transfer to the O2 oxygen of the catalytic base, Asp128beta, is found to be the favored reaction both thermodynamically and kinetically, even though O1 is closer in the reactant complex. Comparison of quantum and classical models of proton transfer allows estimation of the contribution of hydrogen tunneling in lowering the barrier to reaction in the enzyme. A reduction of the activation free energy due to tunneling of 3.1 kcal mol-1 is found, which represents a rate enhancement due to tunneling by 2 orders of magnitude. The calculated KIE of 30 is significantly elevated over the semiclassical limit, in agreement with the experimental observations; a semiclassical value of 6 is obtained when tunneling is omitted. A polarization of the C-H bond to be broken is observed due to the close proximity of the catalytic aspartate and the (formally) positively charged imine nitrogen. A comparison is also made with the related quinoprotein methylamine dehydrogenase (MADH)-the much lower KIE of 11 that we obtain for the MADH/methylamine system is found to arise from a more endothermic potential energy surface for the MADH reaction.     

Significant kinetic solvent isotope effects in folding of the catalytic RNA from the hepatitis delta virus

The exchange of deuterium for hydrogen in water often produces solvent kinetic isotope effects (KSIEs) on the rate constants associated with enzyme reactions, including those catalyzed by RNA. Recently, KSIEs have been used to show that proton transfer occurs in the rate-limiting step of cleavage by the hepatitis delta virus (HDV) ribozyme and other catalytic RNAs. To test the underlying assumption that KSIEs are related to the chemistry step of ribozyme-mediated cleavage reactions, we developed fluorescence resonance energy transfer assays to measure KSIEs on the rate constants of conformational changes associated with substrate binding and dissociation by a trans-acting HDV ribozyme. We observe comparable KSIEs ( approximately 2-2.5-fold) of rate constants of conformational change and cleavage, while proton inventory experiments are consistent with a shift in the ensemble of transition states upon increase of D2O in the solvent. Taken together, these results challenge the common assumption that pL profiles of RNA-catalyzed reactions yielding a pKa and KSIE necessarily provide evidence for an ionization (chemistry) step to be rate-limiting. They also suggest that an unusual proton inventory may provide a signature for a conformational change contributing to the rate-limiting step.     

The Mechanism of the Acid-catalysed Hydrolysis of 1-Aryl-2,2,2-trifluorodiazoethanes

The rates of the acid-catalysed hydrolysis of a series of 1-aryl-2,2,2-trifluorodiazoethanes la-d have been measured in dioxan/water/HClO₄. The results are well correlated with the Hammett equation when σ_p substituent constants are used (?_H = ?1.74 and ?_D = ?1.75). Kinetic solvent isotope effects, about 2.1, and general acid catalysis indicate that proton transfer is rate-determining (A-S_E2 mechanism). Rate measurements have also been made at pressures up to 1400 atm. The evaluated activation volumes, about ?20 cm³/mole, indicate that at least one water molecule is bound in the transition state of protonation. Rate measurements at low water concentrations indicate that no apparent change in mechanism has occured.