A subfamily roadmap of the evolutionarily diverse glycoside hydrolase family 16 (GH16).

Glycoside hydrolase family (GH) 16 comprises a large and taxonomically diverse family of glycosidases and transglycosidases that adopt a common ß-jelly-roll fold and are active on a range of terrestrial and marine polysaccharides. Presently, broadly insightful sequence-function correlations in GH16 are hindered by a lack of a systematic subfamily structure. To fill this gap, we have used a highly scalable protein sequence similarity network analysis to delineate nearly 23,000 GH16 sequences into 23 robust subfamilies, which are strongly supported by hidden Markov model and maximum likelihood molecular phylogenetic analyses. Subsequent evaluation of over 40 experimental three-dimensional structures has highlighted key tertiary structural differences, predominantly manifested in active-site loops, that dictate substrate specificity across the GH16 evolutionary landscape. As for other large GH families (i.e. GH5, GH13, and GH43), this new subfamily classification provides a roadmap for functional glycogenomics that will guide future bioinformatics and experimental structure-function analyses. The GH16 subfamily classification is publicly available in the CAZy database. The sequence similarity network workflow used here, SSNpipe, is freely available from GitHub.

Substrate specificity, regiospecificity, and processivity in glycoside hydrolase family 74.

Glycoside hydrolase family 74 (GH74) is a historically important family of endo-ß-glucanases. On the basis of early reports of detectable activity on cellulose and soluble cellulose derivatives, GH74 was originally considered to be a "cellulase" family, although more recent studies have generally indicated a high specificity toward the ubiquitous plant cell wall matrix glycan xyloglucan. Previous studies have indicated that GH74 xyloglucanases differ in backbone cleavage regiospecificities and can adopt three distinct hydrolytic modes of action: exo, endo-dissociative, and endo-processive. To improve functional predictions within GH74, here we coupled in-depth biochemical characterization of 17 recombinant proteins with structural biology-based investigations in the context of a comprehensive molecular phylogeny, including all previously characterized family members. Elucidation of four new GH74 tertiary structures, as well as one distantly related dual seven-bladed ß-propeller protein from a marine bacterium, highlighted key structure-function relationships along protein evolutionary trajectories. We could define five phylogenetic groups, which delineated the mode of action and the regiospecificity of GH74 members. At the extremes, a major group of enzymes diverged to hydrolyze the backbone of xyloglucan nonspecifically with a dissociative mode of action and relaxed backbone regiospecificity. In contrast, a sister group of GH74 enzymes has evolved a large hydrophobic platform comprising 10 subsites, which facilitates processivity. Overall, the findings of our study refine our understanding of catalysis in GH74, providing a framework for future experimentation as well as for bioinformatics predictions of sequences emerging from (meta)genomic studies.

Discovery of α-l-arabinopyranosidases from human gut microbiome expands the diversity within glycoside hydrolase family 42.

Enzymes of the glycoside hydrolase family 42 (GH42) are widespread in bacteria of the human gut microbiome and play fundamental roles in the decomposition of both milk and plant oligosaccharides. All GH42 enzymes characterized so far have ß-galactosidase activity. Here, we report the existence of a GH42 subfamily that is exclusively specific for α-l-arabinopyranoside and describe the first representative of this subfamily. We found that this enzyme (BlArap42B) from a probiotic Bifidobacterium species cannot hydrolyze ß-galactosides. However, BlArap42B effectively hydrolyzed paeonolide and ginsenoside Rb2, plant glycosides containing an aromatic aglycone conjugated to α-l-arabinopyranosyl-(1,6)-ß-d-glucopyranoside. Paeonolide, a natural glycoside from the roots of the plant genus Paeonia, is not hydrolyzed by classical GH42 ß-galactosidases. X-ray crystallography revealed a unique Trp345-X12-Trp358 sequence motif at the BlArap42B active site, as compared with a Phe-X12-His motif in classical GH42 ß-galactosidases. This analysis also indicated that the C6 position of galactose is blocked by the aromatic side chains, hence allowing accommodation only of Arap lacking this carbon. Automated docking of paeonolide revealed that it can fit into the BlArap42B active site. The Glcp moiety of paeonolide stacks onto the aromatic ring of the Trp252 at subsite +1 and C4-OH is hydrogen bonded with Asp249 Moreover, the aglycone stacks against Phe421 from the neighboring monomer in the BlArap42B trimer, forming a proposed subsite +2. These results further support the notion that evolution of metabolic specialization can be tracked at the structural level in key enzymes facilitating degradation of specific glycans in an ecological niche.

A ß1-6/ß1-3 galactosidase from Bifidobacterium animalis subsp. lactis†Bl-04 gives insight into sub-specificities of ß-galactoside catabolism within Bifidobacterium.

The Bifidobacterium genus harbours several health promoting members of the gut microbiota. Bifidobacteria display metabolic specialization by preferentially utilizing dietary or host-derived ß-galactosides. This study investigates the biochemistry and structure of a glycoside hydrolase family 42 (GH42) ß-galactosidase from the probiotic Bifidobacterium animalis subsp. lactis Bl-04 (BlGal42A). BlGal42A displays a preference for undecorated ß1-6 and ß1-3 linked galactosides and populates a phylogenetic cluster with close bifidobacterial homologues implicated in the utilization of N-acetyl substituted ß1-3 galactosides from human milk and mucin. A long loop containing an invariant tryptophan in GH42, proposed to bind substrate at subsite + 1, is identified here as specificity signature within this clade of bifidobacterial enzymes. Galactose binding at the subsite - 1 of the active site induced conformational changes resulting in an extra polar interaction and the ordering of a flexible loop that narrows the active site. The amino acid sequence of this loop provides an additional specificity signature within this GH42 clade. The phylogenetic relatedness of enzymes targeting ß1-6 and ß1-3 galactosides likely reflects structural differences between these substrates and ß1-4 galactosides, containing an axial galactosidic bond. These data advance our molecular understanding of the evolution of sub-specificities that support metabolic specialization in the gut niche.

Functional characterization of fungal lytic polysaccharide monooxygenases for cellulose surface oxidation.

BACKGROUND: Microbial lytic polysaccharide monooxygenases (LPMOs) cleave diverse biomass polysaccharides, including cellulose and hemicelluloses, by initial oxidation at C1 or C4 of glycan chains. Within the Carbohydrate-Active Enzymes (CAZy) classification, Auxiliary Activity Family 9 (AA9) comprises the first and largest group of fungal LPMOs, which are often also found in tandem with non-catalytic carbohydrate-binding modules (CBMs). LPMOs originally attracted attention for their ability to potentiate complete biomass deconstruction to monosaccharides. More recently, LPMOs have been applied for selective surface modification of insoluble cellulose and chitin. RESULTS: To further explore the catalytic diversity of AA9 LPMOs, over 17,000 sequences were extracted from public databases, filtered, and used to construct a sequence similarity network (SSN) comprising 33 phylogenetically supported clusters. From these, 32 targets were produced successfully in the industrial filamentous fungus Aspergillus niger, 25 of which produced detectable LPMO activity. Detailed biochemical characterization of the eight most highly produced targets revealed individual C1, C4, and mixed C1/C4 regiospecificities of cellulose surface oxidation, different redox co-substrate preferences, and CBM targeting effects. Specifically, the presence of a CBM correlated with increased formation of soluble oxidized products and a more localized pattern of surface oxidation, as indicated by carbonyl-specific fluorescent labeling. On the other hand, LPMOs without native CBMs were associated with minimal release of soluble products and comparatively dispersed oxidation pattern. CONCLUSIONS: This work provides insight into the structural and functional diversity of LPMOs, and highlights the need for further detailed characterization of individual enzymes to identify those best suited for cellulose saccharification versus surface functionalization toward biomaterials applications.

Crystal structure of ß-L-arabinobiosidase belonging to glycoside hydrolase family 121.

Enzymes acting on α-L-arabinofuranosides have been extensively studied; however, the structures and functions of ß-L-arabinofuranosidases are not fully understood. Three enzymes and an ABC transporter in a gene cluster of Bifidobacterium longum JCM 1217 constitute a degradation and import system of ß-L-arabinooligosaccharides on plant hydroxyproline-rich glycoproteins. An extracellular ß-L-arabinobiosidase (HypBA2) belonging to the glycoside hydrolase (GH) family 121 plays a key role in the degradation pathway by releasing ß-1,2-linked arabinofuranose disaccharide (ß-Ara2) for the specific sugar importer. Here, we present the crystal structure of the catalytic region of HypBA2 as the first three-dimensional structure of GH121 at 1.85 Å resolution. The HypBA2 structure consists of a central catalytic (α/α)6 barrel domain and two flanking (N- and C-terminal) ß-sandwich domains. A pocket in the catalytic domain appears to be suitable for accommodating the ß-Ara2 disaccharide. Three acidic residues Glu383, Asp515, and Glu713, located in this pocket, are completely conserved among all members of GH121; site-directed mutagenesis analysis showed that they are essential for catalytic activity. The active site of HypBA2 was compared with those of structural homologs in other GH families: GH63 α-glycosidase, GH94 chitobiose phosphorylase, GH142 ß-L-arabinofuranosidase, GH78 α-L-rhamnosidase, and GH37 α,α-trehalase. Based on these analyses, we concluded that the three conserved residues are essential for catalysis and substrate binding. ß-L-Arabinobiosidase genes in GH121 are mainly found in the genomes of bifidobacteria and Xanthomonas species, suggesting that the cleavage and specific import system for the ß-Ara2 disaccharide on plant hydroxyproline-rich glycoproteins are shared in animal gut symbionts and plant pathogens.

Synergy between Cell Surface Glycosidases and Glycan-Binding Proteins Dictates the Utilization of Specific Beta(1,3)-Glucans by Human Gut Bacteroides.

The human gut microbiota (HGM) has far-reaching impacts on human health and nutrition, which are fueled primarily by the metabolism of otherwise indigestible complex carbohydrates commonly known as dietary fiber. However, the molecular basis of the ability of individual taxa of the HGM to address specific dietary glycan structures remains largely unclear. In particular, the utilization of ß(1,3)-glucans, which are widespread in the human diet as yeast, seaweed, and plant cell walls, had not previously been resolved. Through a systems-based approach, here we show that the symbiont Bacteroides uniformis deploys a single, exemplar polysaccharide utilization locus (PUL) to access yeast ß(1,3)-glucan, brown seaweed ß(1,3)-glucan (laminarin), and cereal mixed-linkage ß(1,3)/ß(1,4)-glucan. Combined biochemical, enzymatic, and structural analysis of PUL-encoded glycoside hydrolases (GHs) and surface glycan-binding proteins (SGBPs) illuminates a concerted molecular system by which B. uniformis recognizes and saccharifies these distinct ß-glucans. Strikingly, the functional characterization of homologous ß(1,3)-glucan utilization loci (1,3GUL) in other Bacteroides further demonstrated that the ability of individual taxa to utilize ß(1,3)-glucan variants and/or ß(1,3)/ß(1,4)-glucans arises combinatorially from the individual specificities of SGBPs and GHs at the cell surface, which feed corresponding signals to periplasmic hybrid two-component sensors (HTCSs) via TonB-dependent transporters (TBDTs). These data reveal the importance of cooperativity in the adaptive evolution of GH and SGBP cohorts to address individual polysaccharide structures. We anticipate that this fine-grained knowledge of PUL function will inform metabolic network analysis and proactive manipulation of the HGM. Indeed, a survey of 2,441 public human metagenomes revealed the international, yet individual-specific, distribution of each 1,3GUL.IMPORTANCEBacteroidetes are a dominant phylum of the human gut microbiota (HGM) that target otherwise indigestible dietary fiber with an arsenal of polysaccharide utilization loci (PULs), each of which is dedicated to the utilization of a specific complex carbohydrate. Here, we provide novel insight into this paradigm through functional characterization of homologous PULs from three autochthonous Bacteroides species, which target the family of dietary ß(1,3)-glucans. Through detailed biochemical and protein structural analysis, we observed an unexpected diversity in the substrate specificity of PUL glycosidases and glycan-binding proteins with regard to ß(1,3)-glucan linkage and branching patterns. In combination, these individual enzyme and protein specificities support taxon-specific growth on individual ß(1,3)-glucans. This detailed metabolic insight, together with a comprehensive survey of individual 1,3GULs across human populations, further expands the fundamental roadmap of the HGM, with potential application to the future development of microbial intervention therapies.

Biochemical and kinetic characterisation of a novel xylooligosaccharide-upregulated GH43 ß-d-xylosidase/α-l-arabinofuranosidase (BXA43) from the probiotic Bifidobacterium animalis subsp. lactis BB-12.

The Bifidobacterium animalis subsp. lactis BB-12 gene BIF_00092, assigned to encode a ß-d-xylosidase (BXA43) of glycoside hydrolase family 43 (GH43), was cloned with a C-terminal His-tag and expressed in Escherichia coli. BXA43 was purified to homogeneity from the cell lysate and found to be a dual-specificity exo-hydrolase active on para-nitrophenyl-ß-d-xylopyranoside (pNPX), para-nitrophenyl-α-L-arabinofuranoside (pNPA), ß-(1 → 4)-xylopyranosyl oligomers (XOS) of degree of polymerisation (DP) 2-4, and birchwood xylan. A phylogenetic tree of the 92 characterised GH43 enzymes displayed five distinct groups (I - V) showing specificity differences. BXA43 belonged to group IV and had an activity ratio for pNPA:pNPX of 1:25. BXA43 was stable below 40°C and at pH 4.0-8.0 and showed maximum activity at pH 5.5 and 50°C. Km and kcat for pNPX were 15.6 ± 4.2 mM and 60.6 ± 10.8 s-1, respectively, and substrate inhibition became apparent above 18 mM pNPX. Similar kinetic parameters and catalytic efficiency values were reported for ß-d-xylosidase (XynB3) from Geobacillus stearothermophilus T‒6 also belonging to group IV. The activity of BXA43 for xylooligosaccharides increased with the size and was 2.3 and 5.6 fold higher, respectively for xylobiose and xylotetraose compared to pNPX. BXA43 showed clearly metal inhibition for Zn2+ and Ag+, which is different to its close homologues. Multiple sequence alignment and homology modelling indicated that Arg505Tyr506 present in BXA43 are probably important for binding to xylotetraose at subsite +3 and occur only in GH43 from the Bifidobacterium genus.