Glyco-Mapper: A Chinese hamster ovary (CHO) genome-specific glycosylation prediction tool

In silico simulation of glycosylation and related pathways

Glycans participate in a vast number of recognition systems in diverse organisms in health and in disease. However, glycans cannot be sequenced because there is no sequencer technology that can fully characterize them. There is no "template" for replicating glycans as there are for amino acids and nucleic acids. Instead, glycans are synthesized by a complicated orchestration of multitudes of glycosyltransferases and glycosidases. Thus glycans can vary greatly in structure, but they are not genetically reproducible and are usually isolated in minute amounts. To characterize (sequence) the glycome (defined as the glycans in a particular organism, tissue, cell, or protein), glycosylation pathway prediction using in silico methods based on glycogene expression data, and glycosylation simulations have been attempted. Since many of the mammalian glycogenes have been identified and cloned, it has become possible to predict the glycan biosynthesis pathway in these systems. By then incorporating systems biology and bioprocessing technologies to these pathway models, given the right enzymatic parameters including enzyme and substrate concentrations and kinetic reaction parameters, it is possible to predict the potentially synthesized glycans in the pathway. This review presents information on the data resources that are currently available to enable in silico simulations of glycosylation and related pathways. Then some of the software tools that have been developed in the past to simulate and analyze glycosylation pathways will be described, followed by a summary and vision for the future developments and research directions in this area.

Cell-culture process optimization via model-based predictions of metabolism and protein glycosylation

In order to meet the rising demand for biologics and become competitive on the developing biosimilar market, there is a need for process intensification of biomanufacturing processes. Process development of biologics has historically relied on extensive experimentation to develop and optimize biopharmaceutical manufacturing. Experimentation to optimize media formulations, feeding schedules, bioreactor operations and bioreactor scale up is expensive, labor intensive and time consuming. Mathematical modeling frameworks have the potential to enable process intensification while reducing the experimental burden. This review focuses on mathematical modeling of cellular metabolism and N-linked glycosylation as applied to upstream manufacturing of biologics. We review developments in the field of modeling cellular metabolism of mammalian cells using kinetic and stoichiometric modeling frameworks along with their applications to simulate, optimize and improve mechanistic understanding of the process. Interest in modeling N-linked glycosylation has led to the creation of various types of parametric and non-parametric models. Most published studies on mammalian cell metabolism have performed experiments in shake flasks where the pH and dissolved oxygen cannot be controlled. Efforts to understand and model the effect of bioreactor-specific parameters such as pH, dissolved oxygen, temperature, and bioreactor heterogeneity are critically reviewed. Most modeling efforts have focused on the Chinese Hamster Ovary (CHO) cells, which are most commonly used to produce monoclonal antibodies (mAbs). However, these modeling approaches can be generalized and applied to any mammalian cell-based manufacturing platform. Current and potential future applications of these models for Vero cell-based vaccine manufacturing, CAR-T cell therapies, and viral vector manufacturing are also discussed. We offer specific recommendations for improving the applicability of these models to industrially relevant processes.

CHOmpact: A reduced metabolic model of Chinese hamster ovary cells with enhanced interpretability

Metabolic modeling has emerged as a key tool for the characterization of biopharmaceutical cell culture processes. Metabolic models have also been instrumental in identifying genetic engineering targets and developing feeding strategies that optimize the growth and productivity of Chinese hamster ovary (CHO) cells. Despite their success, metabolic models of CHO cells still present considerable challenges. Genome-scale metabolic models (GeMs) of CHO cells are very large (>6000 reactions) and are difficult to constrain to yield physiologically consistent flux distributions. The large scale of GeMs also makes the interpretation of their outputs difficult. To address these challenges, we have developed CHOmpact, a reduced metabolic network that encompasses 101 metabolites linked through 144 reactions. Our compact reaction network allows us to deploy robust, nonlinear optimization and ensure that the computed flux distributions are physiologically consistent. Furthermore, our CHOmpact model delivers enhanced interpretability of simulation results and has allowed us to identify the mechanisms governing shifts in the anaplerotic consumption of asparagine and glutamate as well as an important mechanism of ammonia detoxification within mitochondria. CHOmpact, thus, addresses key challenges of large-scale metabolic models and will serve as a platform to develop dynamic metabolic models for the control and optimization of biopharmaceutical cell culture processes.

Preparing glycomics data for robust statistical analysis with GlyCompareCT

GlyCompareCT is a portable command-line tool to facilitate downstream glycomic data analyses, by addressing data inherent sparsity and non-independence. Inputting glycan abundances, users can run GlyCompareCT with one line of code to obtain the abundances of a minimal substructure set, named glycomotif, thereby quantifying hidden biosynthetic relationships between measured glycans. Optional parameters tuning and annotation are supported for personal preference. For complete details on the use and execution of this protocol, please refer to Bao et al. (2021).1.

CHOGlycoNET: Comprehensive glycosylation reaction network for CHO cells

Chinese hamster ovary (CHO) cells are extensively used for the production of glycoprotein therapeutics proteins, for which N-linked glycans are a critical quality attribute due to their influence on activity and immunogenicity. Manipulation of protein glycosylation is commonly achieved through cell or process engineering, which are often guided by mathematical models. However, each study considers a unique glycosylation reaction network that is tailored around the cell line and product at hand. Herein, we use 200 glycan datasets for both recombinantly produced and native proteins from different CHO cell lines to reconstruct a comprehensive reaction network, CHOGlycoNET, based on the individual minimal reaction networks describing each dataset. CHOGlycoNET is used to investigate the distribution of mannosidase and glycosyltransferase enzymes in the Golgi apparatus and identify key network reactions using machine learning and dimensionality reduction techniques. CHOGlycoNET can be used for accelerating glycomodel development and predicting the effect of glycoengineering strategies. Finally, CHOGlycoNET is wrapped in a SBML file to be used as a standalone model or in combination with CHO cell genome scale models.

Analysis of carbohydrates and glycoconjugates by matrix-assisted laser desorption/ionization mass spectrometry: An update for 2017-2018

This review is the tenth update of the original article published in 1999 on the application of matrix-assisted laser desorption/ionization mass spectrometry (MALDI) mass spectrometry to the analysis of carbohydrates and glycoconjugates and brings coverage of the literature to the end of 2018. Also included are papers that describe methods appropriate to glycan and glycoprotein analysis by MALDI, such as sample preparation techniques, even though the ionization method is not MALDI. Topics covered in the first part of the review include general aspects such as theory of the MALDI process, new methods, matrices, derivatization, MALDI imaging, fragmentation and the use of arrays. The second part of the review is devoted to applications to various structural types such as oligo- and poly-saccharides, glycoproteins, glycolipids, glycosides, and biopharmaceuticals. Most of the applications are presented in tabular form. The third part of the review covers medical and industrial applications of the technique, studies of enzyme reactions, and applications to chemical synthesis. The reported work shows increasing use of combined new techniques such as ion mobility and highlights the impact that MALDI imaging is having across a range of diciplines. MALDI is still an ideal technique for carbohydrate analysis and advancements in the technique and the range of applications continue steady progress.

Advanced strategies in glycosylation prediction and control during biopharmaceutical development: Avenues toward Industry 4.0

Glycosylation has been shown to define the safety and efficacy of biopharmaceuticals, thus classified as a critical quality attribute. However, controlling glycan heterogeneity has always been a major challenge owing to the multi-variate factors that govern the glycosylation process. Conventional approaches for controlling glycosylation such as gene editing and metabolic control have succeeded in obtaining desired glycan profiles in accordance with the Quality by Design paradigm. Nonetheless, the development of smart algorithms and omics-enabled complete cell characterization have made it possible to predict glycan profiles beforehand, and manipulate process variables accordingly. This review thus discusses the various approaches available for control and prediction of glycosylation in biopharmaceuticals. Further, the futuristic goal of integrating such technologies is discussed in order to attain an automated and digitized continuous bioprocess for control of glycosylation. Given, control of a process as complex as glycosylation requires intense monitoring intervention, we examine the current technologies that enable automation. Finally, we discuss the challenges and the technological gap that currently limits incorporation of an automated process in routine bio-manufacturing, with a glimpse into the economic bearing. This article is protected by copyright. All rights reserved.

Correcting for sparsity and interdependence in glycomics by accounting for glycan biosynthesis

Glycans are fundamental cellular building blocks, involved in many organismal functions. Advances in glycomics are elucidating the essential roles of glycans. Still, it remains challenging to properly analyze large glycomics datasets, since the abundance of each glycan is dependent on many other glycans that share many intermediate biosynthetic steps. Furthermore, the overlap of measured glycans can be low across samples. We address these challenges with GlyCompare, a glycomic data analysis approach that accounts for shared biosynthetic steps for all measured glycans to correct for sparsity and non-independence in glycomics, which enables direct comparison of different glycoprofiles and increases statistical power. Using GlyCompare, we study diverse N-glycan profiles from glycoengineered erythropoietin. We obtain biologically meaningful clustering of mutant cell glycoprofiles and identify knockout-specific effects of fucosyltransferase mutants on tetra-antennary structures. We further analyze human milk oligosaccharide profiles and find mother’s fucosyltransferase-dependent secretor-status indirectly impact the sialylation. Finally, we apply our method on mucin-type O-glycans, gangliosides, and site-specific compositional glycosylation data to reveal tissues and disease-specific glycan presentations. Our substructure-oriented approach will enable researchers to take full advantage of the growing power and size of glycomics data.

Glycomics can uncover important molecular changes but measured glycans are highly interconnected and incompatible with common statistical methods, introducing pitfalls during analysis. Here, the authors develop an approach to identify glycan dependencies across samples to facilitate comparative glycomics.

An integrated modular framework for modeling the effect of ammonium on the sialylation process of monoclonal antibodies produced by CHO cells

BACKGROUND

Monoclonal antibodies (mABs) have emerged as one of the most important therapeutic recombinant proteins in the pharmaceutical industry. Their immunogenicity and therapeutic efficacy are influenced by post-translational modifications, specifically the glycosylation process. Bioprocess conditions can influence the intracellular process of glycosylation. Among all the process conditions that have been recognized to affect the mAB glycoforms, the detailed mechanism underlying how ammonium could perturb glycosylation remains to be fully understood. It was shown that ammonium induces heterogeneity in protein glycosylation by altering the sialic acid content of glycoproteins. Hence, understanding this mechanism would aid pharmaceutical manufacturers to ensure consistent protein glycosylation.

METHODS

Three different mechanisms have been proposed to explain how ammonium influences the sialylation process. In the first, the inhibition of CMP-sialic acid transporter, which transports CMP-sialic acid (sialylation substrate) into the Golgi, by an increase in UDP-GlcNAc content that is brought about by the augmented incorporation of ammonium into glucosamine formation. In the second, ammonia diffuses into the Golgi and raises its pH, thereby decreasing the sialyltransferase enzyme activity. In the third, the reduction of sialyltransferase enzyme expression level in the presence of ammonium. We employed these mechanisms in a novel integrated modular platform to link dynamic alteration in mAB sialylation process with extracellular ammonium concentration to elucidate how ammonium alters the sialic acid content of glycoproteins.

RESULTS

Our results show that the sialylation reaction rate is insensitive to the first mechanism. At low ammonium concentration, the second mechanism is the controlling mechanism in mAB sialylation and by increasing the ammonium level (< 8 mM) the third mechanism becomes the controlling mechanism. At higher ammonium concentrations (> 8 mM) the second mechanism becomes predominant again.

CONCLUSION

The presented model in this study provides a connection between extracellular ammonium and the monoclonal antibody sialylation process. This computational tool could help scientists to develop and formulate cell culture media. The model illustrated here can assist the researchers to select culture media that ensure consistent mAB sialylation.

Systems glycobiology for discovering drug targets, biomarkers, and rational designs for glyco-immunotherapy

Cancer immunotherapy has revolutionized treatment and led to an unprecedented wave of immuno-oncology research during the past two decades. In 2018, two pioneer immunotherapy innovators, Tasuku Honjo and James P. Allison, were awarded the Nobel Prize for their landmark cancer immunotherapy work regarding “cancer therapy by inhibition of negative immune regulation” –CTLA4 and PD-1 immune checkpoints. However, the challenge in the coming decade is to develop cancer immunotherapies that can more consistently treat various patients and cancer types. Overcoming this challenge requires a systemic understanding of the underlying interactions between immune cells, tumor cells, and immunotherapeutics. The role of aberrant glycosylation in this process, and how it influences tumor immunity and immunotherapy is beginning to emerge. Herein, we review current knowledge of miRNA-mediated regulatory mechanisms of glycosylation machinery, and how these carbohydrate moieties impact immune cell and tumor cell interactions. We discuss these insights in the context of clinical findings and provide an outlook on modulating the regulation of glycosylation to offer new therapeutic opportunities. Finally, in the coming age of systems glycobiology, we highlight how emerging technologies in systems glycobiology are enabling deeper insights into cancer immuno-oncology, helping identify novel drug targets and key biomarkers of cancer, and facilitating the rational design of glyco-immunotherapies. These hold great promise clinically in the immuno-oncology field.

Improvement of the efficiency and quality in developing a new CHO host cell line

Chinese hamster ovary (CHO) cells are a ubiquitous tool for industrial therapeutic recombinant protein production. However, consistently generating high-producing clones remains a major challenge during the cell line development process. The glutamine synthetase (GS) and dihydrofolate reductase (DHFR) selection systems are commonly used CHO expression platforms based on controlling the balance of expression between the transgenic and endogenous GS or DHFR genes. Since the expression of the endogenous selection gene in CHO hosts can interfere with selection, generating a corresponding null CHO cell line is required to improve selection stringency, productivity, and stability. However, the efficiency of generating bi-allelic genetic knockouts using conventional protocols is very low (<5%). This significantly affects clone screening efficiency and reduces the chance of identifying robust knockout host cell lines. In this study, we use the GS expression system as an example to improve the genome editing process with zinc finger nucleases (ZFNs), resulting in improved GS-knockout efficiency of up to 46.8%. Furthermore, we demonstrate a process capable of enriching knockout CHO hosts with robust bioprocess traits. This integrated host development process yields a larger number of GS-knockout hosts with desired growth and recombinant protein expression characteristics.

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.

Global mapping of glycosylation pathways in human-derived cells

Glycans are one of the fundamental classes of macromolecules and are involved in a broad range of biological phenomena. A large variety of glycan structures can be synthesized depending on tissue or cell types and environmental changes. Here, we developed a comprehensive glycosylation mapping tool, termed GlycoMaple, to visualize and estimate glycan structures based on gene expression. We informatically selected 950 genes involved in glycosylation and its regulation. Expression profiles of these genes were mapped onto global glycan metabolic pathways to predict glycan structures, which were confirmed using glycomic analyses. Based on the predictions of N-glycan processing, we constructed 40 knockout HEK293 cell lines and analyzed the effects of gene knockout on glycan structures. Finally, the glycan structures of 64 cell lines, 37 tissues, and primary colon tumor tissues were estimated and compared using publicly available databases. Our systematic approach can accelerate glycan analyses and engineering in mammalian cells.

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.

N-Glycosylation of IgG and IgG-Like Recombinant Therapeutic Proteins: Why Is It Important and How Can We Control It?

Regulatory bodies worldwide consider N-glycosylation to be a critical quality attribute for immunoglobulin G (IgG) and IgG-like therapeutics. This consideration is due to the importance of posttranslational modifications in determining the efficacy, safety, and pharmacokinetic properties of biologics. Given its critical role in protein therapeutic production, we review N-glycosylation beginning with an overview of the myriad interactions of N-glycans with other biological factors. We examine the mechanism and drivers for N-glycosylation during biotherapeutic production and the several competing factors that impact glycan formation, including the abundance of precursor nucleotide sugars, transporters, glycosidases, glycosyltransferases, and process conditions. We explore the role of these factors with a focus on the analytical approaches used to characterize glycosylation and associated processes, followed by the current state of advanced glycosylation modeling techniques. This combination of disciplines allows for a deeper understanding of N-glycosylation and will lead to more rational glycan control.

Harnessing the potential of artificial neural networks for predicting protein glycosylation

Kinetic models offer incomparable insight on cellular mechanisms controlling protein glycosylation. However, their ability to reproduce site-specific glycoform distributions depends on accurate estimation of a large number of protein-specific kinetic parameters and prior knowledge of enzyme and transport protein levels in the Golgi membrane. Herein we propose an artificial neural network (ANN) for protein glycosylation and apply this to four recombinant glycoproteins produced in Chinese hamster ovary (CHO) cells, two monoclonal antibodies and two fusion proteins. We demonstrate that the ANN model accurately predicts site-specific glycoform distributions of up to eighteen glycan species with an average absolute error of 1.1%, correctly reproducing the effect of metabolic perturbations as part of a hybrid, kinetic/ANN, glycosylation model (HyGlycoM), as well as the impact of manganese supplementation and glycosyltransferase knock out experiments as a stand-alone machine learning algorithm. These results showcase the potential of machine learning and hybrid approaches for rapidly developing performance-driven models of protein glycosylation.

A Markov model of glycosylation elucidates isozyme specificity and glycosyltransferase interactions for glycoengineering

Glycosylated biopharmaceuticals are important in the global pharmaceutical market. Despite the importance of their glycan structures, our limited knowledge of the glycosylation machinery still hinders controllability of this critical quality attribute. To facilitate discovery of glycosyltransferase specificity and predict glycoengineering efforts, here we extend the approach to model N-linked protein glycosylation as a Markov process. Our model leverages putative glycosyltransferase (GT) specificity to define the biosynthetic pathways for all measured glycans, and the Markov chain modelling is used to learn glycosyltransferase isoform activities and predict glycosylation following glycosyltransferase knock-in/knockout. We apply our methodology to four different glycoengineered therapeutics (i.e., Rituximab, erythropoietin, Enbrel, and alpha-1 antitrypsin) produced in CHO cells. Our model accurately predicted N-linked glycosylation following glycoengineering and further quantified the impact of glycosyltransferase mutations on reactions catalyzed by other glycosyltransferases. By applying these learned GT-GT interaction rules identified from single glycosyltransferase mutants, our model further predicts the outcome of multi-gene glycosyltransferase mutations on the diverse biotherapeutics. Thus, this modeling approach enables rational glycoengineering and the elucidation of relationships between glycosyltransferases, thereby facilitating biopharmaceutical research and aiding the broader study of glycosylation to elucidate the genetic basis of complex changes in glycosylation.

Model-based optimization of antibody galactosylation in CHO cell culture

Exerting control over the glycan moieties of antibody therapeutics is highly desirable from a product safety and batch-to-batch consistency perspective. Strategies to improve antibody productivity may compromise quality, while interventions for improving glycoform distribution can adversely affect cell growth and productivity. Process design therefore needs to consider the trade-off between preserving cellular health and productivity while enhancing antibody quality. In this work, we present a modeling platform that quantifies the impact of glycosylation precursor feeding - specifically that of galactose and uridine - on cellular growth, metabolism as well as antibody productivity and glycoform distribution. The platform has been parameterized using an initial training data set yielding an accuracy of ±5% with respect to glycoform distribution. It was then used to design an optimized feeding strategy that enhances the final concentration of galactosylated antibody in the supernatant by over 90% compared with the control without compromising the integral of viable cell density or final antibody titer. This work supports the implementation of Quality by Design towards higher-performing bioprocesses.