Catalysis by entropic effects: the action of cytidine deaminase on 5,6-dihydrocytidine

In neutral solution, 5,6-dihydrocytidine undergoes spontaneous deamination (k25 approximately 3.2 x 10(-5) s(-1)) much more rapidly than does cytidine (k25 approximately 3.0 x 10(-10) s(-1)), with a more favorable enthalpy of activation (DeltaDeltaH# = -8.7 kcal/mol) compensated by a less favorable entropy of activation (TDeltaDeltaS# = -1.8 kcal/mol at 25 degrees C). E. coli cytidine deaminase enhances the rate of deamination of 5,6-dihydrocytidine (kcat/k(non) = 4.4 x 10(5)) by enhancing the entropy of activation (DeltaDeltaH# = 0 kcal/mol; TDeltaDeltaS# = +7.6 kcal/mol, at 25 degrees C). Binding of the competitive inhibitor 3,4,5,6-tetrahydrouridine (THU), a stable analogue of 5,6-dihydrocytidine in the transition state for its deamination, is accompanied by a release of enthalpy (DeltaH = -7.1 kcal/mol, TDeltaDeltaS = +2.2 kcal/mol) that approaches the estimated enthalpy of binding of the actual substrate in the transition state for deamination of 5,6-dihydrocytidine (DeltaH = -8.1 kcal/mol, TDeltaDeltaS = +6.0 kcal/mol). Thus, the shortcomings of THU in capturing all of the binding affinity expected of an ideal transition-state analogue reflect a less favorable entropy of association. That difference may arise from the analogue's inability to displace a water molecule from the "leaving group site" at which ammonia is generated in the normal reaction. The effect on binding of removing the 4-OH group from the transition-state analogue THU, to form 3,4,5,6-tetrahydrozebularine (THZ) (DeltaDeltaH = -2.1 kcal/mol, TDeltaDeltaS = -4.4 kcal/mol), is mainly entropic, consistent with the inability of THZ to displace water from the "attacking group site". These results are consistent with earlier indications [Snider, M. J., and Wolfenden, R. (2001) Biochemistry 40, 11364] that site-bound water plays a prominent role in substrate activation and inhibitor binding by cytidine deaminase.     

Insight into the phosphodiesterase mechanism from combined QM/MM free energy simulations

Molecular dynamics simulations employing a combined quantum mechanical and molecular mechanical potential have been carried out to elucidate the reaction mechanism of the hydrolysis of a cyclic nucleotide cAMP substrate by phosphodiesterase 4B (PDE4B). PDE4B is a member of the PDE superfamily of enzymes that play crucial roles in cellular signal transduction. We have determined a two-dimensional potential of mean force (PMF) for the coupled phosphoryl bond cleavage and proton transfer through a general acid catalysis mechanism in PDE4B. The results indicate that the ring-opening process takes place through an S(N)2 reaction mechanism, followed by a proton transfer to stabilize the leaving group. The computed free energy of activation for the PDE4B-catalyzed cAMP hydrolysis is about 13 kcal·mol(-1) and an overall reaction free energy is about -17 kcal·mol(-1), both in accord with experimental results. In comparison with the uncatalyzed reaction in water, the enzyme PDE4B provides a strong stabilization of the transition state, lowering the free energy barrier by 14 kcal·mol(-1). We found that the proton transfer from the general acid residue His234 to the O3' oxyanion of the ribosyl leaving group lags behind the nucleophilic attack, resulting in a shallow minimum on the free energy surface. A key contributing factor to transition state stabilization is the elongation of the distance between the divalent metal ions Zn(2+) and Mg(2+) in the active site as the reaction proceeds from the Michaelis complex to the transition state.     

Powering DNA repair through substrate electrostatic interactions

The reaction catalyzed by the DNA repair enzyme uracil DNA glycosylase (UDG) proceeds through an unprecedented stepwise mechanism involving a positively charged oxacarbenium ion sugar and uracil anion leaving group. Here we use a novel approach to evaluate the catalytic contribution of electrostatic interactions between four essential phosphodiester groups of the DNA substrate and the cationic transition state. Our strategy was to substitute each of these phosphate groups with an uncharged (R)- or (S)-methylphosphonate linkage (MeP). We then compared the damaging effects of these methylphosphonate substitutions on catalysis with their damaging effects on binding of a cationic 1-azadeoxyribose (1-aza-dR(+)) oxacarbenium ion analogue to the UDG-uracil anion binary complex. A plot of log k(cat)/K(m) for the series of MeP-substituted substrates against log K(D) for binding of the 1-aza-dR(+) inhibitors gives a linear correlation of unit slope, confirming that the electronic features of the transition state resemble that of the 1-aza-dR(+), and that the anionic backbone of DNA is used in transition state stabilization. We estimate that all of the combined phosphodiester interactions with the substrate contribute 6-8 kcal/mol toward lowering the activation barrier, a stabilization that is significant compared to the 16 kcal/mol catalytic power of UDG. However, unlike groups of the enzyme that selectively stabilize the charged transition state by an estimated 7 kcal/mol, these phosphodiester groups also interact strongly in the ground state. To our knowledge, these results provide the first experimental evidence for electrostatic stabilization of a charged enzymatic transition state and intermediate using the anionic backbone of DNA.     

Role of enzyme-ribofuranosyl contacts in the ground state and transition state for orotidine 5'-phosphate decarboxylase: a role for substrate destabilization?

The crystal structure of yeast orotidine 5'-monophosphate decarboxylase (ODCase) complexed with the inhibitor 6-hydroxyuridine 5'-phosphate (BMP) reveals the presence of a series of strong interactions between enzyme residues and functional groups of this ligand. Enzyme contacts with the phosphoribofuranosyl moiety of orotidine 5'-phosphate (OMP) have been shown to contribute at least 16.6 kcal/mol of intrinsic binding free energy to the stabilization of the transition state for the reaction catalyzed by yeast ODCase. In addition to these enzyme-ligand contacts, active site residues contributed by both subunits of the dimeric enzyme are positioned to form hydrogen bonds with the 2'- and 3'-OH groups of the ligand's ribosyl moiety. These involve Thr-100 of one subunit and Asp-37 of the opposite subunit, respectively. To evaluate the contributions of these ribofuranosyl contacts to ground state and transition state stabilization, Thr-100 and Asp-37 were each mutated to alanine. Elimination of the enzyme's capacity to contact individual ribosyl OH groups reduced the k(cat)/K(m) value of the T100A enzyme by 60-fold and that of the D37A enzyme by 300-fold. Removal of the 2'-OH group from the substrate OMP decreased the binding affinity by less than a factor of 10, but decreased k(cat) by more that 2 orders of magnitude. Upon removal of the complementary hydroxymethyl group from the enzyme, little further reduction in k(cat)/K(m) for 2'-deoxyOMP was observed. To assess the contribution made by contacts involving both ribosyl hydroxyl groups at once, the ability of the D37A mutant enzyme to decarboxylate 2'-deoxyOMP was measured. The value of k(cat)/K(m) for this enzyme-substrate pair was 170 M(-1) s(-1), representing a decrease of more than 7.6 kcal/mol of binding free energy in the transition state. To the extent that electrostatic repulsion in the ground state can be tested by these simple alterations, the results do not lend obvious support to the view that electrostatic destabilization in the ground state enzyme-substrate complex plays a major role in catalysis.     

Multiple transients in the pre-steady-state of nucleoside hydrolase reveal complex substrate binding, product base release, and two apparent rates of chemistry

We have investigated the transient kinetics of the nucleoside hydrolase from Trypanosoma vivax (TvNH) at low temperatures (5 degrees C). Three novel absorbance transients (termed tau1, tau3, and tau4) were detected during multiple-guanosine turnover stopped-flow absorbance spectroscopy, in addition to a transient (tau2) that had been observed previously at 35 degrees C. At 5 degrees C, TvNH displays full-sites activity and not half-of-the-sites activity as is apparent at 35 degrees C. Both tau1 and tau2 are assigned to chemistry based on rapid-quench results. For tau1, the rate of chemistry is ca. 3000-fold faster than kcat (1-2 orders of magnitude greater than previous estimates). The pH dependencies of the burst amplitudes for tau1 and tau2 indicate that these transients arise from the formation of two different dimeric TvNH.substrate complexes and not from TvNH that contains kinetically asymmetric monomers. The saturating burst rates for tau1 and tau2 are surprisingly pH-independent, given the prominent role of acid-base chemistry in the proposed mechanism for TvNH. tau3 and tau4 are assigned to the substrate binding and base release processes, respectively, and shown to be equivalent to two fluorescence transients (tau3 and tau4, respectively) observed previously by stopped-flow methods at 35 degrees C. The rate of base release is shown to be an apparent rate. Together with steady-state product inhibition results, the data indicate that TvNH follows an ordered uni-bi kinetic mechanism with a TvNH.base dead-end complex, and not the rapid equilibrium random uni-bi mechanism proposed for other NHs. Two applicable kinetic models are presented and their implications for future mechanistic studies discussed.     

Enzyme-substrate interactions in the purine-specific nucleoside hydrolase from Trypanosoma vivax

Nucleoside hydrolases are key enzymes in the purine salvage pathway of Trypanosomatidae and are considered as targets for drug design. We previously reported the first x-ray structure of an inosine-adenosine-guanosine preferring nucleoside hydrolase (IAG-NH) from Trypanosoma vivax (). Here we report the 2.0-A crystal structure of the slow D10A mutant in complex with the inhibitor 3-deaza-adenosine and the 1.6-A crystal structure of the same enzyme in complex with a genuine substrate inosine. The enzyme-substrate complex shows the substrate bound to the enzyme in a different conformation from 3-deaza-adenosine and provides a snapshot along the reaction coordinate of the enzyme-catalyzed reaction. The chemical groups on the substrate important for binding and catalysis are mapped. The 2'-OH, 3'-OH, and 5'-OH contribute 4.6, 7.5, and 5.4 kcal/mol to k(cat)/K(m), respectively. Specific interactions with the exocyclic groups on the purine ring are not required for catalysis. Site-directed mutagenesis indicates that the purine specificity of the IAG-NHs is imposed by a parallel aromatic stacking interaction involving Trp(83) and Trp(260). The pH profiles of k(cat) and k(cat)/K(m) indicate the existence of one or more proton donors, possibly involved in leaving group activation. However, mutagenesis of the active site residues around the nucleoside base and an alanine scan of a flexible loop near the active site fail to identify this general acid. The parallel aromatic stacking seems to provide the most likely alternative mechanism for leaving group activation.     

Temperature dependence of binding and catalysis for the Cdc25B phosphatase

Using a combination of steady-state and single-turnover kinetics, we probe the temperature dependence of substrate association and chemistry for the reaction of Cdc25B phosphatase with its Cdk2-pTpY/CycA protein substrate. The transition state for substrate association is dominated by an enthalpic barrier (DeltaH(++) of 13 kcal/mol) and has a favorable entropic contribution of 4 kcal/mol at 298 K. Phosphate transfer from Cdk2-pTpY/CycA to enzyme (DeltaH(++) of 12 kcal/mol) is enthalpically more favorable than for the small molecule substrate p-nitrophenyl phosphate (DeltaH(++) of 18 kcal/mol), yet entropically less favorable (TDeltaS(++) of 2 vs. -6 kcal/mol at 298 K, respectively). By measuring the temperature dependence of binding and catalysis for several hotspot mutants involved in binding of protein substrate, we determine the enthalpy-entropy compensations for changes in rates of association and phosphate transfer compared to the wild type system. We conclude that the transition state for enzyme-substrate association involves tight and specific contacts at the remote docking site and that phospho-transfer from Cdk2-pTpY/CycA to the pre-organized active site of the enzyme is accompanied by unfavorable entropic rearrangements that promote rapid product dissociation.     

Initial step of B12-dependent enzymatic catalysis: energetic implications regarding involvement of the one-electron-reduced form of adenosylcobalamin cofactor

Density functional theory analysis was performed to elucidate the impact of one-electron reduction upon the initial step of adenosylcobalamin-dependent enzymatic catalysis. The transition state (TS) corresponding to the Co–C bond cleavage and subsequent hydrogen abstraction from the substrate was located. The intrinsic reaction coordinate calculations predicted that the reaction consisting of Co–C5′ bond cleavage in [Co^III(corrin^•)]–Rib (where Rib is ribosyl) and hydrogen-atom abstraction from the CH₃–CH₂–CHO substrate occurs in a concerted fashion. The computed activation energy barrier of the reaction (15.0 kcal/mol) was lowered by approximately 54.5% in comparison with the reaction involving the positively charged cofactor model (Im–[Co^III(corrin)]–Rib⁺, where Im is imidazole; energy barrier = 33.0 kcal/mol). The Im base was detached during the TS search in the reaction involving the one-electron-reduced analogue. Thus, to compare the energetics of the two reactions, the axial Im ligand detachment energy for the Im–[Co^III(corrin^•)]–Rib model was computed [7.6 kcal/mol (gas phase); 4.6 kcal/mol (water)]. Consequently, the effective activation energy barrier for the reaction mediated by the Im-off [Co^III(corrin^•)]–Rib was estimated to be 22.6 kcal/mol, which implied an overall 31.5% reduction in the energetic demands of the reaction. Considering that the lengthened Co–N_axial bond has been observed in X-ray crystal structure studies of B₁₂-dependent mutases, the catalytic impact induced by one-electron reduction of the cofactor is expected to be higher in the presence of the enzymatic environment.     

Substrate-assisted leaving group activation in enzyme-catalyzed N-glycosidic bond cleavage

In enzymatic depurination of nucleosides, the 5'-OH group of the ribose moiety of the substrate is often shown to contribute substantially to catalysis. The purine-specific nucleoside hydrolase from Trypanosoma vivax (TvNH) fixes the 5'-OH group in a gauche,trans orientation about the C4'-C5' bond, enabling the 5'-oxygen to accept an intramolecular hydrogen bond from the C8-atom of the purine leaving group. High level ab initio quantum chemical calculations indicate that this interaction promotes protonation of the purine at N7. Steady state kinetics comprising engineered substrates confirm that a considerable fraction of the catalytic 5'-OH effect can be attributed to leaving group activation.     

Temperature effects on the catalytic efficiency, rate enhancement, and transition state affinity of cytidine deaminase, and the thermodynamic consequences for catalysis of removing a substrate "anchor"

To obtain a clearer understanding of the forces involved in transition state stabilization by Escherichia coli cytidine deaminase, we investigated the thermodynamic changes that accompany substrate binding in the ground state and transition state for substrate hydrolysis. Viscosity studies indicate that the action of cytidine deaminase is not diffusion-limited. Thus, K(m) appears to be a true dissociation constant, and k(cat) describes the chemical reaction of the ES complex, not product release. Enzyme-substrate association is accompanied by a loss of entropy and a somewhat greater release of enthalpy. As the ES complex proceeds to the transition state (ES), there is little further change in entropy, but heat is taken up that almost matches the heat that was released with ES formation. As a result, k(cat)/K(m) (describing the overall conversion of the free substrate to ES is almost invariant with changing temperature. The free energy barrier for the enzyme-catalyzed reaction (k(cat)/K(m)) is much lower than that for the spontaneous reaction (k(non)) (DeltaDeltaG = -21.8 kcal/mol at 25 degrees C). This difference, which also describes the virtual binding affinity of the enzyme for the activated substrate in the transition state (S), is almost entirely enthalpic in origin (DeltaDeltaH = -20.2 kcal/mol), compatible with the formation of hydrogen bonds that stabilize the ES complex. Thus, the transition state affinity of cytidine deaminase increases rapidly with decreasing temperature. When a hydrogen bond between Glu-91 and the 3'-hydroxyl moiety of cytidine is disrupted by truncation of either group, k(cat)/K(m) and transition state affinity are each reduced by a factor of 10(4). This effect of mutation is entirely enthalpic in origin (DeltaDeltaH approximately 7.9 kcal/mol), somewhat offset by a favorable change in the entropy of transition state binding. This increase in entropy is attributed to a loss of constraints on the relative motions of the activated substrate within the ES complex. In an Appendix, some objections to the conventional scheme for transition state binding are discussed.     

Bovine skeletal muscle adenosine deaminase: kinetic and thermodynamic studies

Adenosine deaminase from bovine skeletal muscle catalyzes the hydrolytic deamination of adenosine to inosine and ammonia via an ordered Uni-Bi mechanism, if water is not considered as a true second substrate, as deduced from the inhibition pattern products. The inhibition constants (Ki) obtained for inosine and ammonia were 316 mumol/l and 2 mol/l, respectively. The activation energy of the reaction has been calculated as 10 kcal/mol, delta H* and delta F* as 7.9 and 15.6 kcal/mol, respectively, and delta S* as -23 cal/mol/degrees K.     

Site-bound water and the shortcomings of a less than perfect transition state analogue

Kinetic measurements have shown that substantial enthalpy changes accompany substrate binding by cytidine deaminase, increasing markedly as the reaction proceeds from the ground state (1/K(m), DeltaH = -13 kcal/mol) to the transition state (1/K(tx), DeltaH = -20 kcal/mol) [Snider, M. J., et al. (2000) Biochemistry 39, 9746-9753]. In the present work, we determined the thermodynamic changes associated with the equilibrium binding of inhibitors by cytidine deaminase by isothermal titration calorimetry and van't Hoff analysis of the temperature dependence of their inhibition constants. The results indicate that the binding of the transition state analogue 3,4-dihydrouridine DeltaH = -21 kcal/mol), like that of the transition state itself (DeltaH = -20 kcal/mol), is associated with a large favorable change in enthalpy. The significantly smaller enthalpy change that accompanies the binding of 3,4-dihydrozebularine (DeltaH = -10 kcal/mol), an analogue of 3,4-dihydrouridine in which a hydrogen atom replaces this inhibitor's 4-OH group, is consistent with the view that polar interactions with the substrate at the site of its chemical transformation play a critical role in reducing the enthalpy of activation for substrate hydrolysis. The entropic shortcomings of 3,4-dihydrouridine, in capturing all of the free energy involved in binding the actual transition state, may arise from its inability to displace a water molecule that occupies the binding site normally occupied by product ammonia.     

Thermodynamic studies of the reversible association of Escherichia coli ribosomal subunits.

The kinetics of association of Escherichia coli 30S and 50S ribosomal subunits have been carried out as a function of temperature after a magnesium jump from 1.5 to 3 mM. Turbidimetric recordings combined with a stopped-flow apparatus were used to follow the kinetics. The data show that the rates of formation and dissociation of the 70S particles at 3 mM Mg2+ and +25 degrees C were, respectively: k2 = 10(5) M-1 s-1, k1 = 4,5 X 10(-3) s-1; lowering the temperature decreases the rate constants with activation energies equal to E2 = 7.5 kcal/mol, E1 = 26.5 kcal/mol and enhances the association equilibrium towards the 70S species with an enthalpy change (delta H degrees assoc = -19.9 kcal/mol) dominant over the entropy change (delta S degrees assoc = -33 cal/(deg mol)). These thermodynamic parameters were compared to those obtained from studies on the interactions of codon-anticodon in yeast phenylalanine transfer RNA as well as of ribooligonucleotides. The kinetic and thermodynamic data are shown to be consistent with 16S-23S RNA interaction.     

Kinetics and thermodynamics of mandelate racemase catalysis

Mandelate racemase (EC 5.1.2.2) from Pseudomonas putida catalyzes the interconversion of the two enantiomers of mandelic acid with remarkable proficiency, producing a rate enhancement exceeding 15 orders of magnitude. The rates of the forward and reverse reactions catalyzed by the wild-type enzyme and by a sluggish mutant (N197A) have been studied in the absence and presence of several viscosogenic agents. A partial dependence on relative solvent viscosity was observed for values of kcat and kcat/Km for the wild-type enzyme in sucrose-containing solutions. The value of kcat for the sluggish mutant was unaffected by varying solvent viscosity. However, sucrose did have a slight activating effect on mutant enzyme efficiency. In the presence of the polymeric viscosogens poly(ethylene glycol) and Ficoll, no effect on kcat or kcat/Km for the wild-type enzyme was observed. These results are consistent with both substrate binding and product dissociation being partially rate-determining in both directions. The viscosity variation method was used to estimate the rate constants comprising the steady-state expressions for kcat and kcat/Km. The rate constant for the conversion of bound (R)-mandelate to bound (S)-mandelate (k2) was found to be 889 +/- 40 s(-1) compared with a value of 654 +/- 58 s(-1) for kcat in the same direction. From the temperature dependence of Km (shown to equal K(S)), k2, and the rate constant for the uncatalyzed reaction [Bearne, S. L., and Wolfenden, R. (1997) Biochemistry 36, 1646-1656], we estimated the enthalpic and entropic changes associated with substrate binding (DeltaH = -8.9 +/- 0.8 kcal/mol, TDeltaS = -4.8 +/- 0.8 kcal/mol), the activation barrier for conversion of bound substrate to bound product (DeltaH# = +15.4 +/- 0.4 kcal/mol, TDeltaS# = +2.0 +/- 0.1 kcal/mol), and transition state stabilization (DeltaH(tx) = -22.9 +/- 0.8 kcal/mol, TDeltaS(tx) = +1.8 +/- 0.8 kcal/mol) during mandelate racemase-catalyzed racemization of (R)-mandelate at 25 degrees C. Although the high proficiency of mandelate racemase is achieved principally by enthalpic reduction, there is also a favorable and significant entropic contribution.     

Thermodynamic analysis of catalysis by the dihydroorotases from hamster and Bacillus caldolyticus, as compared with the uncatalyzed reaction

Dihydroorotase (DHOase, EC 3.5.2.3) from the extreme thermophile Bacillus caldolyticus has been subcloned, sequenced, expressed, and purified as a monomer. The catalytic properties of this thermophilic DHOase have been compared with another type I enzyme, the DHOase domain from hamster, to investigate how the thermophilic enzyme is adapted to higher temperatures. B. caldolyticus DHOase has higher Vmax and Ks values than hamster DHOase at the same temperature. The thermodynamic parameters for the binding of L-dihydroorotate were determined at 25 degrees C for hamster DHOase (deltaG = -6.9 kcal/mol, deltaH = -11.5 kcal/mol, TdeltaS = -4.6 kcal/mol) and B. caldolyticus DHOase (deltaG = -5.6 kcal/mol, deltaH = -4.2 kcal/mol, TdeltaS = +1.4 kcal/mol). The smaller enthalpy release and positive entropy for thermophilic DHOase are indicative of a weakly interacting Michaelis complex. Hamster DHOase has an enthalpy of activation of 12.3 kcal/mol, similar to the release of enthalpy upon substrate binding, rendering the kcat/Ks value almost temperature independent. B. caldolyticus DHOase shows a decrease in the enthalpy of activation from 12.2 kcal/mol at temperatures from 30 to 50 degrees C to 5.3 kcal/mol for temperatures of 50-70 degrees C. Vibrational energy at higher temperatures may facilitate the transition ES --> ES(double dagger), making kcat/Ks almost temperature independent. The pseudo-first-order rate constant for water attack on L-dihydroorotate, based on experiments at elevated temperature, is 3.2 x 10(-11) s(-1) at 25 degrees C, with deltaH(double dagger) = 24.7 kcal/mol and TdeltaS(double dagger) = -6.9 kcal/mol. Thus, hamster DHOase enhances the rate of substrate hydrolysis by a factor of 1.6 x 10(14), achieving this rate enhancement almost entirely by lowering the enthalpy of activation (delta deltaH(double dagger) = -19.5 kcal/mol). Both the rate enhancement and transition state affinity of hamster DHOase increase steeply with decreasing temperature, consistent with the development of H-bonds and electrostatic interactions in the transition state that were not present in the enzyme-substrate complex in the ground state.     

Uracil DNA glycosylase: revisiting substrate-assisted catalysis by DNA phosphate anions

Uracil DNA glycosylase (UNG) is a powerful DNA repair enzyme that has been shown to stabilize a glycosyl cation reaction intermediate and a related tight binding inhibitor using electrostatic interactions with the +1 and -1, but not the +2, phosphodiester group of the single-stranded DNA substrate Ap (2+)Ap (1+)Up (1-)ApA. These experimental results differed considerably from computational findings using duplex DNA, where the +2 phosphate was found to stabilize the transition state by approximately 5 kcal/mol, suggesting that UNG uses different catalytic strategies with single-stranded and double-stranded DNA substrates. In addition, the computational studies indicated that the conserved and positively charged His148 (which hydrogen bonds to the +2 phosphate) destabilized the glycosyl cation intermediate by 6-8 kcal/mol through anticatalytic electrostatic interactions. To evaluate these interesting proposals, we measured the kinetic effects of neutral methylphosphonate (MeP) stereoisomers at the +1 and +2 positions of a 12-mer dsDNA substrate and also the catalytic contribution and ionization state of His148. For MeP substitutions at the +1 position, single-turnover kinetic studies showed that the activation barrier was increased by 9.8 and 3.1 kcal/mol, corresponding to a stereoselectivity of nearly 40000-fold for the respective MeP isomers. Identical to the findings with ssDNA, MeP substitutions at the +2 position resulted in only small changes in the activation barrier (+/-0.3 kcal/mol), with little stereoselectivity ( approximately 4-fold). However, the H148A mutation destabilizes both the ground state and transition states by 2.4 and 4.3 kcal/mol, respectively. Thus, His148 is catalytic because it stabilizes the transition state to a greater extent (1.9 kcal/mol) than the ground state. Heteronuclear NMR studies established that His148 was neutral in the free enzyme at neutral pH, and in conformational exchange in a specific DNA complex containing uracil. We conclude that the +1 and +2 phosphate esters play identical catalytic roles in the reactions of single-stranded and double-stranded DNA substrates, and that His148 serves a catalytic role by positioning the substrate and catalytic water, or by an environmental effect.     

Influence of key residues on the reaction mechanism of the cAMP-dependent protein kinase

The reaction mechanism of the catalytic phosphoryl transfer of cAMP-dependent protein kinase (cAPK) was investigated by semi-empirical AM1 molecular orbital computations of an active site model system derived from the crystal structure of the catalytic subunit of the enzyme. The activation barrier is calculated as 20.7 kcal mol(-1) and the reaction itself to be exothermic by 12.2 kcal mol(-1). The active site residue Asp166, which was often proposed to act as a catalytic base, does not accept a proton in any of the reaction steps. Instead, the hydroxyl hydrogen of serine is shifted to the simultaneously transferred phosphate group of ATP. Although the calculated transition state geometry indicates an associative phosphoryl transfer, no concentration of negative charge is found. To study the influence of protein mutations on the reaction mechanism, we compared two-dimensional energy hypersurfaces of the protein kinase wild-type model and a corresponding mutant in which Asp166 was replaced by alanine. Surprisingly, they show similar energy profiles despite the experimentally known decrease of catalytic activity for corresponding mutants. Furthermore, a model structure was examined, where the charged NH3 group of Lys168 was replaced by a neutral methyl group. The energetic hypersurface of this hypothetical mutant shows two possible pathways for phosphoryl transfer, which both require significantly higher activation energies than the other systems investigated, while the energetic stabilization of the reaction product is similar in all systems. As the position of the amino acid side chains and the substrate peptide is virtually unchanged in all model systems, our results suggest that the exchange of Asp166 by other amino acid is less important to the phosphoryl transfer itself, but crucial to maintain the configuration of the active site in vivo. The positively charged side chain of Lys168, however, is necessary to stabilize the intermediate reaction states, particularly the side chain of the substrate peptide.     

Transfer RNA identity contributes to transition state stabilization during aminoacyl-tRNA synthesis.

Sequence-specific interactions between aminoacyl-tRNA synthetases and their cognate tRNAs ensure both accurate RNA recognition and the efficient catalysis of aminoacylation. The effects of tRNA(Trp)variants on the aminoacylation reaction catalyzed by wild-type Escherichia coli tryptophanyl-tRNA synthe-tase (TrpRS) have now been investigated by stopped-flow fluorimetry, which allowed a pre-steady-state analysis to be undertaken. This showed that tRNA(Trp)identity has some effect on the ability of tRNA to bind the reaction intermediate TrpRS-tryptophanyl-adenylate, but predominantly affects the rate at which trypto-phan is transferred from TrpRS-tryptophanyl adenylate to tRNA. Use of the binding ( K (tRNA)) and rate constants ( k (4)) to determine the energetic levels of the various species in the aminoacylation reaction showed a difference of approximately 2 kcal mol(-1)in the barrier to transition state formation compared to wild-type for both tRNA(Trp)A-->C73 and. These results directly show that tRNA identity contributes to the degree of complementarity to the transition state for tRNA charging in the active site of an aminoacyl-tRNA synthetase:aminoacyl-adenylate:tRNA complex.     

Probing the limits of electrostatic catalysis by uracil DNA glycosylase using transition state mimicry and mutagenesis

The DNA repair enzyme uracil DNA glycosylase (UDG) hydrolyzes the glycosidic bond of deoxyuridine in DNA by a remarkable mechanism involving formation of a positively charged oxacarbenium ion-uracil anion intermediate. We have proposed that the positively charged intermediate is stabilized by being sandwiched between the combined negative charges of the anionic uracil leaving group and a conserved aspartate residue that are located on opposite faces of the sugar ring. Here we establish that a duplex DNA oligonucleotide containing a cationic 1-aza-deoxyribose (I) oxacarbenium ion mimic is a potent inhibitor of UDG that binds tightly to the enzyme-uracil anion (EU(-)) product complex (K(D) of EU(-) = 110 pm). The tight binding of I to the EU(-) complex results from its extremely slow off rate (k(off) = 0.0008 s(-1)), which is 25,000-fold slower than substrate analogue DNA. Removal of Asp(64) and His(187), which are involved in stabilization of the cationic sugar and the anionic uracil leaving group, respectively, specifically weakens binding of I to the UDG-uracil complex by 154,000-fold, without significantly affecting substrate or product binding. These results suggest that electrostatic effects can effectively stabilize such an intermediate by at least -7 kcal/mol, without leading to anticatalytic stabilization of the substrate and products.     

Antagonistic substrate binding by a group II intron ribozyme.

In this study, the thermodynamic properties of substrate-ribozyme recognition were explored using a system derived from group II intron ai5gamma. Substrate recognition by group II intron ribozymes is of interest because any nucleic ac?id sequence can be targeted, the recognition sequence can be quite long (>/=13 bp), and reaction can proceed with a very high degree of sequence specificity. Group II introns target their substrates throug?h the formation of base-pairing interactions with two regions of the intron (EBS1 and EBS2), which are usually located far apart in the secondary structure. These structures pair with adjacent, corresponding sites (IBS1 and IBS2) on the substrate. In order to understand the relative energetic contribution of each base-pairing interaction (EBS1-IBS1 or EBS2-IBS2) to substrate binding energy, the free energy of each helix was measured. The individual helices were found to have base-pairing free energies similar to those calculated for regular RNA duplexes of the same sequence, suggesting that each recognition helix derives its binding energy from base-pairing interactions alone and that each helix can form independently. Most interestingly, it was found that the sum of the measured individual free energies (approximately 20 kcal/mol) was much higher than the known free energy for substrate binding (approximately 12 kcal/mol). This indicates that certain group II intron ribozymes can bind their substrates in an antagonistic fashion, paying a net energetic penalty upon binding the full-length substrate. This loss of binding energy is not due to weakening of individual helices, but appears to be linked to ribozyme conformational changes induced by substrate binding. This coupling between substrate binding and ribozyme conformational rearrangement may provide a mechanism for lowering overall substrate binding energy while retaining the full information content of 13 bp, thus resulting in a mechanism for ensuring sequence specificity.