GHL1-GHL15: New families of the hypothetical glycoside hydrolases

The domains of 15 recently discovered families of the hypothetical glycoside hydrolases GHL1-GHL15 were used for iterative screening of the protein database. The evolutionary relationships between these families were revealed, as well as their relationship with the previously known families of protein domains: GH5, GH13, GH13-33, GH17, GH18, GH20, GH27, GH29, GH31, GH35, GH36A, GH36B, GH36C, GH36D, GH36E, GH36F, GH36G, GH36H, GH36J, GH36K, GH39, GH42, GH53, GH66, GH97, GH101, GH107, GH112, GH114, COG1082, COG1306, COG1649, COG2342, DUF3111, and PF00962. The unclassified homologues were grouped in 35 new families of the hypothetical glycoside hydrolases: GHL16-GHL50. The position of the families GHL1-GHL15 in the hierarchical classification of glycoside hydrolases and their homologues is discussed. Several new superfamilies of protein domains are proposed.     

Endo-α-1,4-polygalactosaminidases and their homologs: Structure and evolution

Endo-α-1,4-polygalactosaminidase is a rare enzyme. Its catalytic domain belongs to the GH114 family of glycoside hydrolases. It is shown by phylogenetic analysis that the evolution of the corresponding genes involved duplications, elimination, and horizontal transfer. The domain and secondary structures of endo-α-1,4-polygalactosaminidases are discussed. A hypothesis is put forward as to the structure of the active center of the enzyme. Iterative screening of a protein database reveals evolutionary relationships of the GH114 family with the GH13, GH18, GH20, GH27, GH29, GH31, GH35, GH36, and GH66 families of glycoside hydrolases and with the COG1306, COG1649, COG2342, GHL3, and GHL4 families of proteins with unknown enzymatic functions. Unclassified homologs are grouped into 13 new families of hypothetical glycoside hydrolases: GHL5-GHL15, GH36J, and GH36K.     

Endo-alpha-1-4-polygalactosaminidases and their homologues: structure and evolution

Endo-alpha-1,4-polygalactosaminidase is a rare enzyme. Its catalytic domain belongs to the GH114 family of glycoside hydrolases. Phylogenetic analysis of the family proteins allowed us to show an important role of duplications, eliminations, and horizontal transfer in the evolution of their genes. Domain structure, the secondary structure, and proposed structure of the active center of the endo-alpha-1,4-polygalactosaminidases are discussed. Evolutionary connections of the GH114 family with GH13, GH18, GH20, GH27, GH29, GH31, GH35, GH36, and GH66 families of glycoside hydrolases, as well as, with COG1306, COG1649, COG2342, GHL3, and GHL4 families of enzymatically uncharacterized proteins have been revealed by iterative screening of the protein database. The unclassified homologues have been grouped into 13 new families of hypothetical glycoside hydrolases: GHL5 - GHL15, GH36J, and GH36K.     

Furanosidase superfamily: Search of homologues

The furanosidase superfamily contains the GH32, GH43, GH62, GH68, GH117, DUF377 (GH130), and DUF1861 families of glycoside hydrolases and their homologues. Catalytic domains of these families have five-bladed β-propeller tertiary structure. Iterative screening of the protein database supports of their relationship as well as evolutionary connections with domains from GH33 and GH93 families of glycoside hydrolases. The latter two have the structure of the six-bladed β-propeller. Among detected homologues we found 441 unclassified proteins. These proteins are combined into 39 groups based on homology: FURAN1-FURAN39. FURAN8 and FURAN36 can be considered as separate subfamilies within the GH43 and GH32 families of glycoside hydrolases, respectively. The remaining 37 groups are new families of hypothetical glycoside hydrolases.     

Furanosidase superfamily: search of homologues

The furanosidase superfamily contains GH32, GH43, GH62, GH68, GH117, DUF377, and DUF1861 families of glycoside hydrolases and their homologues. Catalytic domains of these families have five-bladed beta-propeller tertiary structure. Iterative screening of the protein database allowed to support their relationship as well as evolutionary connections with domains from GH33 and GH93 families of glycoside hydrolases. The latter two have structure of the six-bladed beta-propeller. Among revealed homologues we found 441 unclassified proteins. These proteins are combined into 39 groups based on homology: FURAN1-FURAN39. FURAN8 and FURAN36 can be considered as separate subfamilies within GH43 and GH32 families of glycoside hydrolases, respectively. The remaining 37 groups are new families of hypothetical glycoside hydrolases.     

GH10 family of glycoside hydrolases: Structure and evolutionary connections

Evolutionary connections were analyzed for endo-β-xylanases, which possess the GH10 family catalytic domains. A homology search yielded thrice as many proteins as are available from the Carbohydrate-Active Enzymes (CAZy) database. Lateral gene transfer was shown to play an important role in evolution of bacterial proteins of the family, especially in the phyla Acidobacteria, Cyanobacteria, Planctomycetes, Spirochaetes, and Verrucomicrobia. In the case of Verrucomicrobia, 23 lateral transfers from organisms of other phyla were detected. Evolutionary relationships were observed between the GH10 family domains and domains with the TIM-barrel tertiary structure from several other glycosidase families. The GH39 family of glycoside hydrolases showed the closest relationship. Unclassified homologs were grouped into 12 novel families of putative glycoside hydrolases (GHL51–GHL62).     

Structural and evolutionary aspects of two families of non-catalytic domains present in starch and glycogen binding proteins from microbes, plants and animals

Starch-binding domains (SBDs) comprise distinct protein modules that bind starch, glycogen or related carbohydrates and have been classified into different families of carbohydrate-binding modules (CBMs). The present review focuses on SBDs of CBM20 and CBM48 found in amylolytic enzymes from several glycoside hydrolase (GH) families GH13, GH14, GH15, GH31, GH57 and GH77, as well as in a number of regulatory enzymes, e.g., phosphoglucan, water dikinase-3, genethonin-1, laforin, starch-excess protein-4, the β-subunit of AMP-activated protein kinase and its homologues from sucrose non-fermenting-1 protein kinase SNF1 complex, and an adaptor-regulator related to the SNF1/AMPK family, AKINβγ. CBM20s and CBM48s of amylolytic enzymes occur predominantly in the microbial world, whereas the non-amylolytic proteins containing these modules are mostly of plant and animal origin. Comparison of amino acid sequences and tertiary structures of CBM20 and CBM48 reveals the close relatedness of these SBDs and, in some cases, glycogen-binding domains (GBDs). The families CBM20 and CBM48 share both an ancestral form and the mode of starch/glycogen binding at one or two binding sites. Phylogenetic analyses demonstrate that they exhibit independent behaviour, i.e. each family forms its own part in an evolutionary tree, with enzyme specificity (protein function) being well represented within each family. The distinction between CBM20 and CBM48 families is not sharp since there are representatives in both CBM families that possess an intermediate character. These are, for example, CBM20s from hypothetical GH57 amylopullulanase (probably lacking the starch-binding site 2) and CBM48s from the GH13 pullulanase subfamily (probably lacking the starch/glycogen-binding site 1). The knowledge gained concerning the occurrence of these SBDs and GBDs through the range of taxonomy will support future experimental research.     

Multiple lateral transfers and duplications of genes as sources of diversity of α-L-rhamnosidases in Clostridium methylpentosum DSM5476

α-L-Rhamnosidases are an important group of glycoside hydrolases represented in many organisms from various prokaryotic phyla. Based on the homology of catalytic domains, all these proteins are assigned to the GH78 and GH106 families of glycoside hydrolases. However, most prokaryotic genomes contain no genes encoding proteins from these two families. We found that the unique genome of Clostridium methylpentosum DSM5476 contains 83 genes of proteins from these families and undertook investigation of their phylogeny. The absence of homologous genes in most of strains of the genus Clostridium suggests an important ecological role of these genes, in C. methylpentosum in particular. Phylogenetic analysis revealed multiple lateral transfers and duplications of the corresponding genes.     

PSI protein classifier: A new program automating PSI-BLAST search results

A new program, PSI Protein Classifier, generalizing the results of both successive and independent iterations of the PSI-BLAST program was developed. The technical opportunities of the program are described and illustrated by two examples. An iterative screening of the amino acid sequence database detected potential evolutionary relationships between GH5, GH13, GH27, GH31, GH36, GH66, GH101 and GH114 families of glycoside hydrolases. Analysis of the statistically significant sequence similarity (E-value analysis) allowed us to divide the family GH31 into 38 subfamilies.     

Function and structure studies of GH family 31 and 97 α-glycosidases

A huge number of glycoside hydrolases are classified into the glycoside hydrolase family (GH family) based on their amino-acid sequence similarity. The glycoside hydrolases acting on α-glucosidic linkage are in GH family 4, 13, 15, 31, 63, 97, and 122. This review deals mainly with findings on GH family 31 and 97 enzymes. Research on two GH family 31 enzymes is described: clarification of the substrate recognition of Escherichia coli α-xylosidase, and glycosynthase derived from Schizosaccharomyces pombe α-glucosidase. GH family 97 is an aberrant GH family, containing inverting and retaining glycoside hydrolases. The inverting enzyme in GH family 97 displays significant similarity to retaining α-glycosidases, including GH family 97 retaining α-glycosidase, but the inverting enzyme has no catalytic nucleophile residue. It appears that a catalytic nucleophile has been eliminated during the molecular evolution in the same way as a man-made nucleophile mutant enzyme, which catalyzes the inverting reaction, as in glycosynthase and chemical rescue.     

Function and Structure Studies of GH Family 31 and 97 α-Glycosidases

A huge number of glycoside hydrolases are classified into the glycoside hydrolase family (GH family) based on their amino-acid sequence similarity. The glycoside hydrolases acting on α-glucosidic linkage are in GH family 4, 13, 15, 31, 63, 97, and 122. This review deals mainly with findings on GH family 31 and 97 enzymes. Research on two GH family 31 enzymes is described: clarification of the substrate recognition of Escherichia coli α-xylosidase, and glycosynthase derived from Schizosaccharomyces pombe α-glucosidase. GH family 97 is an aberrant GH family, containing inverting and retaining glycoside hydrolases. The inverting enzyme in GH family 97 displays significant similarity to retaining α-glycosidases, including GH family 97 retaining α-glycosidase, but the inverting enzyme has no catalytic nucleophile residue. It appears that a catalytic nucleophile has been eliminated during the molecular evolution in the same way as a man-made nucleophile mutant enzyme, which catalyzes the inverting reaction, as in glycosynthase and chemical rescue.     

Structural and functional analysis of a glycoside hydrolase family 97 enzyme from Bacteroides thetaiotaomicron

SusB, an 84-kDa alpha-glucoside hydrolase involved in the starch utilization system (sus) of Bacteroides thetaiotaomicron, belongs to glycoside hydrolase (GH) family 97. We have determined the enzymatic characteristics and the crystal structures in free and acarbose-bound form at 1.6A resolution. SusB hydrolyzes the alpha-glucosidic linkage, with inversion of anomeric configuration liberating the beta-anomer of glucose as the reaction product. The substrate specificity of SusB, hydrolyzing not only alpha-1,4-glucosidic linkages but also alpha-1,6-, alpha-1,3-, and alpha-1,2-glucosidic linkages, is clearly different from other well known glucoamylases belonging to GH15. The structure of SusB was solved by the single-wavelength anomalous diffraction method with sulfur atoms as anomalous scatterers using an in-house x-ray source. SusB includes three domains as follows: the N-terminal, catalytic, and C-terminal domains. The structure of the SusB-acarbose complex shows a constellation of carboxyl groups at the catalytic center; Glu532 is positioned to provide protonic assistance to leaving group departure, with Glu439 and Glu508 both positioned to provide base-catalyzed assistance for inverting nucleophilic attack by water. A structural comparison with other glycoside hydrolases revealed significant similarity between the catalytic domain of SusB and those of alpha-retaining glycoside hydrolases belonging to GH27, -36, and -31 despite the differences in catalytic mechanism. SusB and the other retaining enzymes appear to have diverged from a common ancestor and individually acquired the functional carboxyl groups during the process of evolution. Furthermore, sequence comparison of the active site based on the structure of SusB indicated that GH97 included both retaining and inverting enzymes.     

β-Mannanase (Man26A) and α-galactosidase (Aga27A) synergism – A key factor for the hydrolysis of galactomannan substrates

This study investigated the behavior of mannan-degrading enzymes, specifically focusing on differences with respect to their substrate specificities and their synergistic associations with enzymes from different glycoside hydrolase (GH) families. Galactosidases from Cyamopsis tetragonolobus seeds (Aga27A, GH27) and Aspergillus niger (AglC, GH36) were evaluated for their abilities to synergistically interact with mannanases from Clostridium cellulovorans (ManA, GH5) and A. niger (Man26A, GH26) in hydrolysis of guar gum and locust bean gum. Among the mannanases, Man26A was more efficient at hydrolyzing both galactomannan substrates, while among the galactosidases; Aga27A was the most effective at removing galactose substituents on both galactomannan substrates and galactose-containing oligosaccharides. An optimal protein mass ratio of glycoside hydrolases required to maximize the release of both reducing sugar and galactose residues was determined. Clear synergistic enhancement of locust bean gum hydrolysis with respect to reducing sugar release was observed when both mannanases at 75% enzyme dosage were supplemented with 25% enzyme protein dosage of Aga27A. At a protein ratio of 75% Man26A to 25% Aga27A, the presence of Man26A significantly enhanced galactose release by 25% Aga27A (2.36 fold) with locust bean gum, compared to when Aga27A was used alone at 100% enzyme protein dosage. A dosage of Aga27A at 75% and ManA at 25% protein content liberated the highest reducing sugar release on guar gum hydrolysis. A dosage of Man26A and Aga27A at 75–25% protein content, respectively, liberated reducing sugar release equivalent to that when Man26A was used alone at 100% protein content. From the findings obtained in this study, it was observed that the GH family classification of an enzyme affects its substrate specificity and synergistic interactions with other glycoside hydrolases from different families (more so than its EC classification). The GH26 Man26A and GH27 Aga27A enzymes appeared to be more promising for applications that involve the hydrolysis of galactomannan containing biomass. This method of screening for maximal compatibility between various GH families can ultimately lead to a more rational development of tailored enzyme cocktails for lignocellulose hydrolysis.     

Comprehensive Analysis of Carbohydrate-Active Enzymes from the Filamentous Fungus <Emphasis Type="Italic">Scytalidium candidum</Emphasis> 3C

Complete enzymatic degradation of plant polysaccharides is a result of combined action of various carbohydrate-active enzymes (CAZymes). In this paper, we demonstrate the potential of the filamentous fungus Scytalidium candidum 3C for processing of plant biomass. Structural annotation of the improved assembly of S. candidum 3C genome and functional annotation of CAZymes revealed putative gene sequences encoding such proteins. A total of 190 CAZyme-encoding genes were identified, including 104 glycoside hydrolases, 52 glycosyltransferases, 28 oxidative enzymes, and 6 carbohydrate esterases. In addition, 14 carbohydrate-binding modules were found. Glycoside hydrolases secreted during the growth of S. candidum 3C in three media were analyzed with a variety of substrates. Mass spectrometry analysis of the fungal culture liquid revealed the presence of peptides identical to 36 glycoside hydrolases, three proteins without known enzymatic function belonging to the same group of families, and 11 oxidative enzymes. The activity of endohemicellulases was determined using specially synthesized substrates in which the glycosidic bond between monosaccharide residues was replaced by a thiolinkage. During analysis of the CAZyme profile of S. candidum 3C, four β-xylanases from the GH10 family and two β-glucanases from the GH7 and GH55 families were detected, partially purified, and identified.     

Discovery and structural characterization of a novel glycosidase family of marine origin

The genomic data on heterotrophic marine bacteria suggest the crucial role that microbes play in the global carbon cycle. However, the massive presence of hypothetical proteins hampers our understanding of the mechanisms by which this carbon cycle is carried out. Moreover, genomic data from marine microorganisms are essentially annotated in the light of the biochemical knowledge accumulated on bacteria and fungi which decompose terrestrial plants. However marine algal polysaccharides clearly differ from their terrestrial counterparts, and their associated enzymes usually constitute novel protein families. In this study, we have applied a combination of bioinformatics, targeted activity screening and structural biology to characterize a hypothetical protein from the marine bacterium Zobellia galactanivorans, which is distantly related to GH43 family. This protein is in fact a 1,3-α-3,6-anhydro-l-galactosidase (AhgA) which catalyses the last step in the degradation pathway of agars, a family of polysaccharides unique to red macroalgae. AhgA adopts a β-propeller fold and displays a zinc-dependent catalytic machinery. This enzyme is the first representative of a new family of glycoside hydrolases, especially abundant in coastal waters. Such genes of marine origin have been transferred to symbiotic microbes associated with marine fishes, but also with some specific human populations.     

Composition and Expression of Genes Encoding Carbohydrate-Active Enzymes in the Straw-Degrading Mushroom Volvariella volvacea

Volvariella volvacea is one of a few commercial cultivated mushrooms mainly using straw as carbon source. In this study, the genome of V. volcacea was sequenced and assembled. A total of 285 genes encoding carbohydrate-active enzymes (CAZymes) in V. volvacea were identified and annotated. Among 15 fungi with sequenced genomes, V. volvacea ranks seventh in the number of genes encoding CAZymes. In addition, the composition of glycoside hydrolases in V. volcacea is dramatically different from other basidiomycetes: it is particularly rich in members of the glycoside hydrolase families GH10 (hemicellulose degradation) and GH43 (hemicellulose and pectin degradation), and the lyase families PL1, PL3 and PL4 (pectin degradation) but lacks families GH5b, GH11, GH26, GH62, GH93, GH115, GH105, GH9, GH53, GH32, GH74 and CE12. Analysis of genome-wide gene expression profiles of 3 strains using 3′-tag digital gene expression (DGE) reveals that 239 CAZyme genes were expressed even in potato destrose broth medium. Our data also showed that the formation of a heterokaryotic strain could dramatically increase the expression of a number of genes which were poorly expressed in its parental homokaryotic strains.     

Hierarchical classification of glycoside hydrolases

This review deals with structural and functional features of glycoside hydrolases, a widespread group of enzymes present in almost all living organisms. Their catalytic domains are grouped into 120 amino acid sequence-based families in the international classification of the carbohydrate-active enzymes (CAZy database). At a higher hierarchical level some of these families are combined in 14 clans. Enzymes of the same clan have common evolutionary origin of their genes and share the most important functional characteristics such as composition of the active center, anomeric configuration of cleaved glycosidic bonds, and molecular mechanism of the catalyzed reaction (either inverting, or retaining). There are now extensive data in the literature concerning the relationship between glycoside hydrolase families belonging to different clans and/or included in none of them, as well as information on phylogenetic protein relationship within particular families. Summarizing these data allows us to propose a multilevel hierarchical classification of glycoside hydrolases and their homologs. It is shown that almost the whole variety of the enzyme catalytic domains can be brought into six main folds, large groups of proteins having the same three-dimensional structure and the supposed common evolutionary origin.     

The crystal structure of a family GH25 lysozyme from Bacillus anthracis implies a neighboring-group catalytic mechanism with retention of anomeric configuration

Lysozymes are found in many of the sequence-based families of glycoside hydrolases (www.cazy.org) where they show considerable structural and mechanistic diversity. Lysozymes from glycoside hydrolase family GH25 adopt a (α/β)₅(β)₃-barrel-like fold with a proposal in the literature that these enzymes act with inversion of anomeric configuration; the lack of a suitable substrate, however, means that no group has successfully demonstrated the configuration of the product. Here we report the 3-D structure of the GH25 enzyme from Bacillus anthracis at 1.4 Å resolution. We show that the active center is extremely similar to those from glycoside hydrolase families GH18, GH20, GH56, GH84, and GH85 implying that, in the absence of evidence to the contrary, GH25 enzymes also act with net retention of anomeric configuration using the neighboring-group catalytic mechanism that is common to this ‘super-family’ of enzymes.     

The amoebapore superfamily

Amoebapores, synthesized by human protozoan parasites, form ion channels in target cells and artificial lipid membranes. The major pathogenic effect of these proteins is due to their cytolytic capability which results in target cell death. They comprise a coherent family and are homologous to other proteins and protein domains found in eight families. These families include in addition to the amoebapores (1) the saposins, (2) the NK-lysins and granulysins, (3) the pulmonary surfactant proteins B, (4) the acid sphingomyelinases, (5) acyloxyacyl hydrolases and (6) the aspartic proteases. These amoebapore homologues have many properties in common including membrane binding and stability. We note for the first time that a new protein, countin, from the cellular slime mold, Dictyostelium discoideum, comprises the eighth family within this superfamily. All currently sequenced members of these eight families are identified, and the structural, functional and phylogenetic properties of these proteins are discussed.     

Several Genes Encoding Enzymes with the Same Activity Are Necessary for Aerobic Fungal Degradation of Cellulose in Nature

The cellulose-degrading fungal enzymes are glycoside hydrolases of the GH families and lytic polysaccharide monooxygenases. The entanglement of glycoside hydrolase families and functions makes it difficult to predict the enzymatic activity of glycoside hydrolases based on their sequence. In the present study we further developed the method Peptide Pattern Recognition to an automatic approach not only to find all genes encoding glycoside hydrolases and lytic polysaccharide monooxygenases in fungal genomes but also to predict the function of the genes. The functional annotation is an important feature as it provides a direct route to predict function from primary sequence. Furthermore, we used Peptide Pattern Recognition to compare the cellulose-degrading enzyme activities encoded by 39 fungal genomes. The results indicated that cellobiohydrolases and AA9 lytic polysaccharide monooxygenases are hallmarks of cellulose-degrading fungi except brown rot fungi. Furthermore, a high number of AA9, endocellulase and β-glucosidase genes were identified, not in what are known to be the strongest, specialized lignocellulose degraders but in saprophytic fungi that can use a wide variety of substrates whereas only few of these genes were found in fungi that have a limited number of natural, lignocellulotic substrates. This correlation suggests that enzymes with different properties are necessary for degradation of cellulose in different complex substrates. Interestingly, clustering of the fungi based on their predicted enzymes indicated that Ascomycota and Basidiomycota use the same enzymatic activities to degrade plant cell walls.