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.     

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 and structural basis for transition-state stabilization in antibody-catalyzed hydrolysis

Catalytic antibodies 6D9 and 9C10, which were induced by immunization with a haptenic transition-state analog (TSA), catalyze the hydrolysis of a nonbioactive chloramphenicol monoester derivative to generate a bioactive chloramphenicol. These antibodies stabilize the transition state to catalyze the hydrolysis reaction, strictly according to the theoretical relationship: for 6D9, k(cat)/k(uncat)=895 and K(S)/K(TSA)=900, and for 9C10, k(cat)/k(uncat)=56 and K(S)/K(TSA)=60. To elucidate the molecular basis of the antibody-catalyzed reaction, the crystal structure of 6D9 was determined, and the binding thermodynamics of 6D9 and 9C10 with both the substrate and the TSA were analyzed using isothermal titration calorimetry. The crystal structure of the unliganded 6D9 Fab was determined at 2.25 A resolution and compared with that of the TSA-liganded 6D9 Fab reported previously, showing that the TSA is bound into the hydrophobic pocket of the antigen-combining site in an "induced fit" manner, especially at the L1 and H3 CDR loops. Thermodynamic analyses showed that 6D9 binds the substrate of the TSA with a positive DeltaS, differing from general thermodynamic characteristics of antigen-antibody interactions. This positive DeltaS could be due to the hydrophobic interactions between 6D9 and the substrate or the TSA mediated by Trp H100i. The difference in DeltaG between substrate and TSA-binding to 6D9 was larger than that to 9C10, which is in good correlation with the larger k(cat) value of 6D9. Interestingly, the DeltaDeltaG was mainly because of the DeltaDeltaH. The correlation between k(cat) and DeltaDeltaH is suggestive of "enthalpic strain" leading to destabilization of antibody-substrate complexes. Together with X-ray structural analyses, the thermodynamic analyses suggest that upon binding the substrate, the antibody alters the conformation of the ester moiety in the substrate from the planar Z form to a thermodynamically unstable twisted conformation, followed by conversion into the transition state. Enthalpic strain also contributes to the transition-state stabilization by destabilizing the ground state, and its degree is much larger for the more efficient catalytic antibody, 6D9.     

Role of active site binding interactions in 4-chlorobenzoyl-coenzyme A dehalogenase catalysis.

4-Chlorobenzoyl-coenzyme A (4-CBA-CoA) dehalogenase catalyzes the hydrolytic dehalogenation of 4-CBA-CoA to 4-hydroxybenzoyl-CoA (4-HBA-CoA) via a multistep mechanism involving initial attack of Asp145 on C(4) of the substrate benzoyl ring to form a Meisenheimer intermediate (EMc), followed by expulsion of the chloride ion to form an arylated enzyme intermediate (EAr) and then ester hydrolysis in the EAr to form product. This study examines the role of binding interactions in dehalogenase catalysis. The enzyme and substrate groups positioned for favorable binding interaction were identified from the X-ray crystal structure of the enzyme-4-HBA-3'-dephospho-CoA complex. These groups were individually modified (via site-directed mutagenesis or chemical synthesis) for the purpose of disrupting the binding interaction. The changes in the Gibbs free energy of the enzyme-substrate complex (DeltaDeltaG(ES)) and enzyme-transition state complex (DeltaDeltaG) brought about by the modification were measured. Cases where DeltaDeltaG exceeds DeltaDeltaG(ES) are indicative of binding interactions used for catalysis. On the basis of this analysis, we show that the H-bond interactions between the Gly114 and Phe64 backbone amide NHs and the substrate benzoyl C=O group contribute an additional 3.1 kcal/mol of stabilization at the rate-limiting transition state. The binding interactions between the enzyme and the substrate CoA nucleotide moiety also intensify in the rate-limiting transition state, reducing the energy barrier to catalysis by an additional 3.3 kcal/mol. Together, these binding interactions contribute approximately 10(6) to the k(cat)/K(m).     

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.     

Hydrophobic nature of the active site of mandelate racemase

Mandelate racemase (EC 5.1.2.2) from Pseudomonas putida catalyzes the interconversion of the two enantiomers of mandelic acid with remarkable proficiency, stabilizing the altered substrate in the transition state by approximately 26 kcal/mol. We have used a series of substrate analogues (glycolates) and intermediate analogues (hydroxamates) to evaluate the contribution of the hydrophobic cavity within the enzyme's active site to ligand binding. Free energy changes accompanying binding of glycolate derivatives correlated well with the hydrophobic substituent constant pi and the van der Waals surface areas of the ligands. The observed dependence of the apparent binding free energy on surface area of the ligand was -30 +/- 5 cal mol(-1) A(-2) at 25 degrees C. Free energy changes accompanying binding of hydroxamate derivatives also correlated well with pi values and the van der Waals surface areas of the ligands, giving a slightly greater free energy dependence equal to -41 +/- 3 cal mol(-1) A(-2) at 25 degrees C. Surprisingly, mandelate racemase exhibited a binding affinity for the intermediate analogue benzohydroxamate that was 2 orders of magnitude greater than that predicted solely on the basis of hydrophobic interactions. This suggests that there are additional specific interactions that stabilize the altered substrate in the transition state. Mandelate racemase was competitively inhibited by (R,S)-1-naphthylglycolate (apparent K(i) = 1.9 +/- 0.1 mM) and (R,S)-2-naphthylglycolate (apparent K(i) = 0.52 +/- 0.03 mM), demonstrating the plasticity of the hydrophobic pocket. Both (R)- (K(m) = 0.46 +/- 0.06 mM, k(cat) = 33 +/- 1 s(-1)) and (S)-2-naphthylglycolate (K(m) = 0.41 +/- 0.03 mM, k(cat) = 25 +/- 1 s(-1)) were shown to be alternative substrates for mandelate racemase. These kinetic results demonstrate that no major steric restrictions are imposed on the binding of this bulkier substrate in the ground state but that steric factors appear to impair transition state/intermediate stabilization. 2-Naphthohydroxamate was identified as a competitive inhibitor of mandelate racemase, binding with an affinity (K(i) = 57 +/- 18 microM) that was reduced relative to that observed for benzohydroxamate and that was in accord with the approximately 10-fold reduction in the value of k(cat)/K(m) for the racemization of 2-naphthylglycolate. These findings indicate that, for mandelate racemase, steric constraints within the hydrophobic cavity of the enzyme-intermediate complex are more stringent than those in the enzyme-substrate complex.     

Effects of temperature and viscosity on R67 dihydrofolate reductase catalysis

R67 dihydrofolate reductase (DHFR) is a novel homotetrameric protein that possesses 222 symmetry and a single, voluminous active site pore. This symmetry poses numerous limitations on catalysis; for example, two dihydrofolate (DHF) molecules or two NADPH molecules, or one substrate plus one cofactor can bind. Only the latter combination leads to catalysis. To garner additional information on how this enzyme facilitates transition-state formation, the temperature dependence of binding and catalysis was monitored. The binding of NADPH and DHF is enthalpy-driven. Previous primary isotope effect studies indicate hydride transfer is at least partially rate-determining. Accordingly, the activation energy associated with transition-state formation was measured and is found to be 6.9 kcal/mol (DeltaH(++)(25) = 6.3 kcal/mol). A large entropic component is also found associated with catalysis, TDeltaS(++)(25) = -11.3 kcal/mol. The poor substrate, dihydropteroate, binds more weakly than dihydrofolate (DeltaDeltaG = 1.4 kcal/mol) and displays a large loss in the binding enthalpy value (DeltaDeltaH = 3.8 kcal/mol). The k(cat) value for dihydropteroate reduction is decreased 1600-fold compared to DHF usage. This effect appears to derive mostly from the DeltaDeltaH difference in binding, demonstrating that the glutamate tail is important for catalysis. This result is surprising, as the para-aminobenzoyl-glutamate tail of DHF has been previously shown to be disordered by both NMR and crystallography studies. Viscosity studies were also performed and confirmed that the hydride transfer rate is not sensitive to sucrose addition. Surprisingly, binding of DHF, by both K(m) and K(d) determination, was found to be sensitive to added viscogens, suggesting a role for water in DHF binding.     

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.     

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.     

Contribution of active-site residues to the function of onconase, a ribonuclease with antitumoral activity

Onconase (ONC), a homologue of ribonuclease A (RNase A), is in clinical trials for the treatment of cancer. ONC possesses a conserved active-site catalytic triad, which is composed of His10, Lys31, and His97. The three-dimensional structure of ONC suggests that two additional residues, Lys9 and an N-terminal lactam formed from a glutamine residue (Pca1), could also contribute to catalysis. To determine the role of Pca1, Lys9, and Lys31 in the function of ONC, site-directed mutagenesis was used to replace each with alanine. Values of k(cat)/K(M) for the variants were determined with a novel fluorogenic substrate, which was designed to match the nucleobase specificity of ONC and gives the highest known k(cat)/K(M) value for the enzyme. The K9A and K31A variants display 10(3)-fold lower k(cat)/K(M) values than the wild-type enzyme, and a K9A/K31A double variant suffers a >10(4)-fold decrease in catalytic activity. In addition, replacing Lys9 or Lys31 eliminates the antitumoral activity of ONC. The side chains of Pca1 and Lys9 form a hydrogen bond in crystalline ONC. Replacing Pca1 with an alanine residue lowers the catalytic activity of ONC by 20-fold. Yet, replacing Pca1 in the K9A variant enzyme does not further reduce catalytic activity, revealing that the function of the N-terminal pyroglutamate residue is to secure Lys9. The thermodynamic cycle derived from k(cat)/K(M) values indicates that the Pca1...Lys9 hydrogen bond contributes 2.0 kcal/mol to the stabilization of the rate-limiting transition state during catalysis. Finally, binding isotherms with a substrate analogue indicate that Lys9 and Lys31 contribute little to substrate binding and that the low intrinsic catalytic activity of ONC originates largely from the low affinity of the enzyme for its substrate. These findings could assist the further development of ONC as a cancer chemotherapeutic.     

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.     

Entropy-driven binding of picomolar transition state analogue inhibitors to human 5'-methylthioadenosine phosphorylase

Human 5'-methylthioadenosine phosphorylase (MTAP) links the polyamine biosynthetic and S-adenosyl-l-methionine salvage pathways and is a target for anticancer drugs. p-Cl-PhT-DADMe-ImmA is a 10 pM, slow-onset tight-binding transition state analogue inhibitor of the enzyme. Titration of homotrimeric MTAP with this inhibitor established equivalent binding and independent catalytic function of the three catalytic sites. Thermodynamic analysis of MTAP with tight-binding inhibitors revealed entropic-driven interactions with small enthalpic penalties. A large negative heat capacity change of -600 cal/(mol K) upon inhibitor binding to MTAP is consistent with altered hydrophobic interactions and release of water. Crystal structures of apo MTAP and MTAP in complex with p-Cl-PhT-DADMe-ImmA were determined at 1.9 and 2.0 ? resolution, respectively. Inhibitor binding caused condensation of the enzyme active site, reorganization at the trimer interfaces, the release of water from the active sites and subunit interfaces, and compaction of the trimeric structure. These structural changes cause the entropy-favored binding of transition state analogues. Homotrimeric human MTAP is contrasted to the structurally related homotrimeric human purine nucleoside phosphorylase. p-Cl-PhT-DADMe-ImmA binding to MTAP involves a favorable entropy term of -17.6 kcal/mol with unfavorable enthalpy of 2.6 kcal/mol. In contrast, binding of an 8.5 pM transition state analogue to human PNP has been shown to exhibit the opposite behavior, with an unfavorable entropy term of 3.5 kcal/mol and a favorable enthalpy of -18.6 kcal/mol. Transition state analogue interactions reflect protein architecture near the transition state, and the profound thermodynamic differences for MTAP and PNP suggest dramatic differences in contributions to catalysis from protein architecture.     

Reaction intermediate analogues for mandelate racemase: interaction between Asn 197 and the alpha-hydroxyl of the substrate promotes catalysis

Mandelate racemase (MR) catalyzes the interconversion of the enantiomers of mandelic acid, stabilizing the altered substrate in the transition state by 26 kcal/mol relative to the substrate in the ground state. To understand the origins of this binding discrimination, carboxylate-, phosphonate-, and hydroxamate-containing substrate and intermediate analogues were examined for their ability to inhibit MR. Comparison of the competitive inhibition constants revealed that an alpha-hydroxyl function is required for recognition of the ligand as an intermediate analogue. Two intermediate analogues, alpha-hydroxybenzylphosphonate (alpha-HBP) and benzohydroxamate, were bound with affinities approximately 100-fold greater than that observed for the substrate. Furthermore, MR bound alpha-HBP enantioselectively, displaying a 35-fold higher affinity for the (S)-enantiomer relative to the (R)-enantiomer. In the X-ray structure of mandelate racemase [Landro, J. A., Gerlt, J. A., Kozarich, J. W., Koo, C. W., Shah, V. J., Kenyon, G. L., Neidhart, D. J., Fujita, J., and Petsko, G. A. (1994) Biochemistry 33, 635-643], the alpha-hydroxyl function of the competitive inhibitor (S)-atrolactate is within hydrogen bonding distance of Asn 197. To demonstrate the importance of the alpha-hydroxyl function in intermediate binding, the N197A mutant was constructed. The values of k(cat) for N197A were reduced 30-fold for (R)-mandelate and 179-fold for (S)-mandelate relative to wild-type MR; the values of k(cat)/K(m) were reduced 208-fold for (R)-mandelate and 556-fold for (S)-mandelate. N197A shows only a 3.5-fold reduction in its affinity for the substrate analogue (R)-atrolactate but a 51- and 18-fold reduction in affinity for alpha-HBP and benzohydroxamate, respectively. Thus, interaction between Asn 197 and the substrate's alpha-hydroxyl function provides approximately 3.5 kcal/mol of transition-state stabilization free energy to differentially stabilize the transition state relative to the ground state.     

Rescue of the orphan enzyme isoguanine deaminase

Cytosine deaminase (CDA) from Escherichia coli was shown to catalyze the deamination of isoguanine (2-oxoadenine) to xanthine. Isoguanine is an oxidation product of adenine in DNA that is mutagenic to the cell. The isoguanine deaminase activity in E. coli was partially purified by ammonium sulfate fractionation, gel filtration, and anion exchange chromatography. The active protein was identified by peptide mass fingerprint analysis as cytosine deaminase. The kinetic constants for the deamination of isoguanine at pH 7.7 are as follows: k(cat) = 49 s(-1), K(m) = 72 μM, and k(cat)/K(m) = 6.7 × 10(5) M(-1) s(-1). The kinetic constants for the deamination of cytosine are as follows: k(cat) = 45 s(-1), K(m) = 302 μM, and k(cat)/K(m) = 1.5 × 10(5) M(-1) s(-1). Under these reaction conditions, isoguanine is the better substrate for cytosine deaminase. The three-dimensional structure of CDA was determined with isoguanine in the active site.     

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.     

Binding energy and catalysis by D-xylose isomerase: kinetic, product, and X-ray crystallographic analysis of enzyme-catalyzed isomerization of (R)-glyceraldehyde

D-Xylose isomerase (XI) and triosephosphate isomerase (TIM) catalyze the aldose-ketose isomerization reactions of D-xylose and d-glyceraldehyde 3-phosphate (DGAP), respectively. D-Glyceraldehyde (DGA) is the triose fragment common to the substrates for XI and TIM. The XI-catalyzed isomerization of DGA to give dihydroxyacetone (DHA) in D(2)O was monitored by (1)H nuclear magnetic resonance spectroscopy, and a k(cat)/K(m) of 0.034 M(-1) s(-1) was determined for this isomerization at pD 7.0. This is similar to the k(cat)/K(m) of 0.017 M(-1) s(-1) for the TIM-catalyzed carbon deprotonation reaction of DGA in D(2)O at pD 7.0 [Amyes, T. L., O'Donoghue, A. C., and Richard, J. P. (2001) J. Am. Chem. Soc. 123, 11325-11326]. The much larger activation barrier for XI-catalyzed isomerization of D-xylose (k(cat)/K(m) = 490 M(-1) s(-1)) versus that for the TIM-catalyzed isomerization of DGAP (k(cat)/K(m) = 9.6 × 10(6) M(-1) s(-1)) is due to (i) the barrier to conversion of cyclic d-xylose to the reactive linear sugar (5.4 kcal/mol) being larger than that for conversion of DGAP hydrate to the free aldehyde (1.7 kcal/mol) and (ii) the intrinsic binding energy [Jencks, W. P. (1975) Adv. Enzymol. Relat. Areas Mol. Biol. 43, 219-410] of the terminal ethylene glycol fragment of D-xylose (9.3 kcal/mol) being smaller than that of the phosphodianion group of DGAP (~12 kcal/mol). The XI-catalyzed isomerization of DGA in D(2)O at pD 7.0 gives a 90% yield of [1-(1)H]DHA and a 10% yield of [1-(2)H]DHA, the product of isomerization with incorporation of deuterium from solvent D(2)O. By comparison, the transfer of (3)H from the labeled hexose substrate to solvent is observed only once in every 10(9) turnovers for the XI-catalyzed isomerization of [2-(3)H]glucose in H(2)O [Allen, K. N., Lavie, A., Farber, G. K., Glasfeld, A., Petsko, G. A., and Ringe, D. (1994) Biochemistry 33, 1481-1487]. We propose that truncation of the terminal ethylene glycol fragment of d-xylose to give DGA results in a large decrease in the rate of XI-catalyzed isomerization with hydride transfer compared with that for proton transfer. An ultra-high-resolution (0.97 ?) X-ray crystal structure was determined for the complex obtained by soaking crystals of XI with 50 mM DGA. The triose binds to XI as the unreactive hydrate, but ligand binding induces metal cofactor movement and conformational changes in active site residues similar to those observed for XI·sugar complexes.     

Role of the flavin domain residues,His689 and Asn732, in the catalytic mechanism of cellobiose dehydrogenase from phanerochaete chrysosporium

Cellobiose dehydrogenase is an extracellular flavocytochrome, which catalyzes the oxidation of cellobiose and other soluble oligosaccharides to their respective lactones, while reducing various one- and two-electron acceptors. Two residues at the active site of the flavin domain, His689 and Asn732, have been proposed to play critical roles in the oxidation of the substrate. To test these proposals, each residue was substituted with either a Gln, Asn, Glu, Asp, Val, Ala, and/or a His residue by site-directed mutagenesis, using a homologous expression system previously developed in our laboratory. This enabled an examination of the functional, stereochemical, and electrostatic constraints for binding and oxidation of the substrate. The steady-state kinetic parameters for the variant proteins were compared using cellobiose and its epimer, lactose, as the substrates. The H689 variants all exhibit >1000-fold lower k(cat) values, while the K(m) values for both substrates in these variants are similar to that of the wild-type enzyme. This supports the proposed role of this His residue as a general base in catalysis. The N732 variants exhibit a range of kinetic parameters: the k(cat) values for oxidation are 5-4000-fold lower than that for the wild-type enzyme, while the K(m) values vary between similar to and 60-fold higher than that for the wild-type. The difference in binding energy between cellobiose and lactose was calculated using the relationship delta(delta G) = -RT ln[(k(cat)/K(m))(lactose)/(k(cat)/K(m))(cellobiose)]. This calculation for the wild-type enzyme suggests that lactose binds considerably more weakly than cellobiose (7.2 kJ/mol difference), which corresponds to one extra (cumulative) hydrogen bond for cellobiose over lactose. Mutations at Asn732 result in a further weakening of lactose binding over cellobiose (2-4 kJ/mol difference). The results support a role for Asn732 in the binding of the substrate.     

Catalytic reaction pathway for the mitogen-activated protein kinase ERK2

The structural, functional, and regulatory properties of the mitogen-activated protein kinases (MAP kinases) have long attracted considerable attention owing to the critical role that these enzymes play in signal transduction. While several MAP kinase X-ray crystal structures currently exist, there is by comparison little mechanistic information available to correlate the structural data with the known biochemical properties of these molecules. We have employed steady-state kinetic and solvent viscosometric techniques to characterize the catalytic reaction pathway of the MAP kinase ERK2 with respect to the phosphorylation of a protein substrate, myelin basic protein (MBP), and a synthetic peptide substrate, ERKtide. A minor viscosity effect on k(cat) with respect to the phosphorylation of MBP was observed (k(cat) = 10 +/- 2 s(-1), k(cat)(eta) = 0.18 +/- 0.05), indicating that substrate processing occurs via slow phosphoryl group transfer (12 +/- 4 s(-1)) followed by the faster release of products (56 +/- 4 s(-1)). At an MBP concentration extrapolated to infinity, no significant viscosity effect on k(cat)/K(m(ATP)) was observed (k(cat)/K(m(ATP)) = 0.2 +/- 0.1 microM(-1) s(-1), k(cat)/K(m(ATP))(eta) = -0.08 +/- 0.04), consistent with rapid-equilibrium binding of the nucleotide. In contrast, at saturating ATP, a full viscosity effect on k(cat)/K(m) for MBP was apparent (k(cat)/K(m(MBP)) = 2.4 +/- 1 microM(-1) s(-1), k(cat)/K(m(MBP))(eta) = 1.0 +/- 0.1), while no viscosity effect was observed on k(cat)/K(m) for the phosphorylation of ERKtide (k(cat)/K(m(ERKtide)) = (4 +/- 2) x 10(-3) microM(-1) s(-1), k(cat)/K(m(ERKtide))(eta) = -0.02 +/- 0.02). This is consistent with the diffusion-limited binding of MBP, in contrast to the rapid-equilibrium binding of ERKtide, to form the ternary Michaelis complex. Calculated values for binding constants show that the estimated value for K(d(MBP)) (/= 1.5 mM). The dramatically higher catalytic efficiency of MBP in comparison to that of ERKtide ( approximately 600-fold difference) is largely attributable to the slow dissociation rate of MBP (/=56 s(-1)), from the ERK2 active site.     

Ground-state destabilization in orotate phosphoribosyltransferases by binding isotope effects

Orotate phosphoribosyltransferases (OPRTs) form and break the N-ribosidic bond to pyrimidines by way of ribocation-like transition states (TSs) and therefore exhibit large α-secondary 1'-(3)H k(cat)/K(m) kinetic isotope effects (KIEs) [Zhang, Y., and Schramm, V. L. (2010) J. Am. Chem. Soc. 132, 8787-8794]. Substrate binding isotope effects (BIEs) with OPRTs report on the degree of ground-state destabilization for these complexes and permit resolution of binding and transition-state effects from the k(cat)/K(m) KIEs. The BIEs for interactions of [1'-(3)H]orotidine 5'-monophosphate (OMP) with the catalytic sites of Plasmodium falciparum and human OPRTs are 1.104 and 1.108, respectively. These large BIEs establish altered sp(3) bond hybridization of C1' toward the sp(2) geometry of the transition states upon OMP binding. Thus, the complexes of these OPRTs distort OMP part of the way toward the transition state. As the [1'-(3)H]OMP k(cat)/K(m) KIEs are approximately 1.20, half of the intrinsic k(cat)/K(m) KIEs originate from BIEs. Orotidine, a slow substrate for these enzymes, binds to the catalytic site with no significant [1'-(3)H]orotidine BIEs. Thus, OPRTs are unable to initiate ground-state destabilization of orotidine by altered C1' hybridization because of the missing 5'-phosphate. However the k(cat)/K(m) KIEs for [1'-(3)H]orotidine are also approximately 1.20. The C1' distortion for OMP happens in two steps, half upon binding and half on going from the Michaelis complex to the TS. With orotidine as the substrate, there is no ground-state destabilization in the Michaelis complexes, but the C1' distortion at the TS is equal to that of OMP. The large single barrier for TS formation with orotidine slows the rate of barrier crossing.     

Solvation and the hidden thermodynamics of a zinc finger probed by nonstandard repair of a protein crevice

The classical Zn finger contains a phenylalanine at the crux of its three architectural elements: a beta-hairpin, an alpha-helix, and a Zn(2+)-binding site. Surprisingly, phenylalanine is not required for high-affinity Zn2+ binding, but instead contributes to the specification of a precise DNA-binding surface. Substitution of phenylalanine by leucine leads to a floppy but native-like structure whose Zn affinity is maintained by marked entropy-enthalpy compensation (DeltaDeltaH -8.3 kcal/mol and -TDeltaDeltaS 7.7 kcal/mol). Phenylalanine and leucine differ in shape, size, and aromaticity. To distinguish which features correlate with dynamic stability, we have investigated a nonstandard finger containing cyclohexanylalanine at this site. The structure of the nonstandard finger is similar to that of the native domain. The cyclohexanyl ring assumes a chair conformation, and conformational fluctuations characteristic of the leucine variant are damped. Although the nonstandard finger exhibits a lower affinity for Zn2+ than does the native domain (DeltaDeltaG -1.2 kcal/mol), leucine-associated perturbations in enthalpy and entropy are almost completely attenuated (DeltaDeltaH -0.7 kcal/mol and -TDeltaDeltaS -0.5 kcal/mol). Strikingly, global changes in entropy (as inferred from calorimetry) are in each case opposite in sign from changes in configurational entropy (as inferred from NMR). This seeming paradox suggests that enthalpy-entropy compensation is dominated by solvent reorganization rather than nominal molecular properties. Together, these results demonstrate that dynamic and thermodynamic perturbations correlate with formation or repair of a solvated packing defect rather than type of physical interaction (aromatic or aliphatic) within the core.