Conformational Itinerary of Sucrose During Hydrolysis by Retaining Amylosucrase

Resveratrol glucosylation by GTF-SI from Streptococcus mutans: computational insights into a GH70 family enzyme

Resveratrol, a polyphenolic compound known for its health benefits but limited by poor water solubility and low bioavailability, represents a valuable substrate for glucosylation by carbohydrate-active enzymes such as glucosyltransferase-SI (GTF-SI). Using quantum mechanics/molecular mechanics (QM/MM) calculations and molecular dynamics simulations, this study reveals the atomic scale dynamics of resveratrol glucosylation by wild-type GTF-SI. This enzyme exhibited an energy barrier of 8.8 kcal mol-1 and an exothermic process, both consistent with experimental data of similar enzymes. We report a concerted and synchronous reaction mechanism for the catalytic step, characterized by an oxocarbenium ion-like transition state, and elucidate a conformational itinerary of the glucosyl moiety (4H3/E3) → [E3]‡ → 4C1, which aligns with the consistent patterns observed across enzymes of the GH13 and GH70 families. A key interaction was observed between Asp477 and the OH group on carbon 6 of the glucosyl moiety, together with a 2.0 kcal mol-1 transition state stabilization by three water molecules within the active site. Comparative insights with the previously studied Q345F SP enzyme system shed light on the unique and common features that govern transglucosylation reactions. Importantly, the calculated activation barriers strongly support the capability of GTF-SI to facilitate resveratrol glucosylation. This study advances our understanding of the transglucosylation reaction and opens up new ways for the glycodiversification of organic compounds such as polyphenols, thus expanding their potential applications in the food, cosmetic, and pharmaceutical industries.

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.

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.

Role of Enzyme and Active Site Conformational Dynamics in the Catalysis by α-Amylase Explored with QM/MM Molecular Dynamics

We assessed enzyme:substrate conformational dynamics and the rate-limiting glycosylation step of a human pancreatic α-amylase:maltopentose complex. Microsecond molecular dynamics simulations suggested that the distance of the catalytic Asp197 nucleophile to the anomeric carbon of the buried glucoside is responsible for most of the enzyme active site fluctuations and that both Asp197 and Asp300 interact the most with the buried glucoside unit. The buried glucoside binds either in a ⁴C₁ chair or ²S_O skew conformations, both of which can change to TS-like conformations characteristic of retaining glucosidases. Starting from four distinct enzyme:substrate complexes, umbrella sampling quantum mechanics/molecular mechanics simulations (converged within less than 1 kcal·mol^-1 within a total simulation time of 1.6 ns) indicated that the reaction occurrs with a Gibbs barrier of 13.9 kcal·mol ^-1, in one asynchronous concerted step encompassing an acid-base reaction with Glu233 followed by a loose S_N2-like nucleophilic substitution by the Asp197. The transition state is characterized by a ²H₃ half-chair conformation of the buried glucoside that quickly changes to the E₃ envelope conformation preceding the attack of the anomeric carbon by the Asp197 nucleophile. Thermodynamic analysis of the reaction supported that a water molecule tightly hydrogen bonded to the glycosidic oxygen of the substrate at the reactant state (∼1.6 Å) forms a short hydrogen bond with Glu233 at the transition state (∼1.7 Å) and lowers the Gibbs barrier in over 5 kcal·mol^-1. The resulting Asp197-glycosyl was mostly found in the ⁴C₁ conformation, although the more endergonic B_3,O conformation was also observed. Altogether, the combination of short distances for the acid-base reaction with the Glu233 and for the nucleophilic attack by the Asp197 nucleophile and the availability of water within hydrogen bonding distance of the glycosidic oxygen provides a reliable criteria to identify reactive conformations of α-amylase complexes.

A Single Hydrogen Bond Controls the Selectivity of Transglycosylation vs Hydrolysis in Family 13 Glycoside Hydrolases

Converting glycoside hydrolases (GHs) from hydrolytic to synthetic enzymes via transglycosylation is a long-standing goal for the biosynthesis of complex carbohydrates. However, the molecular determinants for the selectivity of transglycosylation (T) vs hydrolysis (H) are still not fully unraveled. Herein, we show experimentally that mutation of one active site residue can switch the enzyme activity between hydrolysis and transglycosylation in two highly homologous GHs. Further QM/MM simulations reveal that the mutation modulates the T vs H reaction barriers via the presence/absence of a single H-bond with the nucleophile Asp. Such a H-bond controls the product selectivity via a dual effect: on one hand, it facilitates the breaking of the glycosyl-enzyme intermediate. On the other, it displaces the sugar acceptor, resulting in a reduced affinity and significant steric repulsion for transglycosylation. These findings expand our understanding of the molecular mechanisms that modulate the T/H balance in GHs.

Glucosylation mechanism of resveratrol through the mutant Q345F sucrose phosphorylase from the organism Bifidobacterium adolescentis: a computational study

Mainly due to their great antioxidant, anti-inflammatory and anticancer capacities, natural polyphenolic compounds have many properties with important applications in the food, cosmetic and pharmaceutical industries. Unfortunately, these molecules have very low water solubility and bioavailability. Glucosylation of polyphenols is an excellent alternative to overcome these drawbacks. Specifically, for the natural polyphenol resveratrol this process is very inefficiently performed by the native enzyme sucrose phosphorylase (BaSP) from the organism Bifidobacterium adolescentis (4%). However, the Q345F point mutation of the sucrose phosphorylase (BaSP Q345F) has been shown to achieve 97% monoglucosylation for the same substrate and the mechanism is still unknown. This report presents an analysis of MD simulations performed with the BaSP Q345F and BaSP systems in complex with resveratrol monoglucoside, followed by a study of the transglucosylation reaction of the mutant enzyme BaSP Q345F with resveratrol through the QM/MM hybrid method. With respect to the MD simulations, both protein structures showed greater similarity to the phosphate-binding conformation, and a larger active site and conformational changes in certain structures were found for the mutant system compared with the native enzyme; all this is in agreement with experimental data. With regard to the QM/MM calculations, the structure of an oxocarbenium ion-like transition state and the minimum energy adiabatic path (MEP) that connects the reactants with the products were obtained with a 20.3 kcal mol^-1 energy barrier, which fits within the known experimental range for this type of enzyme. Finally, the analyses performed along the MEP suggest a concerted but asynchronous mechanism. In particular, they show that the interactions involving the residues of the catalytic triad (Asp¹⁹², Glu²³², and Asp²⁹⁰) together with two water molecules at the active site strongly contribute to the stabilization of the transition state. The understanding of this glucosylation mechanism of the polyphenol resveratrol carried out by the mutant sucrose phosphorylase enzyme presented in this work could serve as the basis for subsequent studies on related carbohydrate-active enzymes.

Stereoselective synthesis of a 4-⍺-glucoside of valienamine and its X-ray structure in complex with Streptomyces coelicolor GlgE1-V279S

Glycoside hydrolases (GH) are a large family of hydrolytic enzymes found in all domains of life. As such, they control a plethora of normal and pathogenic biological functions. Thus, understanding selective inhibition of GH enzymes at the atomic level can lead to the identification of new classes of therapeutics. In these studies, we identified a 4-⍺-glucoside of valienamine (8) as an inhibitor of Streptomyces coelicolor (Sco) GlgE1-V279S which belongs to the GH13 Carbohydrate Active EnZyme family. The results obtained from the dose-response experiments show that 8 at a concentration of 1000 µM reduced the enzyme activity of Sco GlgE1-V279S by 65%. The synthetic route to 8 and a closely related 4-⍺-glucoside of validamine (7) was achieved starting from readily available D-maltose. A key step in the synthesis was a chelation-controlled addition of vinylmagnesium bromide to a maltose-derived enone intermediate. X-ray structures of both 7 and 8 in complex with Sco GlgE1-V279S were solved to resolutions of 1.75 and 1.83 Å, respectively. Structural analysis revealed the valienamine derivative 8 binds the enzyme in an E₂ conformation for the cyclohexene fragment. Also, the cyclohexene fragment shows a new hydrogen-bonding contact from the pseudo-diaxial C(3)-OH to the catalytic nucleophile Asp 394 at the enzyme active site. Asp 394, in fact, forms a bidentate interaction with both the C(3)-OH and C(7)-OH of the inhibitor. In contrast, compound 7 disrupts the catalytic sidechain interaction network of Sco GlgE1-V279S via steric interactions resulting in a conformation change in Asp 394. These findings will have implications for the design other aminocarbasugar-based GH13-inhibitors and will be useful for identifying more potent and selective inhibitors.

Glycosyl Oxocarbenium Ions: Structure, Conformation, Reactivity, and Interactions

Carbohydrates (glycans, saccharides, and sugars) are essential molecules in all domains of life. Research on glycoscience spans from chemistry to biomedicine, including material science and biotechnology. Access to pure and well-defined complex glycans using synthetic methods depends on the success of the employed glycosylation reaction. In most cases, the mechanism of the glycosylation reaction is believed to involve the oxocarbenium ion. Understanding the structure, conformation, reactivity, and interactions of this glycosyl cation is essential to predict the outcome of the reaction. In this Account, building on our contributions on this topic, we discuss the theoretical and experimental approaches that have been employed to decipher the key features of glycosyl cations, from their structures to their interactions and reactivity.We also highlight that, from a chemical perspective, the glycosylation reaction can be described as a continuum, from unimolecular SN1 with naked oxocarbenium cations as intermediates to bimolecular SN2-type mechanisms, which involve the key role of counterions and donors. All these factors should be considered and are discussed herein. The importance of dissociative mechanisms (involving contact ion pairs, solvent-separated ion pairs, solvent-equilibrated ion pairs) with bimolecular features in most reactions is also highlighted.The role of theoretical calculations to predict the conformation, dynamics, and reactivity of the oxocarbenium ion is also discussed, highlighting the advances in this field that now allow access to the conformational preferences of a variety of oxocarbenium ions and their reactivities under SN1-like conditions.Specifically, the ground-breaking use of superacids to generate these cations is emphasized, since it has permitted characterization of the structure and conformation of a variety of glycosyl oxocarbenium ions in superacid solution by NMR spectroscopy.We also pay special attention to the reactivity of these glycosyl ions, which depends on the conditions, including the counterions, the possible intra- or intermolecular participation of functional groups that may stabilize the cation and the chemical nature of the acceptor, either weak or strong nucleophile. We discuss recent investigations from different experimental perspectives, which identified the involved ionic intermediates, estimating their lifetimes and reactivities and studying their interactions with other molecules. In this context, we also emphasize the relationship between the chemical methods that can be employed to modulate the sensitivity of glycosyl cations and the way in which glycosyl modifying enzymes (glycosyl hydrolases and transferases) build and cleave glycosidic linkages in nature. This comparison provides inspiration on the use of molecules that regulate the stability and reactivity of glycosyl cations.

Computational modeling of carbohydrate processing enzymes reactions

Carbohydrate processing enzymes are of biocatalytic interest. Glycoside hydrolases and the recently discovered lytic polysaccharide monooxygenase for their use in biomass degradation to obtain biofuels or valued chemical entities. Glycosyltransferases or engineered glycosidases and phosphorylases for the synthesis of carbohydrates and glycosylated products. Quantum mechanics-molecular mechanics (QM/MM) methods are highly contributing to establish their different chemical reaction mechanisms. Other computational methods are also used to study enzyme conformational changes, ligand pathways, and processivity, e.g. for processive glycosidases like cellobiohydrolases. There is still a long road to travel to fully understand the role of conformational dynamics in enzyme activity and also to disclose the variety of reaction mechanisms these enzymes employ. Additionally, computational tools for enzyme engineering are beginning to be applied to evaluate substrate specificity or aid in the design of new biocatalysts with increased thermostability or tailored activity, a growing field where molecular modeling is finding its way.

Activity-Based Protein Profiling of Retaining α-Amylases in Complex Biological Samples

Amylases are key enzymes in the processing of starch in many kingdoms of life. They are important catalysts in industrial biotechnology where they are applied in, among others, food processing and the production of detergents. In man amylases are the first enzymes in the digestion of starch to glucose and arguably also the preferred target in therapeutic strategies aimed at the treatment of type 2 diabetes patients through down-tuning glucose assimilation. Efficient and sensitive assays that report selectively on retaining amylase activities irrespective of the nature and complexity of the biomaterial studied are of great value both in finding new and effective human amylase inhibitors and in the discovery of new microbial amylases with potentially advantageous features for biotechnological application. Activity-based protein profiling (ABPP) of retaining glycosidases is inherently suited for the development of such an assay format. We here report on the design and synthesis of 1,6-epi-cyclophellitol-based pseudodisaccharides equipped with a suite of reporter entities and their use in ABPP of retaining amylases from human saliva, murine tissue as well as secretomes from fungi grown on starch. The activity and efficiency of the inhibitors and probes are substantiated by extensive biochemical analysis, and the selectivity for amylases over related retaining endoglycosidases is validated by structural studies.

Structure-function analysis of silkworm sucrose hydrolase uncovers the mechanism of substrate specificity in GH13 subfamily 17 exo-α-glucosidases

The domestic silkworm Bombyx mori expresses two sucrose-hydrolyzing enzymes, BmSUH and BmSUC1, belonging to glycoside hydrolase family 13 subfamily 17 (GH13_17) and GH32, respectively. BmSUH has little activity on maltooligosaccharides, whereas other insect GH13_17 α-glucosidases are active on sucrose and maltooligosaccharides. Little is currently known about the structural mechanisms and substrate specificity of GH13_17 enzymes. In this study, we examined the crystal structures of BmSUH without ligands; in complexes with substrates, products, and inhibitors; and complexed with its covalent intermediate at 1.60-1.85 Å resolutions. These structures revealed that the conformations of amino acid residues around subsite -1 are notably different at each step of the hydrolytic reaction. Such changes have not been previously reported among GH13 enzymes, including exo- and endo-acting hydrolases, such as α-glucosidases and α-amylases. Amino acid residues at subsite +1 are not conserved in BmSUH and other GH13_17 α-glucosidases, but subsite -1 residues are absolutely conserved. Substitutions in three subsite +1 residues, Gln¹⁹¹, Tyr²⁵¹, and Glu⁴⁴⁰, decreased sucrose hydrolysis and increased maltase activity of BmSUH, indicating that these residues are key for determining its substrate specificity. These results provide detailed insights into structure-function relationships in GH13 enzymes and into the molecular evolution of insect GH13_17 α-glucosidases.

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.