Atomistic insight into the catalytic mechanism of glycosyltransferases by combined quantum mechanics/molecular mechanics (QM/MM) methods

Specificity determinants revealed by the structure of glycosyltransferase Campylobacter concisus PglA

In selected Campylobacter species, the biosynthesis of N-linked glycoconjugates via the pgl pathway is essential for pathogenicity and survival. However, most of the membrane-associated GT-B fold glycosyltransferases responsible for diversifying glycans in this pathway have not been structurally characterized which hinders the understanding of the structural factors that govern substrate specificity and prediction of resulting glycan composition. Herein, we report the 1.8 Å resolution structure of Campylobacter concisus PglA, the glycosyltransferase responsible for the transfer of N-acetylgalatosamine (GalNAc) from uridine 5'-diphospho-N-acetylgalactosamine (UDP-GalNAc) to undecaprenyl-diphospho-N,N'-diacetylbacillosamine (UndPP-diNAcBac) in complex with the sugar donor GalNAc. This study identifies distinguishing characteristics that set PglA apart within the GT4 enzyme family. Computational docking of the structure in the membrane in comparison to homologs points to differences in interactions with the membrane-embedded acceptor and the structural analysis of the complex together with bioinformatics and site-directed mutagenesis identifies donor sugar binding motifs. Notably, E113, conserved solely among PglA enzymes, forms a hydrogen bond with the GalNAc C6″-OH. Mutagenesis of E113 reveals activity consistent with this role in substrate binding, rather than stabilization of the oxocarbenium ion transition state, a function sometimes ascribed to the corresponding residue in GT4 homologs. The bioinformatic analyses reveal a substrate-specificity motif, showing that Pro281 in a substrate binding loop of PglA directs configurational preference for GalNAc over GlcNAc. This proline is replaced by a conformationally flexible glycine, even in distant homologs, which favor substrates with the same stereochemistry at C4, such as glucose. The signature loop is conserved across all Campylobacter PglA enzymes, emphasizing its importance in substrate specificity.

Complete genome analysis of pathogenic Metschnikowia bicuspidata strain MQ2101 isolated from diseased ridgetail white prawn, Exopalaemon carinicauda

BACKGROUND

Metschnikowia bicuspidata is a pathogenic yesst that can cause disease in many different economic aquatic animal species. In recent years, there was a new disease outbreak in ridgetail white prawn (Exopalaemon carinicauda) in coastal areas of Jiangsu Province China that was referred to as zombie disease by local farmers. The pathogen was first isolated and identified as M. bicuspidata. Although the pathogenicity and pathogenesis of this pathogen in other animals have been reported in some previous studies, research on its molecular mechanisms is still very limited. Therefore, a genome-wide study is necessary to better understand the physiological and pathogenic mechanisms of M. bicuspidata.

RESULT

In this study, we obtained a pathogenic strain, MQ2101, of M. bicuspidata from diseased E. carinicauda and sequenced its whole genome. The size of the whole genome was 15.98 Mb, and it was assembled into 5 scaffolds. The genome contained 3934 coding genes, among which 3899 genes with biological functions were annotated in multiple underlying databases. In KOG database, 2627 genes were annotated, which were categorized into 25 classes including general function prediction only, posttranslational modification, protein turnover, chaperones, and signal transduction mechanisms. In KEGG database, 2493 genes were annotated, which were categorized into five classes, including cellular processes, environmental information processing, genetic information processing, metabolism and organismal systems. In GO database, 2893 genes were annotated, which were mainly classified in cell, cell part, cellular processes and metabolic processes. There were 1055 genes annotated in the PHI database, accounting for 26.81% of the total genome, among which 5 genes were directly related to pathogenicity (identity ≥ 50%), including hsp90, PacC, and PHO84. There were also some genes related to the activity of the yeast itself that could be targeted by antiyeast drugs. Analysis based on the DFVF database showed that strain MQ2101 contained 235 potential virulence genes. BLAST searches in the CAZy database showed that strain MQ2101 may have a more complex carbohydrate metabolism system than other yeasts of the same family. In addition, two gene clusters and 168 putative secretory proteins were predicted in strain MQ2101, and functional analysis showed that some of the secretory proteins may be directly involved in the pathogenesis of the strain. Gene family analysis with five other yeasts revealed that strain MQ2101 has 245 unique gene families, including 274 genes involved in pathogenicity that could serve as potential targets.

CONCLUSION

Genome-wide analysis elucidated the pathogenicity-associated genes of M. bicuspidate while also revealing a complex metabolic mechanism and providing putative targets of action for the development of antiyeast drugs for this pathogen. The obtained whole-genome sequencing data provide an important theoretical basis for transcriptomic, proteomic and metabolic studies of M. bicuspidata and lay a foundation for defining its specific mechanism of host infestation.

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.

Hybrid Quantum Mechanics/Molecular Mechanics (QM/MM) Simulation: A Tool for Structure-based Drug Design and Discovery

Quantum mechanics (QM) is physics based theory which explains the physical properties of nature at the level of atoms and sub-atoms. Molecular mechanics (MM) construct molecular systems through the use of classical mechanics. So, hybrid quantum mechanics and molecular mechanics (QM/MM) when combined together can act as computer-based methods which can be used to calculate structure and property data of molecular structures. Hybrid QM/MM combines the strengths of QM with accuracy and MM with speed. QM/MM simulation can also be applied for the study of chemical process in solutions as well as in the proteins, and has a great scope in structure-based drug design (CADD) and discovery. Hybrid QM/MM also applied to HTS, to derive QSAR models and due to availability of many protein crystal structures; it has a great role in computational chemistry, especially in structure- and fragment-based drug design. Fused QM/MM simulations have been developed as a widespread method to explore chemical reactions in condensed phases. In QM/MM simulations, the quantum chemistry theory is used to treat the space in which the chemical reactions occur; however the rest is defined through molecular mechanics force field (MMFF). In this review, we have extensively reviewed recent literature pertaining to the use and applications of hybrid QM/MM simulations for ligand and structure-based computational methods for the design and discovery of therapeutic agents.

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.

Three-Dimensional Structures of Carbohydrates and Where to Find Them

Analysis and systematization of accumulated data on carbohydrate structural diversity is a subject of great interest for structural glycobiology. Despite being a challenging task, development of computational methods for efficient treatment and management of spatial (3D) structural features of carbohydrates breaks new ground in modern glycoscience. This review is dedicated to approaches of chemo- and glyco-informatics towards 3D structural data generation, deposition and processing in regard to carbohydrates and their derivatives. Databases, molecular modeling and experimental data validation services, and structure visualization facilities developed for last five years are reviewed.

Enzymatic Synthesis of Glycans and Glycoconjugates

Glycoconjugates have great potential to improve human health in a multitude of different ways and fields. Prominent examples are human milk oligosaccharides and glycosaminoglycans. The typical choice for the production of homogeneous glycoconjugates is enzymatic synthesis. Through the availability of expression and purification protocols, recombinant Leloir glycosyltransferases are widely applied as catalysts for the synthesis of a wide range of glycoconjugates. Extensive utilization of these enzymes also depends on the availability of activated sugars as building blocks. Multi-enzyme cascades have proven a versatile technique to synthesize and in situ regenerate nucleotide sugar.In this chapter, the functions and mechanisms of Leloir glycosyltransferases are revisited, and the advantage of prokaryotic sources and production systems is discussed. Moreover, in vivo and in vitro pathways for the synthesis of nucleotide sugar are reviewed. In the second part, recent and prominent examples of the application of Leloir glycosyltransferase are given, i.e., the synthesis of glycosaminoglycans, glycoconjugate vaccines, and human milk oligosaccharides as well as the re-glycosylation of biopharmaceuticals, and the status of automated glycan assembly is revisited.

Synthesis, docking study and biological evaluation of ᴅ-fructofuranosyl and ᴅ-tagatofuranosyl sulfones as potential inhibitors of the mycobacterial galactan synthesis targeting the galactofuranosyltransferase GlfT2

A series of ten novel ᴅ-fructofuranosyl and ᴅ-tagatofuranosyl sulfones bearing a 1-O-phosphono moiety and three different substituents at C-2 has been prepared. Due to the structural similarities of these scaffolds to the native substrate of mycobacterial galactofuranosyltransferase GlfT2 in the transition state, we evaluated these compounds by computational methods, as well as in an enzyme assay for the possible inhibition of the mycobacterial galactan biosynthesis. Our data show that despite favorable docking scores to the active site of GlfT2, none of these compounds serve as efficient inhibitors of the enzymes involved in the mycobacterial galactan biosynthesis.

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.

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 function of Asp252 in the catalysis of β-1,4-galactosyltransferase-1 may be the stabilization of a high charge of the by-product UDP³⁻by a substrate-assisted proton transfer reaction.

Leloir Glycosyltransferases in Applied Biocatalysis: A Multidisciplinary Approach

Enzymes are nature’s catalyst of choice for the highly selective and efficient coupling of carbohydrates. Enzymatic sugar coupling is a competitive technology for industrial glycosylation reactions, since chemical synthetic routes require extensive use of laborious protection group manipulations and often lack regio- and stereoselectivity. The application of Leloir glycosyltransferases has received considerable attention in recent years and offers excellent control over the reactivity and selectivity of glycosylation reactions with unprotected carbohydrates, paving the way for previously inaccessible synthetic routes. The development of nucleotide recycling cascades has allowed for the efficient production and reuse of nucleotide sugar donors in robust one-pot multi-enzyme glycosylation cascades. In this way, large glycans and glycoconjugates with complex stereochemistry can be constructed. With recent advances, LeLoir glycosyltransferases are close to being applied industrially in multi-enzyme, programmable cascade glycosylations.

UDP-N-acetyl-α-D-galactosamine:polypeptide N-acetylgalactosaminyltransferase from the snail Biomphalaria glabrata - structural reflections

UDP-GalNAc:polypeptide GalNAc transferase (ppGalNAcT; EC 2.4.1.41) is the initiating enzyme for mucin-type O-glycosylation in animals. Members of this highly conserved glycosyltransferase family catalyse a single glycosidic linkage. They transfer an N-acetylgalactosamine (GalNAc) residue from an activated donor (UDP-GalNAc) to a serine or threonine of an acceptor polypeptide chain. A ppGalNAcT from the freshwater snail Biomphalaria glabrata is the only characterised member of this enzyme family from mollusc origin. In this work, we interpret previously published experimental characterization of this enzyme in the context of in silico models of the enzyme and its acceptor substrates. A homology model of the mollusc ppGalNAcT is created and various substrate peptides are modelled into the active site. We hypothesize about possible molecular interpretations of the available experimental data and offer potential explanations for observed substrate and cofactor specificity. Here, we review and synthesise the current knowledge of Bge-ppGalNAcT, supported by a molecular interpretation of the available data.

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.

Computational insights into active site shaping for substrate specificity and reaction regioselectivity in the EXTL2 retaining glycosyltransferase

QM(DFT)/MM calculations and molecular dynamics simulations on wild-type retaining α1,4-N-acetylhexosaminyltransferase (EXTL2) and Arg293Ala, Asp246Ala, Arg293Ala/Asp246Ala and Asp246Glu mutants are used to understand the role of these two residues.

Different QM/MM Approaches To Elucidate Enzymatic Reactions: Case Study on ppGalNAcT2

Hybrid QM/MM computational studies can provide invaluable insight into the mechanisms of enzymatic reactions that can be exploited for rational drug design. Various approaches are available for such studies. However, their strengths and weaknesses may not be immediately apparent. Using the retaining glycosyltransferase ppGalNAcT2 as a case study, we compare different methodologies used to obtain reaction paths and transition state information. Ab Initio MD using CPMD coupled with the String Method is used to derive the minimum free energy reaction path. The geometrical features of the free energy path, especially around the transition state, agree with the minimum potential energy path obtained by the much less computationally expensive Nudged Elastic Band method. The barrier energy, however, differs by 8 kcal/mol. The free energy surface generated by metadynamics provides a rough overview of the reaction and can confirm the physical relevance of optimized paths or provide an initial guess for path optimization methods. Calculations of enzymatic reactions are usually performed at best at the DFT level of theory. A comparison of widely used functionals with high-level DLPNO-CCSD(T)/CBS data on the potential energy profile serves as a validation of the usability of DFT for this type of enzymatic reaction. The M06-2X meta-hybrid functional in particular matches the DLPNO-CCSD(T)/CBS reference extremely well with errors within 1 kcal/mol. However, even pure-GGA functional OPBE provides sufficient accuracy for this type of study.

Quantum Mechanics/Molecular Mechanics Study of the Sialyltransferase Reaction Mechanism

The sialyltransferase is an enzyme that transfers the sialic acid moiety from cytidine 5'-monophospho-N-acetyl-neuraminic acid (CMP-NeuAc) to the terminal position of glycans. To elucidate the catalytic mechanism of sialyltransferase, we explored the potential energy surface along the sialic acid transfer reaction coordinates by the hybrid quantum mechanics/molecular mechanics method on the basis of the crystal structure of sialyltransferase CstII. Our calculation demonstrated that CstII employed an S_N1-like reaction mechanism via the formation of a short-lived oxocarbenium ion intermediate. The computational barrier height was 19.5 kcal/mol, which reasonably corresponded with the experimental reaction rate. We also found that two tyrosine residues (Tyr156 and Tyr162) played a vital role in stabilizing the intermediate and the transition states by quantum mechanical interaction with CMP.

Insights into the catalytic mechanism of N-acetylglucosaminidase glycoside hydrolase from Bacillus subtilis: a QM/MM study

QM/MM calculations on NagZs fromBacillus subtilisfurther confirm NagZs to be glycoside phosphorylases rather than glycoside hydrolases.

Exploring reaction pathways for O-GlcNAc transferase catalysis. A string method study

The inverting O-GlcNAc glycosyltransferase (OGT) is an important post-translation enzyme, which catalyzes the transfer of N-acetylglucosamine from UDP-N-acetylglucosamine (UDP-GlcNAc) to the hydroxyl group of the Ser/Thr of cytoplasmic, nuclear, and mitochondrial proteins. In the past, three different catalytic bases were proposed for the reaction: His498, α-phosphate, and Asp554. In this study, we used hybrid quantum mechanics/molecular mechanics (QM/MM) Car-Parrinello molecular dynamics to investigate reaction paths using α-phosphate and Asp554 as the catalytic bases. The string method was used to calculate the free-energy reaction profiles of the tested mechanisms. During the investigations, an additional mechanism was observed. In this mechanism, a proton is transferred to α-phosphate via a water molecule. Our calculations show that the mechanism with α-phosphate acting as the base is favorable. This reaction has a rate-limiting free-energy barrier of 23.5 kcal/mol, whereas reactions utilizing Asp554 and water-assisted α-phosphate have barriers of 41.7 and 40.9 kcal/mol, respectively. Our simulations provide a new insight into the catalysis of OGT and may thus guide rational drug design of transition-state analogue inhibitors with potential therapeutic use.

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.

Ab initio modelling of the anomeric and exo anomeric effects in 2-methoxytetrahydropyran and 2-methoxythiane corrected for intramolecular BSSE

Oxygen substitution in 2-methoxytetrahydropyran by sulphur in 2-methoxythiane approximately doubles the anomeric effect, which slows down enzymatic processing of 2-methoxythiane.