Deep evolutionary analysis reveals the design principles of fold A glycosyltransferases

Glycosyltransferases: glycoengineers in human milk oligosaccharide synthesis and manufacturing

Human milk oligosaccharides (HMOs) are a diverse group of complex carbohydrates that play crucial roles in infant health, promoting a beneficial gut microbiota, modulating immune responses, and protecting against pathogens. Central to the synthesis of HMOs are glycosyltransferases, a specialized class of enzymes that catalyse the transfer of sugar moieties to form the complex glycan structures characteristic of HMOs. This review provides an in-depth analysis of glycosyltransferases, beginning with their classification based on structural and functional characteristics. The catalytic activity of these enzymes is explored, highlighting the mechanisms by which they facilitate the precise addition of monosaccharides in HMO biosynthesis. Structural insights into glycosyltransferases are also discussed, shedding light on how their conformational features enable specific glycosidic bond formations. This review maps out the key biosynthetic pathways involved in HMO production, including the synthesis of lactose, and subsequent fucosylation and sialylation processes, all of which are intricately regulated by glycosyltransferases. Industrial methods for HMO synthesis, including chemical, enzymatic, and microbial approaches, are examined, emphasizing the role of glycosyltransferases in these processes. Finally, the review discusses future directions in glycosyltransferase research, particularly in enhancing the efficiency of HMO synthesis and developing advanced analytical techniques to better understand the structural complexity and biological functions of HMOs.

Molecular structure and enzymatic mechanism of the human collagen hydroxylysine galactosyltransferase GLT25D1/COLGALT1

During collagen biosynthesis, lysine residues undergo extensive post-translational modifications through the alternate action of two distinct metal ion-dependent enzyme families (i.e., LH/PLODs and GLT25D/COLGALT), ultimately producing the highly conserved α-(1,2)-glucosyl-β-(1,O)-galactosyl-5-hydroxylysine pattern. Malfunctions in these enzymes are linked to developmental pathologies and extracellular matrix alterations associated to enhanced aggressiveness of solid tumors. Here, we characterized human GLT25D1/COLGALT1, revealing an elongated head-to-head homodimeric assembly. Each monomer encompasses two domains (named GT1 and GT2), both unexpectedly capable of binding metal ion cofactors and UDP-α-galactose donor substrates, resulting in four candidate catalytic sites per dimer. We identify the catalytic site in GT2, featuring an unusual Glu-Asp-Asp motif critical for Mn2+ binding, ruling out direct catalytic roles for the GT1 domain, but showing that in this domain the unexpectedly bound Ca2+ and UDP-α-galactose cofactors are critical for folding stability. Dimerization, albeit not essential for GLT25D1/COLGALT1 activity, provides a critical molecular contact site for multi-enzyme assembly interactions with partner multifunctional LH/PLOD lysyl hydroxylase-glycosyltransferase enzymes.

Cryo-EM SPR structures of Salmonella typhimurium ArnC; the key enzyme in lipid-A modification conferring polymyxin resistance

Polymyxins are last-resort antimicrobial peptides administered clinically against multi-drug resistant bacteria, specifically in the case of Gram-negative species. However, an increasing number of these pathogens employ a defense strategy that involves a relay of enzymes encoded by the pmrE (ugd) loci and the arnBCDTEF operon. The pathway modifies the lipid-A component of the outer membrane (OM) lipopolysaccharide (LPS) by adding a 4-amino-4-deoxy-l-arabinose (L-Ara4N) headgroup, which renders polymyxins ineffective. Here, we report the cryo-EM SPR structures of glycosyltransferase ArnC from Salmonella typhimurium determined in apo and UDP-bound forms at resolutions 2.75 Å and 3.8 Å, respectively. The structure of the ArnC protomer comprises three distinct regions: an N-terminal glycosyltransferase domain, transmembrane region, and the interface helices (IHs). ArnC forms a tetramer with C2 symmetry, where the C-terminal strand inserts into the adjacent protomer. This tetrameric state is further stabilized by two distinct interfaces formed by ArnC that form a network of hydrogen bonds and salt bridges. The binding of UDP induces conformational changes that stabilize the loop between residues H201 to S213, and part of the putative catalytic pocket formed by IH1 and IH2. The surface property analysis revealed a hydrophobic cavity formed by TM1 and TM2 in the apo state, which is disrupted upon UDP binding. The comparison of ArnC structures to their homologs GtrB and DPMS suggests the key residues involved in ArnC catalytic activity.

Structure Analysis and Site-Directed Mutagenesis of the Glycosyltransferase UGT71B8 Leads to Increased Stability and Substrate Activity in Arabidopsis thaliana

The uridine diphosphate-glycosyltransferase (UGT) family catalyses the glucuronidation of the glycosyl group of a nucleotide sugar to an acceptor compound (substrate), it serves as controlling reaction for bioactivity, storage and decrease toxicity of different compounds in living organisms. UGT71B8 belongs to 71B family of UGTs and its donor sugars are UDP glucose, UDP galactose and UDP 5S glucose, respectively. The current study was designed to induce site-directed mutagenesis (SDM) to investigate the activity in UGT71B8 enzyme. During first step, in silico conformational change through 3D structure model was drawn and it was found that all the amino acids of mutation site were found in allowed region. The relative surface accessibility (RSA) and absolute surface accessibility (ASA) of UGT71B8 were found as 0.042-0.037 and 7.424, respectively, which shows that UGT71B8 T138M remains stable after SDM. This prediction model thus led to the efficacious mutation of UGT71B8 enzyme. Mass spectrometric analysis of UGT71B8T138M showed reduced activity with its substrate UDP glucose and kaempherol as acceptor molecule. Moreover, no new substrate activity of UGT71B8 was found. This data would direct future endeavors to engineer more glycosyltransferases of plants to augment its activity with different substrates and provide a basis for more exploration of UGT71B8 as an active compound for potential anti-cancer therapeutics.

Global atlas of predicted functional domains in Legionella pneumophila Dot/Icm translocated effectors

Legionella pneumophila utilizes the Dot/Icm type IVB secretion system to deliver hundreds of effector proteins inside eukaryotic cells to ensure intracellular replication. Our understanding of the molecular functions of the largest pathogenic arsenal known to the bacterial world remains incomplete. By leveraging advancements in 3D protein structure prediction, we provide a comprehensive structural analysis of 368 L. pneumophila effectors, representing a global atlas of predicted functional domains summarized in a database ( https://pathogens3d.org/legionella-pneumophila ). Our analysis identified 157 types of diverse functional domains in 287 effectors, including 159 effectors with no prior functional annotations. Furthermore, we identified 35 cryptic domains in 30 effector models that have no similarity with experimentally structurally characterized proteins, thus, hinting at novel functionalities. Using this analysis, we demonstrate the activity of thirteen functional domains, including three cryptic domains, predicted in L. pneumophila effectors to cause growth defects in the Saccharomyces cerevisiae model system. This illustrates an emerging strategy of exploring synergies between predictions and targeted experimental approaches in elucidating novel effector activities involved in infection.

Structural insights into translocation and tailored synthesis of hyaluronan

Hyaluronan (HA) is an essential component of the vertebrate extracellular matrix. It is a heteropolysaccharide of N-acetylglucosamine (GlcNAc) and glucuronic acid (GlcA) reaching several megadaltons in healthy tissues. HA is synthesized and translocated in a coupled reaction by HA synthase (HAS). Here, structural snapshots of HAS provide insights into HA biosynthesis, from substrate recognition to HA elongation and translocation. We monitor the extension of a GlcNAc primer with GlcA, reveal the coordination of the uridine diphosphate product by a conserved gating loop and capture the opening of a translocation channel to coordinate a translocating HA polymer. Furthermore, we identify channel-lining residues that modulate HA product lengths. Integrating structural and biochemical analyses suggests an avenue for polysaccharide engineering based on finely tuned enzymatic activity and HA coordination.

Evolution of sequence, structural and functional diversity of the ubiquitous DNA/RNA-binding Alba domain

The DNA/RNA-binding Alba domain is prevalent across all kingdoms of life. First discovered in archaea, this protein domain has evolved from RNA- to DNA-binding, with a concomitant expansion in the range of cellular processes that it regulates. Despite its widespread presence, the full extent of its sequence, structural, and functional diversity remains unexplored. In this study, we employed iterative searches in PSI-BLAST to identify 15,161 unique Alba domain-containing proteins from the NCBI non-redundant protein database. Sequence similarity network (SSN) analysis clustered them into 13 distinct subgroups, including the archaeal Alba and eukaryotic Rpp20/Pop7 and Rpp25/Pop6 groups, as well as novel fungal and Plasmodium-specific Albas. Sequence and structural conservation analysis of the subgroups indicated high preservation of the dimer interface, with Alba domains from unicellular eukaryotes notably exhibiting structural deviations towards their C-terminal end. Finally, phylogenetic analysis, while supporting SSN clustering, revealed the evolutionary branchpoint at which the eukaryotic Rpp20- and Rpp25-like clades emerged from archaeal Albas, and the subsequent taxonomic lineage-based divergence within each clade. Taken together, this comprehensive analysis enhances our understanding of the evolutionary history of Alba domain-containing proteins across diverse organisms.

Unveiling the GT114 family: Structural characterization of A075L, a glycosyltransferase from Paramecium bursaria chlorella virus-1 (PBCV-1)

Protein A075L is a β-xylosyltransferase that participates in producing the core of the N-glycans found in VP54, the major viral capsid protein of Paramecium bursaria chlorella virus-1 (PBCV-1). In this study, we present an X-ray crystallographic analysis of the apo form of A075L, along with its complexes with the sugar donor and with a trisaccharide acceptor. The protein structure shows a typical GT-B folding, with two Rossmann-like fold domains, in which the acceptor substrate binds to the N-terminal region, and the nucleotide-sugar donor binds to the C-terminal region. We propose that the catalytic mechanism follows a direct displacement SN2-like reaction, where Asp73 serves as a catalytic base that deprotonates the incoming nucleophile of the acceptor, facilitating direct displacement of the UDP with the inversion of the anomeric configuration of the acceptor without metal ion dependence, while the interactions with side chains of Arg158 and Arg208 stabilize the developing negative charge. Using isothermal titration calorimetry, nuclear magnetic resonance spectroscopy, high-performance liquid chromatography, and molecular dynamics simulations, the catalytic activity and specificity of this enzyme have been unraveled.

Protein representations: Encoding biological information for machine learning in biocatalysis

Enzymes offer a more environmentally friendly and low-impact solution to conventional chemistry, but they often require additional engineering for their application in industrial settings, an endeavour that is challenging and laborious. To address this issue, the power of machine learning can be harnessed to produce predictive models that enable the in silico study and engineering of improved enzymatic properties. Such machine learning models, however, require the conversion of the complex biological information to a numerical input, also called protein representations. These inputs demand special attention to ensure the training of accurate and precise models, and, in this review, we therefore examine the critical step of encoding protein information to numeric representations for use in machine learning. We selected the most important approaches for encoding the three distinct biological protein representations - primary sequence, 3D structure, and dynamics - to explore their requirements for employment and inductive biases. Combined representations of proteins and substrates are also introduced as emergent tools in biocatalysis. We propose the division of fixed representations, a collection of rule-based encoding strategies, and learned representations extracted from the latent spaces of large neural networks. To select the most suitable protein representation, we propose two main factors to consider. The first one is the model setup, which is influenced by the size of the training dataset and the choice of architecture. The second factor is the model objectives such as consideration about the assayed property, the difference between wild-type models and mutant predictors, and requirements for explainability. This review is aimed at serving as a source of information and guidance for properly representing enzymes in future machine learning models for biocatalysis.

Enzymatic Routes to Designer Hemicelluloses for Use in Biobased Materials

Various enzymes can be used to modify the structure of hemicelluloses directly in vivo or following extraction from biomass sources, such as wood and agricultural residues. Generally, these enzymes can contribute to designer hemicelluloses through four main strategies: (1) enzymatic hydrolysis such as selective removal of side groups by glycoside hydrolases (GH) and carbohydrate esterases (CE), (2) enzymatic cross-linking, for instance, the selective addition of side groups by glycosyltransferases (GT) with activated sugars, (3) enzymatic polymerization by glycosynthases (GS) with activated glycosyl donors or transglycosylation, and (4) enzymatic functionalization, particularly via oxidation by carbohydrate oxidoreductases and via amination by amine transaminases. Thus, this Perspective will first highlight enzymes that play a role in regulating the degree of polymerization and side group composition of hemicelluloses, and subsequently, it will explore enzymes that enhance cross-linking capabilities and incorporate novel chemical functionalities into saccharide structures. These enzymatic routes offer a precise way to tailor the properties of hemicelluloses for specific applications in biobased materials, contributing to the development of renewable alternatives to conventional materials derived from fossil fuels.

Computational-guided discovery of UDP-glycosyltransferases for lauryl glucoside production using engineered E. coli

Defining suitable enzymes for reaction steps in novel synthetic pathways is crucial for developing microbial cell factories for non-natural products. Here, we developed a computational workflow to identify C12 alcohol-active UDP-glycosyltransferases. The workflow involved three steps: (1) assembling initial candidates of putative UDP-glycosyltransferases, (2) refining selection by examining conserved regions, and (3) 3D structure prediction and molecular docking. Genomic sequences from Candida, Pichia, Rhizopus, and Thermotoga, known for lauryl glucoside synthesis via whole-cell biocatalysis, were screened. Out of 240 predicted glycosyltransferases, 8 candidates annotated as glycosyltransferases were selected after filtering out those with signal peptides and identifying conserved UDP-glycosyltransferase regions. These proteins underwent 3D structure prediction and molecular docking with 1-dodecanol. RO3G, a candidate from Rhizopus delemar RA 99-880 with a relatively high ChemPLP fitness score, was selected and expressed in Escherichia coli BL21 (DE3). It was further characterized using a feeding experiment with 1-dodecanol. Results confirmed that the RO3G-expressing strain could convert 1-dodecanol to lauryl glucoside, as quantified by HPLC and identified by targeted LC-MS. Monitoring the growth and fermentation profiles of the engineered strain revealed that RO3G expression did not affect cell growth. Interestingly, acetate, a major fermentation product, was reduced in the RO3G-expressing strain compared to the GFP-expressing strain, suggesting a redirection of flux from acetate to other pathways. Overall, this work presents a successful workflow for discovering UDP-glycosyltransferase enzymes with confirmed activity toward 1-dodecanol for lauryl glucoside production.

Structure of a truncated human GlcNAc-1-phosphotransferase variant reveals the basis for its hyperactivity

Mutations that cause loss of function of GlcNAc-1-phosphotransferase (PTase) lead to the lysosomal storage disorder mucolipidosis II. PTase is the key enzyme of the mannose 6-phosphate (M6P) targeting system that is responsible for tagging lysosomal hydrolases with the M6P moiety for their delivery to the lysosome. We had previously generated a truncated hyperactive form of PTase termed S1S3 which was shown to notably increase the phosphorylation level of secreted lysosomal enzymes and enhance their uptake by cells. Here, we report the 3.4 Å cryo-EM structure of soluble S1S3 lacking both transmembrane domains and cytosolic tails. The structure reveals a high degree of conservation of the catalytic core to full-length PTase. In this dimeric structure, the EF-hand of one protomer is observed interacting with the conserved region four of the other. In addition, we present a high-quality EM 3D map of the UDP-GlcNAc bound form of the full-length soluble protein showing the key molecular interactions between the nucleotide sugar donor and side chain amino acids of the protein. Finally, although the domain organization of S1S3 is very similar to that of the Drosophila melanogaster (fruit fly) PTase homolog, we establish that the latter does not act on lysosomal hydrolases.

Virus-encoded glycosyltransferases hypermodify DNA with diverse glycans

Enzymatic modification of DNA nucleobases can coordinate gene expression, nuclease protection, or mutagenesis. We recently discovered a clade of phage-specific cytosine methyltransferase (MT) and 5-methylpyrimidine dioxygenase (5mYOX) enzymes that produce 5-hydroxymethylcytosine (5hmC) as a precursor for enzymatic hypermodifications on viral genomes. Here, we identify phage MT- and 5mYOX-associated glycosyltransferases (GTs) that catalyze linkage of diverse sugars to 5hmC nucleobase substrates. Metavirome mining revealed thousands of biosynthetic gene clusters containing enzymes with predicted roles in cytosine sugar hypermodification. We developed a platform for high-throughput screening of GT-containing pathways, relying on the Escherichia coli metabolome as a substrate pool. We successfully reconstituted several pathways and isolated diverse sugar modifications appended to cytosine, including mono-, di-, or tri-saccharides comprised of hexoses, N-acetylhexosamines, or heptose. These findings expand our knowledge of hypermodifications on nucleic acids and the origins of corresponding sugar-installing enzymes.

Structure, biosynthesis and regulation of the T1 antigen, a phase-variable surface polysaccharide conserved in many Salmonella serovars

The bacterial genus Salmonella includes diverse isolates with multiple variations in the structure of the main polysaccharide component (O antigen) of membrane lipopolysaccharides. In addition, some isolates produce a transient (T) antigen, such as the T1 polysaccharide identified in the 1960s in an isolate of Salmonella enterica Paratyphi B. The structure and biosynthesis of the T1 antigen have remained enigmatic. Here, we use biophysical, biochemical and genetic methods to show that the T1 antigen is a complex linear glycan containing tandem homopolymeric domains of galactofuranose and ribofuranose, linked to lipid A-core, like a typical O antigen. T1 is a phase-variable antigen, regulated by recombinational inversion of the promoter upstream of the T1 genetic locus through a mechanism not observed for other bacterial O antigens. The T1 locus is conserved across many Salmonella isolates, but is mutated or absent in most typhoidal serovars and in serovar Enteritidis.

Klebsiella pneumoniae O-polysaccharide biosynthesis highlights the diverse organization of catalytic modules in ABC transporter-dependent glycan assembly

Klebsiella pneumoniae provides influential prototypes for lipopolysaccharide O antigen (OPS) biosynthesis in Gram-negative bacteria. Sequences of OPS-biosynthesis gene clusters in serotypes O4 and O7 suggest fundamental differences in the organization of required enzyme modules compared to other serotypes. Furthermore, some required activities were not assigned by homology shared with characterized enzymes. The goal of this study was therefore to resolve the serotype O4 and O7 pathways to expand our broader understanding of glycan polymerization and chain termination processes. The O4 and O7 antigens were produced from cloned genetic loci in recombinant Escherichia coli. Systematic in vivo and in vitro approaches were then applied to assign each enzyme in each of the pathways, defining the necessary components for polymerization and chain termination. OPS assembly is accomplished by multiprotein complexes formed by interactions between polymerase components variably distributed in single and multimodule proteins. In each complex, a terminator function is present in a protein containing a characteristic coiled-coil molecular ruler, which determines glycan chain length. In serotype O4, we discovered a CMP-α-3-deoxy-ᴅ-manno-octulosonic acid-dependent chain-terminating glycosyltransferase that is the founding member of a new glycosyltransferase family (GT137) and potentially identifies a new glycosyltransferase fold. The O7 OPS is terminated by a methylphosphate moiety, like the K. pneumoniae O3 antigen, but the methyltransferase-kinase enzyme pairs responsible for termination in these serotypes differ in sequence and predicted structures. Together, the characterization of O4 and O7 has established unique enzyme activities and provided new insight into glycan-assembly strategies that are widely distributed in bacteria.

GASP: A Pan-Specific Predictor of Family 1 Glycosyltransferase Acceptor Specificity Enabled by a Pipeline for Substrate Feature Generation and Large-Scale Experimental Screening

Glycosylation represents a major chemical challenge; while it is one of the most common reactions in Nature, conventional chemistry struggles with stereochemistry, regioselectivity, and solubility issues. In contrast, family 1 glycosyltransferase (GT1) enzymes can glycosylate virtually any given nucleophilic group with perfect control over stereochemistry and regioselectivity. However, the appropriate catalyst for a given reaction needs to be identified among the tens of thousands of available sequences. Here, we present the glycosyltransferase acceptor specificity predictor (GASP) model, a data-driven approach to the identification of reactive GT1:acceptor pairs. We trained a random forest-based acceptor predictor on literature data and validated it on independent in-house generated data on 1001 GT1:acceptor pairs, obtaining an AUROC of 0.79 and a balanced accuracy of 72%. The performance was stable even in the case of completely new GT1s and acceptors not present in the training data set, highlighting the pan-specificity of GASP. Moreover, the model is capable of parsing all known GT1 sequences, as well as all chemicals, the latter through a pipeline for the generation of 153 chemical features for a given molecule taking the CID or SMILES as input (freely available at https://github.com/degnbol/GASP). To investigate the power of GASP, the model prediction probability scores were compared to GT1 substrate conversion yields from a newly published data set, with the top 50% of GASP predictions corresponding to reactions with >50% synthetic yields. The model was also tested in two comparative case studies: glycosylation of the antihelminth drug niclosamide and the plant defensive compound DIBOA. In the first study, the model achieved an 83% hit rate, outperforming a hit rate of 53% from a random selection assay. In the second case study, the hit rate of GASP was 50%, and while being lower than the hit rate of 83% using expert-selected enzymes, it provides a reasonable performance for the cases when an expert opinion is unavailable. The hierarchal importance of the generated chemical features was investigated by negative feature selection, revealing properties related to cyclization and atom hybridization status to be the most important characteristics for accurate prediction. Our study provides a GT1:acceptor predictor which can be trained on other data sets enabled by the automated feature generation pipelines. We also release the new in-house generated data set used for testing of GASP to facilitate the future development of GT1 activity predictors and their robust benchmarking.

Tyrosine-modifying glycosylation by Yersinia effectors

Mono-O-glycosylation of target proteins by bacterial toxins or effector proteins is a well-known mechanism by which bacteria interfere with essential functions of host cells. The respective glycosyltransferases are important virulence factors such as the Clostridioides difficile toxins A and B. Here, we describe two glycosyltransferases of Yersinia species that have a high sequence identity: YeGT from the zoonotic pathogen Yersinia enterocolitica and YkGT from the murine pathogen Yersinia kristensenii. We show that both modify Rho family proteins by attachment of GlcNAc at tyrosine residues (Tyr-34 in RhoA). Notably, the enzymes differed in their target protein specificity. While YeGT modified RhoA, B, and C, YkGT possessed a broader substrate spectrum and glycosylated not only Rho but also Rac and Cdc42 subfamily proteins. Mutagenesis studies indicated that residue 177 is important for this broader target spectrum. We determined the crystal structure of YeGT shortened by 16 residues N terminally (sYeGT) in the ligand-free state and bound to UDP, the product of substrate hydrolysis. The structure assigns sYeGT to the GT-A family. It shares high structural similarity to glycosyltransferase domains from toxins. We also demonstrated that the 16 most N-terminal residues of YeGT and YkGT are important for the mediated translocation into the host cell using the pore-forming protective antigen of anthrax toxin. Mediated introduction into HeLa cells or ectopic expression of YeGT and YkGT caused morphological changes and redistribution of the actin cytoskeleton. The data suggest that YeGT and YkGT are likely bacterial effectors belonging to the family of tyrosine glycosylating bacterial glycosyltransferases.

The sugar donor specificity of plant family 1 glycosyltransferases

Plant family 1 glycosyltransferases (UGTs) represent a formidable tool to produce valuable natural and novel glycosides. Their regio- and stereo-specific one-step glycosylation mechanism along with their inherent wide acceptor scope are desirable traits in biotechnology. However, their donor scope and specificity are not well understood. Since different sugars have different properties in vivo and in vitro, the ability to easily glycodiversify target acceptors is desired, and this depends on our improved understanding of the donor binding site. In the aim to unlock the full potential of UGTs, studies have attempted to elucidate the structure-function relationship governing their donor specificity. These efforts have revealed a complex phenomenon, and general principles valid for multiple enzymes are elusive. Here, we review the studies of UGT donor specificity, and attempt to group the information into key concepts which can help shape future research. We zoom in on the family-defining PSPG motif, on two loop residues reported to interact with the C6 position of the sugar, and on the role of active site arginines in donor specificity. We continue to discuss attempts to alter and expand the donor specificity by enzyme engineering, and finally discuss future research directions.

The underlying mechanisms of arenaviral entry through matriglycan

Matriglycan, a recently characterized linear polysaccharide, is composed of alternating xylose and glucuronic acid subunits bound to the ubiquitously expressed protein α-dystroglycan (α-DG). Pathogenic arenaviruses, like the Lassa virus (LASV), hijack this long linear polysaccharide to gain cellular entry. Until recently, it was unclear through what mechanisms LASV engages its matriglycan receptor to initiate infection. Additionally, how matriglycan is synthesized onto α-DG by the Golgi-resident glycosyltransferase LARGE1 remained enigmatic. Recent structural data for LARGE1 and for the LASV spike complex informs us about the synthesis of matriglycan as well as its usage as an entry receptor by arenaviruses. In this review, we discuss structural insights into the system of matriglycan generation and eventual recognition by pathogenic viruses. We also highlight the unique usage of matriglycan as a high-affinity host receptor compared with other polysaccharides that decorate cells.

The evolving world of milk oligosaccharides: Biochemical diversity understood by computational advances

Milk oligosaccharides, complex carbohydrates unique to mammalian milk, play crucial roles in infant nutrition and immune development. This review explores their biochemical diversity, tracing the evolutionary paths that have led to their variation across different species. We highlight the intersection of nutrition, biology, and chemistry in understanding these compounds. Additionally, we discuss the latest computational methods and analytical techniques that have revolutionized the study of milk oligosaccharides, offering insights into their structural complexity and functional roles. This brief but essential review not only aims to provide a deeper understanding of milk oligosaccharides but also discuss the road toward their potential applications.

Machine-Learning-Assisted Nanozyme Design: Lessons from Materials and Engineered Enzymes

Nanozymes are nanomaterials that exhibit enzyme-like biomimicry. In combination with intrinsic characteristics of nanomaterials, nanozymes have broad applicability in materials science, chemical engineering, bioengineering, biochemistry, and disease theranostics. Recently, the heterogeneity of published results has highlighted the complexity and diversity of nanozymes in terms of consistency of catalytic capacity. Machine learning (ML) shows promising potential for discovering new materials, yet it remains challenging for the design of new nanozymes based on ML approaches. Alternatively, ML is employed to promote optimization of intelligent design and application of catalytic materials and engineered enzymes. Incorporation of the successful ML algorithms used in the intelligent design of catalytic materials and engineered enzymes can concomitantly facilitate the guided development of next-generation nanozymes with desirable properties. Here, recent progress in ML, its utilization in the design of catalytic materials and enzymes, and how emergent ML applications serve as promising strategies to circumvent challenges associated with time-expensive and laborious testing in nanozyme research and development are summarized. The potential applications of successful examples of ML-aided catalytic materials and engineered enzymes in nanozyme design are also highlighted, with special focus on the unified aims in enhancing design and recapitulation of substrate selectivity and catalytic activity.

Glycosyltransferase 8 domain-containing protein 1 (GLT8D1) is a UDP-dependent galactosyltransferase

Glycosyltransferases (GTs) are enzymes that catalyze the formation of glycosidic bonds and hundreds of GTs have been identified so far in humans. Glycosyltransferase 8 domain-containing protein 1 (GLT8D1) has been associated with central nervous system diseases and cancer. However, evidence on its enzymatic properties, including its substrates, has been scarcely described. In this paper, we have produced and purified recombinant secretory GLT8D1. The enzyme was found to be N-glycosylated. Differential scanning fluorimetry was employed to analyze the stabilization of GLT8D1 by Mn2+ and nucleotides, revealing UDP as the most stabilizing nucleotide scaffold. GLT8D1 displayed glycosyltransferase activity from UDP-galactose onto N-acetylgalactosamine but with a low efficiency. Modeling of the structure revealed similarities with other GT-A fold enzymes in CAZy family GT8 and glycosyltransferases in other families with galactosyl-, glucosyl-, and xylosyltransferase activities, each with retaining catalytic mechanisms. Our study provides novel structural and functional insights into the properties of GLT8D1 with implications in pathological processes.

Rewiring the pneumococcal capsule pathway for investigating glycosyltransferase specificity and genetic glycoengineering

Virtually all living cells are covered with glycans. Their structures are primarily controlled by the specificities of glycosyltransferases (GTs). GTs typically adopt one of the three folds, namely, GT-A, GT-B, and GT-C. However, what defines their specificities remain poorly understood. Here, we developed a genetic glycoengineering platform by reprogramming the capsular polysaccharide pathways in Streptococcus pneumoniae to interrogate GT specificity and manipulate glycan structures. Our findings suggest that the central cleft of GT-B enzymes is important for determining acceptor specificity. The constraint of the glycoengineering platform was partially alleviated when the specificity of the precursor transporter was reduced, indicating that the transporter contributes to the overall fidelity of glycan synthesis. We also modified the pneumococcal capsule to produce several medically important mammalian glycans, as well as demonstrated the importance of regiochemistry in a glycosidic linkage on binding lung epithelial cells. Our work provided mechanistic insights into GT specificity and an approach for investigating glycan functions.

Mammalian milk glycomes: Connecting the dots between evolutionary conservation and biosynthetic pathways

Milk oligosaccharides (MOs) are among the most abundant constituents of breast milk and are essential for health and development. Biosynthesized from monosaccharides into complex sequences, MOs differ considerably between taxonomic groups. Even human MO biosynthesis is insufficiently understood, hampering evolutionary and functional analyses. Using a comprehensive resource of all published MOs from >100 mammals, we develop a pipeline for generating and analyzing MO biosynthetic networks. We then use evolutionary relationships and inferred intermediates of these networks to discover (1) systematic glycome biases, (2) biosynthetic restrictions, such as reaction path preference, and (3) conserved biosynthetic modules. This allows us to prune and pinpoint biosynthetic pathways despite missing information. Machine learning and network analysis cluster species by their milk glycome, identifying characteristic sequence relationships and evolutionary gains/losses of motifs, MOs, and biosynthetic modules. These resources and analyses will advance our understanding of glycan biosynthesis and the evolution of breast milk.

Bioinformatics characterization of BcsA-like orphan proteins suggest they form a novel family of pseudomonad cyclic-β-glucan synthases

Bacteria produce a variety of polysaccharides with functional roles in cell surface coating, surface and host interactions, and biofilms. We have identified an 'Orphan' bacterial cellulose synthase catalytic subunit (BcsA)-like protein found in four model pseudomonads, P. aeruginosa PA01, P. fluorescens SBW25, P. putida KT2440 and P. syringae pv. tomato DC3000. Pairwise alignments indicated that the Orphan and BcsA proteins shared less than 41% sequence identity suggesting they may not have the same structural folds or function. We identified 112 Orphans among soil and plant-associated pseudomonads as well as in phytopathogenic and human opportunistic pathogenic strains. The wide distribution of these highly conserved proteins suggest they form a novel family of synthases producing a different polysaccharide. In silico analysis, including sequence comparisons, secondary structure and topology predictions, and protein structural modelling, revealed a two-domain transmembrane ovoid-like structure for the Orphan protein with a periplasmic glycosyl hydrolase family GH17 domain linked via a transmembrane region to a cytoplasmic glycosyltransferase family GT2 domain. We suggest the GT2 domain synthesises β-(1,3)-glucan that is transferred to the GH17 domain where it is cleaved and cyclised to produce cyclic-β-(1,3)-glucan (CβG). Our structural models are consistent with enzymatic characterisation and recent molecular simulations of the PaPA01 and PpKT2440 GH17 domains. It also provides a functional explanation linking PaPAK and PaPA14 Orphan (also known as NdvB) transposon mutants with CβG production and biofilm-associated antibiotic resistance. Importantly, cyclic glucans are also involved in osmoregulation, plant infection and induced systemic suppression, and our findings suggest this novel family of CβG synthases may provide similar range of adaptive responses for pseudomonads.

In-depth analysis of biocatalysts by microfluidics: An emerging source of data for machine learning

Nowadays, the vastly increasing demand for novel biotechnological products is supported by the continuous development of biocatalytic applications which provide sustainable green alternatives to chemical processes. The success of a biocatalytic application is critically dependent on how quickly we can identify and characterize enzyme variants fitting the conditions of industrial processes. While miniaturization and parallelization have dramatically increased the throughput of next-generation sequencing systems, the subsequent characterization of the obtained candidates is still a limiting process in identifying the desired biocatalysts. Only a few commercial microfluidic systems for enzyme analysis are currently available, and the transformation of numerous published prototypes into commercial platforms is still to be streamlined. This review presents the state-of-the-art, recent trends, and perspectives in applying microfluidic tools in the functional and structural analysis of biocatalysts. We discuss the advantages and disadvantages of available technologies, their reproducibility and robustness, and readiness for routine laboratory use. We also highlight the unexplored potential of microfluidics to leverage the power of machine learning for biocatalyst development.

Spindly is a nucleocytosolic O-fucosyltransferase in Dictyostelium and related proteins are widespread in protists and bacteria

O-GlcNAcylation is a prominent modification of nuclear and cytoplasmic proteins in animals and plants and is mediated by a single O-GlcNAc transferase (OGT). Spindly (Spy), a paralog of OGT first discovered in higher plants, has an ortholog in the apicomplexan parasite Toxoplasma gondii, and both enzymes are now recognized as O-fucosyltransferases (OFTs). Here we investigate the evolution of spy-like genes and experimentally confirm OFT activity in the social amoeba Dictyostelium-a protist that is more related to fungi and metazoa. Immunofluorescence probing with the fucose-specific Aleuria aurantia lectin (AAL) and biochemical cell fractionation combined with western blotting suggested the occurrence of nucleocytoplasmic fucosylation. The absence of reactivity in mutants deleted in spy or gmd (unable to synthesize GDP-Fuc) suggested monofucosylation mediated by Spy. Genetic ablation of the modE locus, previously predicted to encode a GDP-fucose transporter, confirmed its necessity for fucosylation in the secretory pathway but not for the nucleocytoplasmic proteins. Affinity capture of these proteins combined with mass spectrometry confirmed monofucosylation of Ser and Thr residues of several known nucleocytoplasmic proteins. As in Toxoplasma, the Spy OFT was required for optimal proliferation of Dictyostelium under laboratory conditions. These findings support a new phylogenetic analysis of OGT and OFT evolution that indicates their occurrence in the last eukaryotic common ancestor but mostly complementary presence in its eukaryotic descendants with the notable exception that both occur in red algae and plants. Their generally exclusive expression, high degree of conservation, and shared monoglycosylation targets suggest overlapping roles in physiological regulation.

Prediction and verification of glycosyltransferase activity by bioinformatics analysis and protein engineering

A significant number of proteins are annotated as functionally uncharacterized proteins. Within this protocol, we describe how to use protein family multiple sequence alignments and structural bioinformatics resources to design loss-of-function mutations of previously uncharacterized proteins within the glycosyltransferase family. We detail approaches to determine target protein active sites using three-dimensional modeling. We generate active site mutants and quantify any changes in enzymatic function by a glycosyltransferase assay. With modifications, this protocol could be applied to other metal-dependent enzymes. For complete details on the use and execution of this protocol, please refer to Ilina et al. (2022).1.

Structural and biochemical insight into a modular β-1,4-galactan synthase in plants

Rhamnogalacturonan I (RGI) is a structurally complex pectic polysaccharide with a backbone of alternating rhamnose and galacturonic acid residues substituted with arabinan and galactan side chains. Galactan synthase 1 (GalS1) transfers galactose and arabinose to either extend or cap the β-1,4-galactan side chains of RGI, respectively. Here we report the structure of GalS1 from Populus trichocarpa, showing a modular protein consisting of an N-terminal domain that represents the founding member of a new family of carbohydrate-binding module, CBM95, and a C-terminal glycosyltransferase family 92 (GT92) catalytic domain that adopts a GT-A fold. GalS1 exists as a dimer in vitro, with stem domains interacting across the chains in a 'handshake' orientation that is essential for maintaining stability and activity. In addition to understanding the enzymatic mechanism of GalS1, we gained insight into the donor and acceptor substrate binding sites using deep evolutionary analysis, molecular simulations and biochemical studies. Combining all the results, a mechanism for GalS1 catalysis and a new model for pectic galactan side-chain addition are proposed.

Integrative Bioinformatics Approaches Indicate a Particular Pattern of Some SARS-CoV-2 and Non-SARS-CoV-2 Proteins

Pattern recognition plays a critical role in integrative bioinformatics to determine the structural patterns of proteins of viruses such as SARS-CoV-2. This study identifies the pattern of SARS-CoV-2 proteins to depict the structure-function relationships of the protein alphabets of SARS-CoV-2 and COVID-19. The assembly enumeration algorithm, Anisotropic Network Model, Gaussian Network Model, Markovian Stochastic Model, and image comparison protein-like alphabets were used. The distance score was the lowest with 22 for "I" and highest with 40 for "9". For post-processing and decision, two protein alphabets "C" (PDB ID: 6XC3) and "S" (PDB ID: 7OYG) were evaluated to understand the structural, functional, and evolutionary relationships, and we found uniqueness in the functionality of proteins. Here, models were constructed using "SARS-CoV-2 proteins" (12 numbers) and "non-SARS-CoV-2 proteins" (14 numbers) to create two words, "SARS-CoV-2" and "COVID-19". Similarly, we developed two slogans: "Vaccinate the world against COVID-19" and "Say no to SARS-CoV-2", which were made with the proteins structure. It might generate vaccine-related interest to broad reader categories. Finally, the evolutionary process appears to enhance the protein structure smoothly to provide suitable functionality shaped by natural selection.

Structure of the human heparan sulfate polymerase complex EXT1-EXT2

Heparan sulfates are complex polysaccharides that mediate the interaction with a broad range of protein ligands at the cell surface. A key step in heparan sulfate biosynthesis is catalyzed by the bi-functional glycosyltransferases EXT1 and EXT2, which generate the glycan backbone consisting of repeating N-acetylglucosamine and glucuronic acid units. The molecular mechanism of heparan sulfate chain polymerization remains, however, unknown. Here, we present the cryo-electron microscopy structure of human EXT1-EXT2, which reveals the formation of a tightly packed hetero-dimeric complex harboring four glycosyltransferase domains. A combination of in vitro and in cellulo mutational studies is used to dissect the functional role of the four catalytic sites. While EXT1 can catalyze both glycosyltransferase reactions, our results indicate that EXT2 might only have N-acetylglucosamine transferase activity. Our findings provide mechanistic insight into heparan sulfate chain elongation as a nonprocessive process and lay the foundation for future studies on EXT1-EXT2 function in health and disease.

Polymerization of the backbone of the pectic polysaccharide rhamnogalacturonan I

Rhamnogalacturonan I (RG-I) is a major plant cell wall pectic polysaccharide defined by its repeating disaccharide backbone structure of [4)-α-D-GalA-(1,2)-α-L-Rha-(1,]. A family of RG-I:Rhamnosyltransferases (RRT) has previously been identified, but synthesis of the RG-I backbone has not been demonstrated in vitro because the identity of Rhamnogalacturonan I:Galaturonosyltransferase (RG-I:GalAT) was unknown. Here a putative glycosyltransferase, At1g28240/MUCI70, is shown to be an RG-I:GalAT. The name RGGAT1 is proposed to reflect the catalytic activity of this enzyme. When incubated together with the rhamnosyltransferase RRT4, the combined activities of RGGAT1 and RRT4 result in elongation of RG-I acceptors in vitro into a polymeric product. RGGAT1 is a member of a new GT family categorized as GT116, which does not group into existing GT-A clades and is phylogenetically distinct from the GALACTURONOSYLTRANSFERASE (GAUT) family of GalA transferases that synthesize the backbone of the pectin homogalacturonan. RGGAT1 has a predicted GT-A fold structure but employs a metal-independent catalytic mechanism that is rare among glycosyltransferases with this fold type. The identification of RGGAT1 and the 8-member Arabidopsis GT116 family provides a new avenue for studying the mechanism of RG-I synthesis and the function of RG-I in plants.

Glycoinformatics in the Artificial Intelligence Era

Artificial intelligence (AI) methods have been and are now being increasingly integrated in prediction software implemented in bioinformatics and its glycoscience branch known as glycoinformatics. AI techniques have evolved in the past decades, and their applications in glycoscience are not yet widespread. This limited use is partly explained by the peculiarities of glyco-data that are notoriously hard to produce and analyze. Nonetheless, as time goes, the accumulation of glycomics, glycoproteomics, and glycan-binding data has reached a point where even the most recent deep learning methods can provide predictors with good performance. We discuss the historical development of the application of various AI methods in the broader field of glycoinformatics. A particular focus is placed on shining a light on challenges in glyco-data handling, contextualized by lessons learnt from related disciplines. Ending on the discussion of state-of-the-art deep learning approaches in glycoinformatics, we also envision the future of glycoinformatics, including development that need to occur in order to truly unleash the capabilities of glycoscience in the systems biology era.

Common and unique features of glycosylation and glycosyltransferases in African trypanosomes

Eukaryotic protein glycosylation is mediated by glycosyl- and oligosaccharyl-transferases. Here, we describe how African trypanosomes exhibit both evolutionary conservation and significant divergence compared with other eukaryotes in how they synthesise their glycoproteins. The kinetoplastid parasites have conserved components of the dolichol-cycle and oligosaccharyltransferases (OSTs) of protein N-glycosylation, and of glycosylphosphatidylinositol (GPI) anchor biosynthesis and transfer to protein. However, some components are missing, and they process and decorate their N-glycans and GPI anchors in unique ways. To do so, they appear to have evolved a distinct and functionally flexible glycosyltransferases (GT) family, the GT67 family, from an ancestral eukaryotic β3GT gene. The expansion and/or loss of GT67 genes appears to be dependent on parasite biology. Some appear to correlate with the obligate passage of parasites through an insect vector, suggesting they were acquired through GT67 gene expansion to assist insect vector (tsetse fly) colonisation. Others appear to have been lost in species that subsequently adopted contaminative transmission. We also highlight the recent discovery of a novel and essential GT11 family of kinetoplastid parasite fucosyltransferases that are uniquely localised to the mitochondria of Trypanosoma brucei and Leishmania major. The origins of these kinetoplastid FUT1 genes, and additional putative mitochondrial GT genes, are discussed.

Oxygen-dependent regulation of E3(SCF)ubiquitin ligases and a Skp1-associated JmjD6 homolog in development of the social amoeba Dictyostelium

E3-SCF (Skp1/cullin-1/F-box protein) polyubiquitin ligases activate the proteasomal degradation of over a thousand proteins, but the evolutionary diversification of the F-box protein (FBP) family of substrate receptor subunits has challenged their elucidation in protists. Here, we expand the FBP candidate list in the social amoeba Dictyostelium and show that the Skp1 interactome is highly remodeled as cells transition from growth to multicellular development. Importantly, a subset of candidate FBPs was less represented when the posttranslational hydroxylation and glycosylation of Skp1 was abrogated by deletion of the O2-sensing Skp1 prolyl hydroxylase PhyA. A role for this Skp1 modification for SCF activity was indicated by partial rescue of development, which normally depends on high O2 and PhyA, of phyA-KO cells by proteasomal inhibitors. Further examination of two FBPs, FbxwD and the Jumonji C protein JcdI, suggested that Skp1 was substituted by other factors in phyA-KO cells. Although a double-KO of jcdI and its paralog jcdH did not affect development, overexpression of JcdI increased its sensitivity to O2. JcdI, a nonheme dioxygenase shown to have physiological O2 dependence, is conserved across protists with its F-box and other domains, and is related to the human oncogene JmjD6. Sensitization of JcdI-overexpression cells to O2 depended on its dioxygenase activity and other domains, but not its F-box, which may however be the mediator of its reduced levels in WT relative to Skp1 modification mutant cells. The findings suggest that activation of JcdI by O2 is tempered by homeostatic downregulation via PhyA and association with Skp1.

Modularity of the hydrophobic core and evolution of functional diversity in fold A glycosyltransferases

Hydrophobic cores are fundamental structural properties of proteins typically associated with protein folding and stability; however, how the hydrophobic core shapes protein evolution and function is poorly understood. Here, we investigated the role of conserved hydrophobic cores in fold-A glycosyltransferases (GT-As), a large superfamily of enzymes that catalyze formation of glycosidic linkages between diverse donor and acceptor substrates through distinct catalytic mechanisms (inverting versus retaining). Using hidden Markov models and protein structural alignments, we identify similarities in the phosphate-binding cassette (PBC) of GT-As and unrelated nucleotide-binding proteins, such as UDP-sugar pyrophosphorylases. We demonstrate that GT-As have diverged from other nucleotide-binding proteins through structural elaboration of the PBC and its unique hydrophobic tethering to the F-helix, which harbors the catalytic base (xED-Asp). While the hydrophobic tethering is conserved across diverse GT-A fold enzymes, some families, such as B3GNT2, display variations in tethering interactions and core packing. We evaluated the structural and functional impact of these core variations through experimental mutational analysis and molecular dynamics simulations and find that some of the core mutations (T336I in B3GNT2) increase catalytic efficiency by modulating the conformational occupancy of the catalytic base between “D-in” and acceptor-accessible “D-out” conformation. Taken together, our studies support a model of evolution in which the GT-A core evolved progressively through elaboration upon an ancient PBC found in diverse nucleotide-binding proteins, and malleability of this core provided the structural framework for evolving new catalytic and substrate-binding functions in extant GT-A fold enzymes.

Structure-based design of UDP-GlcNAc analogs as candidate GnT-V inhibitors

BACKGROUND

N-Glycan branching regulates various functions of glycoproteins. N-Acetylglucosaminyltransferase V (GnT-V) is a GlcNAc transferase that acts on N-glycans and the GnT-V-producing branch is highly related to cancer progression. This indicates that specific GnT-V inhibitors may be drug candidates for cancer treatment. To design novel GnT-V inhibitors, we focused on the unique and weak recognition of the donor substrate UDP-GlcNAc by GnT-V. On the basis of the catalytic pocket structure, we hypothesized that UDP-GlcNAc analogs with increasing hydrophobicity may be GnT-V inhibitors.

METHODS

We chemically synthesized 10 UDP-GlcNAc analogs in which one or two phosphate groups were replaced with hydrophobic groups. To test these compounds, we set up an HPLC-based enzyme assay system for all N-glycan-branching GlcNAc transferases in which GnT-I-V activity was measured using purified truncated enzymes. Using this system, we assessed the inhibitory effects of the synthesized compounds on GnT-V and their specificity.

RESULTS

Several UDP-GlcNAc analogs inhibited GnT-V activity, although the inhibition potency was modest. Compared with other GnTs, these compounds showed a preference for GnT-V, which suggested that GnT-V was relatively tolerant of hydrophobicity in the donor substrate. Docking models of the inhibitory compounds with GnT-V suggested the mechanisms of how these compounds interacted with GnT-V and inhibited its action.

CONCLUSIONS

Chemical modification of the donor substrate may be a promising strategy to develop selective inhibitors of GnT-V.

GENERAL SIGNIFICANCE

Our findings provide new insights into the design of GnT inhibitors and how GnTs recognize the donor substrate.

Successive remodeling of IgG glycans using a solid-phase enzymatic platform

The success of glycoprotein-based drugs in various disease treatments has become widespread. Frequently, therapeutic glycoproteins exhibit a heterogeneous array of glycans that are intended to mimic human glycopatterns. While immunogenic responses to biologic drugs are uncommon, enabling exquisite control of glycosylation with minimized microheterogeneity would improve their safety, efficacy and bioavailability. Therefore, close attention has been drawn to the development of glycoengineering strategies to control the glycan structures. With the accumulation of knowledge about the glycan biosynthesis enzymes, enzymatic glycan remodeling provides a potential strategy to construct highly ordered glycans with improved efficiency and biocompatibility. In this study, we quantitatively evaluate more than 30 enzymes for glycoengineering immobilized immunoglobulin G, an impactful glycoprotein class in the pharmaceutical field. We demonstrate successive glycan remodeling in a solid-phase platform, which enabled IgG glycan harmonization into a series of complex-type N-glycoforms with high yield and efficiency while retaining native IgG binding affinity.

A solid-phase glycan remodeling (SPGR) platform is presented. Over thirty enzymes were evaluated for successive glycoengineering of immobilized antibodies with outstanding performance in several SPGR reactions.

Enzymatic activity of glycosyltransferase GLT8D1 promotes human glioblastoma cell migration

Glioblastoma (GBM) is the most aggressive primary brain tumor characterized by infiltrative growth of malignant glioma cells into the surrounding brain parenchyma. In this study, our analysis of GBM patient cohorts revealed a significantly higher expression of Glycosyltransferase 8 domain containing 1 (GLT8D1) compared to normal brain tissue and could be associated with impaired patient survival. Increased in vitro expression of GLT8D1 significantly enhanced migration of two different sphere-forming GBM cell lines. By in silico analysis we predicted the 3D-structure as well as the active site residues of GLT8D1. The introduction of point mutations in the predicted active site reduced its glycosyltransferase activity in vitro and consequently impaired GBM tumor cell migration. Examination of GLT8D1 interaction partners by LC-MS/MS implied proteins associated with cytoskeleton and intracellular transport as potential substrates. In conclusion, we demonstrated that the enzymatic activity of glycosyltransferase GLT8D1 promotes GBM cell migration.

•

The glycosyltransferase GLT8D1 is enriched in GBM tissue and cells

•

In silico analysis predicts the 3D structure and the active site of GLT8D1

•

Enzymatically active GLT8D1 promotes GBM migration

•

Manipulation of GLT8D1 enzymatic activity decreases GBM migration

Heparosan oligosaccharide synthesis using engineered single-function glycosyltransferases

A 24-fold increase in GlcNAc-transferase ability through KfiA screening and engineering. An approach for heparosan oligosaccharide synthesis relying on single-function glycosyltransferases.

Detecting Glucose Fluctuations in the Campylobacter jejuni N-Glycan Structure

Campylobacter jejuni is a significant cause of human gastroenteritis worldwide, and all strains express an N-glycan that is added to at least 80 different proteins. We characterized 98 C. jejuni isolates from infants from 7 low- and middle-income countries and identified 4 isolates unreactive with our N-glycan-specific antiserum that was raised against the C. jejuni heptasaccharide composed of GalNAc-GalNAc-GalNAc(Glc)-GalNAc-GalNAc-diNAcBac. Mass spectrometric analyses indicated these isolates express a hexasaccharide lacking the glucose branch. Although all 4 strains encode the PglI glucosyltransferase (GlcTF), one aspartate in the DXDD motif was missing, an alteration also present in ∼4% of all available PglI sequences. Deleting this residue from an active PglI resulted in a nonfunctional GlcTF when the protein glycosylation system was reconstituted in E. coli, while replacement with Glu/Ala was not deleterious. Molecular modeling proposed a mechanism for how the DXDD residues and the structure/length beyond the motif influence activity. Mouse vaccination with an E. coli strain expressing the full-length heptasaccharide produced N-glycan-specific antibodies and a corresponding reduction in Campylobacter colonization and weight loss following challenge. However, the antibodies did not recognize the hexasaccharide and were unable to opsonize C. jejuni isolates lacking glucose, suggesting this should be considered when designing N-glycan-based vaccines to prevent campylobacteriosis.

Therapeutic target database update 2022: facilitating drug discovery with enriched comparative data of targeted agents

Drug discovery relies on the knowledge of not only drugs and targets, but also the comparative agents and targets. These include poor binders and non-binders for developing discovery tools, prodrugs for improved therapeutics, co-targets of therapeutic targets for multi-target strategies and off-target investigations, and the collective structure-activity and drug-likeness landscapes of enhanced drug feature. However, such valuable data are inadequately covered by the available databases. In this study, a major update of the Therapeutic Target Database, previously featured in NAR, was therefore introduced. This update includes (a) 34 861 poor binders and 12 683 non-binders of 1308 targets; (b) 534 prodrug-drug pairs for 121 targets; (c) 1127 co-targets of 672 targets regulated by 642 approved and 624 clinical trial drugs; (d) the collective structure-activity landscapes of 427 262 active agents of 1565 targets; (e) the profiles of drug-like properties of 33 598 agents of 1102 targets. Moreover, a variety of additional data and function are provided, which include the cross-links to the target structure in PDB and AlphaFold, 159 and 1658 newly emerged targets and drugs, and the advanced search function for multi-entry target sequences or drug structures. The database is accessible without login requirement at: https://idrblab.org/ttd/.

Mapping the glycosyltransferase fold landscape using interpretable deep learning

Glycosyltransferases (GTs) play fundamental roles in nearly all cellular processes through the biosynthesis of complex carbohydrates and glycosylation of diverse protein and small molecule substrates. The extensive structural and functional diversification of GTs presents a major challenge in mapping the relationships connecting sequence, structure, fold and function using traditional bioinformatics approaches. Here, we present a convolutional neural network with attention (CNN-attention) based deep learning model that leverages simple secondary structure representations generated from primary sequences to provide GT fold prediction with high accuracy. The model learns distinguishing secondary structure features free of primary sequence alignment constraints and is highly interpretable. It delineates sequence and structural features characteristic of individual fold types, while classifying them into distinct clusters that group evolutionarily divergent families based on shared secondary structural features. We further extend our model to classify GT families of unknown folds and variants of known folds. By identifying families that are likely to adopt novel folds such as GT91, GT96 and GT97, our studies expand the GT fold landscape and prioritize targets for future structural studies.

The polyfunctional polysialic acid: A structural view

Polysialic acid (polySia), a homopolymer of α2,8-linked sialic acid residues, modifies a small number of proteins and has central functions in vertebrate signalling. Here, we review the regulatory functions of polySia in signalling processes and the immune system of adult humans, as well as functions based on their chemical properties. The main focus will be on the structure-function relationship of polySia with its interaction partners in humans. Recent studies have indicated that the degree of polymerisation is an important parameter that can guide the regulatory effect of polySia in addition to its binding to target proteins. Therefore, the structures of polySia in solution and bound to interaction partners are compared in order to identify the key factors that define binding specificity.

GTXplorer: A portal to navigate and visualize the evolutionary information encoded in fold a glycosyltransferases

Glycosyltransferases (GTs) play a central role in sustaining all forms of life through the biosynthesis of complex carbohydrates. Despite significant strides made in recent years to establish computational resources, databases, and tools to understand the nature and role of carbohydrates and related glycoenzymes, a data analytics framework that connects the sequence-structure-function relationships to the evolution of GTs is currently lacking. This hinders the characterization of understudied GTs and the synthetic design of GTs for medical and biotechnology applications. Here, we present GTXplorer as an integrated platform that presents evolutionary information of GTs adopting a GT-A fold in an intuitive format enabling in silico investigation through comparative sequence analysis to derive informed hypotheses about their function. The tree view mode provides an overview of the evolutionary relationships of GT-A families and allows users to select phylogenetically relevant families for comparisons. The selected families can then be compared in the alignment view at the residue level using annotated weblogo stacks of the GT-A core specific to the selected clade, family, or subfamily. All data are easily accessible and can be downloaded for further analysis. GTXplorer can be accessed at https://vulcan.cs.uga.edu/gtxplorer/ or from GitHub at https://github.com/esbgkannan/GTxplorer to deploy locally. By packaging multiple data streams into an accessible, user-friendly format, GTXplorer presents the first evolutionary data analytics platform for comparative glycomics.

Peters plus syndrome mutations affect the function and stability of human β1,3-glucosyltransferase

Peters Plus Syndrome (PTRPLS OMIM #261540) is a severe congenital disorder of glycosylation where patients have multiple structural anomalies, including Peters anomaly of the eye (anterior segment dysgenesis), disproportionate short stature, brachydactyly, dysmorphic facial features, developmental delay, and variable additional abnormalities. PTRPLS patients and some Peters Plus-like (PTRPLS-like) patients (who only have a subset of PTRPLS phenotypes, have mutations in the gene encoding β1,3-glucosyltransferase [B3GLCT]). B3GLCT catalyzes the transfer of glucose to O-linked fucose on thrombospondin type-1 repeats. Most B3GLCT substrate proteins belong to the ADAMTS superfamily and play critical roles in extracellular matrix. We sought to determine whether the PTRPLS or PTRPLS-like mutations abrogated B3GLCT activity. B3GLCT has two putative active sites, one in the N-terminal region and the other in the C-terminal glycosyltransferase domain. Using sequence analysis and in vitro activity assays, we demonstrated that the C-terminal domain catalyzes transfer of glucose to O-linked fucose. We also generated a homology model of B3GLCT and identified D421 as the catalytic base. PTRPLS and PTRPLS-like mutations were individually introduced into B3GLCT, and the mutated enzymes were evaluated using in vitro enzyme assays and cell-based functional assays. Our results demonstrated that PTRPLS mutations caused loss of B3GLCT enzymatic activity and/or significantly reduced protein stability. In contrast, B3GLCT with PTRPLS-like mutations retained enzymatic activity, although some showed a minor destabilizing effect. Overall, our data supports the hypothesis that loss of glucose from B3GLCT substrate proteins is responsible for the defects observed in PTRPLS patients, but not for those observed in PTRPLS-like patients.

Biochemical and Functional Characterization of GALT8, an Arabidopsis GT31 β-(1,3)-Galactosyltransferase That Influences Seedling Development

Arabinogalactan-proteins (AGPs) are members of the hydroxyproline-rich glycoprotein (HRGP) superfamily, a group of highly diverse proteoglycans that are present in the cell wall, plasma membrane as well as secretions of almost all plants, with important roles in many developmental processes. The role of GALT8 (At1g22015), a Glycosyltransferase-31 (GT31) family member of the Carbohydrate-Active Enzyme database (CAZy), was examined by biochemical characterization and phenotypic analysis of a galt8 mutant line. To characterize its catalytic function, GALT8 was heterologously expressed in tobacco leaves and its enzymatic activity tested. GALT8 was shown to be a β-(1,3)-galactosyltransferase (GalT) that catalyzes the synthesis of a β-(1,3)-galactan, similar to the in vitro activity of KNS4/UPEX1 (At1g33430), a homologous GT31 member previously shown to have this activity. Liquid chromatography-tandem mass spectrometry (LC-MS/MS) confirmed the products were of 2-6 degree of polymerisation (DP). Previous reporter studies showed that GALT8 is expressed in the central and synergid cells, from whence the micropylar endosperm originates after the fertilization of the central cell of the ovule. Homozygous mutants have multiple seedling phenotypes including significantly shorter hypocotyls and smaller leaf area compared to wild type (WT) that are attributable to defects in female gametophyte and/or endosperm development. KNS4/UPEX1 was shown to partially complement the galt8 mutant phenotypes in genetic complementation assays suggesting a similar but not identical role compared to GALT8 in β-(1,3)-galactan biosynthesis. Taken together, these data add further evidence of the important roles GT31 β-(1,3)-GalTs play in elaborating type II AGs that decorate AGPs and pectins, thereby imparting functional consequences on plant growth and development.

Missing the sweet spot: one of the two N-glycans on human Gb3/CD77 synthase is expendable

N-glycosylation is a ubiquitous posttranslational modification that may influence folding, subcellular localization, secretion, solubility and oligomerization of proteins. In this study, we examined the effects of N-glycans on the activity of human Gb3/CD77 synthase, which catalyzes the synthesis of glycosphingolipids with terminal Galα1 → 4Gal (Gb3 and the P1 antigen) and Galα1 → 4GalNAc disaccharides (the NOR antigen). The human Gb3/CD77 synthase contains two occupied N-glycosylation sites at positions N121 and N203. Intriguingly, we found that while the N-glycan at N203 is essential for activity and correct subcellular localization, the N-glycan at N121 is dispensable and its absence did not reduce, but, surprisingly, even increased the activity of the enzyme. The fully N-glycosylated human Gb3/CD77 synthase and its glycoform missing the N121 glycan correctly localized in the Golgi, whereas a glycoform without the N203 site partially mislocalized in the endoplasmic reticulum. A double mutein missing both N-glycans was inactive and accumulated in the endoplasmic reticulum. Our results suggest that the decreased specific activity of human Gb3/CD77 synthase glycovariants results from their improper subcellular localization and, to a smaller degree, a decrease in enzyme solubility. Taken together, our findings show that the two N-glycans of human Gb3/CD77 synthase have opposing effects on its properties, revealing a dual nature of N-glycosylation and potentially a novel regulatory mechanism controlling the biological activity of proteins.

Recent progress in synthesis of carbohydrates with sugar nucleotide-dependent glycosyltransferases

Sugar nucleotide-dependent glycosyltransferases (GTs) are key enzymes that catalyze the formation of glycosidic bonds in nature. They have been increasingly applied in the synthesis of complex carbohydrates and glycoconjugates with or without in situ generation of sugar nucleotides. Human GTs are becoming more accessible and new bacterial GTs have been identified and characterized. An increasing number of crystal structures elucidated for GTs from mammalian and bacterial sources facilitate structure-based design of mutants as improved catalysts for synthesis. Automated platforms have also been developed for chemoenzymatic synthesis of carbohydrates. Recent progress in applying sugar nucleotide-dependent GTs in enzymatic and chemoenzymatic synthesis of mammalian glycans and glycoconjugates, bacterial surface glycans, and glycosylated natural products from bacteria and plants are reviewed.