The Reaction Coordinate of a Bacterial GH47 α-Mannosidase: A Combined Quantum Mechanical and Structural Approach

Catalytic Reaction Mechanism of Bacterial GH92 α-1,2-Mannosidase: A QM/MM Metadynamics Study

The catalytic mechanism of aC a + 2 ${C{a}^{+2}}$ -dependent family 92 α ${{\rm \alpha }}$ -mannosidase, which is abundantly present in human gut flora and malfunctions leading to the lysosomal storage disease α-mannosidosis, has been investigated using quantum mechanics/molecular mechanics and metadynamics methods. Computational efforts show that the enzyme follows a conformational itinerary of and theC a + 2 ${C{a}^{+2}}$ ion serves a dual purpose, as it not only distorts the sugar ring but also plays a crucial role in orchestrating the arrangement of catalytic residues. This orchestration, in turn, contributes to the facilitation of O S 2 ${{{\rm \ }}^{{\rm O}}{{\rm S}}_{2}}$ conformers for the ensuing reaction. This mechanistic insight is well-aligned with the experimental predictions of the catalytic pathway, and the computed energies are of the same order of magnitude as the experimental estimations. Hence, our results extend the mechanistic understanding of glycosidases.

Trans-cyclosulfamidate mannose-configured cyclitol allows isoform-dependent inhibition of GH47 α-d-mannosidases through a bump-hole strategy

Class I inverting exo-acting α-1,2-mannosidases (CAZY family GH47) display an unusual catalytic itinerary featuring ring-flipped mannosides, 3S1 → 3H4‡ → 1C4. Conformationally locked 1C4 compounds, such as kifunensine, display nanomolar inhibition but large multigene GH47 mannosidase families render specific "isoform-dependent" inhibition impossible. Here we develop a bump-and-hole strategy in which a new mannose-configured 1,6-trans-cyclic sulfamidate inhibits α-d-mannosidases by virtue of its 1C4 conformation. This compound does not inhibit the wild-type GH47 model enzyme by virtue of a steric clash, a "bump", in the active site. An L310S (a conserved residue amongst human GH47 enzymes) mutant of the model Caulobacter GH47 awoke 574 nM inhibition of the previously dormant inhibitor, confirmed by structural analysis of a 0.97 Å structure. Considering that L310 is a conserved residue amongst human GH47 enzymes, this work provides a unique framework for future biotechnological studies on N-glycan maturation and ER associated degradation by isoform-specific GH47 α-d-mannosidase inhibition through a bump-and-hole approach.

A Fluorogenic Disaccharide Substrate for α-Mannosidases Enables High-Throughput Screening and Identification of an Inhibitor of the GH92 Virulence Factor from Streptococcus pneumoniae

Trimming of host glycans is a mechanism that is broadly employed by both commensal and pathogenic microflora to enable colonization. Host glycan trimming by the opportunistic Gram-positive bacterium Streptococcus pneumoniae has been demonstrated to be an important mechanism of virulence. While S. pneumoniae employs a multitude of glycan processing enzymes, the exo-mannosidase SpGH92 has been shown to be an important virulence factor. Accordingly, SpGH92 is hypothesized to be a target for much-needed new treatments of S. pneumoniae infection. Here we report the synthesis of 4-methylumbelliferyl α-d-mannopyranosyl-(1→2)-β-d-mannopyranoside (Manα1,2Manβ-4MU) as a fluorogenic disaccharide substrate and development of an assay for SpGH92 that overcomes its requirement for +1 binding site occupancy. We miniaturize our in vitro assay and apply it to a high-throughput screen of >65 000 compounds, identifying a single inhibitory chemotype, LIPS-343. We further show that Manα1,2Manβ-4MU is also a substrate of the human Golgi-localized α-mannosidase MAN1A1, suggesting that this substrate should be useful for assessing the activity of this and other mammalian α-mannosidases.

Reaction Mechanism of Glycoside Hydrolase Family 116 Utilizes Perpendicular Protonation

Retaining glycoside hydrolases use acid/base catalysis with an enzymatic acid/base protonating the glycosidic bond oxygen to facilitate leaving-group departure alongside attack by a catalytic nucleophile to form a covalent intermediate. Generally, this acid/base protonates the oxygen laterally with respect to the sugar ring, which places the catalytic acid/base and nucleophile carboxylates within about 4.5-6.5 Å of each other. However, in glycoside hydrolase (GH) family 116, including disease-related human acid β-glucosidase 2 (GBA2), the distance between the catalytic acid/base and the nucleophile is around 8 Å (PDB: 5BVU) and the catalytic acid/base appears to be above the plane of the pyranose ring, rather than being lateral to that plane, which could have catalytic consequences. However, no structure of an enzyme-substrate complex is available for this GH family. Here, we report the structures of Thermoanaerobacterium xylanolyticum β-glucosidase (TxGH116) D593N acid/base mutant in complexes with cellobiose and laminaribiose and its catalytic mechanism. We confirm that the amide hydrogen bonding to the glycosidic oxygen is in a perpendicular rather than lateral orientation. Quantum mechanics/molecular mechanics (QM/MM) simulations of the glycosylation half-reaction in wild-type TxGH116 indicate that the substrate binds with the nonreducing glucose residue in an unusual relaxed 4C1 chair at the -1 subsite. Nevertheless, the reaction can still proceed through a 4H3 half-chair transition state, as in classical retaining β-glucosidases, as the catalytic acid D593 protonates the perpendicular electron pair. The glucose C6OH is locked in a gauche, trans orientation with respect to the C5-O5 and C4-C5 bonds to facilitate perpendicular protonation. These data imply a unique protonation trajectory in Clan-O glycoside hydrolases, which has strong implications for the design of inhibitors specific to either lateral protonators, such as human GBA1, or perpendicular protonators, such as human GBA2.

Glycosidase mechanisms: Sugar conformations and reactivity in endo- and exo-acting enzymes

The enzymatic breakdown of carbohydrates plays a critical role in several biological events and enables the development of sustainable processes to obtain bioproducts and biofuels. In this scenario, the design of efficient inhibitors for glycosidases that can act as drug targets and the engineering of carbohydrate-active enzymes with tailored catalytic properties is of remarkable importance. To guide rational approaches, it is necessary to elucidate enzyme molecular mechanisms, in particular understanding how the microenvironment modulates the conformational space explored by the substrate. Computer simulations, especially those based on ab initio methods, have provided a suitable atomic description of carbohydrate conformations and catalytic reactions in several glycosidase families. In this review, we will focus on how the active-site topology (pocket or cleft) and mode of cleavage (endo or exo) can affect the catalytic mechanisms adopted by glycosidases, in particular the substrate conformations along the reaction coordinate.

Impact of Polarization on the Ring Puckering Dynamics of Hexose Monosaccharides

Analysis of crystal structures of hexose monosaccharides α-d-mannose (α-MAN), β-d-mannose (β-MAN), α-d-glucose (α-GLC), β-d-glucose (β-GLC), α-d-galactose (α-GAL), β-d-galactose (β-GAL), α-d-altrose (α-ALT), β-d-altrose (β-ALT), α-d-idose (α-IDO), and β-d-idose (β-IDO) reveals that the monosaccharide ring adopts multiple ring conformations. These ring conformations can be broadly classified as chair, half-chair, envelope, boat, and skew-boat conformations. The ability of the monosaccharide ring to adopt multiple conformations has been closely tied with their bioactivity. However, it has been difficult to capture the dynamic information of these conformations from experimental studies. Even from simulations, capturing these different conformations is challenging because of the energy barriers involved in the transitions between the stable ⁴C₁ and ¹C₄ chair forms. In this study, we analyze the influence of the polarizable force field on the ring dynamics of five major types of unsubstituted aldohexoses─glucose, mannose, galactose, altrose, and idose─and their anomers. We simulate microsecond trajectories to capture the influence of the CHARMM36 additive and polarizable carbohydrate force fields on the ring dynamics. The microsecond trajectories allow us to comment on the issues associated with equilibrium molecular dynamics simulations. Further, we use the extended system adaptive biasing force (eABF) method to compare the conformational sampling efficiencies of the additive and polarizable force fields. Our studies reveal that inclusion of polarization enhances the sampling of ring conformations and lowers the energy barriers between the ⁴C₁ and ¹C₄ conformations. Overall, the CHARMM36 additive force field is observed to be rigid and favor the ⁴C₁ conformations. Although the inclusion of polarizability results in enhancing ring flexibility, we observe sampling that does not agree with experimental results, warranting a revision of the polarizable Drude parameters.

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.

Molecular simulations of complex carbohydrates and glycoconjugates

Complex carbohydrates (glycans) are the most abundant and versatile biopolymers in nature. The broad diversity of biochemical functions that carbohydrates cover is a direct consequence of the variety of 3D architectures they can adopt, displaying branched or linear arrangements, widely ranging in sizes, and with the highest diversity of building blocks of any other natural biopolymer. Despite this unparalleled complexity, a common denominator can be found in the glycans' inherent flexibility, which hinders experimental characterization, but that can be addressed by high-performance computing (HPC)-based molecular simulations. In this short review, I present and discuss the state-of-the-art of molecular simulations of complex carbohydrates and glycoconjugates, highlighting methodological strengths and weaknesses, important insights through emblematic case studies, and suggesting perspectives for future developments.

QM/MM investigation of the catalytic mechanism of processive endoglucanase Cel9G from Clostridium cellulovorans

Carbohydrate degradation catalyzed by glucoside hydrolases (GHs) is a major mechanism in biomass conversion. GH family 9 endoglucanase (Cel9G) from Clostridium cellulovorans, a typical multimodular enzyme, contains a catalytic domain closely linked to a family 3c carbohydrate-binding module (CBM3c). Unlike the conventional behavior proposed for other carbohydrate-binding modules, CBM3c has a direct impact on catalytic activity. In this work, extensive molecular dynamics (MD) simulations were employed to clarify the functional role of CBM3c. Furthermore, the detailed catalytic mechanism of Cel9G was investigated at the atomistic level using the combined quantum mechanical and molecular mechanical (QM/MM) method. Based on these simulations, owing to the rigidity of the peptide linker, CBM3c may affect the enzymatic activity via direct interactions with alpha helix 4 of GH9, especially with the K123 and H125 residues. In addition, using cellohexaose as a substrate, the QM/MM MD simulations confirmed that this enzyme can cleave the β-1,4-glycosidic linkage via an inverting mechanism. An oxocarbenium ion-like transition state was located with a barrier height of 19.6 kcal mol^-1. Furthermore, the G(-1) pyranose unit preferentially adopted a distorted ¹S₅/⁴H₅ conformer in the enzyme-substrate complex. For the cleavage of the glycosidic bond, we were able to identify a plausible route (¹S₅/⁴H₅ → [⁴H₅/⁴E]^# → ⁴C₁) from the reactant to the product at the G(-1) site.

Functional glyco-metagenomics elucidates the role of glycan-related genes in environments

BACKGROUND

Glycan-related genes play a fundamental role in various processes for energy acquisition and homeostasis maintenance while adapting to the environment in which the organism exists; however, their role in the microbiome in the environment is unclear.

METHODS

Sequence alignment was performed between known glycan-related genes and complete genomes of microorganisms, and optimal parameters for identifying glycan-related genes were determined based on the alignments. Using the constructed scheme (> 90% of identity and > 25 aa of alignment length), glycan-related genes in various environments were identified from 198 different metagenome data.

RESULTS

As a result, we identified 86.73 million glycan-related genes from the metagenome data. Among the 12 environments classified in this study, the percentage of glycan-related genes was high in the human-associated environment, suggesting that these environments utilize glycan metabolism better than other environments. On the other hand, the relative abundances of both glycoside hydrolases and glycosyltransferases surprisingly had a coverage of over 80% in all the environments. These glycoside hydrolases and glycosyltransferases were classified into two groups of (1) general enzyme families identified in various environments and (2) specific enzymes found only in certain environments. The general enzyme families were mostly from genes involved in monosaccharide metabolism, and most of the specific enzymes were polysaccharide degrading enzymes.

CONCLUSION

These findings suggest that environmental microorganisms could change the composition of their glycan-related genes to adapt the processes involved in acquiring energy from glycans in their environments. Our functional glyco-metagenomics approach has made it possible to clarify the relationship between the environment and genes from the perspective of carbohydrates, and the existence of glycan-related genes that exist specifically in the environment.

GH47 and Other Glycoside Hydrolases Catalyze Glycosidic Bond Cleavage with the Assistance of Substrate Super-arming at the Transition State

Super-armed glycosyl donors, whose substituents are predominantly held in pseudoaxial positions, exhibit strongly increased reactivity in glycosylation through significant stabilization of oxocarbenium-like transition states. Examination of X-ray crystal structures reveals that the GH47 family of glycoside hydrolases have evolved so as to distort their substrates away from the ground state conformation in such a manner as to present multiple C-O bonds in pseudoaxial positions and so benefit from conformational super-arming of their substrates, thereby enhancing catalysis. Through analysis of literature mutagenic studies, we show that a suitably placed aromatic residue in GHs 6 and 47 sterically enforces super-armed conformations on their substrates. GH families 45, 81, and 134 on the other hand impose conformational super-arming on their substrates, by maintaining the more active ring conformation through hydrogen bonding rather than steric interactions. The recognition of substrate super-arming by select GH families provides a further parallel with synthetic carbohydrate chemistry and nature and opens further avenues for the design of improved glycosidase inhibitors.

Approaches to the Synthesis of Perfluoroalkyl-Modified Carbohydrates and Derivatives: Thiosugars, Iminosugars, and Tetrahydro(thio)pyrans

Fluoroalkyl-substituted carbohydrates play relevant roles in diverse areas such as supramolecular chemistry, glycoconjugation, liquid crystals, and surfactants, with direct applications as wetting, antifreeze, and coating agents. In light of these promising applications, new methodologies for the late-stage incorporation of fluoroalkyl R_F groups into carbohydrates and derivatives are herein presented as they are relevant to the synthetic carbohydrate community. Previously reviewed protocols for the installation of R_F groups onto carbohydrates and derivatives will be succinctly summarized in the light of the new achievements. Fluoroalkyl-substituted iminosugars, on the other hand, are also interesting glycomimetic derivatives with prominent roles as glycosidases and glycosyltransferases inhibitors, as has recently been demonstrated. Also, they positively contribute to the study of sugar-protein interactions and enzyme mechanisms. New advances in the syntheses of fluoroalkyl-substituted iminosugars will also be presented here.

Two distinct catalytic pathways for GH43 xylanolytic enzymes unveiled by X-ray and QM/MM simulations

Xylanolytic enzymes from glycoside hydrolase family 43 (GH43) are involved in the breakdown of hemicellulose, the second most abundant carbohydrate in plants. Here, we kinetically and mechanistically describe the non-reducing-end xylose-releasing exo-oligoxylanase activity and report the crystal structure of a native GH43 Michaelis complex with its substrate prior to hydrolysis. Two distinct calcium-stabilized conformations of the active site xylosyl unit are found, suggesting two alternative catalytic routes. These results are confirmed by QM/MM simulations that unveil the complete hydrolysis mechanism and identify two possible reaction pathways, involving different transition state conformations for the cleavage of xylooligosaccharides. Such catalytic conformational promiscuity in glycosidases is related to the open architecture of the active site and thus might be extended to other exo-acting enzymes. These findings expand the current general model of catalytic mechanism of glycosidases, a main reaction in nature, and impact on our understanding about their interaction with substrates and inhibitors.

Mechanistic Insights into the Hydrolysis of O-GlcNAcylation Catalyzed by Human O-GlcNAcase

O-GlcNAc glycosylation occurs on specific serine/threonine residues of intracellular proteins, which is widely related to various diseases, including type II diabetes, cancer, and Alzheimer's disease. Human O-GlcNAcase (hOGA) is responsible for the removal of O-GlcNAc modification and thus serves as the main target for inhibitor design. In this work, we systematically investigated the mechanism catalyzed by hOGA using the quantum mechanical/molecular mechanical method. Based on calculated free energy profiles, two essential steps named cyclization (Cyc) step and ring opening step are required to generate the final hemiacetal product. The Cyc of the 2-acetamido group, the rate-limiting step, leads to the generation of the intermediate of a bicyclic oxazolinium ion (EI1). Next, the oxazoline ring could be broken via the nucleophilic attack of a water molecule at the C1 position, which generates the final product. Along with this, our simulations clearly suggest the existence of an oxazoline intermediate (EI2), which is produced via proton transfer (PT) from the 2-acetamido group (EI1) to D174. This PT step features a reversible process with a low energy barrier, which could be attributed to a low barrier hydrogen bond between the donor and acceptor. The stabilizing effect of the low barrier hydrogen bond on EI1 is proposed to be very important for accelerating the overall reaction. In fact, the site-directed mutagenesis simulations of D174A and D175A strongly indicate that the catalytic residues mainly affect the observed reaction rate by affecting the stability of the intermediate.

The hydrolysis mechanism of a GH45 cellulase and its potential relation to lytic transglycosylase and expansin function

Family 45 glycoside hydrolases (GH45) are endoglucanases that are integral to cellulolytic secretomes, and their ability to break down cellulose has been successfully exploited in textile and detergent industries. In addition to their industrial relevance, understanding the molecular mechanism of GH45-catalyzed hydrolysis is of fundamental importance because of their structural similarity to cell wall-modifying enzymes such as bacterial lytic transglycosylases (LTs) and expansins present in bacteria, plants, and fungi. Our understanding of the catalytic itinerary of GH45s has been incomplete because a crystal structure with substrate spanning the -1 to +1 subsites is currently lacking. Here we constructed and validated a putative Michaelis complex in silico and used it to elucidate the hydrolytic mechanism in a GH45, Cel45A from the fungus Humicola insolens, via unbiased simulation approaches. These molecular simulations revealed that the solvent-exposed active-site architecture results in lack of coordination for the hydroxymethyl group of the substrate at the -1 subsite. This lack of coordination imparted mobility to the hydroxymethyl group and enabled a crucial hydrogen bond with the catalytic acid during and after the reaction. This suggests the possibility of a nonhydrolytic reaction mechanism when the catalytic base aspartic acid is missing, as is the case in some LTs (murein transglycosylase A) and expansins. We calculated reaction free energies and demonstrate the thermodynamic feasibility of the hydrolytic and nonhydrolytic reaction mechanisms. Our results provide molecular insights into the hydrolysis mechanism in HiCel45A, with possible implications for elucidating the elusive catalytic mechanism in LTs and expansins.

Normal mode analysis calculation for a full-atom model with a smaller number of degrees of freedom for huge protein molecules

The number of degrees of freedom (DOF), N, in normal mode analysis (NMA) calculations of proteins is a crucial problem in huge systems because the eigenvalue problem of an N-by-N matrix must be solved. If it were possible to perform the analysis with a smaller number of DOF for the same system with minimal deterioration in accuracy, this would make a significant impact on the computational study of protein dynamics. We examined two models in which the number of DOF was reduced. Both of them adopted a full-atom model with dihedral angles as independent variables. In one model, side-chain dihedral angles, χ's, and a main-chain dihedral angle, ω, were fixed and only the main-chain dihedral angles, ϕ and ψ, were variable. In another model, the dihedral angles around virtual bonds that connect neighboring Cα atoms were tested. The number of DOFs for the two models was two and one per residue, respectively. The residue-by-residue fluctuation profiles for atoms and dihedral angles were well reproduced in both models. The motion of atoms in the individual lowest-frequency normal modes of the two models was also very similar to those of the original model in which all rotatable dihedral angles were variable. Consequently, these models could predict large-amplitude concerted motion. These results also imply that proteins in a full-atom model can undergo only limited large-scale conformational changes around the native conformation, and consequently, NMA results do not strongly depend on the independent variables adopted.

Mannosidase mechanism: at the intersection of conformation and catalysis

Mannosidases are a diverse group of enzymes that are important in the biological processing of mannose-containing polysaccharides and complex glycoconjugates. They are found in 12 of the >160 sequence-based glycosidase families. We discuss evidence that nature has evolved a small set of common mechanisms that unite almost all of these mannosidase families. Broadly, mannosidases (and the closely related rhamnosidases) perform catalysis through just two conformations of the oxocarbenium ion-like transition state: a B_2,5 (or enantiomeric ^2,5B) boat and a ³H₄ half-chair. This extends to a new family (GT108) of GDPMan-dependent β-1,2-mannosyltransferases/phosphorylases that perform mannosyl transfer through a boat conformation as well as some mannosidases that are metalloenzymes and require divalent cations for catalysis. Yet, among this commonality lies diversity. New evidence shows that one unique family (GH99) of mannosidases use an unusual mechanism involving anchimeric assistance via a 1,2-anhydro sugar (epoxide) intermediate.

Modeling catalytic reaction mechanisms in glycoside hydrolases

Modeling catalysis in carbohydrate-active enzymes is a daunting challenge because of the high flexibility and diversity of both enzymes and carbohydrates. Glycoside hydrolases (GHs) are an illustrative example, where conformational changes and subtle interactions have been shown to be critical for catalysis. GHs have pivotal roles in industry (e.g. biofuel or detergent production) and biomedicine (e.g. targets for cancer and diabetes), and thus, a huge effort is devoted to unveil their molecular mechanisms. Besides experimental techniques, computational methods have served to provide an in-depth understanding of GH mechanisms, capturing complex reaction coordinates and the conformational itineraries that substrates follow during the whole catalytic pathway, providing a framework that ultimately may assist the engineering of these enzymes and the design of new inhibitors.

Distortion of mannoimidazole supports a B2,5 boat transition state for the family GH125 α-1,6-mannosidase from Clostridium perfringens

Enzyme transition-state mimics can act as powerful inhibitors and allow structural studies that report on the conformation of the transition-state.

α-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.

Exploration of Strategies for Mechanism-Based Inhibitor Design for Family GH99 endo-α-1,2-Mannanases

endo-α-1,2-Mannosidases and -mannanases, members of glycoside hydrolase family 99 (GH99), cleave α-Glc/Man-1,3-α-Man-OR structures within mammalian N-linked glycans and fungal α-mannan, respectively. They are proposed to act through a two-step mechanism involving a 1,2-anhydrosugar "epoxide" intermediate incorporating two conserved catalytic carboxylates. In the first step, one carboxylate acts as a general base to deprotonate the 2-hydroxy group adjacent to the fissile glycosidic bond, and the other provides general acid assistance to the departure of the aglycon. We report herein the synthesis of two inhibitors designed to interact with either the general base (α-mannosyl-1,3-(2-aminodeoxymannojirimycin), Man2NH2 DMJ) or the general acid (α-mannosyl-1,3-mannoimidazole, ManManIm). Modest affinities were observed for an endo-α-1,2-mannanase from Bacteroides thetaiotaomicron. Structural studies revealed that Man2NH2 DMJ binds like other iminosugar inhibitors, which suggests that the poor inhibition shown by this compound is not a result of a failure to achieve the expected interaction with the general base, but rather the reduction in basicity of the endocyclic nitrogen caused by introduction of a vicinal, protonated amine at C2. ManManIm binds with the imidazole headgroup distorted downwards, a result of an unfavourable interaction with a conserved active site tyrosine. This study has identified important limitations associated with mechanism-inspired inhibitor design for GH99 enzymes.

Bacteroides thetaiotaomicron generates diverse α-mannosidase activities through subtle evolution of a distal substrate-binding motif

A dominant human gut microbe, the well studied symbiont Bacteroides thetaiotaomicron (Bt), is a glyco-specialist that harbors a large repertoire of genes devoted to carbohydrate processing. Despite strong similarities among them, many of the encoded enzymes have evolved distinct substrate specificities, and through the clustering of cognate genes within operons termed polysaccharide-utilization loci (PULs) enable the fulfilment of complex biological roles. Structural analyses of two glycoside hydrolase family 92 α-mannosidases, BT3130 and BT3965, together with mechanistically relevant complexes at 1.8-2.5 Å resolution reveal conservation of the global enzyme fold and core catalytic apparatus despite different linkage specificities. Structure comparison shows that Bt differentiates the activity of these enzymes through evolution of a highly variable substrate-binding region immediately adjacent to the active site. These observations unveil a genetic/biochemical mechanism through which polysaccharide-processing bacteria can evolve new and specific biochemical activities from otherwise highly similar gene products.

1,6-Cyclophellitol Cyclosulfates: A New Class of Irreversible Glycosidase Inhibitor

The essential biological roles played by glycosidases, coupled to the diverse therapeutic benefits of pharmacologically targeting these enzymes, provide considerable motivation for the development of new inhibitor classes. Cyclophellitol epoxides and aziridines are recently established covalent glycosidase inactivators. Inspired by the application of cyclic sulfates as electrophilic equivalents of epoxides in organic synthesis, we sought to test whether cyclophellitol cyclosulfates would similarly act as irreversible glycosidase inhibitors. Here we present the synthesis, conformational analysis, and application of novel 1,6-cyclophellitol cyclosulfates. We show that 1,6-epi-cyclophellitol cyclosulfate (α-cyclosulfate) is a rapidly reacting α-glucosidase inhibitor whose ⁴C₁ chair conformation matches that adopted by α-glucosidase Michaelis complexes. The 1,6-cyclophellitol cyclosulfate (β-cyclosulfate) reacts more slowly, likely reflecting its conformational restrictions. Selective glycosidase inhibitors are invaluable as mechanistic probes and therapeutic agents, and we propose cyclophellitol cyclosulfates as a valuable new class of carbohydrate mimetics for application in these directions.

Carba-cyclophellitols Are Neutral Retaining-Glucosidase Inhibitors

The conformational analysis of glycosidases affords a route to their specific inhibition through transition-state mimicry. Inspired by the rapid reaction rates of cyclophellitol and cyclophellitol aziridine-both covalent retaining β-glucosidase inhibitors-we postulated that the corresponding carba "cyclopropyl" analogue would be a potent retaining β-glucosidase inhibitor for those enzymes reacting through the 4H3 transition-state conformation. Ab initio metadynamics simulations of the conformational free energy landscape for the cyclopropyl inhibitors show a strong bias for the 4H3 conformation, and carba-cyclophellitol, with an N-(4-azidobutyl)carboxamide moiety, proved to be a potent inhibitor (Ki = 8.2 nM) of the Thermotoga maritima TmGH1 β-glucosidase. 3-D structural analysis and comparison with unreacted epoxides show that this compound indeed binds in the 4H3 conformation, suggesting that conformational strain induced through a cyclopropyl unit may add to the armory of tight-binding inhibitor designs.

Systematic analysis of residual density suggests that a major limitation in well-refined X-ray structures of proteins is the omission of ordered solvent

In the residual electron density map of a fully refined X-ray protein model, there should be no peaks arising from modeling errors or missing atoms. Any residual peaks that do occur should be contributed by random residual intensity differences between the model and the data. If the model is incomplete (i.e., some atoms are missing), there will be more positive peaks than negative ones. On the other hand, if the model includes inappropriately located atoms, there will be an excess of negative peaks. In this study, random residual peaks are quantified using the probability density function P(x), which is defined as the probability for a peak having peak height between x and x + dx. It is found that P(x) is single-exponential and symmetric for both positive and negative peaks. Thus, P(x) can be used to discriminate residual peaks contributed by random noise in complete models from residual peaks being attributable to modeling errors in incomplete models. For a number of representative structures in the PDB it is found that P(x) has far more large (greater than 5 sigma) positive peaks than large negative peaks. This excess of large positive peaks suggests that the main defect in these refined structures is the omission of ordered water molecules.

In silico analysis of interaction pattern switching in ligand⋯receptor binding in Golgi α-mannosidase II induced by the protonated states of inhibitors

Different binding modes for charge-neutral and protonated inhibitor forms in Golgi α-mannosidase II active sites may influence their biological activities.

A β-Mannanase with a Lysozyme-like Fold and a Novel Molecular Catalytic Mechanism

The enzymatic cleavage of β-1,4-mannans is achieved by endo-β-1,4-mannanases, enzymes involved in germination of seeds and microbial hemicellulose degradation, and which have increasing industrial and consumer product applications. β-Mannanases occur in a range of families of the CAZy sequence-based glycoside hydrolase (GH) classification scheme including families 5, 26, and 113. In this work we reveal that β-mannanases of the newly described GH family 134 differ from other mannanase families in both their mechanism and tertiary structure. A representative GH family 134 endo-β-1,4-mannanase from a Streptomyces sp. displays a fold closely related to that of hen egg white lysozyme but acts with inversion of stereochemistry. A Michaelis complex with mannopentaose, and a product complex with mannotriose, reveal ligands with pyranose rings distorted in an unusual inverted chair conformation. Ab initio quantum mechanics/molecular mechanics metadynamics quantified the energetically accessible ring conformations and provided evidence in support of a ¹C₄ → ³H₄^‡ → ³S₁ conformational itinerary along the reaction coordinate. This work, in concert with that on GH family 124 cellulases, reveals how the lysozyme fold can be co-opted to catalyze the hydrolysis of different polysaccharides in a mechanistically distinct manner.

Mannoside recognition and degradation by bacteria

Mannosides constitute a vast group of glycans widely distributed in nature. Produced by almost all organisms, these carbohydrates are involved in numerous cellular processes, such as cell structuration, protein maturation and signalling, mediation of protein-protein interactions and cell recognition. The ubiquitous presence of mannosides in the environment means they are a reliable source of carbon and energy for bacteria, which have developed complex strategies to harvest them. This review focuses on the various mannosides that can be found in nature and details their structure. It underlines their involvement in cellular interactions and finally describes the latest discoveries regarding the catalytic machinery and metabolic pathways that bacteria have developed to metabolize them.

Substrate recognition and catalysis by GH47 α-mannosidases involved in Asn-linked glycan maturation in the mammalian secretory pathway

Maturation of Asn-linked oligosaccharides in the eukaryotic secretory pathway requires the trimming of nascent glycan chains to remove all glucose and several mannose residues before extension into complex-type structures on the cell surface and secreted glycoproteins. Multiple glycoside hydrolase family 47 (GH47) α-mannosidases, including endoplasmic reticulum (ER) α-mannosidase I (ERManI) and Golgi α-mannosidase IA (GMIA), are responsible for cleavage of terminal α1,2-linked mannose residues to produce uniquely trimmed oligomannose isomers that are necessary for ER glycoprotein quality control and glycan maturation. ERManI and GMIA have similar catalytic domain structures, but each enzyme cleaves distinct residues from tribranched oligomannose glycan substrates. The structural basis for branch-specific cleavage by ERManI and GMIA was explored by replacing an essential enzyme-bound Ca²⁺ ion with a lanthanum (La³⁺) ion. This ion swap led to enzyme inactivation while retaining high-affinity substrate interactions. Cocrystallization of La³⁺-bound enzymes with Man₉GlcNAc₂ substrate analogs revealed enzyme-substrate complexes with distinct modes of glycan branch insertion into the respective enzyme active-site clefts. Both enzymes had glycan interactions that extended across the entire glycan structure, but each enzyme engaged a different glycan branch and used different sets of glycan interactions. Additional mutagenesis and time-course studies of glycan cleavage probed the structural basis of enzyme specificity. The results provide insights into the enzyme catalytic mechanisms and reveal structural snapshots of the sequential glycan cleavage events. The data also indicate that full steric access to glycan substrates determines the efficiency of mannose-trimming reactions that control the conversion to complex-type structures in mammalian cells.

Enzymatic Cleavage of Glycosidic Bonds: Strategies on How to Set Up and Control a QM/MM Metadynamics Simulation

Carbohydrates play crucial roles in many biological processes, from cell-cell adhesion to chemical signaling. Their complexity and diversity, related to α/β anomeric configuration, ring substituents, and conformational variations, require a diverse set of enzymes for their processing. Among them, glycoside hydrolases (GHs) are responsible for the hydrolysis of one of the strongest bonds in nature: the glycosidic bond. These highly specialized biological catalysts select particular conformations their carbohydrate substrates to enhance catalysis. The evolution of this conformation during the reaction of glycosidic bond cleavage, known as the conformational catalytic itinerary, is of fundamental interest in glycobiology, with impact on inhibitor and drug design. Here we review some of the aspects and the main strategies one needs to take into account when simulating a reaction in a GH enzyme using QM/MM metadynamics. Several specific aspects are highlighted, from the importance of the distortion of the substrate at the Michaelis complex to the variable control during the metadynamics simulation or the analysis of the reaction mechanism and conformational itinerary. The increasing speed of computer power and methodological advances have added a vital tool to the study of GH mechanisms, as shown here and recent reviews. It is hoped that this chapter will serve as a first guide for those attempting to perform a metadynamics simulation of these relevant and fascinating enzymes.

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.

Evidence for a boat conformation at the transition state of GH76 α-1,6-mannanases--key enzymes in bacterial and fungal mannoprotein metabolism

α-Mannosidases and α-mannanases have attracted attention for the insight they provide into nucleophilic substitution at the hindered anomeric center of α-mannosides, and the potential of mannosidase inhibitors as cellular probes and therapeutic agents. We report the conformational itinerary of the family GH76 α-mannanases studied through structural analysis of the Michaelis complex and synthesis and evaluation of novel aza/imino sugar inhibitors. A Michaelis complex in an (O) S2 conformation, coupled with distortion of an azasugar in an inhibitor complex to a high energy B2,5 conformation are rationalized through ab initio QM/MM metadynamics that show how the enzyme surface restricts the conformational landscape of the substrate, rendering the B2,5 conformation the most energetically stable on-enzyme. We conclude that GH76 enzymes perform catalysis using an itinerary that passes through (O) S2 and B2,5 (≠) conformations, information that should inspire the development of new antifungal agents.

Estimation of the quality of refined protein crystal structures

Crystallographic R work and R free values, which are measures of the ability of the models of macromolecular structures to explain the crystallographic data on which they are based, are often used to assess structure quality. It is widely known, and confirmed here that both are sensitive to the methods used to compute them, and can be manipulated to improve the apparent quality of the model. As an alternative it is proposed here that the quality of crystallographic models should be assessed using a global goodness-of-fit metric RO2A /R work where RO2A is the number of reflections used for refinement divided by the number of nonhydrogen atoms in the structure, and R work is the working R-factor of the refined structure. Also, analysis of structures in the Protein Data Bank suggests that many data sets have been truncated at high resolution, thereby improving the R-factor statistics. To discourage this practice, it is proposed that the resolution of a dataset be defined as the resolution of the shell of data where <I/σI > falls to 1. The proposed goodness-of-fit metric encourages investigators to use all the data available rather than a truncated subset.

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.

Dissecting conformational contributions to glycosidase catalysis and inhibition

•

The conformational itinerary describes the changes in sugar shape during catalysis.

•

Stereoelectronic requirements for glycoside hydrolysis are discussed.

•

Major and emerging approaches to define conformational itineraries are reviewed.

•

New assignments of glycosidase conformational itineraries are summarized.

Glycoside hydrolases (GHs) are classified into >100 sequence-based families. These enzymes process a wide variety of complex carbohydrates with varying stereochemistry at the anomeric and other ring positions. The shapes that these sugars adopt upon binding to their cognate GHs, and the conformational changes that occur along the catalysis reaction coordinate is termed the conformational itinerary. Efforts to define the conformational itineraries of GHs have focussed upon the critical points of the reaction: substrate-bound (Michaelis), transition state, intermediate (if relevant) and product-bound. Recent approaches to defining conformational itineraries that marry X-ray crystallography of enzymes bound to ligands that mimic the critical points, along with advanced computational methods and kinetic isotope effects are discussed.