Molecular Mechanism of the Glycosylation Step Catalyzed by Golgi α-Mannosidase II: A QM/MM Metadynamics Investigation

An atomistic scale simulation study of structural properties in the silk-fibrohexamerin complex

The use of Bombyx mori silk fibroin in composite materials has been extensively explored in many studies, owing to its remarkable mechanical properties. Recently, the N-glycan-engineered P25 protein was utilized to improve the mechanical properties of silk. However, the mechanism by which N-glycan-engineered P25 protein enhances the mechanical properties of silk remains unclear. This study analyzed the interaction between the P25 protein and silkworm silk using quantum mechanics/molecular mechanics multiscale simulations and discovered stronger hydrogen bonding between the amorphous domain and the P25 protein. The results confirmed that glycoengineering of the mannose molecule in N-glycan in orders of three, five, and seven increased the hydrogen bonding of the amorphous structures. However, P25 has fewer binding interactions with the crystalline domain. Silk amino acids and mannose molecules were analyzed using QM simulations, and hydroxyl and charged amino acids in the amorphous domains were found to have relatively higher reactivity with mannose molecules in N-glycans than basic and aliphatic amino acids in the crystalline domain. This study demonstrates how the N-glycan-engineered P25 protein can improve the mechanical properties of silk fibroin and identifies a key factor for N-glycan-engineered proteins.

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.

RNA-seq analysis identifies age-dependent changes in expression of mRNAs - encoding N-glycosylation pathway enzymes in mouse gonadotropes

Follicle-stimulating hormone (FSH) is a glycoprotein that is assembled as a heterodimer of α/β subunits in gonadotropes. Each subunit contains two N-glycan chains. Our previous in vivo genetic studies identified that at least one N-glycan chain must be present on the FSHβ subunit for efficient FSH dimer assembly and secretion. Moreover, macroheterogeneity observed uniquely on human FSHβ results in ratiometric changes in age-specific FSH glycoforms, particularly during menopausal transition. Despite the recognition of many prominent roles of sugars on FSH including dimer assembly and secretion, serum half-life, receptor binding and signal transduction, the N-glycosylation machinery in gonadotropes has never been defined. Here, we used a mouse model in which gonadotropes are GFP-labeled in vivo and achieved rapid purification of GFP+ gonadotropes from pituitaries of female mice at reproductively young, middle, and old ages. We identified by RNA-seq analysis 52 mRNAs encoding N-glycosylation pathway enzymes expressed in 3- and 8-10-month-old mouse gonadotropes. We hierarchically mapped and localized the enzymes to distinct subcellular organelles within the N-glycosylation biosynthetic pathway. Of the 52 mRNAs, we found 27 mRNAs are differentially expressed between the 3- and 8-10-month old mice. We subsequently selected 8 mRNAs which showed varying changes in expression for confirmation of abundance in vivo via qPCR analysis, using more expanded aging time points with distinct 8-month and 14-month age groups. Real time qPCR analysis indicated dynamic changes in expression of N-glycosylation pathway enzyme-encoding mRNAs across the life span. Notably, computational analysis predicted the promoters of genes encoding these 8 mRNAs contain multiple high probability binding sites for estrogen receptor-1 and progesterone receptor. Collectively, our studies define the N-glycome and identify age-specific dynamic changes in mRNAs encoding N-glycosylation pathway enzymes in mouse gonadotropes. Our studies suggest the age-related decline in ovarian steroids may regulate expression of N-glycosylation enzymes in mouse gonadotropes and explain the age-related N-glycosylation shift previously observed on human FSHβ subunit in pituitaries of women.

In silico modelling of the function of disease-related CAZymes

In silico modelling of proteins comprises a diversity of computational tools aimed to obtain structural, electronic, and/or dynamic information about these biomolecules, capturing mechanistic details that are challenging to experimental approaches, such as elusive enzyme-substrate complexes, short-lived intermediates, and reaction transition states (TS). The present article gives the reader insight on the use of in silico modelling techniques to understand complex catalytic reaction mechanisms of carbohydrate-active enzymes (CAZymes), along with the underlying theory and concepts that are important in this field. We start by introducing the significance of carbohydrates in nature and the enzymes that process them, CAZymes, highlighting the conformational flexibility of their carbohydrate substrates. Three commonly used in silico methods (classical molecular dynamics (MD), hybrid quantum mechanics/molecular mechanics (QM/MM), and enhanced sampling techniques) are described for nonexpert readers. Finally, we provide three examples of the application of these methods to unravel the catalytic mechanisms of three disease-related CAZymes: β-galactocerebrosidase (GALC), responsible for Krabbe disease; α-mannoside β-1,6-N-acetylglucosaminyltransferase V (MGAT5), involved in cancer; and O-fucosyltransferase 1 (POFUT1), involved in several human diseases such as leukemia and the Dowling-Degos disease.

1,4-Dideoxy-1,4-imino-D- and L-lyxitol-based inhibitors bind to Golgi α-mannosidase II in different protonation forms

The development of effective inhibitors of Golgi α-mannosidase II (GMII, E.C.3.2.1.114) with minimal off-target effects on phylogenetically-related lysosomal α-mannosidase (LMan, E.C.3.2.1.24) is a complex task due to the complicated structural and chemical properties of their active sites. The pKa values (and also protonation forms in some cases) of several ionizable amino acids, such as Asp, Glu, His or Arg of enzymes, can be changed upon the binding of the inhibitor. Moreover, GMII and LMan work under different pH conditions. The pKa calculations on large enzyme-inhibitor complexes and FMO-PIEDA energy decomposition analysis were performed on the structures of selected inhibitors obtained from docking and hybrid QM/MM calculations. Based on the calculations, the roles of the amino group incorporated in the ring of the imino-D-lyxitol inhibitors and some ionizable amino acids of Golgi-type (Asp270-Asp340-Asp341 of Drosophila melanogaster α-mannosidase dGMII) and lysosomal-type enzymes (His209-Asp267-Asp268 of Canavalia ensiformis α-mannosidase, JBMan) were explained in connection with the observed inhibitory properties. The pyrrolidine ring of the imino-D-lyxitols prefers at the active site of dGMII the neutral form while in JBMan the protonated form, whereas that of imino-L-lyxitols prefers the protonation form in both enzymes. The calculations indicate that the binding mechanism of inhibitors to the active-site of α-mannosidases is dependent on the inhibitor structure and could be used to design new selective inhibitors of GMII. A series of novel synthetic N-substituted imino-D-lyxitols were evaluated with four enzymes from the glycoside hydrolase GH38 family (two of Golgi-type, Drosophila melanogaster GMIIb and Caenorhabditis elegans AMAN-2, and two of lysosomal-type, Drosophila melanogaster LManII and Canavalia ensiformis JBMan, enzymes). The most potent structures [N-9-amidinononyl and N-2-(1-naphthyl)ethyl derivatives] inhibited GMIIb (Ki = 40 nM) and AMAN-2 (Ki = 150 nM) with a weak selectivity index (SI) toward Golgi-type enzymes of IC50(LManII)/IC50(GMIIb) = 35 or IC50(JBMan)/IC50(AMAN-2) = 86. On the other hand, weaker micromolar inhibitors, such as N-2-naphthylmethyl or 4-iodobenzyl derivatives [IC50(GMIIb) = 2.4 μM and IC50 (AMAN-2) = 7.6 μM], showed a significant SI in the range from 111 to 812.

Enzymatic Hydrolysis of Human Milk Oligosaccharides. The Molecular Mechanism of Bifidobacterium Bifidum Lacto-N-biosidase

Bifidobacterium bifidum lacto-N-biosidase (LnbB) is a critical enzyme for the degradation of human milk oligosaccharides in the gut microbiota of breast-fed infants. Guided by recent crystal structures, we unveil its molecular mechanism of catalysis using QM/MM metadynamics. We show that the oligosaccharide substrate follows ¹S₃/^1,4B → [⁴E]^‡ → ⁴C₁/⁴H₅ and ⁴C₁/⁴H₅ → [⁴E/⁴H₅]^‡ → ^1,4B conformational itineraries for the two successive reaction steps, with reaction free energy barriers in agreement with experiments. The simulations also identify a critical histidine (His263) that switches between two orientations to modulate the pK_a of the acid/base residue, facilitating catalysis. The reaction intermediate of LnbB is best depicted as an oxazolinium ion, with a minor population of neutral oxazoline. The present study sheds light on the processing of oligosaccharides of the early life microbiota and will be useful for the engineering of LnbB and similar glycosidases for biocatalysis.

QM/MM Simulations of Enzymatic Hydrolysis of Cellulose: Probing the Viability of an Endocyclic Mechanism for an Inverting Cellulase

Glycoside hydrolases (GH) cleave carbohydrate glycosidic bonds and play pivotal roles in living organisms and in many industrial processes. Unlike acid-catalyzed hydrolysis of carbohydrates in solution, which can occur either via cyclic or acyclic oxocarbenium-like transition states, it is widely accepted that GH-catalyzed hydrolysis proceeds via a general acid mechanism involving a cyclic oxocarbenium-like transition state with protonation of the glycosidic oxygen. The GH45 subfamily C inverting endoglucanase from Phanerochaete chrysosporium (PcCel45A) defies the classical inverting mechanism as its crystal structure conspicuously lacks a general Asp or Glu base residue. Instead, PcCel45A has an Asn residue, a notoriously weak base in solution, as one of its catalytic residues at position 92. Moreover, unlike other inverting GHs, the relative position of the catalytic residues in PcCel45A impairs the proton abstraction from the nucleophilic water that attacks the anomeric carbon, a key step in the classical mechanism. Here, we investigate the viability of an endocyclic mechanism for PcCel45A using hybrid quantum mechanics/molecular mechanics (QM/MM) simulations, with the QM region treated with the self-consistent-charge density-functional tight-binding level of theory. In this mechanism, an acyclic oxocarbenium-like transition state is stabilized leading to the opening of the glucopyranose ring and formation of an unstable acyclic hemiacetal that can be readily decomposed into hydrolysis product. In silico characterization of the Michaelis complex shows that PcCel45A significantly restrains the sugar ring to the ⁴C₁ chair conformation at the −1 subsite of the substrate binding cleft, in contrast to the classical exocyclic mechanism in which ring puckering is critical. We also show that PcCel45A provides an environment where the catalytic Asn92 residue in its standard amide form participates in a cooperative hydrogen bond network resulting in its increased nucleophilicity due to an increased negative charge on the oxygen atom. Our results for PcCel45A suggest that carbohydrate hydrolysis catalyzed by GHs may take an alternative route from the classical mechanism.

Cysteine Nucleophiles in Glycosidase Catalysis: Application of a Covalent β-l-Arabinofuranosidase Inhibitor

The recent discovery of zinc-dependent retaining glycoside hydrolases (GHs), with active sites built around a Zn(Cys)₃ (Glu) coordination complex, has presented unresolved mechanistic questions. In particular, the proposed mechanism, depending on a Zn-coordinated cysteine nucleophile and passing through a thioglycosyl enzyme intermediate, remains controversial. This is primarily due to the expected stability of the intermediate C-S bond. To facilitate the study of this atypical mechanism, we report the synthesis of a cyclophellitol-derived β-l-arabinofuranosidase inhibitor, hypothesised to react with the catalytic nucleophile to form a non-hydrolysable adduct analogous to the mechanistic covalent intermediate. This β-l-arabinofuranosidase inhibitor reacts exclusively with the proposed cysteine thiol catalytic nucleophiles of representatives of GH families 127 and 146. X-ray crystal structures determined for the resulting adducts enable MD and QM/MM simulations, which provide insight into the mechanism of thioglycosyl enzyme intermediate breakdown. Leveraging the unique chemistry of cyclophellitol derivatives, the structures and simulations presented here support the assignment of a zinc-coordinated cysteine as the catalytic nucleophile and illuminate the finely tuned energetics of this remarkable metalloenzyme clan.

N-Linked Surface Glycan Biosynthesis, Composition, Inhibition, and Function in Cnidarian-Dinoflagellate Symbiosis

The success of symbioses between cnidarian hosts (e.g., corals and sea anemones) and micro-algal symbionts hinges on the molecular interactions that govern the establishment and maintenance of intracellular mutualisms. As a fundamental component of innate immunity, glycan-lectin interactions impact the onset of marine endosymbioses, but our understanding of the effects of cell surface glycome composition on symbiosis establishment remains limited. In this study, we examined the canonical N-glycan biosynthesis pathway in the genome of the dinoflagellate symbiont Breviolum minutum (family Symbiodiniaceae) and found it to be conserved with the exception of the transferase GlcNAc-TII (MGAT2). Using coupled liquid chromatography-mass spectrometry (LC-MS/MS), we characterized the cell surface N-glycan content of B. minutum, providing the first insight into the molecular composition of surface glycans in dinoflagellates. We then used the biosynthesis inhibitors kifunensine and swainsonine to alter the glycan composition of B. minutum. Successful high-mannose enrichment via kifunensine treatment resulted in a significant decrease in colonization of the model sea anemone Aiptasia (Exaiptasia pallida) by B. minutum. Hybrid glycan enrichment via swainsonine treatment, however, could not be confirmed and did not impact colonization. We conclude that functional Golgi processing of N-glycans is critical for maintaining appropriate cell surface glycan composition and for ensuring colonization success by B. minutum.

An Epoxide Intermediate in Glycosidase Catalysis

Retaining glycoside hydrolases cleave their substrates through stereochemical retention at the anomeric position. Typically, this involves two-step mechanisms using either an enzymatic nucleophile via a covalent glycosyl enzyme intermediate or neighboring-group participation by a substrate-borne 2-acetamido neighboring group via an oxazoline intermediate; no enzymatic mechanism with participation of the sugar 2-hydroxyl has been reported. Here, we detail structural, computational, and kinetic evidence for neighboring-group participation by a mannose 2-hydroxyl in glycoside hydrolase family 99 endo-α-1,2-mannanases. We present a series of crystallographic snapshots of key species along the reaction coordinate: a Michaelis complex with a tetrasaccharide substrate; complexes with intermediate mimics, a sugar-shaped cyclitol β-1,2-aziridine and β-1,2-epoxide; and a product complex. The 1,2-epoxide intermediate mimic displayed hydrolytic and transfer reactivity analogous to that expected for the 1,2-anhydro sugar intermediate supporting its catalytic equivalence. Quantum mechanics/molecular mechanics modeling of the reaction coordinate predicted a reaction pathway through a 1,2-anhydro sugar via a transition state in an unusual flattened, envelope (E ₃) conformation. Kinetic isotope effects (k _cat/K _M) for anomeric-²H and anomeric-¹³C support an oxocarbenium ion-like transition state, and that for C2-¹⁸O (1.052 ± 0.006) directly implicates nucleophilic participation by the C2-hydroxyl. Collectively, these data substantiate this unprecedented and long-imagined enzymatic mechanism.

Golgi Alpha-Mannosidase II as a Novel Biomarker Predicts Prognosis in Clear Cell Renal Cell Carcinoma

OBJECTIVE

Golgi alpha-mannosidase II (GM II) is one of the crucial enzymes in the process of N-glycan processing. The aim of our study was to examine the clinical significance of GM II in patients with clear cell renal cell carcinoma (ccRCC).

METHODS

Quantitative reverse transcription polymerase chain reaction analysis and immunohistochemical staining were used to analyze GM II expression in patients with ccRCC. The clinical data of 62 patients with ccRCC were collected to analyze the clinical significance of GM II. The clinical significance among GM II expression, clinicopathological staging, and histological grade of ccRCC was explored. Survival analyses were performed to identify the relevance between the expression of GM II and the overall survival of patients with ccRCC. A uni-/multivariate Cox regression model was used to detect risk factors affecting the prognosis of patients with ccRCC. Subsequently, the proliferation and migration of ccRCC cells were detected after transfecting with GM II-short hairpin RNA (shRNA).

RESULTS

The results of these comparisons suggested that GM II expression of ccRCC tissues was dramatically higher than that of para-carcinoma tissues (p < 0.05). GM II expression in the high-differentiation group was lower than that in the median- and low-differentiation groups (p < 0.05). GM II expression in stage I and II tissues was lower than that in stage III and IV tissues (p < 0.05). The expression levels of GM II in the group without lymph node metastasis were lower than those in the group with lymph node metastasis (p < 0.05). Survival analysis indicated that patients with ccRCC with high GM II expression generally had decreased overall survival. Uni-/multivariate Cox model analyses further suggested an association between GM II expression and prognosis of patients with breast cancer. High GM II expression is a potential and independent prognostic biomarker in ccRCC. The inhibition of GM II by transfecting with GM II-shRNA could reduce the proliferation and migration of ccRCC.

CONCLUSION

GM II expression in human ccRCC tissues was upregulated compared with that found in normal human renal tissue, and GM II may promote the progression and migration of ccRCC. Furthermore, the GM II gene may be used as a promising tumor marker for the diagnosis and prognosis of ccRCC.

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.

Identification of a potential allosteric site of Golgi α-mannosidase II using computer-aided drug design

Golgi α-mannosidase II (GMII) is a glycoside hydrolase playing a crucial role in the N-glycosylation pathway. In various tumour cell lines, the distribution of N-linked sugars on the cell surface is modified and correlates with the progression of tumour metastasis. GMII therefore is a possible molecular target for anticancer agents. Here, we describe the identification of a non-competitive GMII inhibitor using computer-aided drug design methods including identification of a possible allosteric binding site, pharmacophore search and virtual screening.

Comparing Hydrolysis and Transglycosylation Reactions Catalyzed by Thermus thermophilus β-Glycosidase. A Combined MD and QM/MM Study

The synthesis of oligosaccharides and other carbohydrate derivatives is of relevance for the advancement of glycosciences both at the fundamental and applied level. For many years, glycosyl hydrolases (GHs) have been explored to catalyze the synthesis of glycosidic bonds. In particular, retaining GHs can catalyze a transglycosylation (T) reaction that competes with hydrolysis (H). This has been done either employing controlled conditions in wild type GHs or by engineering new mutants. The goal, which is to increase the T/H ratio, has been achieved with moderate success in several cases despite the fact that the molecular basis for T/H modulation are unclear. Here we have used QM(DFT)/MM calculations to compare the glycosylation, hydrolysis and transglycosylation steps catalyzed by wild type Thermus thermophilus β-glycosidase (family GH1), a retaining glycosyl hydrolase for which a transglycosylation yield of 36% has been determined experimentally. The three transition states have a strong oxocarbenium character and ring conformations between ⁴H₃ and ⁴E. The atomic charges at the transition states for hydrolysis and transglycosylation are very similar, except for the more negative charge of the oxygen atom of water when compared to that of the acceptor Glc. The glycosylation transition state has a stronger S_N2 character than the deglycosylation ones and the proton transfer is less advanced. At the QM(PBE0/TZVP)/MM level, the TS for transglycosylation has shorter O4_GLC-C1_FUC (forming bond) distance and longer OE2_GLU338-C1_FUC (breaking) distance than the hydrolysis one, although the H_ACC proton is closer to the Glu164 base in the hydrolysis TS. The QM(SCC-DFTB)/MM free energy maxima show the inverted situation, although the hydrolysis TS presents significant structural fluctuations. The 3-OH_GLC group of the acceptor Glc (transglycosylation) and WAT432 (neighbor water in hydrolysis) are identified to stabilize the oxocarbenium transition states through interaction with O5_FUC and O4_FUC. The analysis of interaction suggests that perturbing the Glu392-Fuc interaction could increase the T/H ratio, either by direct mutation of this residue or indirectly as reported experimentally in the Asn390I and Phe401S cases. The molecular understanding of similarities and differences between hydrolysis and transglycosylation steps may be of help in the design of new biocatalysts for glycan synthesis.

Comparative studies of catalytic pathways for Streptococcus pneumoniae sialidases NanA, NanB and NanC

Streptococcus pneumoniae (S. pneumoniae) is a leading human pathogen, which takes large responsibility for severe otitis media, acute meningitis and septicaemia. It encodes up to three distinct sialidases: NanA, NanB and NanC, which are promising drug targets. Recent experimental studies have shown that these three sialidases might work together up to the ultimate step, where NanA and NanB produce N-acetylneuraminic acid (Neu5Ac) and 2,7-anhydro-Neu5Ac following the functions of sialidase and intramolecular trans-sialidase, whilst NanC carries on a ping-pong mechanism that produces or removes 2-deoxy-2,3-didehydro-Neu5AC. It is intriguing that these sialidases have similar active sites but operate via three distinct reaction pathways. To clarify this issue, herein we present the first systematic computational investigation on the catalytic pathways for S. pneumoniae NanA, NanB and NanC based on combined quantum mechanics/molecular mechanics simulations, and propose the most preferred routes for the three S. pneumoniae sialidases. Our findings support the mechanisms of NanA and NanC that were proposed by previous experimental studies, whereas the role of water in NanB was found to differ slightly from our current understandings. The mechanistic insights obtained from this work are expected to assist in the design of potent inhibitors targeting these key enzymes for therapeutic applications.

A QM/MM approach on the structural and stereoelectronic factors governing glycosylation by GTF-SI fromStreptococcus mutans

This manuscript contains novel insights into the reaction mechanism catalyzed by GTF-SI. Structural and electronic features of the system are revealed, such as the strong hydrogen bond depicted above.

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.

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.

Improving the activity of Trichoderma reesei cel7B through stabilizing the transition state

Trichoderma reesei (Tr.) cellulases, which convert cellulose to reducing sugars, are a promising catalyst used in the lignocellulosic biofuel production. Improving Tr. cellulases activity, though very difficult, is highly desired due to the recalcitrance of lignocellulose. Meanwhile, it is preferable to enhance the cellulase's promiscuity so that substrates other than cellulose can also be hydrolyzed. In this work, an attempt is made to improve the catalytic activity of a major endogluanase Tr. Cel7B against xylan which crosslinks with cellulose in lignocellulose. By using quantum mechanics/molecular mechanics (QM/MM) molecular dynamics (MD) simulations, the transition state of the xylo-oligosaccharide hydrolysis is identified. Then, mutations are introduced and their effect on the transition state stabilization is ranked based on the free energy calculations. Seven top ranked mutants are evaluated experimentally. Three mutants A208Q, A222D, and G230R show a higher activity than the wild-type Tr. Cel7B in the hydrolysis of xylan (by up to 47%) as well as filter paper (by up to 50%). The combination of the single mutants can further improve the enzyme activity. Our work demonstrates that the free energy method is effective in engineering the Tr. Cel7B activity against xylan and cellulose, and thus may also be useful for improving the activity of other Tr. cellulases. Biotechnol. Bioeng. 2016;113: 1171-1177. © 2015 Wiley Periodicals, Inc.

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.

Mutational analysis of divalent metal ion binding in the active site of class II α-mannosidase from Sulfolobus solfataricus

Mutational analysis of Sulfolobus solfataricus class II α-mannosidase was focused on side chains that interact with the hydroxyls of the -1 mannosyl of the substrate (Asp-534) or form ligands to the active site divalent metal ion (His-228 and His-533) judged from crystal structures of homologous enzymes. D534A and D534N appeared to be completely inactive. When compared to the wild-type enzyme, the mutant enzymes in general showed only small changes in K(M) for the substrate, p-nitrophenyl-α-mannoside, but elevated activation constants, K(A), for the divalent metal ion (Co²⁺, Zn²⁺, Mn²⁺, or Cd²⁺). Some mutant enzyme forms displayed an altered preference for the metal ion compared to that of the wild type-enzyme. Furthermore, the H228Q, H533E, and H533Q enzymes were inhibited at increasing Zn²⁺ concentrations. The catalytic rate was reduced for all enzymes compared to that of the wild-type enzyme, although less dramatically with some activating metal ions. No major differences in the pH dependence between wild-type and mutant enzymes were found in the presence of different metal ions. The pH optimum was 5, but enzyme instability was observed at pH <4.5; therefore, only the basic limb of the bell-shaped pH profile was analyzed.

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.

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.

Molecular mechanism of a hotdog-fold acyl-CoA thioesterase

Thioesterases are enzymes that hydrolyze thioester bonds between a carbonyl group and a sulfur atom. They catalyze key steps in fatty acid biosynthesis and metabolism, as well as polyketide biosynthesis. The reaction molecular mechanism of most hotdog-fold acyl-CoA thioesterases remains unknown, but several hypotheses have been put forward in structural and biochemical investigations. The reaction of a human thioesterase (hTHEM2), representing a thioesterase family with a hotdog fold where a coenzyme A moiety is cleaved, was simulated by quantum mechanics/molecular mechanics metadynamics techniques to elucidate atomic and electronic details of its mechanism, its transition-state conformation, and the free energy landscape of the process. A single-displacement acid-base-like mechanism, in which a nucleophilic water molecule is activated by an aspartate residue acting as a base, was found, confirming previous experimental proposals. The results provide unambiguous evidence of the formation of a tetrahedral-like transition state. They also explain the roles of other conserved active-site residues during the reaction, especially that of a nearby histidine/serine pair that protonates the thioester sulfur atom, the participation of which could not be elucidated from mutation analyses alone.

Nylon-Oligomer Hydrolase Promoting Cleavage Reactions in Unnatural Amide Compounds

The active site of 6-aminohexanoate-dimer hydrolase, a nylon-6 byproduct-degrading enzyme with a β-lactamase fold, possesses a Ser112/Lys115/Tyr215 catalytic triad similar to the one of penicillin-recognizing family of serine-reactive hydrolases but includes a unique Tyr170 residue. By using a reactive quantum mechanics/molecular mechanics (QM/MM) approach, we work out its catalytic mechanism and related functional/structural specificities. At variance with other peptidases, we show that the involvement of Tyr170 in the enzyme-substrate interactions is responsible for a structural variation in the substrate-binding state. The acylation via a tetrahedral intermediate is the rate-limiting step, with a free-energy barrier of ∼21 kcal/mol, driven by the catalytic triad Ser112, Lys115, and Tyr215, acting as a nucleophile, general base, and general acid, respectively. The functional interaction of Tyr170 with this triad leads to an efficient disruption of the tetrahedral intermediate, promoting a conformational change of the substrate favorable for proton donation from the general acid.

Combined inhibitor free-energy landscape and structural analysis reports on the mannosidase conformational coordinate

Mannosidases catalyze the hydrolysis of a diverse range of polysaccharides and glycoconjugates, and the various sequence-based mannosidase families have evolved ingenious strategies to overcome the stereoelectronic challenges of mannoside chemistry. Using a combination of computational chemistry, inhibitor design and synthesis, and X-ray crystallography of inhibitor/enzyme complexes, it is demonstrated that mannoimidazole-type inhibitors are energetically poised to report faithfully on mannosidase transition-state conformation, and provide direct evidence for the conformational itinerary used by diverse mannosidases, including β-mannanases from families GH26 and GH113. Isofagomine-type inhibitors are poor mimics of transition-state conformation, owing to the high energy barriers that must be crossed to attain mechanistically relevant conformations, however, these sugar-shaped heterocycles allow the acquisition of ternary complexes that span the active site, thus providing valuable insight into active-site residues involved in substrate recognition.

Structural analysis and insights into the glycon specificity of the rice GH1 Os7BGlu26 β-D-mannosidase

Rice Os7BGlu26 is a GH1 family glycoside hydrolase with a threefold higherkcat/Kmvalue for 4-nitrophenyl β-D-mannoside (4NPMan) compared with 4-nitrophenyl β-D-glucoside (4NPGlc). To investigate its selectivity for β-D-mannoside and β-D-glucoside substrates, the structures of apo Os7BGlu26 at a resolution of 2.20 Å and of Os7BGlu26 with mannose at a resolution of 2.45 Å were elucidated from isomorphous crystals in space groupP212121. The (β/α)8-barrel structure is similar to other GH1 family structures, but with a narrower active-site cleft. The Os7BGlu26 structure with D-mannose corresponds to a product complex, with β-D-mannose in the1S5skew-boat conformation. Docking of the1S3,1S5,2SOand3S1pyranose-ring conformations of 4NPMan and 4NPGlc substrates into the active site of Os7BGlu26 indicated that the lowest energies were in the1S5and1S3skew-boat conformations. Comparison of these docked conformers with other rice GH1 structures revealed differences in the residues interacting with the catalytic acid/base between enzymes with and without β-D-mannosidase activity. The mutation of Tyr134 to Trp in Os7BGlu26 resulted in similarkcat/Kmvalues for 4NPMan and 4NPGlc, while mutation of Tyr134 to Phe resulted in a 37-fold higherkcat/Kmfor 4NPMan than 4NPGlc. Mutation of Cys182 to Thr decreased both the activity and the selectivity for β-D-mannoside. It was concluded that interactions with the catalytic acid/base play a significant role in glycon selection.

Formation of a covalent glycosyl-enzyme species in a retaining glycosyltransferase

Elusive glycosyl-enzyme adduct: Using classical MD simulations and QM/MM metadynamics, the long-time sought glycosyl-enzyme covalent intermediate of a retaining glycosyltransferase, with a putative nucleophile residue in the active site, has been trapped (MD=molecular dynamics; QM/MM=quantum mechanics/molecular mechanics).