Structure of the Fusarium oxysporum endoglucanase I with a nonhydrolyzable substrate analogue: substrate distortion gives rise to the preferred axial orientation for the leaving group

Modulation of the catalytic activity and thermostability of a thermostable GH7 endoglucanase by engineering the key loop B3

The loop B3 of glycoside hydrolase family 7 (GH7) endoglucanases is confined into long and short types. TtCel7 is a thermophilic GH7 endoglucanase from Thermothelomyces thermophilus ATCC 42464 with a long loop B3. TtCel7 was distinct for the excellent thermostability (>30 % residual activity after 1-h incubation at 90 °C). The catalytic efficiency was reduced by removing the disulfide bond in loop B3 (C220A) and truncated the loop B3 (B3cut). However, B3cut exhibited improved thermostability, the remaining enzyme activity increased by 39 %-171 % compared toTtCel7 when treated at 70-90 °C for 1-h. Based on the analysis of molecular dynamics simulation, both loops B1 and A3 of B3cut swing toward the catalytic center, which contributed to the reduced cleft-space and increased structure-rigidity. Conversely, the deletion of disulfide bond resulted in a reduction of structural rigidity in C220A. Through structure-directed enzyme modulation, this study has identified two structural elements that are related to the catalysis and thermostability of TtCel7. The loop B3 of TtCel7 possibly stretches the catalytic pocket, thereby increases the openness of the catalytic tunnel and enhancing flexibility for efficient catalysis. Additionally, the disulfide bond within loop B3 serves to enhance structural stability and maintain a heightened level of activity.

Multifaceted Computational Modeling in Glycoscience

Glycoscience assembles all the scientific disciplines involved in studying various molecules and macromolecules containing carbohydrates and complex glycans. Such an ensemble involves one of the most extensive sets of molecules in quantity and occurrence since they occur in all microorganisms and higher organisms. Once the compositions and sequences of these molecules are established, the determination of their three-dimensional structural and dynamical features is a step toward understanding the molecular basis underlying their properties and functions. The range of the relevant computational methods capable of addressing such issues is anchored by the specificity of stereoelectronic effects from quantum chemistry to mesoscale modeling throughout molecular dynamics and mechanics and coarse-grained and docking calculations. The Review leads the reader through the detailed presentations of the applications of computational modeling. The illustrations cover carbohydrate-carbohydrate interactions, glycolipids, and N- and O-linked glycans, emphasizing their role in SARS-CoV-2. The presentation continues with the structure of polysaccharides in solution and solid-state and lipopolysaccharides in membranes. The full range of protein-carbohydrate interactions is presented, as exemplified by carbohydrate-active enzymes, transporters, lectins, antibodies, and glycosaminoglycan binding proteins. A final section features a list of 150 tools and databases to help address the many issues of structural glycobioinformatics.

Development of Non-Hydrolysable Oligosaccharide Activity-Based Inactivators for Endoglycanases: A Case Study on α-1,6 Mannanases

There is a vast genomic resource for enzymes active on carbohydrates. Lagging far behind, however, are functional chemical tools for the rapid characterization of carbohydrate-active enzymes. Activity-based probes (ABPs) offer one chemical solution to these issues with ABPs based upon cyclophellitol epoxide and aziridine covalent and irreversible inhibitors representing a potent and widespread approach. Such inhibitors for enzymes active on polysaccharides are potentially limited by the requirement for several glycosidic bonds, themselves substrates for the enzyme targets. Here, it is shown that non-hydrolysable trisaccharide can be synthesized and applied even to enzymes with challenging subsite requirements. It was found that incorporation of carbasugar moieties, which was accomplished by cuprate-assisted regioselective trans-diaxial epoxide opening of carba-mannal synthesised for this purpose, yields inactivators that act as powerful activity-based inhibitors for α-1,6 endo-mannanases. 3-D structures at 1.35-1.47 Å resolutions confirm the design rationale and binding to the enzymatic nucleophile. Carbasugar oligosaccharide cyclophellitols offer a powerful new approach for the design of robust endoglycosidase inhibitors, while the synthesis procedures presented here should allow adaptation towards activity-based endoglycosidase probes as well as configurational isosteres targeting other endoglycosidase families.

From lignocellulose to plastics: Knowledge transfer on the degradation approaches by fungi

In this review, we argue that there is much to be learned by transferring knowledge from research on lignocellulose degradation to that on plastic. Plastic waste accumulates in the environment to hazardous levels, because it is inherently recalcitrant to biological degradation. Plants evolved lignocellulose to be resistant to degradation, but with time, fungi became capable of utilising it for their nutrition. Examples of how fungal strategies to degrade lignocellulose could be insightful for plastic degradation include how fungi overcome the hydrophobicity of lignin (e.g. production of hydrophobins) and crystallinity of cellulose (e.g. oxidative approaches). In parallel, knowledge of the methods for understanding lignocellulose degradation could be insightful such as advanced microscopy, genomic and post-genomic approaches (e.g. gene expression analysis). The known limitations of biological lignocellulose degradation, such as the necessity for physiochemical pretreatments for biofuel production, can be predictive of potential restrictions of biological plastic degradation. Taking lessons from lignocellulose degradation for plastic degradation is also important for biosafety as engineered plastic-degrading fungi could also have increased plant biomass degrading capabilities. Even though plastics are significantly different from lignocellulose because they lack hydrolysable C-C or C-O bonds and therefore have higher recalcitrance, there are apparent similarities, e.g. both types of compounds are mixtures of hydrophobic polymers with amorphous and crystalline regions, and both require hydrolases and oxidoreductases for their degradation. Thus, many lessons could be learned from fungal lignocellulose degradation.

Molecular Basis for Substrate Recognition and Catalysis by a Marine Bacterial Laminarinase

Laminarin is an abundant algal polysaccharide that serves as carbon storage and fuel to meet the nutrition demands of heterotrophic microbes. Laminarin depolymerization catalyzed by microbial extracellular enzymes initiates remineralization, a key process in ocean biogeochemical cycles. Here, we described a glycoside hydrolase 16 (GH16) family laminarinase from a marine alga-associated Flavobacterium at the biochemical and structural levels. We found that the endolytic enzyme cleaved laminarin with a preference for β-1,3-glycoside linkages and showed transglycosylation activity across a broad range of acceptors. We also solved and compared high-resolution crystal structures of laminarinase in the apo form and in complex with β-1,3-tetrasaccharides, revealing an expanded catalytic cleft formed following substrate binding. Moreover, structure and mutagenesis studies identified multiple specific contacts between the enzyme and glucosyl residues essential for the substrate specificity for β-1,3-glucan. These results provide novel insights into the structural requirements for substrate binding and catalysis of GH16 family laminarinase, enriching our understanding of bacterial utilization of algal laminarin.IMPORTANCE Heterotrophic bacterial communities are key players in marine biogeochemical cycling due to their ability to remineralize organic carbon. Processing of complex organic matter requires heterotrophic bacteria to produce extracellular enzymes with precise specificity to depolymerize substrates to sizes sufficiently small for uptake. Thus, extracellular enzymatic hydrolysis initiates microbe-driven heterotrophic carbon cycling. In this study, based on biochemical and structural analyses, we revealed the depolymerization mechanism of β-1,3-glucan, a carbon reserve in algae, by laminarinase from an alga-associated marine Flavobacterium The findings provide new insights into the substrate recognition and catalysis of bacterial laminarinase and promote a better understanding of how extracellular enzymes are involved in organic matter cycling.

Crystal structures of peanut lectin in the presence of synthetic β-N- and β-S-galactosides disclose evidence for the recognition of different glycomimetic ligands

Carbohydrate-lectin interactions are involved in important cellular recognition processes, including viral and bacterial infections, inflammation and tumor metastasis. Hence, structural studies of lectin-synthetic glycan complexes are essential for understanding lectin-recognition processes and for the further design of promising chemotherapeutics that interfere with sugar-lectin interactions. Plant lectins are excellent models for the study of the molecular-recognition process. Among them, peanut lectin (PNA) is highly relevant in the field of glycobiology because of its specificity for β-galactosides, showing high affinity towards the Thomsen-Friedenreich antigen, a well known tumor-associated carbohydrate antigen. Given this specificity, PNA is one of the most frequently used molecular probes for the recognition of tumor cell-surface O-glycans. Thus, it has been extensively used in glycobiology for inhibition studies with a variety of β-galactoside and β-lactoside ligands. Here, crystal structures of PNA are reported in complex with six novel synthetic hydrolytically stable β-N- and β-S-galactosides. These complexes disclosed key molecular-binding interactions of the different sugars with PNA at the atomic level, revealing the roles of specific water molecules in protein-ligand recognition. Furthermore, binding-affinity studies by isothermal titration calorimetry showed dissociation-constant values in the micromolar range, as well as a positive multivalency effect in terms of affinity in the case of the divalent compounds. Taken together, this work provides a qualitative structural rationale for the upcoming synthesis of optimized glycoclusters designed for the study of lectin-mediated biological processes. The understanding of the recognition of β-N- and β-S-galactosides by PNA represents a benchmark in protein-carbohydrate interactions since they are novel synthetic ligands that do not belong to the family of O-linked glycosides.

α-d-Gal-cyclophellitol cyclosulfamidate is a Michaelis complex analog that stabilizes therapeutic lysosomal α-galactosidase A in Fabry disease

α-d-Gal-cyclophellitol cyclosulfamidate is a new class of neutral, conformationally-constrained competitive glycosidase inhibitor that stabilizes α-gal A and prevents its degradation both in vitro and in cellulo by mimicry of the Michaelis complex conformation.

A unique fungal strain collection from Vietnam characterized for high performance degraders of bioecological important biopolymers and lipids

Fungal strains are abundantly used throughout all areas of biotechnology and many of them are adapted to degrade complex biopolymers like chitin or lignocellulose. We therefore assembled a collection of 295 fungi from nine different habitats in Vietnam, known for its rich biodiversity, and investigated their cellulase, chitinase, xylanase and lipase activity. The collection consists of 70 isolates from wood, 55 from soil, 44 from rice straw, 3 found on fruits, 24 from oil environments (butchery), 12 from hot springs, 47 from insects as well as 27 from shrimp shells and 13 strains from crab shells. These strains were cultivated and selected by growth differences to enrich phenotypes, resulting in 211 visually different fungi. DNA isolation of 183 isolates and phylogenetic analysis was performed and 164 species were identified. All were subjected to enzyme activity assays, yielding high activities for every investigated enzyme set. In general, enzyme activity corresponded with the environment of which the strain was isolated from. Therefore, highest cellulase activity strains were isolated from wood substrates, rice straw and soil and similar substrate effects were observed for chitinase and lipase activity. Xylanase activity was similarly distributed as cellulase activity, but substantial activity was also found from fungi isolated from insects and shrimp shells. Seven strains displayed significant activities against three of the four tested substrates, while three degraded all four investigated carbon sources. The collection will be an important source for further studies. Therefore representative strains were made available to the scientific community and deposited in the public collection of the Leibniz-Institute DSMZ-German Collection of Microorganisms and Cell Cultures, Braunschweig.

An atypical interaction explains the high-affinity of a non-hydrolyzable S-linked 1,6-α-mannanase inhibitor

The non-hydrolyzable S-linked azasugar 1,6-α-mannobiosylthioisofagomine effects potent inhibition ofBacillus circulansfamily 76endo-1,6-α-mannanase through an atypical interaction involving the acid/base residue of the enzyme.

Biochemical and Structural Characterizations of Two Dictyostelium Cellobiohydrolases from the Amoebozoa Kingdom Reveal a High Level of Conservation between Distant Phylogenetic Trees of Life

Glycoside hydrolase family 7 (GH7) cellobiohydrolases (CBHs) are enzymes commonly employed in plant cell wall degradation across eukaryotic kingdoms of life, as they provide significant hydrolytic potential in cellulose turnover. To date, many fungal GH7 CBHs have been examined, yet many questions regarding structure-activity relationships in these important natural and commercial enzymes remain. Here, we present the crystal structures and a biochemical analysis of two GH7 CBHs from social amoeba: Dictyostelium discoideum Cel7A (DdiCel7A) and Dictyostelium purpureum Cel7A (DpuCel7A). DdiCel7A and DpuCel7A natively consist of a catalytic domain and do not exhibit a carbohydrate-binding module (CBM). The structures of DdiCel7A and DpuCel7A, resolved to 2.1 Å and 2.7 Å, respectively, are homologous to those of other GH7 CBHs with an enclosed active-site tunnel. Two primary differences between the Dictyostelium CBHs and the archetypal model GH7 CBH, Trichoderma reesei Cel7A (TreCel7A), occur near the hydrolytic active site and the product-binding sites. To compare the activities of these enzymes with the activity of TreCel7A, the family 1 TreCel7A CBM and linker were added to the C terminus of each of the Dictyostelium enzymes, creating DdiCel7ACBM and DpuCel7ACBM, which were recombinantly expressed in T. reesei DdiCel7ACBM and DpuCel7ACBM hydrolyzed Avicel, pretreated corn stover, and phosphoric acid-swollen cellulose as efficiently as TreCel7A when hydrolysis was compared at their temperature optima. The Ki of cellobiose was significantly higher for DdiCel7ACBM and DpuCel7ACBM than for TreCel7A: 205, 130, and 29 μM, respectively. Taken together, the present study highlights the remarkable degree of conservation of the activity of these key natural and industrial enzymes across quite distant phylogenetic trees of life.

IMPORTANCE

GH7 CBHs are among the most important cellulolytic enzymes both in nature and for emerging industrial applications for cellulose breakdown. Understanding the diversity of these key industrial enzymes is critical to engineering them for higher levels of activity and greater stability. The present work demonstrates that two GH7 CBHs from social amoeba are surprisingly quite similar in structure and activity to the canonical GH7 CBH from the model biomass-degrading fungus T. reesei when tested under equivalent conditions (with added CBM-linker domains) on an industrially relevant substrate.

Carbohydrate-active enzymes: sequences, shapes, contortions and cells

The enzyme-catalysed degradation of oligo and polysaccharides is of considerable interest in many fields ranging from the fundamental–understanding the intrinsic chemical beauty–through to the applied, including diverse practical applications in medicine and biotechnology. Carbohydrates are the most stereochemically-complex biopolymer, and myriad different natural polysaccharides have led to evolution of multifaceted enzyme consortia for their degradation. The glycosidic bonds that link sugar monomers are among the most chemically-stable, yet enzymatically-labile, bonds in the biosphere. That glycoside hydrolases can achieve a rate enhancement (kcat/kuncat) >1017-fold provides testament to their remarkable proficiency and the sophistication of their catalysis reaction mechanisms. The last two decades have seen significant advances in the discovery of new glycosidase sequences, sequence-based classification into families and clans, 3D structures and reaction mechanisms, providing new insights into enzymatic catalysis. New impetus to these studies has been provided by the challenges inherent in plant and microbial polysaccharide degradation, both in the context of environmentally-sustainable routes to foods and biofuels, and increasingly in human nutrition. Study of the reaction mechanism of glycoside hydrolases has also inspired the development of enzyme inhibitors, both as mechanistic probes and increasingly as therapeutic agents. We are on the cusp of a new era where we are learning how to dovetail powerful computational techniques with structural and kinetic data to provide an unprecedented view of conformational details of enzyme action.

Carbohydrate-protein interactions: molecular modeling insights

The article reviews the significant contributions to, and the present status of, applications of computational methods for the characterization and prediction of protein-carbohydrate interactions. After a presentation of the specific features of carbohydrate modeling, along with a brief description of the experimental data and general features of carbohydrate-protein interactions, the survey provides a thorough coverage of the available computational methods and tools. At the quantum-mechanical level, the use of both molecular orbitals and density-functional theory is critically assessed. These are followed by a presentation and critical evaluation of the applications of semiempirical and empirical methods: QM/MM, molecular dynamics, free-energy calculations, metadynamics, molecular robotics, and others. The usefulness of molecular docking in structural glycobiology is evaluated by considering recent docking- validation studies on a range of protein targets. The range of applications of these theoretical methods provides insights into the structural, energetic, and mechanistic facets that occur in the course of the recognition processes. Selected examples are provided to exemplify the usefulness and the present limitations of these computational methods in their ability to assist in elucidation of the structural basis underlying the diverse function and biological roles of carbohydrates in their dialogue with proteins. These test cases cover the field of both carbohydrate biosynthesis and glycosyltransferases, as well as glycoside hydrolases. The phenomenon of (macro)molecular recognition is illustrated for the interactions of carbohydrates with such proteins as lectins, monoclonal antibodies, GAG-binding proteins, porins, and viruses.

Genomics review of holocellulose deconstruction by aspergilli

SUMMARY

Biomass is constructed of dense recalcitrant polymeric materials: proteins, lignin, and holocellulose, a fraction constituting fibrous cellulose wrapped in hemicellulose-pectin. Bacteria and fungi are abundant in soil and forest floors, actively recycling biomass mainly by extracting sugars from holocellulose degradation. Here we review the genome-wide contents of seven Aspergillus species and unravel hundreds of gene models encoding holocellulose-degrading enzymes. Numerous apparent gene duplications followed functional evolution, grouping similar genes into smaller coherent functional families according to specialized structural features, domain organization, biochemical activity, and genus genome distribution. Aspergilli contain about 37 cellulase gene models, clustered in two mechanistic categories: 27 hydrolyze and 10 oxidize glycosidic bonds. Within the oxidative enzymes, we found two cellobiose dehydrogenases that produce oxygen radicals utilized by eight lytic polysaccharide monooxygenases that oxidize glycosidic linkages, breaking crystalline cellulose chains and making them accessible to hydrolytic enzymes. Among the hydrolases, six cellobiohydrolases with a tunnel-like structural fold embrace single crystalline cellulose chains and cooperate at nonreducing or reducing end termini, splitting off cellobiose. Five endoglucanases group into four structural families and interact randomly and internally with cellulose through an open cleft catalytic domain, and finally, seven extracellular β-glucosidases cleave cellobiose and related oligomers into glucose. Aspergilli contain, on average, 30 hemicellulase and 7 accessory gene models, distributed among 9 distinct functional categories: the backbone-attacking enzymes xylanase, mannosidase, arabinase, and xyloglucanase, the short-side-chain-removing enzymes xylan α-1,2-glucuronidase, arabinofuranosidase, and xylosidase, and the accessory enzymes acetyl xylan and feruloyl esterases.

Characterization of the Dielectric Constant in the Trichoderma reesei Cel7B Active Site

An attempt is made to evaluate the dielectric constant of the Trichoderma reesei Cel7B active site. Through kinetic measurements, the pKa value of the catalytic acid E201 is determined. Mutations (away from E201) with net charge changes are introduced to perturb the E201 pKa. It is shown that the mutation with a +1 charge change (including G225R, G230R, and A335R) decreases the pKa of E201, whereas the mutation with a -1 charge change (including Q149E, A222D, G225D, and G230D) increases the pKa. This effect is consistent with the electrostatic interaction between the changed charge and the E201 side chain. The fitting of the experimental data yields an apparent dielectric constant of 25-80. Molecular dynamics simulations with explicit water molecules indicate that the high solvent accessibility of the active site contributes largely to the high dielectric constant. ONIOM calculations show that high dielectric constant benefits the catalysis through decreasing the energy of the transition state relative to that of the enzyme substrate complex.

An allolactose trapped at the lacZ β-galactosidase active site with its galactosyl moiety in a (4)H3 conformation provides insights into the formation, conformation, and stabilization of the transition state

When lactose was incubated with G794A-β-galactosidase (a variant with a "closed" active site loop that binds transition state analogs well) an allolactose was trapped with its Gal moiety in a (4)H3 conformation, similar to the oxocarbenium ion-like conformation expected of the transition state. The numerous interactions formed between the (4)H3 structure and β-galactosidase indicate that this structure is representative of the transition state. This conformation is also very similar to that of d-galactono-1,5-lactone, a good transition state analog. Evidence indicates that substrates take up the (4)H3 conformation during migration from the shallow to the deep mode. Steric forces utilizing His418 and other residues are important for positioning the O1 leaving group into a quasi-axial position. An electrostatic interaction between the O5 of the distorted Gal and Tyr503 as well as C-H-π bonds with Trp568 are also significant. Computational studies of the energy of sugar ring distortion show that the β-galactosidase reaction itinerary is driven by energetic considerations in utilization of a (4)H3 transition state with a novel (4)C1-(4)H3-(4)C1 conformation itinerary. To our knowledge, this is the first X-ray crystallographic structural demonstration that the transition state of a natural substrate of a glycosidase has a (4)H3 conformation.

Glycosynthesis in a waterworld: new insight into the molecular basis of transglycosylation in retaining glycoside hydrolases

Carbohydrates are ubiquitous in Nature and play vital roles in many biological systems. Therefore the synthesis of carbohydrate-based compounds is of considerable interest for both research and commercial purposes. However, carbohydrates are challenging, due to the large number of sugar subunits and the multiple ways in which these can be linked together. Therefore, to tackle the challenge of glycosynthesis, chemists are increasingly turning their attention towards enzymes, which are exquisitely adapted to the intricacy of these biomolecules. In Nature, glycosidic linkages are mainly synthesized by Leloir glycosyltransferases, but can result from the action of non-Leloir transglycosylases or phosphorylases. Advantageously for chemists, non-Leloir transglycosylases are glycoside hydrolases, enzymes that are readily available and exhibit a wide range of substrate specificities. Nevertheless, non-Leloir transglycosylases are unusual glycoside hydrolases in as much that they efficiently catalyse the formation of glycosidic bonds, whereas most glycoside hydrolases favour the mechanistically related hydrolysis reaction. Unfortunately, because non-Leloir transglycosylases are almost indistinguishable from their hydrolytic counterparts, it is unclear how these enzymes overcome the ubiquity of water, thus avoiding the hydrolytic reaction. Without this knowledge, it is impossible to rationally design non-Leloir transglycosylases using the vast diversity of glycoside hydrolases as protein templates. In this critical review, a careful analysis of literature data describing non-Leloir transglycosylases and their relationship to glycoside hydrolase counterparts is used to clarify the state of the art knowledge and to establish a new rational basis for the engineering of glycoside hydrolases.

Quantifying bond distortions in transient enzyme species by a combination of density functional theory calculations and time-resolved infrared difference spectroscopy. Implications for the mechanism of dephosphorylation of the sarcoplasmic reticulum Ca(2+)-ATPase (SERCA1a)

The sarcoplasmic Ca(2+)-ATPase (SERCA1a) forms two phosphoenzyme intermediates during Ca(2+) pumping. The second intermediate E2P hydrolyzes rapidly, which is essential for the rapid removal of Ca(2+) from the cytosol of muscle cells. The present work studies whether a weakening of the scissile PO bond in the E2P ground state facilitates dephosphorylation. To this end, the experimentally known vibrational spectrum of the E2P phosphate group was calculated with density functional theory (DFT) using structural models at two levels of structural complexity: (i) Models of acetyl phosphate in simple environments and (ii) ~150 atom models of the catalytic site. It was found that the enzyme environment distorts the structure of the phosphate group: one of the terminal PO bonds is shorter in the catalytic site indicating weaker interactions than in water. However, the bond that bridges phosphate and Asp351 is unaffected. This indicates that the scissile PO bond is not weakened by the enzyme environment of E2P. A second finding was that the catalytic site of the E2P state in aqueous solution appears to adopt a structure as in the crystals with BeF3(-), where the ATPase is in a non-reactive conformation. The reactant state of the dephosphorylation reaction differs from the E2P ground state: Glu183 faces Asp351 and positions the attacking water molecule. This state has a 0.04Å longer, and thus weaker, bridging PO bond. The reactant state is not detected in our experiments, indicating that its energy is at least 1kcal/mol higher than that of the E2P ground state.

The complete conformational free energy landscape of β-xylose reveals a two-fold catalytic itinerary for β-xylanases

Unraveling the conformational catalytic itinerary of glycoside hydrolases (GHs) is a growing topic of interest in glycobiology, with major impact in the design of GH inhibitors. β-xylanases are responsible for the hydrolysis of glycosidic bonds in β-xylans, a group of hemicelluloses of high biotechnological interest that are found in plant cell walls. The precise conformations followed by the substrate during catalysis in β-xylanases have not been unambiguously resolved, with three different pathways being proposed from structural analyses. In this work, we compute the conformational free energy landscape (FEL) of β-xylose to predict the most likely catalytic itineraries followed by β-xylanases. The calculations are performed by means of ab initio metadynamics, using the Cremer-Pople puckering coordinates as collective variables. The computed FEL supports only two of the previously proposed itineraries, 2SO → [2,5B]ǂ → 5S1 and 1S3 → [4H3]ǂ → 4C1, which clearly appear in low energy regions of the FEL. Consistently, 2SO and 1S3 are conformations preactivated for catalysis in terms of free energy/anomeric charge and bond distances. The results however exclude the OE → [OS2]ǂ → B2,5 itinerary that has been recently proposed for a family 11 xylanase. Classical and ab initio QM/MM molecular dynamics simulations reveal that, in this case, the observed OE conformation has been enforced by enzyme mutation. These results add a word of caution on using modified enzymes to inform on catalytic conformational itineraries of glycoside hydrolases.

AM1/d-CB1: A Semiempirical Model for QM/MM Simulations of Chemical Glycobiology Systems

A semiempirical method based on the AM1/d Hamiltonian is introduced to model chemical glycobiological systems. We included in the parameter training set glycans and the chemical environment often found about them in glycoenzymes. Starting with RM1 and AM1/d-PhoT models we optimized H, C, N, O, and P atomic parameters targeting the best performing molecular properties that contribute to enzyme catalyzed glycan reaction mechanisms. The training set comprising glycans, amino acids, phosphates and small organic model systems was used to derive parameters that reproduce experimental data or high-level density functional results for carbohydrate, phosphate and amino acid heats of formation, amino acid proton affinities, amino acid and monosaccharide dipole moments, amino acid ionization potentials, water-phosphate interaction energies, and carbohydrate ring pucker relaxation times. The result is the AM1/d-Chemical Biology 1 or AM1/d-CB1 model that is considerably more accurate than existing NDDO methods modeling carbohydrates and the amino acids often present in the catalytic domains of glycoenzymes as well as the binding sites of lectins. Moreover, AM1/d-CB1 computed proton affinities, dipole moments, ionization potentials and heats of formation for transition state puckered carbohydrate ring conformations, observed along glycoenzyme catalyzed reaction paths, are close to values computed using DFT M06-2X. AM1/d-CB1 provides a platform from which to accurately model reactions important in chemical glycobiology.

Carbohydrate-protein interactions that drive processive polysaccharide translocation in enzymes revealed from a computational study of cellobiohydrolase processivity

Translocation of carbohydrate polymers through protein tunnels and clefts is a ubiquitous biochemical phenomenon in proteins such as polysaccharide synthases, glycoside hydrolases, and carbohydrate-binding modules. Although static snapshots of carbohydrate polymer binding in proteins have long been studied via crystallography and spectroscopy, the molecular details of polysaccharide chain processivity have not been elucidated. Here, we employ simulation to examine how a cellulose chain translocates by a disaccharide unit during the processive cycle of a glycoside hydrolase family 7 cellobiohydrolase. Our results demonstrate that these biologically and industrially important enzymes employ a two-step mechanism for chain threading to form a Michaelis complex and that the free energy barrier to chain threading is significantly lower than the hydrolysis barrier. Taken with previous studies, our findings suggest that the rate-limiting step in enzymatic cellulose degradation is the glycosylation reaction, not chain processivity. Based on the simulations, we find that strong electrostatic interactions with polar residues that are conserved in GH7 cellobiohydrolases, but not in GH7 endoglucanases, at the leading glucosyl ring provide the thermodynamic driving force for polysaccharide chain translocation. Also, we consider the role of aromatic-carbohydrate interactions, which are widespread in carbohydrate-active enzymes and have long been associated with processivity. Our analysis suggests that the primary role for these aromatic residues is to provide tunnel shape and guide the carbohydrate chain to the active site. More broadly, this work elucidates the role of common protein motifs found in carbohydrate-active enzymes that synthesize or depolymerize polysaccharides by chain translocation mechanisms coupled to catalysis.

Structural snapshots illustrate the catalytic cycle of β-galactocerebrosidase, the defective enzyme in Krabbe disease

Glycosphingolipids are ubiquitous components of mammalian cell membranes, and defects in their catabolism by lysosomal enzymes cause a diverse array of diseases. Deficiencies in the enzyme β-galactocerebrosidase (GALC) cause Krabbe disease, a devastating genetic disorder characterized by widespread demyelination and rapid, fatal neurodegeneration. Here, we present a series of high-resolution crystal structures that illustrate key steps in the catalytic cycle of GALC. We have captured a snapshot of the short-lived enzyme-substrate complex illustrating how wild-type GALC binds a bona fide substrate. We have extensively characterized the enzyme kinetics of GALC with this substrate and shown that the enzyme is active in crystallo by determining the structure of the enzyme-product complex following extended soaking of the crystals with this same substrate. We have also determined the structure of a covalent intermediate that, together with the enzyme-substrate and enzyme-product complexes, reveals conformational changes accompanying the catalytic steps and provides key mechanistic insights, laying the foundation for future design of pharmacological chaperones.

The β-glucanase ZgLamA from Zobellia galactanivorans evolved a bent active site adapted for efficient degradation of algal laminarin

Laminarinase is commonly used to describe β-1,3-glucanases widespread throughout Archaea, bacteria, and several eukaryotic lineages. Some β-1,3-glucanases have already been structurally and biochemically characterized, but very few from organisms that are in contact with genuine laminarin, the storage polysaccharide of brown algae. Here we report the heterologous expression and subsequent biochemical and structural characterization of ZgLamAGH16 from Zobellia galactanivorans, the first GH16 laminarinase from a marine bacterium associated with seaweeds. ZgLamAGH16 contains a unique additional loop, compared with other GH16 laminarinases, which is composed of 17 amino acids and gives a bent shape to the active site cleft of the enzyme. This particular topology is perfectly adapted to the U-shaped conformation of laminarin chains in solution and thus explains the predominant specificity of ZgLamAGH16 for this substrate. The three-dimensional structure of the enzyme and two enzyme-substrate complexes, one with laminaritetraose and the other with a trisaccharide of 1,3-1,4-β-d-glucan, have been determined at 1.5, 1.35, and 1.13 Å resolution, respectively. The structural comparison of substrate recognition pattern between these complexes allows the proposition that ZgLamAGH16 likely diverged from an ancestral broad specificity GH16 β-glucanase and evolved toward a bent active site topology adapted to efficient degradation of algal laminarin.

Synthesis and evaluation of 3-deoxy and 3-deoxy-3-fluoro derivatives of gluco- and manno-configured tetrahydropyridoimidazole glycosidase inhibitors

Three tetrahydropyridoimidazole-type glycosidase inhibitors have been synthesized with the 3-deoxy ribo- and arabino-, and 3-deoxy-3-fluoro gluco-configurations and two of them screened for activity against α- and β-gluco- and mannosidase enzymes. Only one substance, the 3-deoxy-3-fluoro-derivative of the gluco-configured tetrahydropyridoimidazole was found to have any activity against a single enzyme, sweet almond β-glucosidase, and even then at a level 100-fold lower than that of the corresponding simple gluco-configured tetrahydropyridoimidazole thereby underlining the importance of the 3-hydroxy group in the key substrate-enzyme interactions.

Active site plasticity within the glycoside hydrolase NagZ underlies a dynamic mechanism of substrate distortion

NagZ is a glycoside hydrolase that participates in peptidoglycan (PG) recycling by removing β-N-acetylglucosamine from PG fragments that are excised from the bacterial cell wall during growth. Notably, the products formed by NagZ, 1,6-anhydroMurNAc-peptides, activate β-lactam resistance in many Gram-negative bacteria, making this enzyme of interest as a potential therapeutic target. Crystal structure determinations of NagZ from Salmonella typhimurium and Bacillus subtilis in complex with natural substrate, trapped as a glycosyl-enzyme intermediate, and bound to product, define the reaction coordinate of the NagZ family of enzymes. The structures, combined with kinetic studies, reveal an uncommon degree of structural plasticity within the active site of a glycoside hydrolase, and unveil how NagZ drives substrate distortion using a highly mobile loop that contains a conserved histidine that has been proposed as the general acid/base.

Computational investigation of the pH dependence of loop flexibility and catalytic function in glycoside hydrolases

Cellulase enzymes cleave glycosidic bonds in cellulose to produce cellobiose via either retaining or inverting hydrolysis mechanisms, which are significantly pH-dependent. Many fungal cellulases function optimally at pH ~5, and their activities decrease dramatically at higher or lower pH. To understand the molecular-level implications of pH in cellulase structure, we use a hybrid, solvent-based, constant pH molecular dynamics method combined with pH-based replica exchange to determine the pK(a) values of titratable residues of a glycoside hydrolase (GH) family 6 cellobiohydrolase (Cel6A) and a GH family 7 cellobiohydrolase (Cel7A) from the fungus Hypocrea jecorina. For both enzymes, we demonstrate that a bound substrate significantly affects the pKa values of the acid residues at the catalytic center. The calculated pK(a) values of catalytic residues confirm their proposed roles from structural studies and are consistent with the experimentally measured apparent pKa values. Additionally, GHs are known to impart a strained pucker conformation in carbohydrate substrates in active sites for catalysis, and results from free energy calculations combined with constant pH molecular dynamics suggest that the correct ring pucker is stable near the optimal pH for both Cel6A and Cel7A. Much longer molecular dynamics simulations of Cel6A and Cel7A with fixed protonation states based on the calculated pK(a) values suggest that pH affects the flexibility of tunnel loops, which likely affects processivity and substrate complexation. Taken together, this work demonstrates several molecular-level effects of pH on GH enzymes important for cellulose turnover in the biosphere and relevant to biomass conversion processes.

Structural, biochemical, and computational characterization of the glycoside hydrolase family 7 cellobiohydrolase of the tree-killing fungus Heterobasidion irregulare

Root rot fungi of the Heterobasidion annosum complex are the most damaging pathogens in temperate forests, and the recently sequenced Heterobasidion irregulare genome revealed over 280 carbohydrate-active enzymes. Here, H. irregulare was grown on biomass, and the most abundant protein in the culture filtrate was identified as the only family 7 glycoside hydrolase in the genome, which consists of a single catalytic domain, lacking a linker and carbohydrate-binding module. The enzyme, HirCel7A, was characterized biochemically to determine the optimal conditions for activity. HirCel7A was crystallized and the structure, refined at 1.7 Å resolution, confirms that HirCel7A is a cellobiohydrolase rather than an endoglucanase, with a cellulose-binding tunnel that is more closed than Phanerochaete chrysosporium Cel7D and more open than Hypocrea jecorina Cel7A, suggesting intermediate enzyme properties. Molecular simulations were conducted to ascertain differences in enzyme-ligand interactions, ligand solvation, and loop flexibility between the family 7 glycoside hydrolase cellobiohydrolases from H. irregulare, H. jecorina, and P. chrysosporium. The structural comparisons and simulations suggest significant differences in enzyme-ligand interactions at the tunnel entrance in the -7 to -4 binding sites and suggest that a tyrosine residue at the tunnel entrance of HirCel7A may serve as an additional ligand-binding site. Additionally, the loops over the active site in H. jecorina Cel7A are more closed than loops in the other two enzymes, which has implications for the degree of processivity, endo-initiation, and substrate dissociation. Overall, this study highlights molecular level features important to understanding this biologically and industrially important family of glycoside hydrolases.

The crystallization and structural analysis of cellulases (and other glycoside hydrolases): strategies and tactics

The three-dimensional (3-D) structures of cellulases, and other glycoside hydrolases, are a central feature of research in carbohydrate chemistry and biochemistry. 3-D structure is used to inform protein engineering campaigns, both academic and industrial, which are typically used to improve the stability or activity of an enzyme. Examples of classical protein engineering goals include higher thermal stability, reduced metal-ion dependency, detergent and protease resistance, decreased product inhibition, and altered specificity. 3-D structure may also be used to interpret the behavior of enzyme variants that are derived from screening or random mutagenesis approaches, with a view to establishing an iterative design process. In other areas, 3-D structure is used as one of the many tools to probe enzymatic catalysis, typically dovetailing with physical organic chemistry approaches to provide complete reaction mechanisms for enzymes by visualizing catalytic site interactions at different stages of the reaction. Such mechanistic insight is not only fundamentally important, impacting on inhibitor and drug design approaches with ramifications way beyond cellulose hydrolysis, but also provides the framework for the design of enzyme variants to use as biocatalysts for the synthesis of bespoke oligosaccharides. Here we review some of the strategies and tactics that may be applied to the X-ray structure solution of cellulases (and other carbohydrate-active enzymes). The general approach is first to decide why you are doing the work, then to establish correct domain boundaries for truncated constructs (typically the catalytic domain only), and finally to pursue crystallization of pure, homogeneous, and monodisperse protein with appropriate ligand and additive combinations. Cellulase-specific strategies are important for the delineation of domain boundaries, while glycoside hydrolases generally also present challenges and opportunities for the selection and optimization of ligands to both aid crystallization, and also provide structural and mechanistic insight. As the many roles for plant cell wall degrading enzymes increase, so does the need for rapid high-quality structure determination to provide a sound structural foundation for understanding mechanism and specificity, and for future protein engineering strategies.

Analysis of keystone enzyme in Agar hydrolysis provides insight into the degradation (of a polysaccharide from) red seaweeds

Agars are abundant polysaccharides from marine red algae, and their chemical structure consists of alternating D-galactose and 3,6-anhydro-L-galactose residues, the latter of which are presumed to make the polymer recalcitrant to degradation by most terrestrial bacteria. Here we study a family 117 glycoside hydrolase (BpGH117) encoded within a recently discovered locus from the human gut bacterium Bacteroides plebeius. Consistent with this locus being involved in agarocolloid degradation, we show that BpGH117 is an exo-acting 3,6-anhydro-α-(1,3)-L-galactosidase that removes the 3,6-anhydrogalactose from the non-reducing end of neoagaro-oligosaccharides. A Michaelis complex of BpGH117 with neoagarobiose reveals the distortion of the constrained 3,6-anhydro-L-galactose into a conformation that favors catalysis. Furthermore, this complex, supported by analysis of site-directed mutants, provides evidence for an organization of the active site and positioning of the catalytic residues that are consistent with an inverting mechanism of catalysis and suggests that a histidine residue acts as the general acid. This latter feature differs from the vast majority of glycoside hydrolases, which use a carboxylic acid, highlighting the alternative strategies that enzymes may utilize in catalyzing the cleavage of glycosidic bonds.