Structural Snapshots and Loop Dynamics along the Catalytic Cycle of Glycosyltransferase GpgS

In Vitro and In Silico Explorations of the Protein Conformational Changes of Corynebacterium glutamicum MshA, a Model Retaining GT-B Glycosyltransferase

MshA is a GT-B glycosyltransferase catalyzing the first step in the biosynthesis of mycothiol. While many GT-B enzymes undergo an open-to-closed transition, MshA is unique because its 97° rotation is beyond the usual range of 10-25°. Molecular dynamics (MD) simulations were carried out for MshA in both ligand bound and unbound states to investigate the effect of ligand binding on localized protein dynamics and its conformational free energy landscape. Simulations showed that both the unliganded "opened" and liganded "closed" forms of the enzyme sample a wide degree of dihedral angles and interdomain distances with relatively low overlapping populations. Calculation of the free energy surface using replica exchange MD for the apo "opened" and an artificial generated apo "closed" structure revealed overlaps in the geometries sampled, allowing calculation of a barrier of 2 kcal/mol for the open-to-closed transition in the absence of ligands. MD simulations of fully liganded MshA revealed a smaller sampling of the dihedral angles. The localized protein fluctuation changes suggest that UDP-GlcNAc binding activates the motions of loops in the 1-l-myo-inositol-1-phosphate (I1P)-binding site despite little change in the interactions with UDP-GlcNAc. Circular dichroism, intrinsic fluorescence spectroscopy, and mutagenesis studies were used to confirm the ligand-induced structural changes in MshA. The results support a proposed mechanism where UDP-GlcNAc binds with rigid interactions to the C-terminal domain of MshA and activates flexible loops in the N-terminal domain for binding and positioning of I1P. This model can be used for future structure-based drug development of inhibitors of the mycothiol biosynthetic pathway.

HDX-MS Reveals Substrate-Dependent, Localized EX1 Conformational Dynamics in the Retaining GT-B Glycosyltransferase, MshA

Glycosyltransferases (GTs) are well-characterized with respect to static 3D structures and molecular dynamics simulations, but there is a lack of reports on in-solution dynamics on time scales relevant to turnover. Here, backbone amide hydrogen/deuterium exchange followed by mass spectrometry (HDX-MS) was used to investigate the in-solution dynamics of the model retaining GT MshA from Corynebacterium glutamicum (CgMshA). CgMshA has a GT-B fold and catalyzes the transfer of N-acetyl-glucosamine (GlcNAc) from UDP-GlcNAc to l-myo-inositol-1-phosphate in the first step in mycothiol biosynthesis. HDX-MS results identify several key regions of conformational changes in response to UDP-GlcNAc binding, including residues 159-198 in the N-terminal domain and residues 323-354 in the C-terminal domain. These regions also exhibited substrate-dependent EX1 exchange kinetics consistent with conformational tension on the milliseconds to seconds time scale. A potential source of this conformational change is the flexible β4/α5 loop in the C-terminal domain, which sits at the interface of the two domains and likely interacts with the GlcNAc ring of UDP-GlcNAc. In contrast to UDP-GlcNAc, the UDP-CgMshA product complex exhibited severe decreases in deuterium incorporation, suggesting a less dynamic conformation. The HDX-MS results are complemented by solvent viscosity effects of 1.8-2.3 on the CgMshA kcat value, which are consistent with product release as a rate-determining step and possibly a direct role for protein dynamics in catalysis. The identification of in-solution dynamics that are sensitive to substrate binding allows for the proposal of a more detailed mechanism in CgMshA including conformation tension between the donor sugar and the flexible C-terminal domain β4/α5 loop providing sufficient conformational sampling for substrate-assisted catalysis to occur.

Subfunctionalization of a monolignol to a phytoalexin glucosyltransferase is accompanied by substrate inhibition

Uridine diphosphate-dependent glycosyltransferases (UGTs) mediate the glycosylation of plant metabolites, thereby altering their physicochemical properties and bioactivities. Plants possess numerous UGT genes, with the encoded enzymes often glycosylating multiple substrates and some exhibiting substrate inhibition kinetics, but the biological function and molecular basis of these phenomena are not fully understood. The promiscuous monolignol/phytoalexin glycosyltransferase NbUGT72AY1 exhibits substrate inhibition (K_i) at 4 μM scopoletin, whereas the highly homologous monolignol StUGT72AY2 is inhibited at 190 μM. We therefore used hydrogen/deuterium exchange mass spectrometry and structure-based mutational analyses of both proteins and introduced NbUGT72AY1 residues into StUGT72AY2 and vice versa to study promiscuity and substrate inhibition of UGTs. A single F87I and chimeric mutant of NbUGT72AY1 showed significantly reduced scopoletin substrate inhibition, whereas its monolignol glycosylation activity was almost unaffected. Reverse mutations in StUGT72AY2 resulted in increased scopoletin glycosylation, leading to enhanced promiscuity, which was accompanied by substrate inhibition. Studies of 3D structures identified open and closed UGT conformers, allowing visualization of the dynamics of conformational changes that occur during catalysis. Previously postulated substrate access tunnels likely serve as drainage channels. The results suggest a two-site model in which the second substrate molecule binds near the catalytic site and blocks product release. Mutational studies showed that minor changes in amino acid sequence can enhance the promiscuity of the enzyme and add new capabilities such as substrate inhibition without affecting existing functions. The proposed subfunctionalization mechanism of expanded promiscuity may play a role in enzyme evolution and highlights the importance of promiscuous enzymes in providing new functions.

Docking-guided rational engineering of a macrolide glycosyltransferase glycodiversifies epothilone B

Glycosyltransferases typically display acceptor substrate flexibility but more stringent donor specificity. BsGT-1 is a highly effective glycosyltransferase to glycosylate macrolides, including epothilones, promising antitumor compounds. Here, we show that BsGT-1 has three major regions significantly influencing the glycodiversification of epothilone B based on structural molecular docking, "hot spots" alanine scanning, and site saturation mutagenesis. Mutations in the PSPG-like motif region and the C2 loop region are more likely to expand donor preference; mutations of the flexible N3 loop region located at the mouth of the substrate-binding cavity produce novel epothilone oligosaccharides. These "hot spots" also functioned in homologues of BsGT-1. The glycosides showed significantly enhanced water solubility and decreased cytotoxicity, although the glycosyl appendages of epothilone B also reduced drug permeability and attenuated antitumor efficacy. This study laid a foundation for the rational engineering of other GTs to synthesize valuable small molecules.

Ligand-Induced Conformational and Dynamical Changes in a GT-B Glycosyltransferase: Molecular Dynamics Simulations of Heptosyltransferase I Complexes

Understanding the dynamical motions and ligand recognition motifs of heptosyltransferase I (HepI) can be critical to discerning the behavior of other glycosyltransferase (GT) enzymes. Prior studies in our lab have demonstrated that GTs in the GT-B structural class, which are characterized by their connection of two Rossman-like domains by a linker region, have conserved structural fold and dynamical motions, despite low sequence homology, therefore making discoveries found in HepI transferable to other GT-B enzymes. Through molecular dynamics simulations and ligand binding free energy analysis of HepI in the apo and bound complexes (for all kinetically relevant combinations of the native substrates/products), we have determined the energetically favored enzymatic pathway for ligand binding and release. Our principal component, dynamic cross correlation, and network analyses of the simulations have revealed correlated motions involving residues within the N-terminal domain communicating with C-terminal domain residues via both proximal amino acid residues and also functional groups of the bound substrates. Analyses of the structural changes, energetics of substrate/product binding, and changes in pK_a have elucidated a variety of inter and intradomain interactions that are critical for enzyme catalysis. These data corroborate our experimental observations of protein conformational changes observed in both presteady state kinetic and circular dichroism analyses of HepI. These simulations provided invaluable structural insights into the regions involved in HepI conformational rearrangement upon ligand binding. Understanding the specific interactions governing conformational changes is likely to enhance our efforts to develop novel dynamics disrupting inhibitors against GT-B structural enzymes in the future.

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.

Conserved Conformational Hierarchy across Functionally Divergent Glycosyltransferases of the GT-B Structural Superfamily as Determined from Microsecond Molecular Dynamics

It has long been understood that some proteins undergo conformational transitions en route to the Michaelis Complex to allow chemistry. Examination of crystal structures of glycosyltransferase enzymes in the GT-B structural class reveals that the presence of ligand in the active site triggers an open-to-closed conformation transition, necessary for their catalytic functions. Herein, we describe microsecond molecular dynamics simulations of two distantly related glycosyltransferases that are part of the GT-B structural superfamily, HepI and GtfA. Simulations were performed using the open and closed conformations of these unbound proteins, respectively, and we sought to identify the major dynamical modes and communication networks that interconnect the open and closed structures. We provide the first reported evidence within the scope of our simulation parameters that the interconversion between open and closed conformations is a hierarchical multistep process which can be a conserved feature of enzymes of the same structural superfamily. Each of these motions involves of a collection of smaller molecular reorientations distributed across both domains, highlighting the complexities of protein dynamic involved in the interconversion process. Additionally, dynamic cross-correlation analysis was employed to explore the potential effect of distal residues on the catalytic efficiency of HepI. Multiple distal nonionizable residues of the C-terminal domain exhibit motions anticorrelated to positively charged residues in the active site in the N-terminal domain involved in substrate binding. Mutations of these residues resulted in a reduction in negatively correlated motions and an altered enzymatic efficiency that is dominated by lower K_m values with k_cat effectively unchanged. The findings suggest that residues with opposing conformational motions involved in the opening and closing of the bidomain HepI protein can allosterically alter the population and conformation of the “closed” state, essential to the formation of the Michaelis complex. The stabilization effects of these mutations likely equally influence the energetics of both the ground state and the transition state of the catalytic reaction, leading to the unaltered k_cat. Our study provides new insights into the role of conformational dynamics in glycosyltransferase’s function and new modality to modulate enzymatic efficiency.

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.

Dissecting the Structural and Chemical Determinants of the "Open-to-Closed" Motion in the Mannosyltransferase PimA from Mycobacteria

The phosphatidyl-myo-inositol mannosyltransferase A (PimA) is an essential peripheral membrane glycosyltransferase that initiates the biosynthetic pathway of phosphatidyl-myo-inositol mannosides (PIMs), key structural elements and virulence factors of Mycobacterium tuberculosis. PimA undergoes functionally important conformational changes, including (i) α-helix-to-β-strand and β-strand-to-α-helix transitions and (ii) an "open-to-closed" motion between the two Rossmann-fold domains, a conformational change that is necessary to generate a catalytically competent active site. In previous work, we established that GDP-Man and GDP stabilize the enzyme and facilitate the switch to a more compact active state. To determine the structural contribution of the mannose ring in such an activation mechanism, we analyzed a series of chemical derivatives, including mannose phosphate (Man-P) and mannose pyrophosphate-ribose (Man-PP-RIB), and additional GDP derivatives, such as pyrophosphate ribose (PP-RIB) and GMP, by the combined use of X-ray crystallography, limited proteolysis, circular dichroism, isothermal titration calorimetry, and small angle X-ray scattering methods. Although the β-phosphate is present, we found that the mannose ring, covalently attached to neither phosphate (Man-P) nor PP-RIB (Man-PP-RIB), does promote the switch to the active compact form of the enzyme. Therefore, the nucleotide moiety of GDP-Man, and not the sugar ring, facilitates the "open-to-closed" motion, with the β-phosphate group providing the high-affinity binding to PimA. Altogether, the experimental data contribute to a better understanding of the structural determinants involved in the "open-to-closed" motion not only observed in PimA but also visualized and/or predicted in other glycosyltransfeases. In addition, the experimental data might prove to be useful for the discovery and/or development of PimA and/or glycosyltransferase inhibitors.

Deep evolutionary analysis reveals the design principles of fold A glycosyltransferases

Glycosyltransferases (GTs) are prevalent across the tree of life and regulate nearly all aspects of cellular functions. The evolutionary basis for their complex and diverse modes of catalytic functions remain enigmatic. Here, based on deep mining of over half million GT-A fold sequences, we define a minimal core component shared among functionally diverse enzymes. We find that variations in the common core and emergence of hypervariable loops extending from the core contributed to GT-A diversity. We provide a phylogenetic framework relating diverse GT-A fold families for the first time and show that inverting and retaining mechanisms emerged multiple times independently during evolution. Using evolutionary information encoded in primary sequences, we trained a machine learning classifier to predict donor specificity with nearly 90% accuracy and deployed it for the annotation of understudied GTs. Our studies provide an evolutionary framework for investigating complex relationships connecting GT-A fold sequence, structure, function and regulation.

Carbohydrates are one of the major groups of large biological molecules that regulate nearly all aspects of life. Yet, unlike DNA or proteins, carbohydrates are made without a template to follow. Instead, these molecules are built from a set of sugar-based building blocks by the intricate activities of a large and diverse family of enzymes known as glycosyltransferases.

An incomplete understanding of how glycosyltransferases recognize and build diverse carbohydrates presents a major bottleneck in developing therapeutic strategies for diseases associated with abnormalities in these enzymes. It also limits efforts to engineer these enzymes for biotechnology applications and biofuel production.

Taujale et al. have now used evolutionary approaches to map the evolution of a major subset of glycosyltransferases from species across the tree of life to understand how these enzymes evolved such precise mechanisms to build diverse carbohydrates. First, a minimal structural unit was defined based on being shared among a group of over half a million unique glycosyltransferase enzymes with different activities. Further analysis then showed that the diverse activities of these enzymes evolved through the accumulation of mutations within this structural unit, as well as in much more variable regions in the enzyme that extend from the minimal unit.

Taujale et al. then built an extended family tree for this collection of glycosyltransferases and details of the evolutionary relationships between the enzymes helped them to create a machine learning framework that could predict which sugar-containing molecules were the raw materials for a given glycosyltransferase. This framework could make predictions with nearly 90% accuracy based only on information that can be deciphered from the gene for that enzyme.

These findings will provide scientists with new hypotheses for investigating the complex relationships connecting the genetic information about glycosyltransferases with their structures and activities. Further refinement of the machine learning framework may eventually enable the design of enzymes with properties that are desirable for applications in biotechnology.

Metadynamics to Enhance Sampling in Biomolecular Simulations

Molecular dynamics is a powerful simulation method to provide detailed atomic-scale insight into a range of biological processes including protein folding, biochemical reactions, ligand binding, and many others. Over the last several decades, enhanced sampling methods have been developed to address the large separation in time scales between a molecular dynamics simulation (usually microseconds or shorter) and the time scales of biological processes (often orders of magnitude longer). This chapter specifically focuses on the metadynamics family of methods, which achieves enhanced sampling through the introduction of a history-dependent bias potential that is based on one or more slow degrees of freedom, called collective variables. We introduce the method and its recent variants related to biomolecular studies and then discuss frontier areas of the method. A large part of this chapter is devoted to helping new users of the method understand how to choose metadynamics parameters properly and apply the method to their system of interest.

The Mycobacterium tuberculosis capsule: a cell structure with key implications in pathogenesis

Bacterial capsules have evolved to be at the forefront of the cell envelope, making them an essential element of bacterial biology. Efforts to understand the Mycobacterium tuberculosis (Mtb) capsule began more than 60 years ago, but the relatively recent development of mycobacterial genetics combined with improved chemical and immunological tools have revealed a more refined view of capsule molecular composition. A glycogen-like α-glucan is the major constituent of the capsule, with lower amounts of arabinomannan and mannan, proteins and lipids. The major Mtb capsular components mediate interactions with phagocytes that favor bacterial survival. Vaccination approaches targeting the mycobacterial capsule have proven successful in controlling bacterial replication. Although the Mtb capsule is composed of polysaccharides of relatively low complexity, the concept of antigenic variability associated with this structure has been suggested by some studies. Understanding how Mtb shapes its envelope during its life cycle is key to developing anti-infective strategies targeting this structure at the host-pathogen interface.

Quick-soaking of crystals reveals unprecedented insights into the catalytic mechanism of glycosyltransferases

Glycosyltransferases (GTs) catalyze the transfer of a sugar moiety from nucleotide-sugar or lipid-phospho-sugar donors to a wide range of acceptor substrates, generating a remarkable amount of structural diversity in biological systems. Glycosyl transfer reactions can proceed with either inversion or retention of the anomeric configuration with respect to the sugar donor substrate. In this chapter, we discuss the application of a quick soaking method of substrates and products into protein crystals to visualize native ternary complexes of retaining glycosyltransferases. The crystal structures provide different snapshots of the catalytic cycle, including the Michaelis complex. During this sequence of events, we visualize how the enzyme guides the substrates into the reaction center where the glycosyl transfer reaction takes place, and unveil the mechanism of product release, involving multiple conformational changes not only in the substrates and products but also in the enzyme. The methodology described here provides unprecedented insights into the catalytic mechanism of glycosyltransferases at the molecular level, and can be applied to the study a myriad of enzymatic mediated reactions.