Hybrid Quantum Mechanical/Molecular Mechanical Investigation of the β-1,4-Galactosyltransferase-I Mechanism

Utilization of glycosyltransferases as a seamless tool for synthesis and modification of the oligosaccharides-A review

Glycosyltransferases (GTs) catalyze the transfer of active monosaccharide donors to carbohydrates to create a wide range of oligosaccharide structures. GTs display strong regioselectivity and stereoselectivity in producing glycosidic bonds, making them extremely valuable in the in vitro synthesis of oligosaccharides. The synthesis of oligosaccharides by GTs often gives high yields; however, the enzyme activity may experience product inhibition. Additionally, the higher cost of nucleotide sugars limits the usage of GTs for oligosaccharide synthesis. In this review, we comprehensively discussed the structure and mechanism of GTs based on recent literature and the CAZY website data. To provide innovative ideas for the functional studies of GTs, we summarized several remarkable characteristics of GTs, including folding, substrate specificity, regioselectivity, donor sugar nucleotides, catalytic reversibility, and differences between GTs and GHs. In particular, we highlighted the recent advancements in multi-enzyme cascade reactions and co-immobilization of GTs, focusing on overcoming problems with product inhibition and cost issues. Finally, we presented various types of GT that have been successfully used for oligosaccharide synthesis. We concluded that there is still an opportunity for improvement in enzymatically produced oligosaccharide yield, and future research should focus on improving the yield and reducing the production cost.

UDP-Glycosyltransferases in Edible Fungi: Function, Structure, and Catalytic Mechanism

UDP-glycosyltransferases (UGTs) are the most studied glycosyltransferases, and belong to large GT1 family performing the key roles in antibiotic synthesis, the development of bacterial glycosyltransferase inhibitors, and in animal inflammation. They transfer the glycosyl groups from nucleotide UDP-sugars (UDP-glucose, UDP-galactose, UDP-xylose, and UDP-rhamnose) to the acceptors including saccharides, proteins, lipids, and secondary metabolites. The present review summarized the recent of UDP-glycosyltransferases, including their structures, functions, and catalytic mechanism, especially in edible fungi. The future perspectives and new challenges were also summarized to understand of their structure–function relationships in the future. The outputs in this field could provide a reference to recognize function, structure, and catalytic mechanism of UDP-glycosyltransferases for understanding the biosynthesis pathways of secondary metabolites, such as hydrocarbons, monoterpenes, sesquiterpene, and polysaccharides in edible fungi.

Multifaceted Computational Modeling in Glycoscience

Glycoscience assembles all the scientific disciplines involved in studying various molecules and macromolecules containing carbohydrates and complex glycans. Such an ensemble involves one of the most extensive sets of molecules in quantity and occurrence since they occur in all microorganisms and higher organisms. Once the compositions and sequences of these molecules are established, the determination of their three-dimensional structural and dynamical features is a step toward understanding the molecular basis underlying their properties and functions. The range of the relevant computational methods capable of addressing such issues is anchored by the specificity of stereoelectronic effects from quantum chemistry to mesoscale modeling throughout molecular dynamics and mechanics and coarse-grained and docking calculations. The Review leads the reader through the detailed presentations of the applications of computational modeling. The illustrations cover carbohydrate-carbohydrate interactions, glycolipids, and N- and O-linked glycans, emphasizing their role in SARS-CoV-2. The presentation continues with the structure of polysaccharides in solution and solid-state and lipopolysaccharides in membranes. The full range of protein-carbohydrate interactions is presented, as exemplified by carbohydrate-active enzymes, transporters, lectins, antibodies, and glycosaminoglycan binding proteins. A final section features a list of 150 tools and databases to help address the many issues of structural glycobioinformatics.

Assays for Evaluation of Substrates for and Inhibitors of β-1,4-Galactosyltransferase 7

β-1,4-Galactosyltransferase 7 (β4GalT7) is a key enzyme in the synthesis of two classes of glycosaminoglycans (GAG), i.e., heparan sulfate (HS) and chondroitin/dermatan sulfate (CS/DS). GAG chains are linear polysaccharides of alternating hexuronic acid and N-acetylhexosamine residues, commonly linked to core proteins to form proteoglycans with important roles in the regulation of a range of biological processes. The biosynthesis of GAGs is initiated by xylosylation of a serine residue of the core protein followed by galactosylation, catalyzed by β4GalT7. The biosynthesis can also be initiated by xylosides carrying hydrophobic aglycons, such as 2-naphthyl β-D-xylopyranoside. We have cloned and expressed β4GalT7, and designed a cell-free assay to measure the activity of this enzyme. The assay employs a 96-well plate format for high throughput. In this chapter, we describe the cloning, expression, and purification of β4GalT7, as well as assays proposed for development of substrates for GAG priming and for investigating inhibitors of β4GalT7.

Ab Initio Molecular Dynamics Simulations of the SN1/SN2 Mechanistic Continuum in Glycosylation Reactions

We report a computational approach to evaluate the reaction mechanisms of glycosylation using ab initio molecular dynamics (AIMD) simulations in explicit solvent. The reaction pathways are simulated via free energy calculations based on metadynamics and trajectory simulations using Born-Oppenheimer molecular dynamics. We applied this approach to investigate the mechanisms of the glycosylation of glucosyl α-trichloroacetimidate with three acceptors (EtOH, i-PrOH, and t-BuOH) in three solvents (ACN, DCM, and MTBE). The reactants and the solvents are treated explicitly using density functional theory. We show that the profile of the free energy surface, the synchronicity of the transition state structure, and the time gap between leaving group dissociation and nucleophile association can be used as three complementary indicators to describe the glycosylation mechanism within the SN1/SN2 continuum for a given reaction. This approach provides a reliable means to rationalize and predict reaction mechanisms and to estimate lifetimes of oxocarbenium intermediates and their dependence on the glycosyl donor, acceptor, and solvent environment.

How inverting β-1,4-galactosyltransferase-1 can quench a high charge of the by-product UDP3- in catalysis: a QM/MM study of enzymatic reaction with native and UDP-5' thio galactose substrates

The catalysis of inverting glycosyltransferases consists of several biophysical and biochemical processes during which the transfer of a sugar residue from the purine phosphate donor substrate to an acceptor substrate occurs with stereo-inversion of the anomeric C1 center at a product. During catalysis a highly charged phosphate by-product (UDP3-) is formed and a mechanism of how the enzyme stabilizes it back to the UDP2- form is not known. Using methods of molecular modeling (hybrid DFT-QM/MM calculations) we proposed and validated a catalytic mechanism of bovine inverting β-1,4-galactosyltransferase-1 (β4Gal-T1) with native (UDP-galactose) and thio donor substrates (UDP-5' thio galactose). We focused on three aspects of the mechanism not yet investigated: (i) the formation of an oxocarbenium ion intermediate, which was only found for the retaining glycosyltransferases for the time being; (ii) the mechanism of stabilization of a highly charged phosphate by-product (UDP3-) back to its standard in vivo form (UDP2-); (iii) explanation for why in experimental measurements the rate of catalysis with the thio donor substrate is only 8% of the rate of that with the natural substrate. To understand the differences in the interaction patterns between the complexes enzyme : UDP-Gal and enzyme : UDP-5S-Gal, fragmented molecular orbital (FMO) decomposition energy analysis was carried out at the DFT level.

Computerized Molecular Modeling of Carbohydrates

Computerized molecular modeling continues to increase in capability and applicability to carbohydrates. This chapter covers nomenclature and conformational aspects of carbohydrates, perhaps of greater use to computational chemists who do not have a strong background in carbohydrates, and its comments on various methods and studies might be of more use to carbohydrate chemists who are inexperienced with computation. Work on the intrinsic variability of glucose, an overall theme, is described. Other areas of the authors' emphasis, including evaluation of hydrogen bonding by the atoms-in-molecules approach, and validation of modeling methods with crystallographic results are also presented.

Selectins-The Two Dr. Jekyll and Mr. Hyde Faces of Adhesion Molecules-A Review

Selectins belong to a group of adhesion molecules that fulfill an essential role in immune and inflammatory responses and tissue healing. Selectins are glycoproteins that decode the information carried by glycan structures, and non-covalent interactions of selectins with these glycan structures mediate biological processes. The sialylated and fucosylated tetrasaccharide sLe^x is an essential glycan recognized by selectins. Several glycosyltransferases are responsible for the biosynthesis of the sLe^x tetrasaccharide. Selectins are involved in a sequence of interactions of circulated leukocytes with endothelial cells in the blood called the adhesion cascade. Recently, it has become evident that cancer cells utilize a similar adhesion cascade to promote metastases. However, like Dr. Jekyll and Mr. Hyde’s two faces, selectins also contribute to tissue destruction during some infections and inflammatory diseases. The most prominent function of selectins is associated with the initial stage of the leukocyte adhesion cascade, in which selectin binding enables tethering and rolling. The first adhesive event occurs through specific non-covalent interactions between selectins and their ligands, with glycans functioning as an interface between leukocytes or cancer cells and the endothelium. Targeting these interactions remains a principal strategy aimed at developing new therapies for the treatment of immune and inflammatory disorders and cancer. In this review, we will survey the significant contributions to and the current status of the understanding of the structure of selectins and the role of selectins in various biological processes. The potential of selectins and their ligands as therapeutic targets in chronic and acute inflammatory diseases and cancer will also be discussed. We will emphasize the structural characteristic of selectins and the catalytic mechanisms of glycosyltransferases involved in the biosynthesis of glycan recognition determinants. Furthermore, recent achievements in the synthesis of selectin inhibitors will be reviewed with a focus on the various strategies used for the development of glycosyltransferase inhibitors, including substrate analog inhibitors and transition state analog inhibitors, which are based on knowledge of the catalytic mechanism.

Theoretical Insights into the Reaction and Inhibition Mechanism of Metal-Independent Retaining Glycosyltransferase Responsible for Mycothiol Biosynthesis

Understanding enzymatic reactions with atomic resolution has proven in recent years to be of tremendous interest for biochemical research, and thus, the use of QM/MM methods for the study of reaction mechanisms is experiencing a continuous growth. Glycosyltransferases (GTs) catalyze the formation of glycosidic bonds, and are important for many biotechnological purposes, including drug targeting. Their reaction product may result with only one of the two possible stereochemical outcomes for the reacting anomeric center, and therefore, they are classified as either inverting or retaining GTs. While the inverting GT reaction mechanism has been widely studied, the retaining GT mechanism has always been controversial and several questions remain open to this day. In this work, we take advantage of our recent GPU implementation of a pure QM(DFT-PBE)/MM approach to explore the reaction and inhibition mechanism of MshA, a key retaining GT responsible for the first step of mycothiol biosynthesis, a low weight thiol compound found in pathogens like Mycobacterium tuberculosis that is essential for its survival under oxidative stress conditions. Our results show that the reaction proceeds via a front-side S_Ni-like concerted reaction mechanism (D_NA_N in IUPAC nomenclature) and has a 17.5 kcal/mol free energy barrier, which is in remarkable agreement with experimental data. Detailed analysis shows that the key reaction step is the diphosphate leaving group dissociation, leading to an oxocarbenium-ion-like transition state. In contrast, fluorinated substrate analogues increase the reaction barrier significantly, rendering the enzyme effectively inactive. Detailed analysis of the electronic structure along the reaction suggests that this particular inhibition mechanism is associated with fluorine's high electronegative nature, which hinders phosphate release and proper stabilization of the transition state.

Structure and Mechanism of Staphylococcus aureus TarS, the Wall Teichoic Acid β-glycosyltransferase Involved in Methicillin Resistance

In recent years, there has been a growing interest in teichoic acids as targets for antibiotic drug design against major clinical pathogens such as Staphylococcus aureus, reflecting the disquieting increase in antibiotic resistance and the historical success of bacterial cell wall components as drug targets. It is now becoming clear that β-O-GlcNAcylation of S. aureus wall teichoic acids plays a major role in both pathogenicity and antibiotic resistance. Here we present the first structure of S. aureus TarS, the enzyme responsible for polyribitol phosphate β-O-GlcNAcylation. Using a divide and conquer strategy, we obtained crystal structures of various TarS constructs, mapping high resolution overlapping N-terminal and C-terminal structures onto a lower resolution full-length structure that resulted in a high resolution view of the entire enzyme. Using the N-terminal structure that encapsulates the catalytic domain, we furthermore captured several snapshots of TarS, including the native structure, the UDP-GlcNAc donor complex, and the UDP product complex. These structures along with structure-guided mutants allowed us to elucidate various catalytic features and identify key active site residues and catalytic loop rearrangements that provide a valuable platform for anti-MRSA drug design. We furthermore observed for the first time the presence of a trimerization domain composed of stacked carbohydrate binding modules, commonly observed in starch active enzymes, but adapted here for a poly sugar-phosphate glycosyltransferase.

Historically, β-lactam class antibiotics such as methicillin have been very successful in the treatment of bacterial infections, effectively destroying bacteria by rupturing their cell walls while posing little harm to the human organism. In recent years, however, the alarming emergence of Methicillin Resistant S. aureus or MRSA has resulted in a world-wide health crisis, calling on new strategies to combat pathogenesis and antibiotic resistance. As such, understanding the pathways and players that orchestrate resistance is important for overcoming these mechanisms and restoring our powerful β-lactam antibiotic arsenal. In this article we describe the crystal structure of TarS, an enzyme responsible for the glycosylation of wall teichoic acid polymers of the S. aureus cell wall, a process that has been shown to be specifically responsible for methicillin resistance in MRSA. TarS is therefore a promising drug target whose inhibition in combinational therapies would result in MRSA re-sensitization to β-lactam antibiotics. Here we present the first structure of TarS together with several snap-shots of its substrate/product complexes, and elucidate important catalytic features that are valuable for rational drug design efforts to combat resistance in MRSA.

Who's on base? Revealing the catalytic mechanism of inverting family 6 glycoside hydrolases

In several important classes of inverting carbohydrate-active enzymes, the identity of the catalytic base remains elusive, including in family 6 Glycoside Hydrolase (GH6) enzymes, which are key components of cellulase cocktails for cellulose depolymerization. Despite many structural and kinetic studies with both wild-type and mutant enzymes, especially on the Trichoderma reesei (Hypocrea jecorina) GH6 cellulase (TrCel6A), the catalytic base in the single displacement inverting mechanism has not been definitively identified in the GH6 family. Here, we employ transition path sampling to gain insight into the catalytic mechanism, which provides unbiased atomic-level understanding of key order parameters involved in cleaving the strong glycosidic bond. Our hybrid quantum mechanics and molecular mechanics (QM/MM) simulations reveal a network of hydrogen bonding that aligns two active site water molecules that play key roles in hydrolysis: one water molecule drives the reaction by nucleophilic attack on the substrate and a second shuttles a proton to the putative base (D175) via a short water wire. We also investigated the case where the putative base is mutated to an alanine, an enzyme that is experimentally still partially active. The simulations predict that proton hopping along a water wire via a Grotthuss mechanism provides a mechanism of catalytic rescue. Further simulations reveal that substrate processive motion is 'driven' by strong electrostatic interactions with the protein at the product sites and that the -1 sugar adopts a 2SO ring configuration as it reaches its binding site. This work thus elucidates previously elusive steps in the processive catalytic mechanism of this important class of enzymes.

Simulations of cellulose translocation in the bacterial cellulose synthase suggest a regulatory mechanism for the dimeric structure of cellulose

In addition to suggesting a mechanism for regulating cellulose structure, molecular simulations indicate translocation is not rate-limiting for cellulose biosynthesis.

The processive cycle of the bacterial cellulose synthase (Bcs) includes the addition of a single glucose moiety to the end of a growing cellulose chain followed by the translocation of the nascent chain across the plasma membrane. The mechanism of this translocation and its precise location within the processive cycle are not well understood. In particular, the molecular details of how a polymer (cellulose) whose basic structural unit is a dimer (cellobiose) can be constructed by adding one monomer (glucose) at a time are yet to be elucidated. Here, we have utilized molecular dynamics simulations and free energy calculations to the shed light on these questions. We find that translocation forward by one glucose unit is quite favorable energetically, giving a free energy stabilization of greater than 10 kcal mol^–1. In addition, there is only a small barrier to translocation, implying that translocation is not rate limiting within the Bcs processive cycle (given experimental rates for cellulose synthesis in vitro). Perhaps most significantly, our results also indicate that steric constraints at the transmembrane tunnel entrance regulate the dimeric structure of cellulose. Namely, when a glucose molecule is added to the cellulose chain in the same orientation as the acceptor glucose, the terminal glucose freely rotates upon forward motion, thus suggesting a regulatory mechanism for the dimeric structure of cellulose. We characterize both the conserved and non-conserved enzyme–polysaccharide interactions that drive translocation, and find that 20 of the 25 residues that strongly interact with the translocating cellulose chain in the simulations are well conserved, mostly with polar or aromatic side chains. Our results also allow for a dynamical analysis of the role of the so-called ‘finger helix’ in cellulose translocation that has been observed structurally. Taken together, these findings aid in the elucidation of the translocation steps of the Bcs processive cycle and may be widely relevant to polysaccharide synthesizing or degrading enzymes that couple catalysis with chain translocation.

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.

How Cellulose Elongates--A QM/MM Study of the Molecular Mechanism of Cellulose Polymerization in Bacterial CESA

The catalytic mechanism of bacterial cellulose synthase was investigated by using a hybrid quantum mechanics and molecular mechanics (QM/MM) approach. The Michaelis complex model was built based on the X-ray crystal structure of the cellulose synthase subunits BcsA and BcsB containing a uridine diphosphate molecule and a translocating glucan. Our study identified an SN2-type transition structure corresponding to the nucleophilic attack of the nonreducing end O4 on the anomeric carbon C1, the breaking of the glycosidic bond C1-O1, and the transfer of proton from the nonreducing end O4 to the general base D343. The activation barrier found for this SN2-type transition state is 68 kJ/mol. The rate constant of polymerization is estimated to be ∼8.0 s(-1) via transition state theory. A similar SN2-type transition structure was also identified for a second glucose molecule added to the growing polysaccharide chain, which aligned with the polymer 180° rotated compared to the initially added unit. This study provides detailed insights into how cellulose is extended by one glucose molecule at a time and how the individual glucose units align into cellobiose repeating units.

A molecular description of cellulose biosynthesis

Cellulose is the most abundant biopolymer on Earth, and certain organisms from bacteria to plants and animals synthesize cellulose as an extracellular polymer for various biological functions. Humans have used cellulose for millennia as a material and an energy source, and the advent of a lignocellulosic fuel industry will elevate it to the primary carbon source for the burgeoning renewable energy sector. Despite the biological and societal importance of cellulose, the molecular mechanism by which it is synthesized is now only beginning to emerge. On the basis of recent advances in structural and molecular biology on bacterial cellulose synthases, we review emerging concepts of how the enzymes polymerize glucose molecules, how the nascent polymer is transported across the plasma membrane, and how bacterial cellulose biosynthesis is regulated during biofilm formation. Additionally, we review evolutionary commonalities and differences between cellulose synthases that modulate the nature of the cellulose product formed.

Glycosyltransferases: mechanisms and applications in natural product development

Glycosylation reactions mainly catalyzed by glycosyltransferases (Gts) occur almost everywhere in the biosphere, and always play crucial roles in vital processes.

Rules for priming and inhibition of glycosaminoglycan biosynthesis; probing the β4GalT7 active site

Xylose is the optimal substrate for β4GalT7, an essential enzyme in GAG biosynthesis, but analogs act as effective inhibitors.

Crystal structures of β-1,4-galactosyltransferase 7 enzyme reveal conformational changes and substrate binding

The β-1,4-galactosyltransferase 7 (β4GalT7) enzyme is involved in proteoglycan synthesis. In the presence of a manganese ion, it transfers galactose from UDP-galactose to xylose on a proteoglycan acceptor substrate. We present here the crystal structures of human β4GalT7 in open and closed conformations. A comparison of these crystal structures shows that, upon manganese and UDP or UDP-Gal binding, the enzyme undergoes conformational changes involving a small and a long loop. We also present the crystal structures of Drosophila wild-type β4GalT7 and D211N β4GalT7 mutant enzymes in the closed conformation in the presence of the acceptor substrate xylobiose and the donor substrate UDP-Gal, respectively. To understand the catalytic mechanism, we have crystallized the ternary complex of D211N β4GalT7 mutant enzyme in the presence of manganese with the donor and the acceptor substrates together in the same crystal structure. The galactose moiety of the bound UDP-Gal molecule forms seven hydrogen bonds with the protein molecule. The nonreducing end of the xylose moiety of xylobiose binds to the hydrophobic acceptor sugar binding pocket created by the conformational changes, whereas its extended xylose moiety forms hydrophobic interactions with a Tyr residue. In the ternary complex crystal structure, the nucleophile O4 oxygen atom of the xylose molecule is found in close proximity to the C1 and O5 atoms of the galactose moiety. This is the first time that a Michaelis complex of a glycosyltransferase has been described, and it clearly suggests an SN2 type catalytic mechanism for the β4GalT7 enzyme.

Selecting DFT methods for use in optimizations of enzyme active sites: applications to ONIOM treatments of DNA glycosylases

When using a hybrid methodology to treat an enzymatic reaction, many factors contribute to selecting the method for the high-level region, which can be complicated by the presence of dispersion-driven interactions such as π–π stacking. In addition, the proper treatment of the reaction center often requires a large number of heavy atoms to be included in the high-level region, precluding the use of ab initio methods such as MP2 as well as large basis sets, in the optimization step. In the present work, popular DFT methods were tested to identify an appropriate functional for treating the high-level region in ONIOM optimizations of reactions catalyzed by nonmetalloenzymes. Eight different DFT methods (B3LYP, B97-2, MPW1K, MPWB1K, BB1K, B1B95, M06-2X, and ωB97X-D) in combination with four double-ζ quality Pople basis sets were tested for their ability to optimize noncovalent interactions (hydrogen bonding and π–π) and characterize reactions (proton transfer, SN2 hydrolysis, and unimolecular cleavage). Although the primary focus of this study is accurate structure determination, energetics were also examined at both the optimization level of theory, and with triple-ζ quality basis set and select (M06-2X or ωB97X-D) methods. If dispersion-driven interactions exist within the active site, then MPWB1K/6-31G(d,p) or M06-2X/6-31+G(d,p) are recommended for the optimization step with subsequent triple-ζ quality single-point energies. However, since dispersion-corrected functionals (M06-2X and ωB97X-D) generally require diffuse functions to yield appropriate geometries, the possible size of the high-level region is greatly limited with these methods. In contrast, if the model is large enough to recover steric constraints on π–π interactions, then B3LYP with a small basis set performs comparatively well for the optimization step and is significantly less computationally expensive. Interestingly, the functionals that afford the best geometries often do not yield the best energetics, which emphasizes the importance of structural benchmark studies.

A QM/MM investigation of the catalytic mechanism of metal-ion-independent core 2 β1,6-N-acetylglucosaminyltransferase

β1,6-GlcNAc-transferase (C2GnT) is an important controlling factor of biological functions for many glycoproteins and its activity has been found to be altered in breast, colon, and lung cancer cells, in leukemia cells, in the lymhomonocytes of multiple sclerosis patients, leukocytes from diabetes patients, and in conditions causing an immune deficiency. The result of the action of C2GnT is the core 2 structure that is essential for the further elongation of the carbohydrate chains of O-glycans. The catalytic mechanism of this metal-ion-independent glycosyltransferase is of paramount importance and is investigated here by using quantum mechanical (QM) (density functional theory (DFT))/molecular modeling (MM) methods with different levels of theory. The structural model of the reaction site used in this report is based on the crystal structures of C2GnT. The entire enzyme-substrate system was subdivided into two different subsystems: the QM subsystem containing 206 atoms and the MM region containing 5914 atoms. Three predefined reaction coordinates were employed to investigate the catalytic mechanism. The calculated potential energy surfaces discovered the existence of a concerted SN 2-like mechanism. In this mechanism, a nucleophilic attack by O6 facilitated by proton transfer to the catalytic base and the separation of the leaving group all occur almost simultaneously. The transition state for the proposed reaction mechanism at the M06-2X/6-31G** (with diffuse functions on the O1', O5', OGlu , and O6 atoms) level was located at C1-O6=1.74 Å and C1-O1=2.86 Å. The activation energy for this mechanism was estimated to be between 20 and 29 kcal mol⁻¹, depending on the method used. These calculations also identified a low-barrier hydrogen bond between the nucleophile O6H and the catalytic base Glu320, and a hydrogen bond between the N-acetamino group and the glycosidic oxygen of the donor in the TS. It is proposed that these interactions contribute to a stabilization of TS and participate in the catalytic mechanism.

“Lost sugars” — reality of their biological and medical applications

The glycan chains attached to cell surfaces or to single proteins are highly dynamic structures with various functions. The glycan chains of mammals and of some microorganisms often terminate in sialic acids or α-1,3-galactose. Although these two sugars are completely distinct, there are several similarities in their biological and medical importance. First, one type of sialic acid, N-glycolylneuraminic acid, and the galactose bound by an α-1,3-linkage to LacNAc, that forms an α-gal epitope, were both eliminated in human evolution, resulting in the production of antibodies to these sugars. Both of these evolutionary events have consequences connected with the consumption of foods of mammalian origin, causing medical complications of varying severity. In terms of ageing, sialic acids prevent the clearance of glycoproteins and circulating blood cells, whereas cryptic α-gal epitopes on senescent red blood cells contribute to their removal from circulation. The efficiency of therapeutic proteins can be increased by sialylation. Another common feature is the connection with microorganisms since sialic acids and α-gal epitopes serve as receptors on host cells and can also be expressed on the surfaces of some microorganisms. Whereas, the sialylation of IgG antibodies may help to treat inflammation, the expression of the α-gal epitope on microbial antigens increases the immunogenicity of the corresponding vaccines. Finally, sialic acids and the α-gal epitope have applications in cancer immunotherapy. N-glycolylneuraminic acid is a powerful target for cancer immunotherapy, and the α-gal epitope increases the efficiency of cancer vaccines. The final section of this article contains a brief overview of the methods for oligosaccharide chain synthesis and the characteristics of sialyltransferases and α-1,3-galactosyltransferase.

Substrate-assisted catalytic mechanism of O-GlcNAc transferase discovered by quantum mechanics/molecular mechanics investigation

In higher eukaryotes, a variety of proteins are post-translationally modified by adding O-linked N-acetylglucosamine (GlcNAc) residue to serine or threonine residues. Misregulation of O-GlcNAcylation is linked to a wide variety of diseases, such as diabetes, cancer, and neurodegenerative diseases, including Alzheimer's disease. GlcNAc transfer is catalyzed by an inverting glycosyltransferase O-GlcNAc transferase (uridine diphospho-N-acetylglucosamine:polypeptide β-N-acetylaminyltransferase, OGT) that belongs to the GT-B superfamily. The catalytic mechanism of this metal-independent glycosyltransferase is of primary importance and is investigated here using QM(DFT)/MM methods. The structural model of the reaction site used in this paper is based on the crystal structures of OGT. The entire enzyme-substrate system was partitioned into two different subsystems: the QM subsystem containing 198 atoms, and the MM region containing 11,326 atoms. The catalytic mechanism was monitored by means of three two-dimensional potential energy maps calculated as a function of three predefined reaction coordinates at different levels of theory. These potential energy surfaces revealed the existence of a concerted S(N)2-like mechanism, in which a nucleophilic attack by O(Ser), facilitated by proton transfer to the catalytic base, and the dissociation of the leaving group occur almost simultaneously. The transition state for the proposed reaction mechanism at the MPW1K level was located at C1-O(Ser) = 1.92 Å and C1-O1 = 3.11 Å. The activation energy for this passage was estimated to be ~20 kcal mol(-1). These calculations also identified, for the first time for glycosyltransferases, the substrate-assisted mechanism in which the N-acetamino group of the donor participates in the catalytic mechanism.

Computerized molecular modeling of carbohydrates

Computerized molecular modeling continues to increase in capability and applicability to carbohydrates. This chapter covers nomenclature and conformational aspects of carbohydrates, perhaps of greater use to carbohydrate-inexperienced computational chemists. Its comments on various methods and studies might be of more use to computation-inexperienced carbohydrate chemists. New work on intrinsic variability of glucose, an overall theme, is described.

Current trends in the structure-activity relationships of sialyltransferases

Sialyltransferases (STs) represent an important group of enzymes that transfer N-acetylneuraminic acid (Neu5Ac) from cytidine monophosphate-Neu5Ac to various acceptor substrates. In higher animals, sialylated oligosaccharide structures play crucial roles in many biological processes but also in diseases, notably in microbial infection and cancer. Cell surface sialic acids have also been found in a few microorganisms, mainly pathogenic bacteria, and their presence is often associated with virulence. STs are distributed into five different families in the CAZy database (http://www.cazy.org/). On the basis of crystallographic data available for three ST families and fold recognition analysis for the two other families, STs can be grouped into two structural superfamilies that represent variations of the canonical glycosyltransferase (GT-A and GT-B) folds. These two superfamilies differ in the nature of their active site residues, notably the catalytic base (a histidine or an aspartate residue). The observed structural and functional differences strongly suggest that these two structural superfamilies have evolved independently.

Bioinformatics and molecular modeling in glycobiology

The field of glycobiology is concerned with the study of the structure, properties, and biological functions of the family of biomolecules called carbohydrates. Bioinformatics for glycobiology is a particularly challenging field, because carbohydrates exhibit a high structural diversity and their chains are often branched. Significant improvements in experimental analytical methods over recent years have led to a tremendous increase in the amount of carbohydrate structure data generated. Consequently, the availability of databases and tools to store, retrieve and analyze these data in an efficient way is of fundamental importance to progress in glycobiology. In this review, the various graphical representations and sequence formats of carbohydrates are introduced, and an overview of newly developed databases, the latest developments in sequence alignment and data mining, and tools to support experimental glycan analysis are presented. Finally, the field of structural glycoinformatics and molecular modeling of carbohydrates, glycoproteins, and protein-carbohydrate interaction are reviewed.

Crystal structure of the catalytic domain of Drosophila beta1,4-Galactosyltransferase-7

The beta1,4-galactosyltransferase-7 (beta4Gal-T7) enzyme, one of seven members of the beta4Gal-T family, transfers in the presence of manganese Gal from UDP-Gal to an acceptor sugar (xylose) that is attached to a side chain hydroxyl group of Ser/Thr residues of proteoglycan proteins. It exhibits the least protein sequence similarity with the other family members, including the well studied family member beta4Gal-T1, which, in the presence of manganese, transfers Gal from UDP-Gal to GlcNAc. We report here the crystal structure of the catalytic domain of beta4Gal-T7 from Drosophila in the presence of manganese and UDP at 1.81 A resolution. In the crystal structure, a new manganese ion-binding motif (HXH) has been observed. Superposition of the crystal structures of beta4Gal-T7 and beta4Gal-T1 shows that the catalytic pocket and the substrate-binding sites in these proteins are similar. Compared with GlcNAc, xylose has a hydroxyl group (instead of an N-acetyl group) at C2 and lacks the CH(2)OH group at C5; thus, these protein structures show significant differences in their acceptor-binding site. Modeling of xylose in the acceptor-binding site of the beta4Gal-T7 crystal structure shows that the aromatic side chain of Tyr(177) interacts strongly with the C5 atom of xylose, causing steric hindrance to any additional group at C5. Because Drosophila Cd7 has a 73% protein sequence similarity to human Cd7, the present crystal structure offers a structure-based explanation for the mutations in human Cd7 that have been linked to Ehlers-Danlos syndrome.