Lluís Raich

A Trapped Covalent Intermediate of a Glycoside Hydrolase on the Pathway to Transglycosylation. Insights from Experiments and Quantum Mechanics/Molecular Mechanics Simulations

The conversion of glycoside hydrolases (GHs) into transglycosylases (TGs), i.e., from enzymes that hydrolyze carbohydrates to enzymes that synthesize them, represents a promising solution for the large-scale synthesis of complex carbohydrates for biotechnological purposes. However, the lack of knowledge about the molecular details of transglycosylation hampers the rational design of TGs. Here we present the first crystallographic structure of a natural glycosyl-enzyme intermediate (GEI) of Saccharomyces cerevisiae Gas2 in complex with an acceptor substrate and demonstrate, by means of quantum mechanics/molecular mechanics metadynamics simulations, that it is tuned for transglycosylation (ΔG(⧧) = 12 kcal/mol). The 2-OH···nucleophile interaction is found to be essential for catalysis: its removal raises the free energy barrier significantly (11 and 16 kcal/mol for glycosylation and transglycosylation, respectively) and alters the conformational itinerary of the substrate (from (4)C1 → [(4)E](⧧) → (1,4)B/(4)E to (4)C1 → [(4)H3](⧧) → (4)C1). Our results suggest that changes in the interactions involving the 2-position could have an impact on the transglycosylation activity of several GHs.

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 1C4 → 3H4‡ → 3S1 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.

Contribution of Shape and Charge to the Inhibition of a Family GH99 endo-α-1,2-Mannanase

Inhibitor design incorporating features of the reaction coordinate and transition-state structure has emerged as a powerful approach for the development of enzyme inhibitors. Such inhibitors find use as mechanistic probes, chemical biology tools, and therapeutics. Endo-α-1,2-mannosidases and endo-α-1,2-mannanases, members of glycoside hydrolase family 99 (GH99), are interesting targets for inhibitor development as they play key roles in N-glycan maturation and microbiotal yeast mannan degradation, respectively. These enzymes are proposed to act via a 1,2-anhydrosugar “epoxide” mechanism that proceeds through an unusual conformational itinerary. Here, we explore how shape and charge contribute to binding of diverse inhibitors of these enzymes. We report the synthesis of neutral dideoxy, glucal and cyclohexenyl disaccharide inhibitors, their binding to GH99 endo-α-1,2-mannanases, and their structural analysis by X-ray crystallography. Quantum mechanical calculations of the free energy landscapes reveal how the neutral inhibitors provide shape but not charge mimicry of the proposed intermediate and transition state structures. Building upon the knowledge of shape and charge contributions to inhibition of family GH99 enzymes, we design and synthesize α-Man-1,3-noeuromycin, which is revealed to be the most potent inhibitor (KD 13 nM for Bacteroides xylanisolvens GH99 enzyme) of these enzymes yet reported. This work reveals how shape and charge mimicry of transition state features can enable the rational design of potent inhibitors.

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 maritimaTmGH1 β-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.

Precise Probing of Residue Roles by Post-Translational β,γ-C,N Aza-Michael Mutagenesis in Enzyme Active Sites

Biomimicry valuably allows the understanding of the essential chemical components required to recapitulate biological function, yet direct strategies for evaluating the roles of amino acids in proteins can be limited by access to suitable, subtly-altered unnatural variants. Here we describe a strategy for dissecting the role of histidine residues in enzyme active sites using unprecedented, chemical, post-translational side-chain-β,γ C–N bond formation. Installation of dehydroalanine (as a “tag”) allowed the testing of nitrogen conjugate nucleophiles in “aza-Michael”-1,4-additions (to “modify”). This allowed the creation of a regioisomer of His (iso-His, Hisiso) linked instead through its pros-Nπ atom rather than naturally linked via C4, as well as an aza-altered variant aza-Hisiso. The site-selective generation of these unnatural amino acids was successfully applied to probe the contributing roles (e.g., size, H-bonding) of His residues toward activity in the model enzymes subtilisin protease from Bacillus lentus and Mycobacterium tuberculosis pantothenate synthetase.

Conformational analysis of the mannosidase inhibitor kifunensine: a quantum mechanical and structural approach

The varied yet family‐specific conformational pathways used by individual glycoside hydrolases (GHs) offer a tantalising prospect for the design of tightly binding and specific enzyme inhibitors. A cardinal example of a GH‐family‐specific inhibitor, and one that finds widespread practical use, is the natural product kifunensine, which is a low‐nanomolar inhibitor that is selective for GH family 47 inverting α‐mannosidases. Here we show, through quantum‐mechanical approaches, that kifunensine is restrained to a “ring‐flipped” 1C4 conformation with another accessible, but higher‐energy, region around the 1,4B conformation. The conformations of kifunensine in complex with a range of GH47 enzymes—including an atomic‐level resolution (1 Å) structure of kifunensine with Caulobacter sp. CkGH47 reported herein and with GH family 38 and 92 α‐mannosidases—were mapped onto the kifunensine free‐energy landscape. These studies revealed that kifunensine has the ability to mimic the product state of GH47 enzymes but cannot mimic any conformational states relevant to the reaction coordinate of mannosidases from other families.

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 4C1 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.

The molecular mechanism of the ligand exchange reaction of an antibody against a glutathione-coated gold cluster

The labeling of proteins with heavy atom clusters is of paramount importance in biomedical research, but its detailed molecular mechanism remains unknown. Here we uncover it for the particular case of the anti-influenza N9 neuraminidase NC10 antibody against a glutathione-coated gold cluster by means of ab initio QM/MM calculations. We show that the labeling reaction follows an associative double SN2-like reaction mechanism, involving a proton transfer, with low activation barriers only if one of the two distinct peptide/peptidic ligands (the one that occupies the side position) is substituted. Positively charged residues in the vicinity of the incoming thiol result in strong interactions between the antibody and the AuMPC, favoring the ligand exchange reaction for suitable protein mutants. These results pave the way for future investigations aimed at engineering biomolecules to increase their reactivity towards a desired gold atom cluster.

The Molecular Mechanism of Substrate Recognition and Catalysis of the Membrane Acyltransferase PatA from Mycobacteria

Glycolipids play a central role in a variety of important biological processes in all living organisms. PatA is a membrane acyltransferase involved in the biosynthesis of phosphatidyl-myo-inositol mannosides (PIMs), key structural elements, and virulence factors of Mycobacterium tuberculosis. PatA catalyzes the transfer of a palmitoyl moiety from palmitoyl-CoA to the 6-position of the mannose ring linked to the 2-position of inositol in PIM1/PIM2. We report here the crystal structure of PatA in the presence of 6-O-palmitoyl-α-d-mannopyranoside, unraveling the acceptor binding mechanism. The acceptor mannose ring localizes in a cavity at the end of a surface-exposed long groove where the active site is located, whereas the palmitate moiety accommodates into a hydrophobic pocket deeply buried in the α/β core of the protein. Both fatty acyl chains of the PIM2 acceptor are essential for the reaction to take place, highlighting their critical role in the generation of a competent active site. By the use of combined structural and quantum-mechanics/molecular-mechanics (QM/MM) metadynamics, we unravel the catalytic mechanism of PatA at the atomic-electronic level. Our study provides a detailed structural rationale for a stepwise reaction, with the generation of a tetrahedral transition state for the rate-determining step. Finally, the crystal structure of PatA in the presence of β-d-mannopyranose and palmitate suggests an inhibitory mechanism for the enzyme, providing exciting possibilities for inhibitor design and the discovery of chemotherapeutic agents against this major human pathogen.

Enzymatic Cleavage of Glycosidic Bonds

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.