QM/MM free-energy simulations of reaction in Serratia marcescens Chitinase B reveal the protonation state of Asp142 and the critical role of Tyr214

A trapped covalent intermediate as a key catalytic element in the hydrolysis of a GH3 β-glucosidase: An X-ray crystallographic and biochemical study

Glycoside hydrolases (GHs) are industrially important enzymes that hydrolyze glycosidic bonds in glycoconjugates. In this study, we found a GH3 β-glucosidase (CcBgl3B) from Cellulosimicrobium cellulans sp. 21 was able to selectively hydrolyze the β-1,6-glucosidic bond linked glucose of ginsenosides. X-ray crystallographic studies of the ligand complex ginsenoside-specific β-glucosidase provided a novel finding that support the catalytic mechanism of GH3. The substrate was clearly identified within the catalytic center of wild-type CcBgl3B, revealing that the C1 atom of the glucose was covalently bound to the Oδ1 group of the conserved catalytic nucleophile Asp264 as an enzyme-glycosyl intermediate. The glycosylated Asp264 could be identified by mass spectrometry. Through site-directed mutagenesis studies with Asp264, it was found that the covalent intermediate state formed by Asp264 and the substrate was critical for catalysis. In addition, Glu525 variants (E525A, E525Q and E525D) showed no or marginal activity against pNPβGlc; thus, this residue could supply a proton for the reaction. Overall, our study provides an insight into the catalytic mechanism of the GH3 enzyme CcBgl3B.

Multiscale QM/MM Simulations Identify the Roles of Asp239 and 1-OH···Nucleophile in Transition State Stabilization in Arabidopsis thaliana Cell-Wall Invertase 1

Arabidopsis thaliana cell-wall invertase 1 (AtCWIN1), a key enzyme in sucrose metabolism in plants, catalyzes the hydrolysis of sucrose into fructose and glucose. AtCWIN1 belongs to the glycoside hydrolase GH-J clan, where two carboxylate residues (Asp23 and Glu203 in AtCWIN1) are well documented as a nucleophile and an acid/base catalyst. However, details at the atomic level about the role of neighboring residues and enzyme-substrate interactions during catalysis are not fully understood. Here, quantum mechanical/molecular mechanical (QM/MM) free-energy simulations were carried out to clarify the origin of the observed decreased rates in Asp239Ala, Asp239Asn, and Asp239Phe in AtCWIN1 compared to the wild type and delineate the role of Asp239 in catalysis. The glycosylation and deglycosylation steps were considered in both wild type and mutants. Deglycosylation is predicted to be the rate-determining step in the reaction, with a calculated overall free-energy barrier of 15.9 kcal/mol, consistent with the experimental barrier (15.3 kcal/mol). During the reaction, the -1 furanosyl ring underwent a conformational change corresponding to 3E ↔ [E2]⧧ ↔ 1E according to the nomenclature of saccharide structures along the full catalytic reaction. Asp239 was found to stabilize not only the transition state but also the fructosyl-enzyme intermediate, which explains findings from previous structural and mutagenesis experiments. The 1-OH···nucleophile interaction has been found to provide an important contribution to the transition state stabilization, with a contribution of ∼7 kcal/mol, and affected glycosylation more significantly than deglycosylation. This study provides molecular insights that improve the current understanding of sucrose binding and hydrolysis in members of clan GH-J, which may benefit protein engineering research. Finally, a rationale on the sucrose inhibitor configuration in chicory 1-FEH IIa, proposed a long time ago in the literature, is also provided based on the QM/MM calculations.

Synergic chitin degradation by Streptomyces sp. SCUT-3 chitinases and their applications in chitinous waste recycling and pathogenic fungi biocontrol

The genus Streptomyces comprises the most important chitin decomposers in soil and revealing their chitinolytic machinery is beneficial for the conversion of chitinous wastes. Streptomyces sp. SCUT-3, a chitin-hydrolyzing and a robust feather-degrading bacterium, was isolated previously. The potential chitin-degrading enzymes produced by SCUT-3 were analyzed in the present study. Among these enzymes, three chitinases were successfully expressed in Pichia pastoris at comparatively high yields of 4.8 U/mL (SsExoChi18A), 11.2 U/mL (SsExoChi18B), and 17.8 U/mL (SsEndoChi19). Conserved motifs and constructive 3D structures of these three exo- and endochitinases were also analyzed. These chitinases hydrolyzed colloidal chitin to chitin oligomers. SsExoChi18A showed apparent synergic effects with SsEndoChi19 in colloidal chitin and shrimp shell hydrolysis, with an improvement of 29.3 % and 124.9 %, respectively. Compared with SsExoChi18B and SsEndoChi19, SsExoChi18A exhibited the strongest antifungal effects against four plant pathogens by inhibiting mycelial growth and spore germination. This study provided good candidates for chitinous waste-processing enzymes and antifungal biocontrol agents. These synergic chitin-degrading enzymes of SCUT-3 are good targets for its further genetical modification to construct super chitinous waste-degrading bacteria with strong abilities to hydrolyze both protein and chitin, thereby providing a direction for the future path of the chitinous waste recycling industry.

Mechanistic Modeling of Monoglyceride Lipase Covalent Modification Elucidates the Role of Leaving Group Expulsion and Discriminates Inhibitors with High and Low Potency

Inhibition of monoglyceride lipase (MGL), also known as monoacylglycerol lipase (MAGL), has emerged as a promising approach for treating neurological diseases. To gain useful insights in the design of agents with balanced potency and reactivity, we investigated the mechanism of MGL carbamoylation by the reference triazole urea SAR629 (IC₅₀ = 0.2 nM) and two recently described inhibitors featuring a pyrazole (IC₅₀ = 1800 nM) or a 4-cyanopyrazole (IC₅₀ = 8 nM) leaving group (LG), using a hybrid quantum mechanics/molecular mechanics (QM/MM) approach. Opposite to what was found for substrate 2-arachidonoyl-sn-glycerol (2-AG), covalent modification of MGL by azole ureas is controlled by LG expulsion. Simulations indicated that changes in the electronic structure of the LG greatly affect reaction energetics with triazole and 4-cyanopyrazole inhibitors following a more accessible carbamoylation path compared to the unsubstituted pyrazole derivative. The computational protocol provided reaction barriers able to discriminate between MGL inhibitors with different potencies. These results highlight how QM/MM simulations can contribute to elucidating structure–activity relationships and provide insights for the design of covalent inhibitors.

A QM/MM Evaluation of the Missing Step in the Reduction Mechanism of HMG-CoA by Human HMG-CoA Reductase

Statins are important drugs in the regulation of cholesterol levels in the human body that have as a primary target the enzyme β-hydroxy-β-methylglutaryl-CoA reductase (HMGR). This enzyme plays a crucial role in the mevalonate pathway, catalyzing the four-electron reduction of HMG-CoA to mevalonate. A second reduction step of this reaction mechanism has been the subject of much speculation in the literature, with different conflicting theories persisting to the present day. In this study, the different mechanistic hypotheses were evaluated with atomic-level detail through a combination of molecular dynamics simulations (MD) and quantum mechanics/molecular mechanics (QM/MM) calculations. The obtained Gibbs free activation and Gibbs free reaction energy (15 kcal mol−1 and −40 kcal mol−1) show that this hydride step takes place with the involvement of a cationic His405 and Lys639, and a neutral Glu98, while Asp715 remains in an anionic state. The results provide an atomic-level portrait of this step, clearly demonstrating the nature and protonation state of the amino acid residues involved, the energetics associated, and the structure and charge of the key participating atoms in the several intermediate states, finally elucidating this missing step.

Active site dynamics and catalytic mechanism in arabinan hydrolysis catalyzed by GH43 endo-arabinanase from QM/MM molecular dynamics simulation and potential energy surface

The endo-1,5-α-L-arabinanases, belonging to glycoside hydrolase family 43 (GH43), catalyse the hydrolysis of α-1,5-arabinofuranosidic bonds in arabinose-containing polysaccharides. These enzymes are proposed targets for industrial and medical applications. Here, molecular dynamics (MD), potential energy surface and free energy (potential of mean force) simulations are undertaken using hybrid quantum mechanical/molecular mechanical (QM/MM) potentials to understand the active site dynamics, catalytic mechanism and the electrostatic influence of active site residues of the GH43 endo-arabinanase from G. stearothermophilus. The calculated results give support to the single-displacement mechanism proposed for the inverting GH43 enzymes: first a proton is transferred from the general acid E201 to the substrate, followed by a nucleophilic attack by water, activated by the general base D27, on the anomer carbon. A conformational change (²E ↔E₃ ↔ ⁴E) in the -1 sugar ring is observed involving a transition state featuring an oxocarbenium ion character. Residues D87, K106, H271 are highlighted as potential targets for future mutation experiments in order to increase the efficiency of the reaction. To our knowledge, this is the first QM/MM study providing molecular insights into the glycosidic bond hydrolysis of a furanoside substrate by an inverting GH in solution.Communicated by Ramaswamy H. Sarma.

Functional characterization of a novel thermophilic exo-arabinanase from Thermothielavioides terrestris

Arabinanases from glycoside hydrolase family GH93 are enzymes with exo-activity that hydrolyze the α-1,5 bonds between arabinose residues present on arabinan. Currently, several initiatives aiming to use byproducts rich in arabinan such as pectin and sugar beet pulp as raw material to produce various compounds of interest are being developed. However, it is necessary to use robust enzymes that have an optimal performance under pH and temperature conditions used in the industrial processes. In this work, the first GH93 from the thermophilic fungus Thermothielavioides terrestris (Abn93T) was heterologously expressed in Aspergillus nidulans, purified and biochemically characterized. The enzyme is a thermophilic glycoprotein (optimum activity at 70 °C) with prolonged stability in acid pHs (4.0 to 6.5). The presence of glycosylation affected slightly the hydrolytic capacity of the enzyme, which was further increased by 34% in the presence of 1 mM CoCl₂. Small-angle X-ray scattering results show that Abn93T is a globular-like-shaped protein with a slight bulge at one end. The hydrolytic mechanism of the enzyme was elucidated using capillary zone electrophoresis and molecular docking calculations. Abn93T has an ability to produce (in synergism with arabinofuranosidases) arabinose and arabinobiose from sugar beet arabinan, which can be explored as fermentable sugars and prebiotics. KEY POINTS: • Thermophilic exo-arabinanase from family GH93 • Molecular basis of arabinan depolymerization.

Current Status of Carbohydrates Information in the Protein Data Bank

Carbohydrates are well known for their physicochemical, biological, functional, and therapeutic characteristics. Unfortunately, their chemical nature imposes severe challenges for the structural elucidation of these phenomena, impairing not only the depth of our understanding of carbohydrates but also the development of new biotechnological and therapeutic applications based on these molecules. In the recent past, the amount of structural information, obtained mainly from X-ray crystallography, has increased progressively, as well as its quality. In this context, the current work presents a global analysis of the carbohydrate information available in the Protein Data Bank (PDB). From high quality structures, it is clear that most of the data are highly concentrated on a few sets of residue types, on their monosaccharidic forms, and connected by a small diversity of glycosidic linkages. The geometries of these linkages can be mostly associated with the types of linkages instead of residues, while the level of puckering distortion was characterized, quantified, and located in a pseudorotational equilibrium landscape, not only to local minima but also to transitional states. These qualitative and quantitative analyses offer a global picture of the carbohydrate structural content in the PDB, potentially supporting the building of new models for carbohydrate-related biological phenomena at the atomistic level, including new developments on force field parameters.

Design and SAR Analysis of Covalent Inhibitors Driven by Hybrid QM/MM Simulations

Quantum mechanics/molecular mechanics (QM/MM) hybrid technique is emerging as a reliable computational method to investigate and characterize chemical reactions occurring in enzymes. From a drug discovery perspective, a thorough understanding of enzyme catalysis appears pivotal to assist the design of inhibitors able to covalently bind one of the residues belonging to the enzyme catalytic machinery. Thanks to the current advances in computer power, and the availability of more efficient algorithms for QM-based simulations, the use of QM/MM methodology is becoming a viable option in the field of covalent inhibitor design. In the present review, we summarized our experience in the field of QM/MM simulations applied to drug design problems which involved the optimization of agents working on two well-known drug targets, namely fatty acid amide hydrolase (FAAH) and epidermal growth factor receptor (EGFR). In this context, QM/MM simulations gave valuable information in terms of geometry (i.e., of transition states and metastable intermediates) and reaction energetics that allowed to correctly predict inhibitor binding orientation and substituent effect on enzyme inhibition. What is more, enzyme reaction modelling with QM/MM provided insights that were translated into the synthesis of new covalent inhibitor featured by a unique combination of intrinsic reactivity, on-target activity, and selectivity.

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 novel ring-shaped reaction pathway with interconvertible intermediates in chitinase A as revealed by QM/MM simulation combined with a one-dimensional projection technique

Efficient sampling achieved by the use of a one-dimensional projection technique reveals the catalytic mechanism of chitinase A from Serratia marcescens.

Mechanism of the Escherichia coli MltE lytic transglycosylase, the cell-wall-penetrating enzyme for Type VI secretion system assembly

Lytic transglycosylases (LTs) catalyze the non-hydrolytic cleavage of the bacterial cell wall by an intramolecular transacetalization reaction. This reaction is critically and broadly important in modifications of the bacterial cell wall in the course of its biosynthesis, recycling, manifestation of virulence, insertion of structural entities such as the flagellum and the pili, among others. The first QM/MM analysis of the mechanism of reaction of an LT, that for the Escherichia coli MltE, is undertaken. The study reveals a conformational itinerary consistent with an oxocarbenium-like transition state, characterized by a pivotal role for the active-site glutamic acid in proton transfer. Notably, an oxazolinium intermediate, as a potential intermediate, is absent. Rather, substrate-assisted catalysis is observed through a favorable dipole provided by the N-acetyl carbonyl group of MurNAc saccharide. This interaction stabilizes the incipient positive charge development in the transition state. This mechanism coincides with near-synchronous acetal cleavage and acetal formation.

Data characterizing the energetics of enzyme-catalyzed hydrolysis and transglycosylation reactions by DFT cluster model calculations

The data presented in this paper are related to the research article entitled “QM/MM modeling of the hydrolysis and transfructosylation reactions of fructosyltransferase from Aspergillus japonicas, an enzyme that produces prebiotic fructooligosaccharide” (Jitonnom et al., 2018) [1]. This paper presents the procedure and data for characterizing the whole relative energy profiles of hydrolysis and transglycosylation reactions whose elementary steps differ in chemical composition. The data also reflects the choices of the QM cluster model, the functional/basis set method and the equations in determining the reaction energetics.

Deciphering the chemoselectivity of nickel-dependent quercetin 2,4-dioxygenase

QM/MM calculations were performed to elucidate the reaction mechanism and chemoselectivity of 2,4-QueD. The protonation state of the first-shell ligand Glu74 plays an important role in dictating the selectivity.

Quantum Mechanics/Molecular Mechanics Simulations Identify the Ring-Opening Mechanism of Creatininase

Creatininase catalyzes the conversion of creatinine (a biosensor for kidney function) to creatine via a two-step mechanism: water addition followed by ring opening. Water addition is common to other known cyclic amidohydrolases, but the precise mechanism for ring opening is still under debate. The proton donor in this step is either His178 or a water molecule bound to one of the metal ions, and the roles of His178 and Glu122 are unclear. Here, the two possible reaction pathways have been fully examined by means of combined quantum mechanics/molecular mechanics simulations at the SCC-DFTB/CHARMM22 level of theory. The results indicate that His178 is the main catalytic residue for the whole reaction and explain its role as proton shuttle during the ring-opening step. In the first step, His178 provides electrostatic stabilization to the gem-diolate tetrahedral intermediate. In the second step, His178 abstracts the hydroxyl proton of the intermediate and delivers it to the cyclic amide nitrogen, leading to ring opening. The latter is the rate-limiting step with a free energy barrier of 18.5 kcal/mol, in agreement with the experiment. We find that Glu122 must be protonated during the enzyme reaction, so that it can form a stable hydrogen bond with its neighboring water molecule. Simulations of the E122Q mutant showed that this replacement disrupts the H-bond network formed by three conserved residues (Glu34, Ser78, and Glu122) and water, increasing the energy barrier. Our computational studies provide a comprehensive explanation for previous structural and kinetic observations, including why the H178A mutation causes a complete loss of activity but the E122Q mutation does not.

Identification and characterization of a core fucosidase from the bacterium Elizabethkingia meningoseptica

All reported α-l-fucosidases catalyze the removal of nonreducing terminal l-fucoses from oligosaccharides or their conjugates, while having no capacity to hydrolyze core fucoses in glycoproteins directly. Here, we identified an α-fucosidase from the bacterium Elizabethkingia meningoseptica with catalytic activity against core α-1,3-fucosylated substrates, and we named it core fucosidase I (cFase I). Using site-specific mutational analysis, we found that three acidic residues (Asp-242, Glu-302, and Glu-315) in the predicted active pocket are critical for cFase I activity, with Asp-242 and Glu-315 acting as a pair of classic nucleophile and acid/base residues and Glu-302 acting in an as yet undefined role. These findings suggest a catalytic mechanism for cFase I that is different from known α-fucosidase catalytic models. In summary, cFase I exhibits glycosidase activity that removes core α-1,3-fucoses from substrates, suggesting cFase I as a new tool for glycobiology, especially for studies of proteins with core fucosylation.

Engineered N-acetylhexosamine-active enzymes in glycoscience

BACKGROUND

In recent years, enzymes modifying N-acetylhexosamine substrates have emerged in numerous theoretical studies as well as practical applications from biology, biomedicine, and biotechnology. Advanced enzyme engineering techniques converted them into potent synthetic instruments affording a variety of valuable glycosides.

SCOPE OF REVIEW

This review presents the diversity of engineered enzymes active with N-acetylhexosamine carbohydrates: from popular glycoside hydrolases and glycosyltransferases to less known oxidases, epimerases, kinases, sulfotransferases, and acetylases. Though hydrolases in natura, engineered chitinases, β-N-acetylhexosaminidases, and endo-β-N-acetylglucosaminidases were successfully employed in the synthesis of defined natural and derivatized chitooligomers and in the remodeling of N-glycosylation patterns of therapeutic antibodies. The genes of various N-acetylhexosaminyltransferases were cloned into metabolically engineered microorganisms for producing human milk oligosaccharides, Lewis X structures, and human-like glycoproteins. Moreover, mutant N-acetylhexosamine-active glycosyltransferases were applied, e.g., in the construction of glycomimetics and complex glycostructures, industrial production of low-lactose milk, and metabolic labeling of glycans. In the synthesis of biotechnologically important compounds, several innovative glycoengineered systems are presented for an efficient bioproduction of GlcNAc, UDP-GlcNAc, N-acetylneuraminic acid, and of defined glycosaminoglycans.

MAJOR CONCLUSIONS

The above examples demonstrate that engineering of N-acetylhexosamine-active enzymes was able to solve complex issues such as synthesis of tailored human-like glycoproteins or industrial-scale production of desired oligosaccharides. Due to the specific catalytic mechanism, mutagenesis of these catalysts was often realized through rational solutions.

GENERAL SIGNIFICANCE

Specific N-acetylhexosamine glycosylation is crucial in biological, biomedical and biotechnological applications and a good understanding of its details opens new possibilities in this fast developing area of glycoscience.

Nanoscale Engineering of Designer Cellulosomes

Biocatalysts showcase the upper limit obtainable for high-speed molecular processing and transformation. Efforts to engineer functionality in synthetic nanostructured materials are guided by the increasing knowledge of evolving architectures, which enable controlled molecular motion and precise molecular recognition. The cellulosome is a biological nanomachine, which, as a fundamental component of the plant-digestion machinery from bacterial cells, has a key potential role in the successful development of environmentally-friendly processes to produce biofuels and fine chemicals from the breakdown of biomass waste. Here, the progress toward so-called "designer cellulosomes", which provide an elegant alternative to enzyme cocktails for lignocellulose breakdown, is reviewed. Particular attention is paid to rational design via computational modeling coupled with nanoscale characterization and engineering tools. Remaining challenges and potential routes to industrial application are put forward.

Discerning the catalytic mechanism of Staphylococcus aureus sortase A with QM/MM free energy calculations

Sortases are key virulence factors in Gram-positive bacteria. These enzymes embed surface proteins in the cell wall through a transpeptidation reaction that involves recognizing a penta-peptide "sorting signal" in a target protein, cleaving it, and covalently attaching it to a second substrate that is later inserted into the cell wall. Although well studied, several aspects of the mechanism by which sortases perform these functions remains unclear. In particular, experiments have revealed two potential sorting signal binding motifs: a "Threonine-Out" (Thr-Out) structure in which the catalytically critical threonine residues protrudes into solution, and a "Threonine-In" (Thr-In) configuration in which this residue inserts into the binding site. To determine which of these is the biologically relevant state, we have performed a series of conventional and hybrid quantum mechanics/molecular mechanics (QM/MM) molecular dynamics simulations of the Staphylococcus aureus sortase A (SrtA) enzyme bound to a sorting signal substrate. Through the use of multi-dimensional metadynamics, our simulations were able to both map the acylation mechanism of SrtA in the Thr-In and Thr-Out states, as well as determine the free energy minima and barriers along these reactions. Results indicate that in both states the catalytic mechanisms are similar, however the free energy barriers are lower in the Thr-In configuration, suggesting that Thr-In is the catalytically relevant state. This has important implications for advancing our understanding of the mechanisms of sortase enzymes, as well we for future structure based drug design efforts aimed at inhibiting sortase function in vivo.

Aromatic-Mediated Carbohydrate Recognition in Processive Serratia marcescens Chitinases

Microorganisms use a host of enzymes, including processive glycoside hydrolases, to deconstruct recalcitrant polysaccharides to sugars. Processive glycoside hydrolases closely associate with polymer chains and repeatedly cleave glycosidic linkages without dissociating from the crystalline surface after each hydrolytic step; they are typically the most abundant enzymes in both natural secretomes and industrial cocktails by virtue of their significant hydrolytic potential. The ubiquity of aromatic residues lining the enzyme catalytic tunnels and clefts is a notable feature of processive glycoside hydrolases. We hypothesized that these aromatic residues have uniquely defined roles, such as substrate chain acquisition and binding in the catalytic tunnel, that are defined by their local environment and position relative to the substrate and the catalytic center. Here, we investigated this hypothesis with variants of Serratia marcescens family 18 processive chitinases ChiA and ChiB. We applied molecular simulation and free energy calculations to assess active site dynamics and ligand binding free energies. Isothermal titration calorimetry provided further insight into enthalpic and entropic contributions to ligand binding free energy. Thus, the roles of six aromatic residues, Trp-167, Trp-275, and Phe-396 in ChiA, and Trp-97, Trp-220, and Phe-190 in ChiB, have been examined. We observed that point mutation of the tryptophan residues to alanine results in unfavorable changes in the free energy of binding relative to wild-type. The most drastic effects were observed for residues positioned at the "entrances" of the deep substrate-binding clefts and known to be important for processivity. Interestingly, phenylalanine mutations in ChiA and ChiB had little to no effect on chito-oligomer binding, in accordance with the limited effects of their removal on chitinase functionality.

QM/MM study of the mechanism of reduction of 3-hydroxy-3-methylglutaryl coenzyme A catalyzed by human HMG-CoA reductase

Detailing with atomistic resolution the reaction mechanism of human HMG-CoA reductase (HMG-CoA-R) might provide valuable insights for the development of new cholesterol-lowering drugs.

New insights into the enzymatic mechanism of human chitotriosidase (CHIT1) catalytic domain by atomic resolution X-ray diffraction and hybrid QM/MM

Chitotriosidase (CHIT1) is a human chitinase belonging to the highly conserved glycosyl hydrolase family 18 (GH18). GH18 enzymes hydrolyze chitin, an N-acetylglucosamine polymer synthesized by lower organisms for structural purposes. Recently, CHIT1 has attracted attention owing to its upregulation in immune-system disorders and as a marker of Gaucher disease. The 39 kDa catalytic domain shows a conserved cluster of three acidic residues, Glu140, Asp138 and Asp136, involved in the hydrolysis reaction. Under an excess concentration of substrate, CHIT1 and other homologues perform an additional activity, transglycosylation. To understand the catalytic mechanism of GH18 chitinases and the dual enzymatic activity, the structure and mechanism of CHIT1 were analyzed in detail. The resolution of the crystals of the catalytic domain was improved from 1.65 Å (PDB entry 1waw) to 0.95-1.10 Å for the apo and pseudo-apo forms and the complex with chitobiose, allowing the determination of the protonation states within the active site. This information was extended by hybrid quantum mechanics/molecular mechanics (QM/MM) calculations. The results suggest a new mechanism involving changes in the conformation and protonation state of the catalytic triad, as well as a new role for Tyr27, providing new insights into the hydrolysis and transglycosylation activities.