The role of amino-acid residues Q39 and E451 in the determination of substrate specificity of the Spodoptera frugiperda beta-glycosidase.

The design of beta-glycosidases with planed substrate specificity for biotechnological application has received little attention. This is mostly a consequence of the lack of data on the molecular basis of the beta-glycosidase specificity, namely data on the energy of the noncovalent interactions in the enzyme-transition state complex. In an attempt to fill this gap, site-directed mutagenesis and enzyme steady-state kinetic experiments with different substrates were conducted, using as model a digestive beta-glycosidase (glycoside hydrolase family 1) from Spodoptera frugiperda (Lepidoptera) (Sfbetagly50). The active site of this enzyme was modeled based on its sequence and on crystallographic data of similar enzymes. Energy of noncovalent interactions in transition state between Sfbetagly50 amino acids and glycone hydroxyls was determined. Sfbetagly50 residue E451 seems to be a key residue in determining beta-glycosidase preference for glucosides vs. galactosides based on the following data: (a) The energy of the noncovalent interaction between glycone equatorial hydroxyl 4 with E451 in the transition state is about 60% higher than its interaction with Q39. (b) The energy of the E451-hydroxyl 4 interaction decreases more than the Q39-hydroxyl 4 interaction when hydroxyl 4 is changed from equatorial to axial position. (c) A Sfbetagly50 mutant, where E451 was substituted by A, hydrolyzes galactosides faster than glucosides. It was also shown that glycone hydroxyl 6 interacts favorably with Q39, but not with E451, probably due to steric hindrance. These interactions result in the beta-glycosidase hydrolyzing fucosides (6-deoxygalactosides) faster than glucosides and galactosides.     

Biochemical and structural assessment of the 1-N-azasugar GalNAc-isofagomine as a potent family 20 beta-N-acetylhexosaminidase inhibitor

Azasugar inhibitors of the isofagomine class are potent competitive inhibitors of configuration-retaining beta-glycosidases. This potency results from the formation of a strong electrostatic interaction between a protonated endocyclic nitrogen at the "anomeric" center of the inhibitor and the catalytic nucleophile of the enzyme. Although the majority of retaining beta-glycosidases use a mechanism involving a carboxylate residue as a nucleophile, Streptomyces plicatus beta-N-acetylhexos-aminidase (SpHEX) and related family 20 glycosidases lack such a catalytic residue and use instead the carbonyl oxygen of the 2-acetamido group of the substrate as a nucleophile to "attack" the anomeric center. Thus, a strong electrostatic interaction between the inhibitor and enzyme is not expected to occur; nonetheless, the 1-N-azasugar (2R,3R,4S,5R)-2-acetamido-3,4-dihydroxy-5-hydroxymethyl-piperidinium hydrochloride (GalNAc-isofagomine.HCl), which was synthesized and assayed for its ability to inhibit SpHEX, was found to be a potent competitive inhibitor of the enzyme (K(i) = 2.7 microm). A crystallographic complex of GalNAc-isofagomine bound to SpHEX was solved and refined to 1.75 A and revealed that the lack of a strong electrostatic interaction between the "anomeric" center of GalNAc-isofagomine and SpHEX is compensated for by a novel 2.8-A hydrogen bond formed between the equatorial proton of the endocyclic nitrogen of the azasugar ring and the carboxylate of the general acid-base residue Glu-314 of SpHEX. This interaction appears to contribute to the unexpected potency of GalNAc-isofagomine toward SpHEX.     

Acidophilic adaptation of family 11 endo-beta-1,4-xylanases: modeling and mutational analysis

Xyl1 from Streptomyces sp. S38 belongs to the low molecular mass family 11 of endo-beta-1,4-xylanases. Its three-dimensional structure has been solved at 2.0 A and its optimum temperature and pH for enzymatic activity are 60 degrees C and 6.0, respectively. Aspergillus kawachii xylanase XynC belongs to the same family but is an acidophilic enzyme with an optimum pH of 2.0. Structural comparison of Xyl1 and XynC showed differences in residues surrounding the two glutamic acid side chains involved in the catalysis that could be responsible for the acidophilic adaptation of XynC. Mutations W20Y, N48D, A134E, and Y193W were introduced by site-directed mutagenesis and combined in multiple mutants. Trp 20 and Tyr 193 are involved in substrate binding. The Y193W mutation inactivated Xyl1 whereas W20Y decreased the optimum pH of Xyl1 to 5.0 and slightly increased its specific activity. The N48D mutation also decreased the optimum pH of Xyl1 by one unit. The A134E substitution did not induce any change, but when combined with N48D, a synergistic effect was observed with a 1.4 unit decrease in the optimum pH. Modeling showed that the orientations of residue 193 and of the fully conserved Arg 131 are different in acidophilic and "alkaline" xylanases whereas the introduced Tyr 20 probably modifies the pKa of the acid-base catalyst via residue Asn 48. Docking of a substrate analog in the catalytic site highlighted striking differences between Xyl1 and XynC in substrate binding. Hydrophobicity calculations showed a correlation between acidophilic adaptation and a decreased hydrophobicity around the two glutamic acid side chains involved in catalysis.     

Improving the alkalophilic performances of the Xyl1 xylanase from Streptomyces sp. S38: structural comparison and mutational analysis

Endo-beta-1,4-xylanases of the family 11 glycosyl-hydrolases are catalytically active over a wide range of pH. Xyl1 from Streptomyces sp. S38 belongs to this family, and its optimum pH for enzymatic activity is 6. Xyn11 from Bacillus agaradhaerens and XylJ from Bacillus sp. 41M-1 share 85% sequence identity and have been described as highly alkalophilic enzymes. In an attempt to better understand the alkalophilic adaptation of xylanases, the three-dimensional structures of Xyn11 and Xyl1 were compared. This comparison highlighted an increased number of salt-bridges and the presence of more charged residues in the catalytic cleft as well as an eight-residue-longer loop in the alkalophilic xylanase Xyn11. Some of these charges were introduced in the structure of Xyl1 by site-directed mutagenesis with substitutions Y16D, S18E, G50R, N92D, A135Q, E139K, and Y186E. Furthermore, the eight additional loop residues of Xyn11 were introduced in the homologous loop of Xyl1. In addition, the coding sequence of the XylJ catalytic domain was synthesized by recursive PCR, expressed in a Streptomyces host, purified, and characterized together with the Xyl1 mutants. The Y186E substitution inactivated Xyl1, but the activity was restored when this mutation was combined with the G50R or S18E substitutions. Interestingly, the E139K mutation raised the optimum pH of Xyl1 from 6 to 7.5 but had no effect when combined with the N92D substitution. Modeling studies identified the possible formation of an interaction between the introduced lysine and the substrate, which could be eliminated by the formation of a putative salt-bridge in the N92D/E139K mutant.     

Functional analysis of active site residues of Bacillus thuringiensis WB7 chitinase by site-directed mutagenesis

To investigate the roles of the active site residues in the catalysis of Bacillus thuringiensis WB7 chitinase, twelve mutants, F201L, F201Y, G203A, G203D, D205E, D205N, D207E, D207N, W208C, W208R, E209D and E209Q were constructed by site-directed mutagenesis. The results showed that the mutants F201L, G203D, D205N, D207E, D207N, W208C and E209D were devoid of activity, and the loss of the enzymatic activities for F201Y, G203A, D205E, W208R and E209Q were 72, 70, 48, 31 and 29%, respectively. The pH-activity profiles indicated that the optimum pH for the mutants as well as for the wildtype enzyme was 8.0. E209Q exhibited a broader active pH range while D205E, G203A and F201Y resulted in a narrower active pH range. The pH range of activity reduced 1 unit for D205E, and 2 units for G203A and F201Y. The temperature-activity profiles showed that the optimum temperature for other mutants as well as wildtype enzyme was 60°C, but 50°C for G203A, which suggested that G203A resulted in a reduction of thermostability. The study indicated that the six active site residues involving in mutagenesis played an important part in WB7 chitinase. In addition, the catalytic mechanisms of the six active site residues in WB7 chitinase were discussed.     

Arginine-82 regulates the pKa of the group responsible for the light-driven proton release in bacteriorhodopsin.

In wild-type bacteriorhodopsin light-induced proton release occurs before uptake at neutral pH. In contrast, in mutants in which R82 is replaced by a neutral residue (as in R82A and R82Q), only a small fraction of the protons is released before proton uptake at neutral pH; the major fraction is released after uptake. In R82Q the relative amounts of the two types of proton release, "early" (preceding proton uptake) and "late" (following proton uptake), are pH dependent. The main conclusions are that 1) R82 is not the normal light-driven proton release group; early proton release can be observed in the R82Q mutant at higher pH values, suggesting that the proton release group has not been eliminated. 2) R82 affects the pKa of the proton release group both in the unphotolyzed state of the pigment and during the photocycle. In the wild type (in 150 mM salt) the pKa of this group decreases from approximately 9.5 in the unphotolyzed pigment to approximately 5.8 in the M intermediate, leading to early proton release at neutral pH. In the R82 mutants the respective values of pKa of the proton release group in the unphotolyzed pigment and in M are approximately 8 and 7.5 in R82Q (in 1 M salt) and approximately 8 and 6.5 in R82K (in 150 mM KCl). Thus in R82Q the pKa of the proton release group does not decrease enough in the photocycle to allow early proton release from this group at neutral pH. 3) Early proton release in R82Q can be detected as a photocurrent signal that is kinetically distinct from those photocurrents that are due to proton movements from the Schiff base to D85 during M formation and from D96 to the Schiff base during the M-->N transition. 4) In R82Q, at neutral pH, proton uptake from the medium occurs during the formation of O. The proton is released during the O-->bacteriorhodopsin transition, probably from D85 because the normal proton release group cannot deprotonate at this pH. 5) The time constant of early proton release is increased from 85 microseconds in the wild type to 1 ms in R82Q (in 150 mM salt). This can be directly attributed to the increase in the pKa of the proton release group and also explains the uncoupling of proton release from M formation. 6) In the E204Q mutant only late proton release is observed at both neutral and alkaline pH, consistent with the idea that E204 is the proton release group. The proton release is concurrent with the O-->bacteriorhodopsin transition, as in R82Q at neutral pH.     

Electrostatic switches that mediate the pH-dependent conformational change of "short" recombinant human pseudocathepsin D

Human cathepsin D (hCatD) is an aspartic peptidase with a low pH optimum. X-ray crystal structures have been solved for an active, low pH (pH 5.1) form (CatD(lo)) [Baldwin, E. T., Bhat, T. N., Gulnik, S., Hosur, M. V., Sowder, R. C., Cachau, R. E., Collins, J., Silva, A. M., and Erickson, J. W. (1993) Proc. Natl. Acad. Sci. U.S.A. 90, 6796-6800] and an inactive, high pH (pH 7.5) form (CatD(hi)) [Lee, A. Y., Gulnik, S. V., and Erickson, J. W. (1998) Nat. Struct. Biol. 5, 866-871]. It has been suggested that ionizable switches involving the carboxylate side chains of E5, E180, and D187 may mediate the reversible interconversion between CatD(hi) and CatD(lo) and that Y10 stabilizes CatD(hi) [Lee, A. Y., Gulnik, S. V., and Erickson, J. W. (1998) Nat. Struct. Biol. 5, 866-871]. To test these hypotheses, we generated single point mutants in "short" recombinant human pseudocathepsin D (srCatD), a model kinetically similar to hCatD [Beyer, B. M., and Dunn, B. M. (1996) J. Biol. Chem. 271, 15590-15596]. E180Q, Y10F, and D187N exhibit significantly higher kcat/Km values (2-, 3-, and 6-fold, respectively) at pH 3.7 and 4.75 compared to srCatD, indicating that these residues are important in stabilizing the CatD(hi). E5Q exhibits a 2-fold lower kcat/Km compared to srCatD at both pH values, indicating the importance of E5 in stabilizing the CatD(lo). Accordingly, full time-course "pH-jump" (pH 5.5-4.75) studies of substrate hydrolysis indicate that E180Q, D187N, and Y10F have shorter kinetic lag phases that represent the change from CatD(hi) to CatD(lo) compared to srCatD and E5Q. Intrinsic tryptophan fluorescence reveals that the variants have a native-like structure over the pH range of our assays. The results indicate that E180 and D187 participate as an electrostatic switch that initiates the conformational change of CatD(lo) to CatD(hi) and Y10 stabilizes CatD(hi) by hydrogen bonding to the catalytic Asp 33. E5 appears to play a less significant role as an ionic switch that stabilizes CatD(lo).     

The role of carboxyl, guanidine and imidazole groups in catalysis by a midgut trehalase purified from an insect larvae

A trehalase (EC 3.2.1.28) of 67 kDa was purified to homogeneity from the midgut of Spodoptera frugiperda (Lepidoptera) larvae. The enzyme is inhibited by toxic beta-glucosides produced by plants (amygdalin, prunasin, salicin and phlorezin) and by their aglycones (mandelonitrile, phloretin). From kcat and Km values determined in different pHs, the pKa values of catalytic essential groups were calculated (pKa = 4.5 and pKa = 8.0). These pKa values agree with the ones determined from enzyme chemical in activation with carbodiimide and phenyl glyoxal, respectively, indicating that the enzyme has a carboxyl group that act as a nucleophile and a guanidine group that is the proton donor during the catalytic cycle. The enzyme has two putative subsites for glucose binding. Based on the protection afforded by ligands against chemical modification, the roles of the subsites were inferred. Thus, the one that binds the competitive inhibitors, methyl alpha-glucoside (MalphaGlu) and mandelonitrile, contains the catalytic carboxyl, whereas the other having the catalytic Arg residue binds the competitive inhibitor Tris. Diethyl pyrocarbonate is ineffective except in the presence of MalphaGlu, when it decreases trehalase activity and changes the pKa value of the catalytic Arg residue. This suggests that the pKa value of the Arg residue is modulated by a His residue located near the active site. This also indicates that the enzyme molecule changes its conformation when the subsite containing the carboxyl group is occupied. The increase in trehalase inactivation by phenyl glyoxal in the presence of MalphaGlu agrees with the last observation.     

Identification of active-site residues of the pro-metastatic endoglycosidase heparanase

Heparanase is a beta-D-endoglucuronidase that cleaves heparan sulfate (HS) and has been implicated in many important physiological and pathological processes, including tumor cell metastasis, angiogenesis, and leukocyte migration. We report herein the identification of active-site residues of human heparanase. Using PSI-BLAST and PHI-BLAST searches of sequence databases, similarities were identified between heparanase and members of several of the glycosyl hydrolase families (10, 39, and 51) from glycosyl hydrolase clan A (GH-A), including strong local identities to regions containing the critical active-site catalytic proton donor and nucleophile residues that are conserved in this clan of enzymes. Furthermore, secondary structure predictions suggested that heparanase is likely to contain an (alpha/beta)(8) TIM-barrel fold, which is common to the GH-A families. On the basis of sequence alignments with a number of glycosyl hydrolases from GH-A, Glu(225) and Glu(343) of human heparanase were identified as the likely proton donor and nucleophile residues, respectively. The substitution of these residues with alanine and the subsequent expression of the mutant heparanases in COS-7 cells demonstrated that the HS-degrading capacity of both was abolished. In contrast, the alanine substitution of two other glutamic acid residues (Glu(378) and Glu(396)), both predicted to be outside the active site, did not affect heparanase activity. These data suggest that heparanase is a member of the clan A glycosyl hydrolases and has a common catalytic mechanism that involves two conserved acidic residues, a putative proton donor at Glu(225) and a nucleophile at Glu(343).     

Biochemical characterization and identification of the catalytic residues of a family 43 beta-D-xylosidase from Geobacillus stearothermophilus T-6

Beta-D-xylosidases are hemilcellulases that hydrolyze short xylooligosaccharides into xylose units. Here, we describe the characterization and kinetic analysis of a family 43 beta-xylosidase from Geobacillus stearothermophilus T-6 (XynB3). Enzymes in this family use an inverting single-displacement mechanism with two conserved carboxylic acids, a general acid, and a general base. XynB3 was most active at 65 degrees C and pH 6.5, with clear preference to xylose-based substrates. Products analysis indicated that XynB3 is an exoglycosidase that cleaves single xylose units from the nonreducing end of xylooligomers. On the basis of sequence homology, amino acids Asp15 and Glu187 were suggested to act as the general-base and general-acid catalytic residues, respectively. Kinetic analysis with substrates bearing different leaving groups showed that, for the wild-type enzyme, the k(cat) and k(cat)/K(m) values were only marginally affected by the leaving-group reactivity, whereas for the E187G mutant, both values exhibited significantly greater dependency on the pK(a) of the leaving group. The pH-dependence activity profile of the putative general-acid mutant (E187G) revealed that the protonated catalytic residue was removed. Addition of the exogenous nucleophile azide did not affect the activities of the wild type or the E187G mutant but rescued the activity of the D15G mutant. On the basis of thin-layer chromatography and (1)H NMR analyses, xylose and not xylose azide was the only product of the accelerated reaction, suggesting that the azide ion does not attack the anomeric carbon directly but presumably activates a water molecule. Together, these results confirm the suggested catalytic role of Glu187 and Asp15 in XynB3 and provide the first unequivocal evidence regarding the exact roles of the catalytic residues in an inverting GH43 glycosidase.     

Engineering of the pH optimum of Bacillus cereus beta-amylase: conversion of the pH optimum from a bacterial type to a higher-plant type

The optimum pH of Bacillus cereus beta-amylase (BCB, pH 6.7) differs from that of soybean beta-amylase (SBA, pH 5.4) due to the substitution of a few amino acid residues near the catalytic base residue (Glu 380 in SBA and Glu 367 in BCB). To explore the mechanism for controlling the optimum pH of beta-amylase, five mutants of BCB (Y164E, Y164F, Y164H, Y164Q, and Y164Q/T47M/Y164E/T328N) were constructed and characterized with respect to enzymatic properties and X-ray structural crystal analysis. The optimum pH of the four single mutants shifted to 4.2-4.8, approximately 2 pH units and approximately 1 pH unit lower than those of BCB and SBA, respectively, and their k(cat) values decreased to 41-3% of that of the wild-type enzyme. The X-ray crystal analysis of the enzyme-maltose complexes showed that Glu 367 of the wild type is surrounded by two water molecules (W1 and W2) that are not found in SBA. W1 is hydrogen-bonded to both side chains of Glu 367 and Tyr 164. The mutation of Tyr 164 to Glu and Phe resulted in the disruption of the hydrogen bond between Tyr 164 Oeta and W1 and the introduction of two additional water molecules near position 164. In contrast, the triple mutant of BCB with a slightly decreased pH optimum at pH 6.0 has no water molecules (W1 and W2) around Glu 367. These results suggested that a water-mediated hydrogen bond network (Glu 367...W1...Tyr 164...Thr 328) is the primary requisite for the increased pH optimum of wild-type BCB. This strategy is completely different from that of SBA, in which a hydrogen bond network (Glu 380...Thr 340...Glu 178) reduces the optimum pH in a hydrophobic environment.     

Contribution of tyrosine 6 to the catalytic mechanism of isoenzyme 3-3 of glutathione S-transferase.

The role of the hydroxyl group of tyrosine 6 in the catalytic mechanism of isoenzyme 3-3 of rat glutathione S-transferase has been examined by x-ray crystallography and site-specific replacement of the residue with phenylalanine and evaluation of the catalytic properties of the mutant enzyme. This particuar tyrosine residue is conserved in the sequences of all of the cytosolic enzymes and is found, in crystal structures of both isoenzyme 3-3 from the mu-gene class and an isoenzyme from the pi-gene class, to be proximal to the sulfur of glutathione (GSH) or glutathione sulfonate bound at the active site. The 2.2-A structure of the binary complex of isoenzyme 3-3 and GSH indicates that the hydroxyl group of Tyr6 is located 3.2-3.5 A from the sulfur of GSH, well within hydrogen bonding distance. Removal of the hydroxyl group of Tyr6 has essentially no effect on the dissociation constant (22 +/- 3 microM) for GSH. Nevertheless the Y6F mutant exhibits a turnover number which is only about 1% that of the native enzyme when assayed at pH 6.5 with either 1-chloro-2,4-dinitrobenzene (CDNB) or 4-phenyl-3-buten-2-one. UV difference spectra of the binary enzyme-GSH complexes suggest that the predominant ionization state of GSH in the active site of the Y6F mutant is the neutral thiol (e.g. EY6F.GSH) which is in contrast to the native enzyme in which the thiol is substantially deprotonated (e.g. E.GS-). Spectrophotometric titration suggests that the pKa of the thiol is 6.9 +/- 0.3 in the E.GSH complex and greater than or equal to 8 in the EY6F.GSH binary complex. In addition, the pH dependence of kcat/KmCDNB reveals that the reactions catalyzed by the native enzyme and the Y6F mutant are dependent on a single ionization in the E.GSH and EY6F.GSH complexes with pKa = 6.2 +/- 0.1 and 7.8 +/- 0.3, respectively. The results suggest that the hydrogen bond between Tyr6 and the enzyme-bound nucleophile helps to lower the pKa of GSH in the binary enzyme-substrate complex.     

The proton release group of bacteriorhodopsin controls the rate of the final step of its photocycle at low pH

The factors determining the pH dependence of the formation and decay of the O photointermediate of the bacteriorhodopsin (bR) photocycle were investigated in the wild-type (WT) pigment and in the mutants of Glu-194 and Glu-204, key residues of the proton release group (PRG) in bR. We have found that in the WT the rate constant of O --> bR transition decreases 30-fold upon decreasing the pH from 6 to 3 with a pKa of about 4.3. D2O slows the rise and decay of the O intermediate in the WT at pH 3.5 by a factor of 5.5. We suggest that the rate of the O --> bR transition (which reflects the rate of deprotonation of the primary proton acceptor Asp-85) at low pH is controlled by the deprotonation of the PRG. To test this hypothesis, we studied the E194D mutant. We show that the pKa of the PRG in the ground state of the E194D mutant, when Asp-85 is protonated, is increased by 1.2 pK units compared to that of the WT. We found a similar increase in the pKa of the rate constant of the O --> bR transition in E194D. This provides further evidence that the rate of the O --> bR transition is controlled by the PRG. In a further test, the E194Q mutation, which disables the PRG and slows proton release, almost completely eliminates the pH dependence of O decay at pHs below 6. A second phenomenon we investigated was that in the WT at neutral and alkaline pH the fraction of the O intermediate decreases with pKa 7.5. A similar pH dependence is observed in the mutants in which the PRG is disabled, E194Q and E204Q, suggesting that the decrease in the fraction of the O intermediate with pKa ca. 7.5 is not controlled by the PRG. We propose that the group with pKa 7.5 is Asp-96. The slowing of the reprotonation of Asp-96 at high pH is the cause of the decrease in the rate of the N --> O transition, leading to the decrease in the fraction of O.     

Purification, molecular cloning, and properties of a beta-glycosidase isolated from midgut lumen of Tenebrio molitor (Coleoptera) larvae

Two beta-glycosidases (M(r) 59k) were purified from midgut contents of larvae of the yellow mealworm, Tenebrio molitor (Coleoptera: Tenebrionidae). The two enzymes (betaGly1 and betaGly2) have identical kinetic properties, but differ in hydrophobicity. The two glycosidases were cloned and their sequences differ by only four amino acids. The T. molitor glycosidases are family 1 glycoside hydrolases and have the E379 (nucleophile) and E169 (proton donor) as catalytic amino acids based on sequence alignments. The enzymes share high homology and similarity with other insect, mammalian and plant beta-glycosidases. The two enzymes may hydrolyze several substrates, such as disaccharides, arylglucosides, natural occurring plant glucosides, alkylglucosides, oligocellodextrins and the polymer laminarin. The enzymes have only one catalytic site, as inferred from experiments of competition between substrates and sequence alignments. The observed inhibition by high concentrations of the plant glucoside amygdalin, used as substrate, is an artifact generated by transglucosylation. The active site of each purified beta-glycosidase has four subsites, of which subsites +1 and +2 bind glucose with more affinity. Subsite +2 has more affinity for hydrophobic groups, binding with increasing affinities: glucose, mandelonitrile and nitrophenyl moieties. Subsite +3 has more affinity for glucose than butylene moieties. The intrinsic catalytic constant calculated for hydrolysis of the glucose beta-1,4-glucosidic bond is 21.2 s(-1) x M(-1). The putative physiological role of these enzymes is the digestion of di- and oligosaccharides derived from hemicelluloses.     

Generation of the O630 photointermediate of bacteriorhodopsin is controlled by the state of protonation of several protein residues

The last stages of the photocycle of the photosynthetic pigment all-trans bacteriorhodopsin (bR570), as well as its proton pump mechanism, are markedly pH dependent. We have measured the relative amount of the accumulated O630 intermediate (Phir), as well as its rise and decay rate constants (kr and kd, respectively), over a wide pH range. The experiments were carried out in deionized membrane suspensions to which varying concentrations of metal cations and of large organic cations were added. The observed pH dependencies, s-shaped curves in the case of Phir and bell-shaped curves for kr and kd, are interpreted in terms of the titration of three protein residues denoted as R1, R2, and R3. The R1 titration is responsible for the increase in Phir, kr, and kd upon lowering the pH from pH approximately 9.5 to 7. At low pH Phir exhibits a secondary rise which is attributed to the titration of a low pKa group, R2. After reaching a maximum at pH approximately 7, kr and kd undergo a decrease upon decreasing the pH, which is attributed to the titration of R3. All three titrations exhibit pKa values which decrease upon increasing the salt concentration. As in the case of the Purple (bR570) if Blue (bR605) equilibrium, divalent cations are substantially more effective than monovalent cations in shifting the pKa values. Moreover, bulky organic cations are as effective as small metal cations. It is concluded that analogously to the Purple if Blue equilibrium, the salt binding sites which control the pKa values of R1, R2, and R3 are located on, or close to, the membrane surface. Possible identifications of the three protein residues are considered. Experiments with the E204Q mutant show that the mutation has markedly affected the R2 (Phir) titration, suggesting that R2 should be identified with Glu-204 or with a group whose pKa is affected by Glu-204. The relation between the R1, R2 and R3 titrations and the proton pump mechanism is discussed. It is evident that the pH dependence of Phir is unrelated to the measured pKa of the group (XH) which releases the proton to the extracellular medium during the photocycle. However, since the same residue may exhibit different pKa values at different stages of the photocycle, it cannot be excluded that R2 or R3 may be identified with XH.     

Improved catalytic properties of halohydrin dehalogenase by modification of the halide-binding site

Halohydrin dehalogenase (HheC) from Agrobacterium radiobacter AD1 catalyzes the dehalogenation of vicinal haloalcohols by an intramolecular substitution reaction, resulting in the formation of the corresponding epoxide, a halide ion, and a proton. Halide release is rate-limiting during the catalytic cycle of the conversion of (R)-p-nitro-2-bromo-1-phenylethanol by the enzyme. The recent elucidation of the X-ray structure of HheC showed that hydrogen bonds between the OH group of Tyr187 and between the Odelta1 atom of Asn176 and Nepsilon1 atom of Trp249 could play a role in stabilizing the conformation of the halide-binding site. The possibility that these hydrogen bonds are important for halide binding and release was studied using site-directed mutagenesis. Steady-state kinetic studies revealed that mutant Y187F, which has lost both hydrogen bonds, has a higher catalytic activity (k(cat)) with two of the three tested substrates compared to the wild-type enzyme. Mutant W249F also shows an enhanced k(cat) value with these two substrates, as well as a remarkable increase in enantiopreference for (R)-p-nitro-2-bromo-1-phenylethanol. In case of a mutation at position 176 (N176A and N176D), a 1000-fold lower catalytic efficiency (k(cat)/K(m)) was obtained, which is mainly due to an increase of the K(m) value of the enzyme. Pre-steady-state kinetic studies showed that a burst of product formation precedes the steady state, indicating that halide release is still rate-limiting for mutants Y187F and W249F. Stopped-flow fluorescence experiments revealed that the rate of halide release is 5.6-fold higher for the Y187F mutant than for the wild-type enzyme and even higher for the W249F enzyme. Taken together, these results show that the disruption of two hydrogen bonds around the halide-binding site increases the rate of halide release and can enhance the overall catalytic activity of HheC.     

Hydrogen bonding and catalysis: a novel explanation for how a single amino acid substitution can change the pH optimum of a glycosidase

The pH optima of family 11 xylanases are well correlated with the nature of the residue adjacent to the acid/base catalyst. In xylanases that function optimally under acidic conditions, this residue is aspartic acid, whereas it is asparagine in those that function under more alkaline conditions. Previous studies of wild-type (WT) Bacillus circulans xylanase (BCX), with an asparagine residue at position 35, demonstrated that its pH-dependent activity follows the ionization states of the nucleophile Glu78 (pKa 4.6) and the acid/base catalyst Glu172 (pKa 6.7). As predicted from sequence comparisons, substitution of this asparagine residue with an aspartic acid residue (N35D BCX) shifts its pH optimum from 5.7 to 4.6, with an approximately 20% increase in activity. The bell-shaped pH-activity profile of this mutant enzyme follows apparent pKa values of 3.5 and 5.8. Based on 13C-NMR titrations, the predominant pKa values of its active-site carboxyl groups are 3.7 (Asp35), 5.7 (Glu78) and 8.4 (Glu172). Thus, in contrast to the WT enzyme, the pH-activity profile of N35D BCX appears to be set by Asp35 and Glu78. Mutational, kinetic, and structural studies of N35D BCX, both in its native and covalently modified 2-fluoro-xylobiosyl glycosyl-enzyme intermediate states, reveal that the xylanase still follows a double-displacement mechanism with Glu78 serving as the nucleophile. We therefore propose that Asp35 and Glu172 function together as the general acid/base catalyst, and that N35D BCX exhibits a "reverse protonation" mechanism in which it is catalytically active when Asp35, with the lower pKa, is protonated, while Glu78, with the higher pKa, is deprotonated. This implies that the mutant enzyme must have an inherent catalytic efficiency at least 100-fold higher than that of the parental WT, because only approximately 1% of its population is in the correct ionization state for catalysis at its pH optimum. The increased efficiency of N35D BCX, and by inference all "acidic" family 11 xylanases, is attributed to the formation of a short (2.7 A) hydrogen bond between Asp35 and Glu172, observed in the crystal structure of the glycosyl-enzyme intermediate of this enzyme, that will substantially stabilize the transition state for glycosyl transfer. Such a mechanism may be much more commonly employed than is generally realized, necessitating careful analysis of the pH-dependence of enzymatic catalysis.     

Active site of epoxide hydrolases revisited: a noncanonical residue in potato StEH1 promotes both formation and breakdown of the alkylenzyme intermediate

The carboxylate of Glu35 in the active site of potato epoxide hydrolase StEH1 interacts with the catalytic water molecule and is the first link in a chain of hydrogen bonds connecting the active site with bulk solvent. To probe its importance to catalysis, the carboxylate was replaced with an amide through an E35Q mutation. Comparing enzyme activities using the two trans-stilbene oxide (TSO) enantiomers as substrates revealed the reaction with R,R-TSO to be the one more severely affected by the E35Q mutation, as judged by determined kinetic parameters describing the pre-steady states or the steady states of the catalyzed reactions. The hydrolysis of S,S-TSO afforded by the E35Q mutant was comparable with that of the wild-type enzyme, with only a minor decrease in activity, or a change in pH dependencies of kcat, and the rate of alkylenzyme hydrolysis, k3. The pH dependence of E35Q-catalyzed hydrolysis of R,R-TSO, however, exhibited an inverted titration curve as compared to that of the wild-type enzyme, with a minimal catalytic rate at pH values where the wild-type enzyme exhibited maximum rates. To simulate the pH dependence of the E35Q mutant, a shift in the acidity of the alkylenzyme had to be invoked. The proposed decrease in the pKa of His300 in the E35Q mutant was supported by computer simulations of the active site electrostatics. Hence, Glu35 participates in activation of the Asp nucleophile, presumably by facilitating channeling of protons out of the active site, and during the hydrolysis half-reaction by orienting the catalytic water for optimal hydrogen bonding, to fine-tune the acid-base characteristics of the general base His300.     

Calculation of HyHel10-lysozyme binding free energy changes: Effect of ten point mutations

The change in free energy of binding of hen egg white lysozyme (HEL) to the antibody HyHel-10 arising from ten point mutations in HEL (D101K, D101G, K96M, K97D, K97G, K97G, R21E, R21K, W62Y, and W63Y) was calculated using a combination of the finite difference Poisson-Boltzmann method for the electrostatic contribution, a solvent accessible surface area term for the non-polar contribution, and rotamer counting for the sidechain entropy contribution. Comparison of experimental and calculated results indicate that because of pKa shifts in some of the mutated residues, primarily those involving Aspartate or Glutamate, proton uptake or release occurs in binding. When this effect was incorporated into the binding free energy calculations, the agreement with experiment improved significantly, and resulted in a mean error of about 1.9 kcal/mole. Thus these calculations predict that there should be a significant pH dependence to the change in binding caused by these mutations. The other major contributions to binding energy changes comes from solvation and charge charge interactions, which tend to oppose each other. Smaller contributions come from nonpolar interactions and sidechain entropy changes. The structures of the HyHel-10-HEL complexes with mutant HEL were obtained by modeling, and the effect of the modeled structure on the calculations was also examined. “Knowledge based” modeling and automatic generation of models using molecular mechanics produced comparable results. Proteins 33:39–48, 1998. © 1998 Wiley-Liss, Inc.     

表面电荷突变与胰蛋白酶底物专一性的改造

分别利用酶切重组和“３＋１”引物ＰＣＲ定点突变的方法构建了三个胰蛋白酶表面电荷双突变体：Ｒ６２Ｄ＋Ｋ９７Ｅ、Ｒ６２Ｄ＋Ｋ１７５Ｅ和Ｋ９７Ｅ＋Ｋ１７５Ｅ．对三者在Ｅ．ｃｏｌｉＸ－９０菌中的表达产物进行了动力学测定,分别得到了三种双突变体在两种ｐＨ条件下,水解ＴＡＭＥ、ＴＬＭＥ两种底物的动力学数据．结果表明,Ｒ６２Ｄ＋Ｋ１７５Ｅ和Ｋ９７Ｅ＋Ｋ１７５Ｅ在ｐＨ６．８５时,对两种底物的催化活性与野生型相比下降了２～３个数量级,当ｐＨ升高至８．８５时,它们的活性基本丧失;双突变体Ｒ６２Ｄ＋Ｋ９７Ｅ虽然催化活性也有所降低,但随着ｐＨ的升高,它对Ｌｙｓ底物的特异性（选择性系数２５倍于Ａｒｇ底物）转变为对Ａｒｇ底物略高的特异性,基本符合分子设计．实验结果还表明,各种双突变体催化活性的降低主要是由于酶和底物的亲和力降低引起的．