The active site of human glucocerebrosidase: structural predictions and experimental validations

Gaucher disease is a lysosomal storage disorder caused by a deficiency in glucocerebrosidase which cleaves the beta-glucosidic linkage of glucosylceramide, a normal intermediate in glycolipid metabolism. Glucocerebrosidase belongs to the clan GH-A of glycoside hydrolases, a large group of enzymes which function with retention of the anomeric configuration at the hydrolysis site. Accurate three-dimensional (3D) structure data for glucocerebrosidase should help to better understand the molecular bases of Gaucher disease. As such 3D structure data were not available, we used the two-dimensional hydrophobic cluster analysis (HCA) method to make structure predictions for the catalytic domains of clan GH-A glycoside hydrolases. We found that all the enzymes of clan GH-A may share a similar catalytic domain consisting of an (alpha/beta)8 barrel with the critical acid/base and nucleophile residues located at the C-terminal ends of strands beta 4 and beta 7, respectively. In the case of glucocerebrosidase, Glu 235 was predicted to be the putative acid/base catalyst whereas the nucleophile was located at Glu 340. Next, in order to obtain experimental evidence supporting these HCA-based predictions, we used retroviral vectors to express, in murine null cells, E235A and E340A mutant proteins, in which alanine residues unable to participate in the enzymatic reaction replace the presumed critical glutamic acid residues. Both mutants were found to be catalytically inactive although they were correctly folded/processed and sorted to the lysosome. Thus, Glu 235 and Glu 340 do indeed play key roles in the active site of human glucocerebrosidase as predicted by the HCA analysis. In a broader perspective, our work points out that bioinformatics approaches may be highly useful for generating structure-function predictions based on sequence-structure interrelationships, especially in the context of a rapid increase in protein sequence information through genome sequencing.     

Investigation of the functional relevance of the catalytically important Glu(28) in family 51 arabinosidases

The alpha-L-arabinofuranosidase (AbfD3) from Thermobacillus xylanilyticus is a family 51 glycosyl hydrolase. According to classification hierarchy, family 51 belongs to clan GH-A. While the major GH-A motifs, the catalytic acid-base and nucleophile, are conserved in AbfD3, a third catalytically important residue (Glu(28)) does not appear to be analogous to any known GH-A motif. To evaluate the importance of Glu(28), bioinformatics analyses and site-saturation mutagenesis were performed. The results indicate that Glu(28) forms part of a family 51 arabinosidase motif which might be functionally homologous to a conserved N-terminal motif found in exo-acting enzymes from families 1 and 5. Importantly, the data reveal that Glu(28) is a key determinant of substrate recognition in the -1 subsite, where it may also play an important role in water-mediated deglycosylation of the glycosyl-enzyme covalent intermediate.     

Structure modeling, ligand binding, and binding affinity calculation (LR-MM-PBSA) of human heparanase for inhibition and drug design

Human heparanase is an endo-beta-D-glycosidase that cleaves heparan sulphate (HS) chains in the extracellular matrix and basement membrane. It is known that the cleavage of HS by heparanase results in cell invasion and metastasis of cancer. Therefore, heparanase is considered an important target for cancer drug development. The three-dimensional structure of heparanase would be useful in the rational design of inhibitors targeted to the enzyme; however, the three-dimensional structure has not yet been determined. In our effort to design inhibitors, we developed a three-dimensional structure of heparanase using a homology-modeling approach. The homology-built structure is consistent to previous bioinformatics and site-mutation experimental results. The heparanase features a (alpha/beta)(8) TIM-barrel fold with two glutamate residues (Glu225 and Glu343) located in the active-site cleft. This feature supports the putative mechanism of proton donor and nucleophilic sites. Docking simulations yielded 41 complex structures, which indicate that the bound inhibitor could block ligand binding into the catalytic site. A free energy of binding model was established for 25 heparanase inhibitors with a training set of 25 heparanase inhibitors using the linear response MM-PBSA approach (LR-MM-PBSA). The correlation between calculated and experimental activity was 0.79 and the reliability of the model was validated with leave-one-out cross-validation method. Its predictive capability was further validated using a test set of 16 inhibitors similar to the training set of inhibitors. The correlation between the predicted and observed activities is significantly improved by the protein "induced-fit" that accounts for the flexibility of the receptor. These interaction and pharmacophore elements provide a unique insight to the rational design of new ligands targeted to the enzyme.     

Structures of chitobiase mutants complexed with the substrate Di-N-acetyl-d-glucosamine: the catalytic role of the conserved acidic pair, aspartate 539 and glutamate 540

The catalytic domain of chitobiase (beta-N-1-4 acetylhexosaminidase) from Serratia marcescens, is an alpha/beta TIM-barrel. This enzyme belongs to family 20 of glycosyl hydrolases in which a conserved amino acid pair, aspartate-glutamate, is present (Asp539-Glu540). It was proposed that catalysis by this enzyme family is carried out by glutamate 540 acting as a proton donor and by the acetamido group of the substrate as a nucleophile. We investigated the role of Asp539 and Glu540 by site-directed mutagenesis, biochemical characterization and by structural analyses of chitobiase -substrate co-crystals. We found that both residues are essential for chitobiase activity. The mutations, however, led to subtle changes in the catalytic site. Our results support the model that Glu540 acts as the proton donor and that Asp539 acts in several different ways. Asp539 restrains the acetamido group of the substrate in a specific orientation by forming a hydrogen bond with N2 of the non-reduced (-1) sugar. In addition, this residue participates in substrate binding. It is also required for the correct positioning of Glu540 and may provide additional negative charge at the active site. Thus, these biochemical and structural studies provide a molecular explanation for the functional importance and conservation of these residues.     

Crystal structure of mannanase 26A from Pseudomonas cellulosa and analysis of residues involved in substrate binding

The crystal structure of Pseudomonas cellulosa mannanase 26A has been solved by multiple isomorphous replacement and refined at 1.85 A resolution to an R-factor of 0.182 (R-free = 0.211). The enzyme comprises (beta/alpha)(8)-barrel architecture with two catalytic glutamates at the ends of beta-strands 4 and 7 in precisely the same location as the corresponding glutamates in other 4/7-superfamily glycoside hydrolase enzymes (clan GH-A glycoside hydrolases). The family 26 glycoside hydrolases are therefore members of clan GH-A. Functional analyses of mannanase 26A, informed by the crystal structure of the enzyme, provided important insights into the role of residues close to the catalytic glutamates. These data showed that Trp-360 played a critical role in binding substrate at the -1 subsite, whereas Tyr-285 was important to the function of the nucleophile catalyst. His-211 in mannanase 26A does not have the same function as the equivalent asparagine in the other GH-A enzymes. The data also suggest that Trp-217 and Trp-162 are important for the activity of mannanase 26A against mannooligosaccharides but are less important for activity against polysaccharides.     

Biochemical and structural characterization of the intracellular mannanase AaManA of Alicyclobacillus acidocaldarius reveals a novel glycoside hydrolase family belonging to clan GH-A

An intracellular mannanase was identified from the thermoacidophile Alicyclobacillus acidocaldarius Tc-12-31. This enzyme is particularly interesting, because it shows no significant sequence similarity to any known glycoside hydrolase. Gene cloning, biochemical characterization, and structural studies of this novel mannanase are reported in this paper. The gene consists of 963 bp and encodes a 320-amino acid protein, AaManA. Based on its substrate specificity and product profile, AaManA is classified as an endo-beta-1,4-mannanase that is capable of transglycosylation. Kinetic analysis studies revealed that the enzyme required at least five subsites for efficient hydrolysis. The crystal structure at 1.9 angstroms resolution showed that AaManA adopted a (beta/alpha)8-barrel fold. Two catalytic residues were identified: Glu151 at the C terminus of beta-stand beta4 and Glu231 at the C terminus of beta7. Based on the structure of the enzyme and evidence of its transglycosylation activity, AaManA is placed in clan GH-A. Superpositioning of its structure with that of other clan GH-A enzymes revealed that six of the eight GH-A key residues were functionally conserved in AaManA, with the exceptions being residues Thr95 and Cys150. We propose a model of substrate binding in AaManA in which Glu282 interacts with the axial OH-C(2) in-2 subsites. Based on sequence comparisons, the enzyme was assigned to a new glycoside hydrolase family (GH113) that belongs to clan GH-A.     

Identification of a catalytic residue of Clostridium paraputrificum N-acetyl-beta-D-glucosaminidase Nag3A by site-directed mutagenesis

Clostridium paraputrificum M-21 beta-N-acetylglucosaminidase 3A (Nag3A) is an enzyme classified in family 3 of the glycoside hydrolases. To identify catalytic residues of this enzyme, mutations were introduced into highly conserved Glu and Asp residues. Replacement of Asp175 with Ala abolished the catalytic activity without change in the circular dichroism spectrum, strongly suggesting that this residue is a catalytic residue, a nucleophile/base or a proton donor. Since the K(m) values of mutant enzymes D119N, D229N, D229A and D274N increased 17 to 41 times as compared with that of wild-type enzyme, Asp119, Asp229, and Asp274 appear to be involved in substrate recognition and binding. Taking previous studies into consideration, we presume that Asp303 is the catalytic nucleophile and Asp175 is the proton donor of C. paraputrificum Nag3A.     

Active-site motifs of lysosomal acid hydrolases: invariant features of clan GH-A glycosyl hydrolases deduced from hydrophobic cluster analysis

The clan GH-A is a group of more than 200 proteins representingnine established families of glycosyl hydrolases that act ona large variety of substrates. This clan includes five enzymesimplicated in lysosomal storage diseases: ß-glucuronidase(Sly disease), ß-glucocerebrosidase (Gau-cher disease),ß-galactosidase (Landing disease and Morquio typeB disease), ß-mannosidase (mannosidosis) and     

The crystal structure of human cytosolic beta-glucosidase unravels the substrate aglycone specificity of a family 1 glycoside hydrolase

Human cytosolic beta-glucosidase (hCBG) is a xenobiotic-metabolizing enzyme that hydrolyses certain flavonoid glucosides, with specificity depending on the aglycone moiety, the type of sugar and the linkage between them. In this study, the substrate preference of this enzyme was investigated by mutational analysis, X-ray crystallography and homology modelling. The crystal structure of hCBG was solved by the molecular replacement method and refined at 2.7 A resolution. The main-chain fold of the enzyme belongs to the (beta/alpha)(8) barrel structure, which is common to family 1 glycoside hydrolases. The active site is located at the bottom of a pocket (about 16 A deep) formed by large surface loops, surrounding the C termini of the barrel of beta-strands. As for all the clan of GH-A enzymes, the two catalytic glutamate residues are located on strand 4 (the acid/base Glu165) and on strand 7 (the nucleophile Glu373). Although many features of hCBG were shown to be very similar to previously described enzymes from this family, crucial differences were observed in the surface loops surrounding the aglycone binding site, and these are likely to strongly influence the substrate specificity. The positioning of a substrate molecule (quercetin-4'-glucoside) by homology modelling revealed that hydrophobic interactions dominate the binding of the aglycone moiety. In particular, Val168, Trp345, Phe225, Phe179, Phe334 and Phe433 were identified as likely to be important in determining substrate specificity in hCBG, and site-directed mutagenesis supported a key role for some of these residues.     

Bioinformatics of the glycoside hydrolase family 57 and identification of catalytic residues in amylopullulanase from Thermococcus hydrothermalis.

Fifty-nine amino acid sequences belonging to family 57 (GH-57) of the glycoside hydrolases were collected using the CAZy server, Pfam database and blast tools. Owing to the sequence heterogeneity of the GH-57 members, sequence alignments were performed using mainly manual methods. Likewise, five conserved regions were identified, which are postulated to be GH-57 consensus motifs. In the 659 amino acid-long 4-alpha-glucanotransferase from Thermococcus litoralis, these motifs correspond to 13_HQP (region I), 76_GQLEIV (region II), 120_WLTERV (region III), 212_HDDGEKFGVW (region IV), and 350_AQCNDAYWH (region V). The third and fourth conserved regions contain the catalytic nucleophile and the proton donor, respectively. Based on our sequence alignment, residues Glu291 and Asp394 were proposed as the nucleophile and proton donor, respectively, in a GH-57 amylopullulanase from Thermococcus hydrothermalis. To validate this prediction, site-directed mutagenesis was performed. The results of this work reveal that both residues are critical for the pullulanolytic and amylolytic activities of the amylopullulanase. Therefore, these data support the prediction and strongly suggest that the bifunctionality of the amylopullulanase is determined by a single catalytic centre. Despite this positive validation, our alignment also reveals that certain GH-57 members do not possess the Glu and Asp corresponding to the predicted GH-57 catalytic residues. However, the sequences concerned by this anomaly encode putative proteins for which no biochemical or enzymatic data are yet available. Finally, the evolutionary trees generated for GH-57 reveal that the entire family can be divided into several subfamilies that may reflect the different enzyme specificities.     

Glutamic acid 160 is the acid-base catalyst of beta-xylosidase from Bacillus stearothermophilus T-6: a family 39 glycoside hydrolase

A beta-xylosidase from Bacillus stearothermophilus T-6 was cloned, overexpressed in Escherichia coli and purified to homogeneity. Based on sequence alignment, the enzyme belongs to family 39 glycoside hydrolases, which itself forms part of the wider GH-A clan. The conserved Glu160 was proposed as the acid-base catalyst. An E160A mutant was constructed and subjected to steady state and pre-steady state kinetic analysis together with azide rescue and pH activity profiles. The observed results support the assignment of Glu160 as the acid-base catalytic residue.     

Characterization and evolution of a myrosinase from the cabbage aphid Brevicoryne brassicae.

The aphid myrosinase gene has been elucidated using Rapid Amplification of cDNA Ends-PCR. Sequencing has shown that aphid myrosinase has significant sequence similarity (35%) to plant myrosinases and other members of glycosyl hydrolase family 1 (GHF1). The residues acting as proton donor and nucleophile, in the hydrolysis of glucosinolates by aphid myrosinase, are identified as Glu 167 and Glu 374 respectively. The equivalent residues in plant myrosinase are Gln 187 and Glu 409 and for the cyanogenic beta-glucosidase Glu 183 and Glu 397. Thus it would appear that the absence of a proton donor is not necessary for the hydrolysis of glucosinolates as was thought to be the case for the plant myrosinases. Aphid myrosinase appears to be more similar to animal beta-O-glucosidases than to plant myrosinases, as assessed by sequence similarity and phylogenetic techniques. These results strongly suggest that myrosinase activity has twice arisen from beta-O-glucosidases in plants and animals. Comparison of aphid myrosinase with plant myrosinase has highlighted Lys 173 and Arg 312 as possibly playing a crucial role in the hydrolysis of glucosinolates by aphid myrosinase.     

Evidence that the methylesterase of bacterial chemotaxis may be a serine hydrolase.

CheB, the methylesterase of chemotactic bacteria, catalyzes the hydrolysis of glutamyl-methyl esters in bacterial chemoreceptor proteins. The two cysteines predicted by the amino acid sequence of CheB were replaced by alanine residues. The resulting mutants, Cys207-Ala, Cys309-Ala and a double cysteine mutant Cys207-Ala/Cys309-Ala, retained methylesterase activity, indicating that sulfhydryls are not crucial for CheB mediated catalysis. A homology search revealed a conserved serine active-site region between residues 162 and 166 which is homologous to the active-site region of acetylcholine esterases, suggesting that Ser164 of CheB is the active-site nucleophile. Oligonucleotide-directed mutagenesis was used to change the serine to a cysteine. This Ser164-Cys mutant had less than 2% of the wild-type activity. Unlike the serine proteinases which utilize a 'catalytic triad' mechanism, CheB does not have the conserved histidine and aspartic acid residues located in positions N-terminal to the active-site serine. In addition, CheB is not labeled with di-isopropylfluorophosphate, a potent inhibitor of other serine hydrolases. A novel mechanism is proposed for CheB involving substrate-assisted catalysis to account for these apparent anomalies.     

pKa modulation of the acid/base catalyst within GH32 and GH68: a role in substrate/inhibitor specificity?

Glycoside hydrolases of families 32 (GH32) and 68 (GH68) belong to clan GH-J, containing hydrolytic enzymes (sucrose/fructans as donor substrates) and fructosyltransferases (sucrose/fructans as donor and acceptor substrates). In GH32 members, some of the sugar substrates can also function as inhibitors, this regulatory aspect further adding to the complexity in enzyme functionalities within this family. Although 3D structural information becomes increasingly available within this clan and huge progress has been made on structure-function relationships, it is not clear why some sugars bind as inhibitors without being catalyzed. Conserved aspartate and glutamate residues are well known to act as nucleophile and acid/bases within this clan. Based on the available 3D structures of enzymes and enzyme-ligand complexes as well as docking simulations, we calculated the pKa of the acid-base before and after substrate binding. The obtained results strongly suggest that most GH-J members show an acid-base catalyst that is not sufficiently protonated before ligand entrance, while the acid-base can be fully protonated when a substrate, but not an inhibitor, enters the catalytic pocket. This provides a new mechanistic insight aiming at understanding the complex substrate and inhibitor specificities observed within the GH-J clan. Moreover, besides the effect of substrate entrance on its own, we strongly suggest that a highly conserved arginine residue (in the RDP motif) rather than the previously proposed Tyr motif (not conserved) provides the proton to increase the pKa of the acid-base catalyst.     

Identification by mutagenesis of a conserved glutamate (Glu487) residue important for catalytic activity in rat liver carnitine palmitoyltransferase II

Mammalian mitochondrial membranes express two active but distinct carnitine palmitoyltransferases: carnitine palmitoyltransferase I (CPTI), which is malonyl coA-sensitive and detergent-labile; and carnitine palmitoyltransferase II (CPTII), which is malonyl coA-insensitive and detergent-stable. To determine the role of the highly conserved C-terminal acidic residues glutamate 487 (Glu(487)) and glutamate 500 (Glu(500)) on catalytic activity in rat liver CPTII, we separately mutated these residues to alanine, aspartate, or lysine, and the effect of the mutations on CPTII activity was determined in the Escherichia coli-expressed mutants. Substitution of Glu(487) with alanine, aspartate, or lysine resulted in almost complete loss in CPTII activity. Because a conservative substitution mutation of this residue, Glu(487) with aspartate (E487D), resulted in a 97% loss in activity, we predicted that Glu(487) would be at the active-site pocket of CPTII. The substantial loss in CPTII activity observed with the E487K mutant, along with the previously reported loss in activity observed in a child with a CPTII deficiency disease, establishes that Glu(487) is crucial for maintaining the configuration of the liver isoform of the CPTII active site. Substitution of the conserved Glu(500) in CPTII with alanine or aspartate reduced the V(max) for both substrates, suggesting that Glu(500) may be important in stabilization of the enzyme-substrate complex. A conservative substitution of Glu(500) to aspartate resulted in a significant decrease in the V(max) for the substrates. Thus, Glu(500) may play a role in substrate binding and catalysis. Our site-directed mutagenesis studies demonstrate that Glu(487) in the liver isoform of CPTII is essential for catalysis.     

Crystal structures of beta-galactosidase from Penicillium sp. and its complex with galactose

Beta-galactosidases catalyze the hydrolysis of beta(1-3) and beta(1-4) galactosyl bonds in oligosaccharides as well as the inverse reaction of enzymatic condensation and transglycosylation. Here we report the crystallographic structures of Penicillium sp. beta-galactosidase and its complex with galactose solved by the SIRAS quick cryo-soaking technique at 1.90 A and 2.10 A resolution, respectively. The amino acid sequence of this 120 kDa protein was first assigned putatively on the basis of inspection of the experimental electron density maps and then determined by nucleotide sequence analysis. Primary structure alignments reveal that Penicillium sp. beta-galactosidase belongs to family 35 of glycosyl hydrolases (GHF-35). This model is the first 3D structure for a member of GHF-35. Five distinct domains which comprise the structure are assembled in a way previously unobserved for beta-galactosidases. Superposition of this complex with other beta-galactosidase complexes from several hydrolase families allowed the identification of residue Glu200 as the proton donor and residue Glu299 as the nucleophile involved in catalysis. Penicillium sp. beta-galactosidase is a glycoprotein containing seven N-linked oligosaccharide chains and is the only structure of a glycosylated beta-galactosidase described to date.     

Mannose foraging by Bacteroides thetaiotaomicron: structure and specificity of the beta-mannosidase, BtMan2A

The human colonic bacterium Bacteroides thetaiotaomicron, which plays an important role in maintaining human health, produces an extensive array of exo-acting glycoside hydrolases (GH), including 32 family GH2 glycoside hydrolases. Although it is likely that these enzymes enable the organism to utilize dietary and host glycans as major nutrient sources, the biochemical properties of these GH2 glycoside hydrolases are currently unclear. Here we report the biochemical properties and crystal structure of the GH2 B. thetaiotaomicron enzyme BtMan2A. Kinetic analysis demonstrates that BtMan2A is a beta-mannosidase in which substrate binding energy is provided principally by the glycone binding site, whereas aglycone recognition is highly plastic. The three-dimensional structure, determined to a resolution of 1.7 A, reveals a five-domain structure that is globally similar to the Escherichia coli LacZ beta-galactosidase. The catalytic center is housed mainly within a (beta/alpha)8 barrel although the N-terminal domain also contributes to the active site topology. The nature of the substrate-binding residues is quite distinct from other GH2 enzymes of known structure, instead they are similar to other clan GH-A enzymes specific for manno-configured substrates. Mutagenesis studies, informed by the crystal structure, identified a WDW motif in the N-terminal domain that makes a significant contribution to catalytic activity. The observation that this motif is invariant in GH2 mannosidases points to a generic role for these residues in this enzyme class. The identification of GH-A clan and GH2 specific residues in the active site of BtMan2A explains why this enzyme is able to harness substrate binding at the proximal glycone binding site more efficiently than mannan-hydrolyzing glycoside hydrolases in related enzyme families. The catalytic properties of BtMan2A are consistent with the flexible nutrient acquisition displayed by the colonic bacterium.     

In silico identification of catalytic residues and domain fold of the family GH119 sharing the catalytic machinery with the α-amylase family GH57

The glycoside hydrolase family 119 (GH119) contains the α-amylase from Bacillus circulans and five other hypothetical proteins. Until now, nothing has been reported on the catalytic residues and catalytic-domain fold of GH119. Based on a detailed in silico analysis involving sequence comparison in combination with BLAST searches and structural modelling, an unambiguous relationship was revealed between the families GH119 and GH57. This includes sharing the catalytic residues, i.e. Glu231 and Asp373 as catalytic nucleophile and proton donor, respectively, in the predicted catalytic (β/α)₇-barrel domain of GH119 B. circulans α-amylase. The GH57 and GH119 families may thus define a new CAZy clan.     

Molecular dynamics simulations of a branched tetradecasaccharide substrate in the active site of a xyloglucan endo-transglycosylase

Molecular dynamics simulations of the tetradecasaccharide XXXGXXXG in complex with the hybrid aspen xyloglucan endo-transglycosylase PttXET16-34 have been performed and analysed with respect to structure, dynamics, flexibility and ligand interactions. Notably, the charge state of the so-called ‘helper residue’ aspartate 87 (Asp87), which lies between the catalytic nucleophile [glutamate 85 (Glu85)] and general acid/base (Glu89) residues on the same beta strand, had a significant effect on PttXET16-34 active site structure. When Asp87 was deprotonated, electrostatic repulsion forced the nucleophile away from C1 of the sugar ring in subsite ? 1 and the proton–donating ability of Glu89 was also weakened due to the formation of a hydrogen bond with Asp87, whereas the protonation of Asp87 resulted in the formation of a hydrogen bond with the catalytic nucleophile and correct positioning of the catalytic machinery. The results suggest that catalysis in glycoside hydrolase family 16, and by extension clan GH-B enzymes, is optimal when the catalytic nucleophile is deprotonated for nucleophilic attack on the substrate, whereas the ‘helper residue’ and general acid/base residue are both in their conjugate acid forms to align the nucleophile and deliver a proton to the departing sugar, respectively.     

Crystal structure of Cel44A, a glycoside hydrolase family 44 endoglucanase from Clostridium thermocellum

The crystal structure of Cel44A, which is one of the enzymatic components of the cellulosome of Clostridium thermocellum, was solved at a resolution of 0.96 A. This enzyme belongs to glycoside hydrolase family (GH family) 44. The structure reveals that Cel44A consists of a TIM-like barrel domain and a beta-sandwich domain. The wild-type and the E186Q mutant structures complexed with substrates suggest that two glutamic acid residues, Glu(186) and Glu(359), are the active residues of the enzyme. Biochemical experiments were performed to confirm this idea. The structural features indicate that GH family 44 belongs to clan GH-A and that the reaction catalyzed by Cel44A is retaining type hydrolysis. The stereochemical course of hydrolysis was confirmed by a (1)H NMR experiment using the reduced cellooligosaccharide as a substrate.