Histidine 90 function in 4-chlorobenzoyl-coenzyme a dehalogenase catalysis.

4-chlorobenzoyl-coenzyme A (4-CBA-CoA) dehalogenase catalyzes the hydrolytic dehalogenation of 4-CBA-CoA by attack of Asp145 on the C4 of the substrate benzoyl ring to form a Meisenheimer intermediate (EMc), followed by expulsion of chloride ion to form an arylated enzyme intermediate (EAr) and, finally, ester hydrolysis in EAr to form 4-hydroxybenzoyl-CoA (4-HBA-CoA). This study examines the contribution of the active site His90 to catalysis of this reaction pathway. The His90 residue was replaced with glutamine by site-directed mutagenesis. X-ray crystallographic analysis of H90Q dehalogenase complexed with 4-HBA-CoA revealed that the positions of the catalytic groups are unchanged from those observed in the structure of the 4-HBA-CoA-wild-type dehalogenase complex. The one exception is the Gln90 side chain, which is rotated away from the position of the His90 side chain. The vacated His90 site is occupied by two water molecules. Kinetic techniques were used to evaluate ligand binding and catalytic turnover rates in the wild-type and H90Q mutant dehalogenases. The rate constants for 4-CBA-CoA (both 7 microM(-1) x s(-1)) and 4-HBA-CoA (33 and 11 microM(-1) x s(-1)) binding to the two dehalogenases are similar in value. For wild-type dehalogenase, the rate constant for a single turnover is 2.3 s(-1) while that for multiple turnovers is 0.7 s(-1). For H90Q dehalogenase, these rate constants are 1.6 x 10(-2) and 2 x 10(-4) s(-1). The rate constants for EMc formation in wild-type and mutant dehalogenase are approximately 200 s(-1) while the rate constants for EAr formation are 40 and 0.3 s(-1), respectively. The rate constant for hydrolysis of EAr in wild-type dehalogenase is 20 s(-1) and in the H90Q mutant, 0.13 s(-1). The 133-fold reduction in the rate of EAr formation in the mutant may be the result of active site hydration, while the 154-fold reduction in the rate EAr hydrolysis may be the result of lost general base catalysis. Substitution of the His90 with Gln also introduces a rate-limiting step which follows catalysis, and may involve renewing the catalytic site through a slow conformational change.     

The strength of dehalogenase-substrate hydrogen bonding correlates with the rate of Meisenheimer intermediate formation

4-Chlorobenzoyl-coenzyme A (4-CBA-CoA) dehalogenase catalyzes the hydrolytic dehalogenation of 4-CBA-CoA to 4-hydroxybenzoyl-CoA by using an active site aspartate as the nucleophile. Formation of the corresponding Meisenheimer complex (EMc) is followed by chloride ion expulsion which forms the arylated intermediate (EAr). This is then hydrolyzed to the product. In this paper, we explore the relationship between active site polarizing forces acting on the benzoyl carbonyl and the rate of formation of the Meisenheimer complex. The polarizing forces at the C[double bond]O group were modulated by introducing site-selected mutations (A112V, Y65D, G113A, G113S, G113N, and F64P), near the C[double bond]O binding site. Using either the substrate, 4-CBA-CoA, or the substrate analogue, 4-methylbenzoyl-CoA (4-MBA-CoA), Raman difference spectroscopy provided the position of the C[double bond]O stretching frequency (nu(C)[double bond](O)) for a total of 10 enzyme-ligand complexes. In turn, the values of the C[double bond]O frequencies could be converted to differences in effective hydrogen bonding strengths between members of the series, based on earlier model studies [Clarkson, J., Tonge, P. J., Taylor, K. L., Dunaway-Mariano, D., and Carey, P. (1997) Biochemistry 36, 10192-10199]. Catalysis in the F64P, G113A, G113S, and G113N dehalogenase mutants was very slow with k(cat) values ranging from 8 x 10(-3) to 7.6 x 10(-6) s(-1). The EAr intermediate did not accumulate to a detectable level on these enzymes during a single turnover. Catalysis in the Y65D and A112V dehalogenase mutants were almost as efficient as catalysis in wild-type dehalogenase with k(cat) values of 0.1-0.6 s(-1). In wild-type dehalogenase, 22% of the bound substrate accumulated as the EAr intermediate during a single turnover (k(obs) for EAr formation = 24 s(-(1)); in the Y65D mutant, the level of accumulation is 17% (k(obs) for EAr formation = 3 s(-1)), and in the A112V mutant, the level is 23% (k(obs) for EAr formation = 17 s(-1)). The k(obs) for EAr formation in wild-type dehalogenase and the more active dehalogenase mutants (Y65D and A112V) was taken to be an estimate of the k for EMc formation, and the k(obs) for EP formation in a single turnover was taken to be an estimate of the k for EMc formation in the severely impaired mutants (F64P, G113A, G113S, and G113N). A plot of the log k(obs) for EMc formation versus the C[double bond]O stretching frequency of bound 4-CBA-CoA (or 4-MBA-CoA) is a straight line (R(2) = 0.9584). Throughout the series, nu(C)[double bond](O) varied by 61 cm(-1), corresponding to the change in hydrogen bonding enthalpy of 67 kJ/mol. The results show that changes in polarizing forces at the benzoyl carbonyl are transmitted to the benzoyl (4) position and correlate with the rate of aromatic nucleophilic addition five chemical bonds away. Interestingly, the relationship between effective polarizing forces and reactivity seen here for dehalogenase is similar to that reported for the addition-elimination reaction involving the hydrolysis of a series of acyl serine proteases.     

Raman evidence for Meisenheimer complex formation in the hydrolysis reactions of 4-fluorobenzoyl- and 4-nitrobenzoyl-coenzyme A catalyzed by 4-chlorobenzoyl-coenzyme A dehalogenase

4-Chlorobenzoyl-coenzyme A (4-CBA-CoA) dehalogenase catalyzes the hydrolytic dehalogenation of 4-CBA-CoA to 4-hydroxybenzoyl-CoA by using an active site Asp145 carboxylate as the nucleophile. Formation of the corresponding Meisenheimer complex (EMc) is followed by chloride ion expulsion to form arylated enzyme (EAr). The EAr is then hydrolyzed to product. In this paper, we report the kinetics for dehalogenase-catalyzed 4-fluorobenzoyl-CoA (4-FBA-CoA) and 4-nitrobenzoyl-CoA (4-NBA-CoA) hydrolysis and provide Raman spectral evidence for the accumulation of EMc in these reactions. The 4-FBA-CoA and 4-NBA-CoA substrate analogues were selected for the poor leaving group ability of their C(4) substituents. Thus, the formation of the EAr from EMc should be hindered, giving rise to a quasi-steady-state equilibrium between EMc and the Michaelis complex. Detailed kinetic studies were carried out to quantitate the composition of the reaction mixtures. Quench experiments demonstrated that significant populations of EAr do not exist in reaction mixtures involving the 4-F- or 4-N-substrates. A kinetic model enabled us to estimate that approximately 10-20% of the enzyme-substrate complexes in the reaction mixtures are present as EMc. Raman difference spectra of 4-NBA-CoA and 4-FBA-CoA bound to WT and H90Q mutant dehalogenase have broad features near 1500 and 1220 cm(-1) that are absent in the free ligand. Crucially, these features are also absent in the Raman spectra of the complexes involving the D145A dehalogenase mutant that are unable to form an EMc. Quantum mechanical calculations, at the DFT level, provide strong support for assigning the novel 1500 and 1220 cm(-1) features to an EMc.     

Contributions of long-range electrostatic interactions to 4-chlorobenzoyl-CoA dehalogenase catalysis: a combined theoretical and experimental study

It is well established that electrostatic interactions play a vital role in enzyme catalysis. In this work, we report theory-guided mutation experiments that identified strong electrostatic contributions of a remote residue, namely, Glu232 located on the adjacent subunit, to 4-chlorobenzoyl-CoA dehalogenase catalysis. The Glu232Asp mutant was found to bind the substrate analogue 4-methylbenzoyl-CoA more tightly than does the wild-type dehalogenase. In contrast, the kcat for 4-chlorobenzoyl-CoA conversion to product was reduced 10000-fold in the mutant. UV difference spectra measured for the respective enzyme-ligand complexes revealed an approximately 3-fold shift in the equilibrium of the two active site conformers away from that inducing strong pi-electron polarization in the ligand benzoyl ring. Increased substrate binding, decreased ring polarization, and decreased catalytic efficiency indicated that the repositioning of the point charge in the Glu232Asp mutant might affect the orientation of the Asp145 carboxylate with respect to the substrate aromatic ring. The time course for formation and reaction of the arylated enzyme intermediate during a single turnover was measured for wild-type and Glu232Asp mutant dehalogenases. The accumulation of arylated enzyme in the wild-type dehalogenase was not observed in the mutant. This indicates that the reduced turnover rate in the mutant is the result of a slow arylation of Asp145, owing to decreased efficiency in substrate nucleophilic attack by Asp145. To rationalize the experimental observations, a theoretical model is proposed, which computes the potential of mean force for the nucleophilic aromatic substitution step using a hybrid quantum mechanical/molecular mechanical method. To this end, the removal or reorientation of the side chain charge of residue 232, modeled respectively by the Glu232Gln and Glu232Asp mutants, is shown to increase the rate-limiting energy barrier. The calculated 23.1 kcal/mol free energy barrier for formation of the Meisenheimer intermediate in the Glu232Asp mutant represents an increase of 6 kcal/mol relative to that of the wild-type enzyme, consistent with the 5.6 kcal/mol increase calculated from the difference in experimentally determined rate constants. On the basis of the combination of the experimental and theoretical evidence, we hypothesize that the Glu232(B) residue contributes to catalysis by providing an electrostatic force that acts on the Asp145 nucleophile.     

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.     

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.     

Linkage between substrate recognition and catalysis during cleavage of sarcin/ricin loop RNA by restrictocin

Restrictocin is a site-specific endoribonuclease that inactivates ribosomes by cleaving the sarcin/ricin loop (SRL) of 23S-28S rRNA. Here we present a kinetic and thermodynamic analysis of the SRL cleavage reaction based on monitoring the cleavage of RNA oligonucleotides (2-27-mers). Restrictocin binds to a 27-mer SRL model substrate (designated wild-type SRL) via electrostatic interactions to form a nonspecific ground state complex E:S. At pH 6.7, physical steps govern the reaction rate: the wild-type substrate reacts at a partially diffusion-limited rate, and a faster-reacting SRL, containing a 3'-sulfur atom at the scissile phosphate, reacts at a fully diffusion-limited rate (k2/K1/2 = 1.1 x 10(9) M-1 s-1). At pH 7.4, the chemical step apparently limits the SRL cleavage rate. After the nonspecific binding step, restrictocin recognizes the SRL structure, which imparts 4.3 kcal/mol transition state stabilization relative to a single-stranded RNA. The two conserved SRL modules, bulged-G motif and GAGA tetraloop, contribute at least 2.4 and 1.9 kcal/mol, respectively, to the recognition. These findings suggest a model of SRL recognition in which restrictocin contacts the GAGA tetraloop and the bulged guanosine of the bulged-G motif to progress from the nonspecific ground state complex (E:S) to the higher-energy-specific complex (E.S) en route to the chemical transition state. Comparison of restrictocin with other ribonucleases revealed that restrictocin exhibits a 10(3)-10(6)-fold smaller ribonuclease activity against single-stranded RNA than do the restrictocin homologues, non-structure-specific ribonucleases T1 and U2. Together, these findings show how structural features of the SRL substrate facilitate catalysis and provide a mechanism for distinguishing between cognate and noncognate RNA.     

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.     

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.     

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.     

Role of the minor energetic determinants of chicken egg white lysozyme (HEWL) to the stability of the HEWL.antibody scFv-10 complex

Seven of the 13 non-glycine contact amino acids in the hen (chicken) egg white lysozyme (HEWL) epitope for antibody Fab-10 each contribute < or =0.3 kcal/mol to the change in free energy (DeltaDeltaG(D)) from wild type (WT) when replaced by alanine (nullspots), and three others each give (0.7 < DeltaDeltaG(D) < or = 1. 0) kcal/mol (warm spots) (Rajpal et al. Protein Sci 1998;7:1868-1874). The low DeltaDeltaG(D) values introduced by alanine mutations present an opportunity to explore accurately their cumulative effects, as the sum of the combined DeltaDeltaG(D) values is not so large as to destabilize the complex beyond the range of accurate measurement. Substitution of six of the seven null spot residues by alanine leads to a cumulative DeltaDeltaG(D) = 2.25 +/- 0.04 kcal/mol, whereas the sum of the six individual changes is only -0.36 +/- 0.32 kcal/mol. The triple warm spot mutation generates a DeltaDeltaG(D) = 5.11 +/- 0.06 kcal/mol versus DeltaDeltaG(D) = 2.52 +/- 0.22 kcal/mol for the sum of the three individuals. The non-additivity in the individual DeltaDeltaG(D) values for the alanine mutations may indicate that these residues provide a conformationally stabilizing effect on the hot spot residues, each of which exhibits DeltaDeltaG(D) > 4.0 kcal/mol on alanine substitution.     

Computational studies of the farnesyltransferase ternary complex part II: the conformational activation of farnesyldiphosphate

Studies aimed at elucidating the reaction mechanism of farnesyltransferase (FTase), which catalyzes the prenylation of many cellular signaling proteins including Ras, has been an active area of research. Much is known regarding substrate binding and the impact of various catalytic site residues on catalysis. However, the molecular level details regarding the conformational rearrangement of farnesyldiphosphate (FPP), which has been proposed via structural analysis and mutagenesis studies to occur prior to the chemical step, is still poorly understood. Following on our previous computational characterization of the resting state of the FTase ternary complex, the thermodynamics of the conformational rearrangement step in the absence of magnesium was investigated for the wild type FTase and the Y300Fbeta mutant complexed with the peptide CVIM. In addition, we also explored the target dependence of the conformational activation step by perturbing isoleucine into a leucine (CVLM). The calculated free energy profiles of the proposed conformational transition confirm the presence of a stable intermediate state, which was identified only when the diphosphate is monoprotonated (FPP2-). The farnesyl group in the computed intermediate state assumes a conformation similar to that of the product complex, particularly for the first two isoprene units. We found that Y300beta can readily form hydrogen bonds with either of the phosphates of FPP. Removing the hydroxyl group on Y300beta does not significantly alter the thermodynamics of the conformational transition, but shifts the location of the intermediate farther away from the nucleophile by 0.5 A, which suggests that Y300beta facilitate the reaction by stabilizing the chemical step. Our results also showed an increased transition barrier height for CVLM (1.5 kcal/mol higher than that of CVIM). Although qualitatively consistent with the findings from the recent kinetic isotope experiments by Fierke and co-workers, the magnitude is not large enough to affect the rate-limiting step.     

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.     

Evaluation of direct and cooperative contributions towards the strength of buried hydrogen bonds and salt bridges

An experimental approach to evaluate the net binding free energy of buried hydrogen bonds and salt bridges is presented. The approach, which involves a modified multiple-mutant cycle protocol, was applied to selected interactions between TEM-1-beta-lactamase and its protein inhibitor, BLIP. The selected interactions (two salt bridges and two hydrogen bonds) all involving BLIP-D49, define a distinct binding unit. The penta mutant, where all side-chains constructing the binding unit were mutated to Ala, was used as a reference state to which combinations of side-chains were introduced. At first, pairs of interacting residues were added allowing the determination of interaction energies in the absence of neighbors, using double mutant cycles. Addition of neighboring residues allowed the evaluation of their cooperative effects on the interaction. The two isolated salt bridges were either neutral or repulsive whereas the two hydrogen bonds contribute 0.3 kcal mol(-1 )each. Conversely, a double mutant cycle analysis of these interactions in their native environment showed that they all stabilize the complex by 1-1.5 kcal mol(-1). Examination of the effects of neighboring residues on each of the interactions revealed that the formation of a salt bridge triad, which involves two connected salt bridges, had a strong cooperative effect on stabilizing the complex independent of the presence or absence of additional neighbors. These results demonstrate the importance of forming net-works of buried salt bridges. We present theoretical electrostatic calculations which predict the observed mode of cooperativity, and suggest that the cooperative networking effect results from the favorable contribution of the protein to the interaction. Furthermore, a good correlation between calculated and experimentally determined interaction energies for the two salt bridges, and to a lesser extent for the two hydrogen bonds, is shown. The data analysis was performed on values of DeltaDeltaG(double dagger)K(d) which reflect the strength of short range interactions, while DeltaDeltaG(o)K(D) values which include the effects of long range electrostatic forces that alter specifically DeltaDeltaG(double dagger)k(a) were treated separately.     

A binding free energy hot spot in the ankyrin repeat protein GABPbeta mediated protein-protein interaction

The frequently observed ankyrin repeat motif represents a structural scaffold evolved for mediating protein-protein interactions. As such, these repeats modulate a diverse range of cellular functions. We thermodynamically characterized the heterodimeric GA-binding protein (GABP) alphabeta complex and focused specifically on the interaction mediated by the ankyrin repeat domain of the GABPbeta. Our isothermal titration calorimetric analysis of the interaction between the GABP subunits determined an association constant (K(A)) of 6.0 x 10(8) M(-1) and that the association is favorably driven by a significant change in enthalpy (DeltaH) and a minor change in entropy (-TDeltaS). A total of 16 GABPbeta interface residues were chosen for alanine scanning mutagenesis. The calorimetrically measured differences in the free energy of binding were compared to computationally calculated values resulting in a correlation coefficient r = 0.71. We identified three spatially contiguous hydrophobic and aromatic residues that form a binding free energy hot spot (DeltaDeltaG > 2.0 kcal/mol). One residue provides structural support to the hot spot residues. Three non-hot spot residues are intermediate contributors (DeltaDeltaG approximately 1.0 kcal/mol) and create a canopy-like structure over the hot spot residues to possibly occlude solvent and orientate the subunits. The remaining interface residues are located peripherally and have weak contributions. Finally, our mutational analysis revealed a significant entropy-enthalpy compensation for this interaction.     

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.     

Energetic contribution of non-essential 5' sequence to catalysis in a hepatitis delta virus ribozyme

Hepatitis delta virus (HDV) ribozymes employ multiple catalytic strategies to achieve overall rate enhancement of RNA cleavage. These strategies include general acid-base catalysis by a cytosine side chain and involvement of divalent metal ions. Here we used a trans-acting form of the antigenomic ribozyme to examine the contribution of the 5' sequence in the substrate to HDV ribozyme catalysis. The cleavage rate constants increased for substrates with 5' sequence alterations that reduced ground-state binding to the ribozyme. Quantitatively, a plot of activation free energy of chemical conversion versus Gibb's free energy of substrate binding revealed a linear relationship with a slope of -1. This relationship is consistent with a model in which components of the substrate immediately 5' to the cleavage site in the HDV ribozyme-substrate complex destabilize ground-state binding. The intrinsic binding energy derived from the ground-state destabilization could contribute up to 2 kcal/mol toward the total 8.5 kcal/mol reduction in activation free energy for RNA cleavage catalyzed by the HDV ribozyme.     

Structure of mandelate racemase with bound intermediate analogues benzohydroxamate and cupferron

Mandelate racemase (MR, EC 5.1.2.2) from Pseudomonas putida catalyzes the Mg(2+)-dependent interconversion of the enantiomers of mandelate, 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, we determined the X-ray crystal structures of wild-type MR complexed with two analogues of the putative aci-carboxylate intermediate, benzohydroxamate and Cupferron, to 2.2-? resolution. Benzohydroxamate is shown to be a reasonable mimic of the transition state and/or intermediate because its binding affinity for 21 MR variants correlates well with changes in the free energy of transition state stabilization afforded by these variants. Both benzohydroxamate and Cupferron chelate the active site divalent metal ion and are bound in a conformation with the phenyl ring coplanar with the hydroxamate and diazeniumdiolate moieties, respectively. Structural overlays of MR complexed with benzohydroxamate, Cupferron, and the ground state analogue (S)-atrolactate reveal that the para carbon of the substrate phenyl ring moves by 0.8-1.2 ? between the ground state and intermediate state, consistent with the proposal that the phenyl ring moves during MR catalysis while the polar groups remain relatively fixed. Although the overall protein structure of MR with bound intermediate analogues is very similar to that of MR with bound (S)-atrolactate, the intermediate-Mg(2+) distance becomes shorter, suggesting a tighter complex with the catalytic Mg(2+). In addition, Tyr 54 moves closer to the phenyl ring of the bound intermediate analogues, contributing to an overall constriction of the active site cavity. However, site-directed mutagenesis experiments revealed that the role of Tyr 54 in MR catalysis is relatively minor, suggesting that alterations in enzyme structure that contribute to discrimination between the altered substrate in the transition state and the ground state by this proficient enzyme are extremely subtle.     

Modulating electron density in the bound product, 4-hydroxybenzoyl-CoA, by mutations in 4-chlorobenzoyl-CoA dehalogenase near the 4-hydroxy group

The enzyme 4-chlorobenzoyl-CoA dehalogenase hydrolyzes 4-chlorobenzoyl-CoA (4-CBA-CoA) to 4-hydroxybenzoyl-CoA (4-HBA-CoA). Biochemical and crystallographic studies have identified a critical role for the dehalogenase residue Asp 145 in close proximity to the ligand's 4-hydroxy group in the structure of the product-enzyme complex. In the present study the effects of site selective mutations at Asp 145 on the product complex are explored by Raman spectroscopy. The spectral signatures of the WT-product complex, the large red shift in lambdamax, and the complete reorganization of the benzoyl ring modes in Raman data are absent for the D145E complex. The major spectral perturbations in the WT complex are brought about by strong electron "pull" at the benzoyl carbonyl and electron "push" by the side chain of Asp 145 near the 4-OH group. Acting in concert, these factors polarize the benzoyl's pi-electrons. Since the Raman data show that very strong electron pull occurs at the benzoyl's carbonyl in the D145E complex, it is apparent that the needed electron push near the benzoyl's 4-OH group is missing. Thus, very precise positioning of Asp 145's side chain near the benzoyl's 4-position is needed to bring about the dramatic electron reorganization seen in the WT complex, and this criterion cannot be met by the glutamate side chain with its additional CH2 group. For two other Asp145 mutants D145A and D145S that lack catalytic activity, Raman difference spectroscopic data for product complexes demonstrate the presence of a population of ionized product (i.e., 4-O-) in the active sites. The presence of the ionized phenolate form explains the observation that these complexes have highly red-shifted absorbance maxima with lambdamaxs near 400 nm. For the WT complex only the 4-OH form is seen, ionization being energetically expensive with the presence of the proximal negative charge on the Asp 145 side chain. Semiquantitative estimates of the pKa for the bound product in D145S and D145A indicate that this ionization lies in the pH 6.5-7.0 range. This is approximately 2 pH units below the pKa for the free product. The Raman spectrum of 4-dimethylaminobenzoyl-CoA undergoes major changes upon binding to dehalogenase. The bound form has two features near 1562 and 1529 cm-1 and therefore closely resembles the spectrum of product bound to wild-type enzyme, which underlines the quinonoid nature in these complexes. The use of a newly developed Raman system allowed us to obtain normal (nonresonance) Raman data for the dehalogenase complexes in the 100-300 microM range and heralds an important advance in the application of Raman spectroscopy to dilute solutions of macromolecules.     

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.