Influence of protein/glycan interaction on site-specific glycan heterogeneity

A hierarchical structure in the N-glycosylation process governs the N-glycosylation output: prolonged cultivation induces glycoenzymes expression variations that are reflected in the cellular N-glycome but not in the protein and site-specific glycoprofile of CHO cells

N-glycosylation is a central component in the modification of secretory proteins. One characteristic of this process is a heterogeneous output. The heterogeneity is the result of both structural constraints of the glycoprotein as well as the composition of the cellular glycosylation machinery. Empirical data addressing correlations between glycosylation output and glycosylation machinery composition are seldom due to the low abundance of glycoenzymes. We assessed how differences in the glycoenzyme expression affected the N-glycosylation output at a cellular as well as at a protein-specific level. Our results showed that cellular N-glycome changes could be correlated with the variation of glycoenzyme expression, whereas at the protein level differential responses to glycoenzymes alterations were observed. We therefore identified a hierarchical structure in the N-glycosylation process: the enzyme levels in this complex pathway determine its capacity (reflected in the N-glycome), while protein-specific parameters determine the glycosite-specificity. What emerges is a highly variable and adaptable protein modification system that represents a hallmark of eukaryotic cells.

Decoding glycosylation potential from protein structure across human glycoproteins with a multi-view recurrent neural network

Glycosylation is described as a non-templated biosynthesis. Yet, the template-free premise is antithetical to the observation that different N-glycans are consistently placed at specific sites. It has been proposed that glycosite-proximal protein structures could constrain glycosylation and explain the observed microheterogeneity. Using site-specific glycosylation data, we trained a hybrid neural network to parse glycosites (recurrent neural network) and match them to feasible N-glycosylation events (graph neural network). From glycosite-flanking sequences, the algorithm predicts most human N-glycosylation events documented in the GlyConnect database and proposed structures corresponding to observed monosaccharide composition of the glycans at these sites. The algorithm also recapitulated glycosylation in Enhanced Aromatic Sequons, SARS-CoV-2 spike, and IgG3 variants, thus demonstrating the ability of the algorithm to predict both glycan structure and abundance. Thus, protein structure constrains glycosylation, and the neural network enables predictive in silico glycosylation of uncharacterized or novel protein sequences and genetic variants.

SLC3A2 N-glycosylation and Golgi remodeling regulate SLC7A amino acid exchangers and stress mitigation

Proteostasis requires oxidative metabolism (ATP) and mitigation of the associated damage by glutathione, in an increasingly dysfunctional relationship with aging. SLC3A2 (4F2hc, CD98) plays a role as a disulfide-linked adaptor to the SLC7A5 and SLC7A11 exchangers which import essential amino acids and cystine while exporting Gln and Glu, respectively. The positions of N-glycosylation sites on SLC3A2 have evolved with the emergence of primates, presumably in synchrony with metabolism. Herein, we report that each of the four sites in SLC3A2 has distinct profiles of Golgi-modified N-glycans. N-glycans at the primate-derived site N381 stabilized SLC3A2 in the galectin-3 lattice against coated-pit endocytosis, while N365, the site nearest the membrane promoted glycolipid-galectin-3 (GL-Lect)-driven endocytosis. Our results indicate that surface retention and endocytosis are precisely balanced by the number, position, and remodeling of N-glycans on SLC3A2. Furthermore, proteomics and functional assays revealed an N-glycan-dependent clustering of the SLC3A2∗SLC7A5 heterodimer with amino-acid/Na+ symporters (SLC1A4, SLC1A5) that balances branched-chain amino acids and Gln levels, at the expense of ATP to maintain the Na+/K+ gradient. In replete conditions, SLC3A2 interactions require Golgi-modified N-glycans at N365D and N381D, whereas reducing N-glycosylation in the endoplasmic reticulum by fluvastatin treatment promoted the recruitment of CD44 and transporters needed to mitigate stress. Thus, SLC3A2 N-glycosylation and Golgi remodeling of the N-glycans have distinct roles in amino acids import for growth, maintenance, and metabolic stresses.

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.

Sequential in vitro enzymatic N-glycoprotein modification reveals site-specific rates of glycoenzyme processing

N-glycosylation is an essential eukaryotic posttranslational modification that affects various glycoprotein properties, including folding, solubility, protein–protein interactions, and half-life. N-glycans are processed in the secretory pathway to form varied ensembles of structures, and diversity at a single site on a glycoprotein is termed ‘microheterogeneity’. To understand the factors that influence glycan microheterogeneity, we hypothesized that local steric and electrostatic factors surrounding each site influence glycan availability for enzymatic modification. We tested this hypothesis via expression of reporter N-linked glycoproteins in N-acetylglucosaminyltransferase MGAT1-null HEK293 cells to produce immature Man₅GlcNAc₂ glycoforms (38 glycan sites total). These glycoproteins were then sequentially modified in vitro from high mannose to hybrid and on to biantennary, core-fucosylated, complex structures by a panel of N-glycosylation enzymes, and each reaction time course was quantified by LC-MS/MS. Substantial differences in rates of in vitro enzymatic modification were observed between glycan sites on the same protein, and differences in modification rates varied depending on the glycoenzyme being evaluated. In comparison, proteolytic digestion of the reporters prior to N-glycan processing eliminated differences in in vitro enzymatic modification. Furthermore, comparison of in vitro rates of enzymatic modification with the glycan structures found on the mature reporters expressed in WT cells correlated well with the enzymatic bottlenecks observed in vivo. These data suggest higher order local structures surrounding each glycosylation site contribute to the efficiency of modification both in vitro and in vivo to establish the spectrum of microheterogeneity in N-linked glycoproteins.

Decoding human-macaque interspecies differences in Fc-effector functions: The structural basis for CD16-dependent effector function in Rhesus macaques

Fc mediated effector functions of antibodies play important roles in immunotherapies and vaccine efficacy but assessing those functions in animal models can be challenging due to species differences. Rhesus macaques, Macaca mulatta (Mm) share approximately 93% sequence identity with humans but display important differences in their adaptive immune system that complicates their use in validating therapeutics and vaccines that rely on Fc effector functions. In contrast to humans, macaques only have one low affinity FcγRIII receptor, CD16, which shares a polymorphism at position 158 with human FcγRIIIa with Ile158 and Val158 variants. Here we describe structure-function relationships of the Ile/Val158 polymorphism in Mm FcγRIII. Our data indicate that the affinity of the allelic variants of Mm FcγRIII for the macaque IgG subclasses vary greatly with changes in glycan composition both on the Fc and the receptor. However, unlike the human Phe/Val158 polymorphism in FcγRIIIa, the higher affinity variant corresponds to the larger, more hydrophobic side chain, Ile, even though it is not directly involved in the binding interface. Instead, this side chain appears to modulate glycan-glycan interactions at the Fc/FcγRIII interface. Furthermore, changes in glycan composition on the receptor have a greater effect for the Val158 variant such that with oligomannose type glycans and with glycans only on Asn45 and Asn162, Val158 becomes the variant with higher affinity to Fc. These results have implications not only for the better interpretation of nonhuman primate studies but also for studies performed with human effector cells carrying different FcγRIIIa alleles.

Production of microhomogeneous glycopeptide by a mutated NGT according FuncLib with unique sugar as substrate

N-glycosylation is one of the most important post-translational modifications of proteins. Cytoplasmic soluble N-glycosyltransferase (NGT) exists in bacteria, which is able to transfer different monosaccharide from sugar nucleotide to the NXS/T(X ≠ Pro) consensus sequence in a polypeptide. At present, the NGT enzymes reported could transfer a variety of different sugars to protein, which will lead to the heterogeneity of the sugar chain and the complexity and instability of the structure and function of glycopeptides. According to the FuncLib algorithm, we obtained mutant ApNGT-P1 from ApNGT (the NGT from Actinobacillus pleuropneumoniae) with increased substrate specificity. Compared with the wild-type ApNGT, mutant ApNGT-P1 could only utilize UDP-Glc as sugar donors. The optimum temperature of ApNGT-P1 was about 40 °C and the optimum pH was 7.5-8.0 in PBS buffer. ApNGT-P1 exhibited better tolerance for K⁺, Mn²⁺, Ca²⁺, and Mg²⁺, but was strongly inhibited by Na⁺, Cu²⁺ and Zn²⁺. The mutant can be applied to the efficient production of glycosylated peptides or proteins with uniform glucose at their glycosylation sites. Besides, this work provided a feasible pathway for further studies on the improving donor substrates selectivity of NGTs.

A graphical representation of glycan heterogeneity

A substantial shortcoming of large scale datasets is often the inability to easily represent and visualize key features. This problem becomes acute when considering the increasing technical ability to profile large numbers of glycopeptides and glycans in recent studies. Here, we describe a simple, concise graphical representation intended to capture the microheterogeneity associated with glycan modification at specific sites. We illustrate this method by showing visual representations of the glycans and glycopeptides from a variety of species. The graphical representation presented allows one to easily discern the compositions of all glycans, similarities and differences of modifications found in different samples, and, in the case of N-linked glycans, the initial steps in the biosynthetic pathway.

N-Glycosylation Enhances Conformational Flexibility of Protein Disulfide Isomerase Revealed by Microsecond Molecular Dynamics and Markov State Modeling

Secreted proteins of eukaryotes are decorated with branched carbohydrate oligomers called glycans. This fact is only starting to be considered for in silico investigations of protein dynamics. Using all-atom molecular dynamics (MD) simulations and Markov state modeling (MSM), we unveil the influence of glycans on the conformational flexibility of the multidomain protein disulfide isomerase (PDI), which is a ubiquitous chaperone in the endoplasmic reticulum (ER). Yeast PDI (yPDI) from Saccharomyces cerevisiae is glycosylated at asparagine side chains and the knowledge of its five modified sites enables a realistic computational modeling. We compare simulations of glycosylated and unglycosylated yPDI and find that the presence of glycan-glycan and glycan-protein interactions influences the flexibility of PDI in different ways. For example, glycosylation reduces interdomain interactions, shifting the conformational ensemble toward more open, extended structures. In addition, we compare our results on yPDI with structural information of homologous proteins such as human PDI (hPDI), which is natively unglycosylated. Interestingly, hPDI lacks a surface recess that is present in yPDI. We find that glycosylation of yPDI facilitates its catalytic site to reach close to this surface recess. Hence, this might point to a possible functional relevance of glycosylation in yeast to act on substrates, while glycosylation seems redundant for the human homologous protein. We conclude that glycosylation is fundamental for protein dynamics, making it a necessity for a truthful representation of the flexibility and function in in silico studies of glycoproteins.

Towards rational glyco-engineering in CHO: from data to predictive models

Metabolic modelling strives to develop modelling approaches that are robust and highly predictive. To achieve this, various modelling designs, including hybrid models, and parameter estimation methods that define the type and number of parameters used in the model, are adapted. Accurate input data play an important role so that the selection of experimental methods that provide input data of the required precision with low measurement errors is crucial. For the biopharmaceutically relevant protein glycosylation, the most prominent available models are kinetic models which are able to capture the dynamic nature of protein N-glycosylation. In this review we focus on how to choose the most suitable model for a specific research question, as well as on parameters and considerations to take into account before planning relevant experiments.

Glycoengineering Chinese hamster ovary cells: a short history

Biotherapeutic glycoproteins have revolutionised the field of pharmaceuticals, with new discoveries and continuous improvements underpinning the rapid growth of this industry. N-glycosylation is a critical quality attribute of biotherapeutic glycoproteins that influences the efficacy, half-life and immunogenicity of these drugs. This review will focus on the advances and future directions of remodelling N-glycosylation in Chinese hamster ovary (CHO) cells, which are the workhorse of recombinant biotherapeutic production, with particular emphasis on antibody products, using strategies such as cell line and protein backbone engineering.

Glycan-protein interactions determine kinetics of N-glycan remodeling

A hallmark of N-linked glycosylation in the secretory compartments of eukaryotic cells is the sequential remodeling of an initially uniform oligosaccharide to a site-specific, heterogeneous ensemble of glycostructures on mature proteins. To understand site-specific processing, we used protein disulfide isomerase (PDI), a model protein with five glycosylation sites, for molecular dynamics (MD) simulations and compared the result to a biochemical in vitro analysis with four different glycan processing enzymes. As predicted by an analysis of the accessibility of the N-glycans for their processing enzymes derived from the MD simulations, N-glycans at different glycosylation sites showed different kinetic properties for the processing enzymes. In addition, altering the tertiary structure of the glycoprotein PDI affected its N-glycan remodeling in a site-specific way. We propose that the observed differential N-glycan reactivities depend on the surrounding protein tertiary structure and lead to different glycan structures in the same protein through kinetically controlled processing pathways.

Atomistic glycoprotein simulations reveal a site-specific availability of glycan substrates in time-resolved mass spectrometry of maturating enzyme kinetics.

Big-Data Glycomics: Tools to Connect Glycan Biosynthesis to Extracellular Communication

Characteristically, cells must sense and respond to environmental cues. Despite the importance of cell-cell communication, our understanding remains limited and often lacks glycans. Glycans decorate proteins and cell membranes at the cell-environment interface, and modulate intercellular communication, from development to pathogenesis. Providing further challenges, glycan biosynthesis and cellular behavior are co-regulating systems. Here, we discuss how glycosylation contributes to extracellular responses and signaling. We further organize approaches for disentangling the roles of glycans in multicellular interactions using newly available datasets and tools, including glycan biosynthesis models, omics datasets, and systems-level analyses. Thus, emerging tools in big data analytics and systems biology are facilitating novel insights on glycans and their relationship with multicellular behavior.

Synthetic Glycobiology: Parts, Systems, and Applications

Protein glycosylation, the attachment of sugars to amino acid side chains, can endow proteins with a wide variety of properties of great interest to the engineering biology community. However, natural glycosylation systems are limited in the diversity of glycoproteins they can synthesize, the scale at which they can be harnessed for biotechnology, and the homogeneity of glycoprotein structures they can produce. Here we provide an overview of the emerging field of synthetic glycobiology, the application of synthetic biology tools and design principles to better understand and engineer glycosylation. Specifically, we focus on how the biosynthetic and analytical tools of synthetic biology have been used to redesign glycosylation systems to obtain defined glycosylation structures on proteins for diverse applications in medicine, materials, and diagnostics. We review the key biological parts available to synthetic biologists interested in engineering glycoproteins to solve compelling problems in glycoscience, describe recent efforts to construct synthetic glycoprotein synthesis systems, and outline exemplary applications as well as new opportunities in this emerging space.

Modeling Glycan Processing Reveals Golgi-Enzyme Homeostasis upon Trafficking Defects and Cellular Differentiation

The decoration of proteins by carbohydrates is essential for eukaryotic life yet heterogeneous due to a lack of biosynthetic templates. This complex carbohydrate mixture—the glycan profile—is generated in the compartmentalized Golgi, in which level and localization of glycosylation enzymes are key determinants. Here, we develop and validate a computational model for glycan biosynthesis to probe how the biosynthetic machinery creates different glycan profiles. We combined stochastic modeling with Bayesian fitting that enables rigorous comparison to experimental data despite starting with uncertain initial parameters. This is an important development in the field of glycan modeling, which revealed biological insights about the glycosylation machinery in altered cellular states. We experimentally validated changes in N-linked glycan-modifying enzymes in cells with perturbed intra-Golgi-enzyme sorting and the predicted glycan-branching activity during osteogenesis. Our model can provide detailed information on altered biosynthetic paths, with potential for advancing treatments for glycosylation-related diseases and glyco-engineering of cells.

•

Developed a stochastic model of N-glycosylation coupled with Bayesian fitting

•

Validated predicted changes of Golgi organization in trafficking mutants

•

Model pinpointed functionally relevant glycan alterations in osteogenesis

How glycosylation affects glycosylation: the role of N-glycans in glycosyltransferase activity

N-glycosylation is one of the most important posttranslational modifications of proteins. It plays important roles in the biogenesis and functions of proteins by influencing their folding, intracellular localization, stability and solubility. N-glycans are synthesized by glycosyltransferases, a complex group of ubiquitous enzymes that occur in most kingdoms of life. A growing body of evidence shows that N-glycans may influence processing and functions of glycosyltransferases, including their secretion, stability and substrate/acceptor affinity. Changes in these properties may have a profound impact on glycosyltransferase activity. Indeed, some glycosyltransferases have to be glycosylated themselves for full activity. N-glycans and glycosyltransferases play roles in the pathogenesis of many diseases (including cancers), so studies on glycosyltransferases may contribute to the development of new therapy methods and novel glycoengineered enzymes with improved properties. In this review, we focus on the role of N-glycosylation in the activity of glycosyltransferases and attempt to summarize all available data about this phenomenon.

Sequential Glycosylation of Proteins with Substrate-Specific N-Glycosyltransferases

Protein glycosylation is a common post-translational modification that influences the functions and properties of proteins. Despite advances in methods to produce defined glycoproteins by chemoenzymatic elaboration of monosaccharides, the understanding and engineering of glycoproteins remain challenging, in part, due to the difficulty of site-specifically controlling glycosylation at each of several positions within a protein. Here, we address this limitation by discovering and exploiting the unique, conditionally orthogonal peptide acceptor specificities of N-glycosyltransferases (NGTs). We used cell-free protein synthesis and mass spectrometry of self-assembled monolayers to rapidly screen 41 putative NGTs and rigorously characterize the unique acceptor sequence preferences of four NGT variants using 1254 acceptor peptides and 8306 reaction conditions. We then used the optimized NGT-acceptor sequence pairs to sequentially install monosaccharides at four sites within one target protein. This strategy to site-specifically control the installation of N-linked monosaccharides for elaboration to a variety of functional N-glycans overcomes a major limitation in synthesizing defined glycoproteins for research and therapeutic applications.

Revisiting the Language of Glycoscience: Readers, Writers and Erasers in Carbohydrate Biochemistry

The roles of carbohydrates in nature are many and varied. However, the lack of template encoding in glycoscience distances carbohydrate structure, and hence function, from gene sequence. This challenging situation is compounded by descriptors of carbohydrate structure and function that have tended to emphasise their complexity. Herein, we suggest that revising the language of glycoscience could make interdisciplinary discourse more accessible to all interested parties.

Translation of genome to glycome: role of the Golgi apparatus

Glycans are one of the four biopolymers of the cell and they play important roles in cellular and organismal physiology. They consist of both linear and branched structures and are synthesized in a nontemplated manner in the secretory pathway of mammalian cells with the Golgi apparatus playing a key role in the process. In spite of the absence of a template, the glycans synthesized by a cell are not a random collection of possible glycan structures but a distribution of specific glycans in defined quantities that is unique to each cell type (Cell type here refers to distinct cell forms present in an organism that can be distinguished based on morphological, phenotypic and/or molecular criteria.) While information to produce cell type-specific glycans is encoded in the genome, how this information is translated into cell type-specific glycome (Glycome refers to the quantitative distribution of all glycan structures present in a given cell type.) is not completely understood. We summarize here the factors that are known to influence the fidelity of glycan biosynthesis and integrate them into known glycosylation pathways so as to rationalize the translation of genetic information to cell type-specific glycome.

Generation of site-distinct N-glycan variants for in vitro bioactivity testing

Glycosylation, a critical product quality attribute, may affect the efficacy and safety of therapeutic proteins in vivo. Chinese hamster ovary fed-batch cell culture batches yielded consistent glycoprofiles of a Fc-fusion antibody comprizing three different N-glycosylation sites. By adding media supplements at specific concentrations in cell culture and applying enzymatic glycoengineering, a diverse N-glycan variant population was generated, including high mannose, afucosylated, fucosylated, agalactosylated, galactosylated, asialylated, and sialylated forms. Site-specific glycosylation profiles were elucidated by glycopeptide mapping and the effect of the glycosylation variants on the FcγRIIIa receptor binding affinity and the biological activity (cell-based and surface plasmon resonance) was assessed. The two fusion body glycosylation sites were characterized by a high degree of sialic acid, more complex N-glycan structures, a higher degree of antennarity, and a site-specific behavior in the presence of a media supplement. On the other hand, the media supplements affected the Fc-site glycosylation heterogeneity similarly to the various studies described in the literature with classical monoclonal antibodies. Enzymatic glycoengineering solely managed to generate high levels of galactosylation at the fusion body sites. Variants with low core fucosylation, and to a lower extent, high mannose glycans exhibited increased FcγRIIIa receptor binding affinity. All N-glycan variants exhibited weak effects on the biological activity of the fusion body. Both media supplementation and enzymatic glycoengineering are suitable to generate sufficient diversity to assess the effect of glycostructures on the biological activity.

The contributions of individual galactosyltransferases to protein specific N-glycan processing in Chinese Hamster Ovary cells

Galactosylation as part of N-glycan processing is conducted by a set of beta-1,4-galactosyltransferases (B4GALTs), with B4GALT1 as the dominant isoenzyme for this reaction. Nevertheless, the exact contributions of this key-player as well as of the other isoenzymes involved in N-glycosylation, B4GALT2, B4GALT3 and B4GALT4, have not been studied in-depth. To increase the understanding of the protein- and site-specific activities of individual galactosyltransferases in Chinese Hamster Ovary cells, a panel of triple deletion cell lines was generated that expressed only one isoform of B4GALT each. Two model proteins were selected for this study to cover a large spectrum of possible N-glycan structures: erythropoietin and deamine-oxidase. They were expressed as Fc-fusion constructs (EPO-Fc and Fc-DAO) and their N-glycan processing status was analyzed by site-specific mass spectrometry. The sole activity of B4GALT1 resulted in a decrease of 15-21 % of fully galactosylated structures for erythropoietin, emphasizing the involvement of other isoenzymes. Interestingly, the contributions of B4GALT2 and B4GALT3 differed for the two model proteins. Unexpectedly, removal of galactosyltransferases influenced the overall process of N-glycan maturation, with the result of a higher occurrence of poorly processed oligosaccharides. In the context of high productivity cell lines, which can push N-glycan maturation towards incomplete galactosylation, galactosyltransferases are potential targets to ensure stable product quality. In view of our results, specifically engineered "designer" cell lines may be required for different proteins.