A mathematical model of N-linked glycoform biosynthesis

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.

Improving recombinant protein production by yeast through genome-scale modeling using proteome constraints

Eukaryotic cells are used as cell factories to produce and secrete multitudes of recombinant pharmaceutical proteins, including several of the current top-selling drugs. Due to the essential role and complexity of the secretory pathway, improvement for recombinant protein production through metabolic engineering has traditionally been relatively ad-hoc; and a more systematic approach is required to generate novel design principles. Here, we present the proteome-constrained genome-scale protein secretory model of yeast Saccharomyces cerevisiae (pcSecYeast), which enables us to simulate and explain phenotypes caused by limited secretory capacity. We further apply the pcSecYeast model to predict overexpression targets for the production of several recombinant proteins. We experimentally validate many of the predicted targets for α-amylase production to demonstrate pcSecYeast application as a computational tool in guiding yeast engineering and improving recombinant protein production.

Due to the complexity of the protein secretory pathway, strategy suitable for the production of a certain recombination protein cannot be generalized. Here, the authors construct a proteome-constrained genome-scale protein secretory model for yeast and show its application in the production of different misfolded or recombinant proteins.

Computational Modeling of O-Linked Glycan Biosynthesis in CHO Cells

Glycan biosynthesis simulation research has progressed remarkably since 1997, when the first mathematical model for N-glycan biosynthesis was proposed. An O-glycan model has also been developed to predict O-glycan biosynthesis pathways in both forward and reverse directions. In this work, we started with a set of O-glycan profiles of CHO cells transiently transfected with various combinations of glycosyltransferases. The aim was to develop a model that encapsulated all the enzymes in the CHO transfected cell lines. Due to computational power restrictions, we were forced to focus on a smaller set of glycan profiles, where we were able to propose an optimized set of kinetics parameters for each enzyme in the model. Using this optimized model we showed that the abundance of more processed glycans could be simulated compared to observed abundance, while predicting the abundance of glycans earlier in the pathway was less accurate. The data generated show that for the accurate prediction of O-linked glycosylation, additional factors need to be incorporated into the model to better reflect the experimental conditions.

Glycan processing in the Golgi: optimal information coding and constraints on cisternal number and enzyme specificity

Many proteins that undergo sequential enzymatic modification in the Golgi cisternae are displayed at the plasma membrane as cell identity markers. The modified proteins, called glycans, represent a molecular code. The fidelity of this glycan code is measured by how accurately the glycan synthesis machinery realises the desired target glycan distribution for a particular cell type and niche. In this paper, we construct a simplified chemical synthesis model to quantitatively analyse the tradeoffs between the number of cisternae, and the number and specificity of enzymes, required to synthesize a prescribed target glycan distribution of a certain complexity to within a given fidelity. We find that to synthesize complex distributions, such as those observed in real cells, one needs to have multiple cisternae and precise enzyme partitioning in the Golgi. Additionally, for fixed number of enzymes and cisternae, there is an optimal level of specificity (promiscuity) of enzymes that achieves the target distribution with high fidelity. The geometry of the fidelity landscape in the multidimensional space of the number and specificity of enzymes, inter-cisternal transfer rates, and number of cisternae, provides a measure for robustness and identifies stiff and sloppy directions. Our results show how the complexity of the target glycan distribution and number of glycosylation enzymes places functional constraints on the Golgi cisternal number and enzyme specificity.

Computational Modeling of Glycan Processing in the Golgi for Investigating Changes in the Arrangements of Biosynthetic Enzymes

Modeling glycan biosynthesis is becoming increasingly important due to the far-reaching implications that glycosylation can exhibit, from pathologies to biopharmaceutical manufacturing. Here we describe a stochastic simulation approach, to overcome the deterministic nature of previous models, that aims to simulate the action of glycan modifying enzymes to produce a glycan profile. This is then coupled with an approximate Bayesian computation methodology to systematically fit to empirical data in order to determine which set of parameters adequately describes the organization of enzymes within the Golgi. The model is described in detail along with a proof of concept and therapeutic applications.

Functions of Glycosylation and Related Web Resources for Its Prediction

Glycosylation involves the attachment of carbohydrate sugar chains, or glycans, onto an amino acid residue of a protein. These glycans are often branched structures and serve to modulate the function of proteins. Glycans are synthesized through a complex process of enzymatic reactions that occur in the Golgi apparatus in mammalian systems. Because there is currently no sequencer for glycans, technologies such as mass spectrometry is used to characterize glycans in a biological sample to ascertain its glycome. This is a tedious process that requires high levels of expertise and equipment. Thus, the enzymes that work on glycans, called glycogenes or glycoenzymes, have been studied to better understand glycan function. With the development of glycan-related databases and a glycan repository, bioinformatics approaches have attempted to predict the glycosylation pathway and the glycosylation sites on proteins. This chapter introduces these methods and related Web resources for understanding glycan function.

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.

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.

Simulating the enzymes of ganglioside biosynthesis with Glycologue

Gangliosides are an important class of sialylated glycosphingolipids linked to ceramide that are a component of the mammalian cell surface, especially those of the central nervous system, where they function in intercellular recognition and communication. We describe an in silico method for determining the metabolic pathways leading to the most common gangliosides, based on the known enzymes of their biosynthesis. A network of 41 glycolipids is produced by the actions of the 10 enzymes included in the model. The different ganglioside nomenclature systems in common use are compared and a systematic variant of the widely used Svennerholm nomenclature is described. Knockouts of specific enzyme activities are used to simulate congenital defects in ganglioside biosynthesis, and altered ganglioside status in cancer, and the effects on network structure are predicted. The simulator is available at the Glycologue website, https://glycologue.org/.

Insights into Bioinformatic Applications for Glycosylation: Instigating an Awakening towards Applying Glycoinformatic Resources for Cancer Diagnosis and Therapy

Glycosylation plays a crucial role in various diseases and their etiology. This has led to a clear understanding on the functions of carbohydrates in cell communication, which eventually will result in novel therapeutic approaches for treatment of various disease. Glycomics has now become one among the top ten technologies that will change the future. The direct implication of glycosylation as a hallmark of cancer and for cancer therapy is well established. As in proteomics, where bioinformatics tools have led to revolutionary achievements, bioinformatics resources for glycosylation have improved its practical implication. Bioinformatics tools, algorithms and databases are a mandatory requirement to manage and successfully analyze large amount of glycobiological data generated from glycosylation studies. This review consolidates all the available tools and their applications in glycosylation research. The achievements made through the use of bioinformatics into glycosylation studies are also presented. The importance of glycosylation in cancer diagnosis and therapy is discussed and the gap in the application of widely available glyco-informatic tools for cancer research is highlighted. This review is expected to bring an awakening amongst glyco-informaticians as well as cancer biologists to bridge this gap, to exploit the available glyco-informatic tools for cancer.

Metabolic pathway engineering: Perspectives and applications

BACKGROUND

Metabolic engineering aims at contriving microbes as biocatalysts for enhanced and cost-effective production of countless secondary metabolites. These secondary metabolites can be treated as the resources of industrial chemicals, pharmaceuticals and fuels. Plants are also crucial targets for metabolic engineers to produce necessary secondary metabolites. Metabolic engineering of both microorganism and plants also contributes towards drug discovery. In order to implement advanced metabolic engineering techniques efficiently, metabolic engineers should have detailed knowledge about cell physiology and metabolism. Principle behind methodologies: Genome-scale mathematical models of integrated metabolic, signal transduction, gene regulatory and protein-protein interaction networks along with experimental validation can provide such knowledge in this context. Incorporation of omics data into these models is crucial in the case of drug discovery. Inverse metabolic engineering and metabolic control analysis (MCA) can help in developing such models. Artificial intelligence methodology can also be applied for efficient and accurate metabolic engineering.

CONCLUSION

In this review, we discuss, at the beginning, the perspectives of metabolic engineering and its application on microorganism and plant leading to drug discovery. At the end, we elaborate why inverse metabolic engineering and MCA are closely related to modern metabolic engineering. In addition, some crucial steps ensuring efficient and optimal metabolic engineering strategies have been discussed. Moreover, we explore the use of genomics data for the activation of silent metabolic clusters and how it can be integrated with metabolic engineering. Finally, we exhibit a few applications of artificial intelligence to metabolic engineering.

Cell-Free Synthetic Glycobiology: Designing and Engineering Glycomolecules Outside of Living Cells

Glycans and glycosylated biomolecules are directly involved in almost every biological process as well as the etiology of most major diseases. Hence, glycoscience knowledge is essential to efforts aimed at addressing fundamental challenges in understanding and improving human health, protecting the environment and enhancing energy security, and developing renewable and sustainable resources that can serve as the source of next-generation materials. While much progress has been made, there remains an urgent need for new tools that can overexpress structurally uniform glycans and glycoconjugates in the quantities needed for characterization and that can be used to mechanistically dissect the enzymatic reactions and multi-enzyme assembly lines that promote their construction. To address this technology gap, cell-free synthetic glycobiology has emerged as a simplified and highly modular framework to investigate, prototype, and engineer pathways for glycan biosynthesis and biomolecule glycosylation outside the confines of living cells. From nucleotide sugars to complex glycoproteins, we summarize here recent efforts that harness the power of cell-free approaches to design, build, test, and utilize glyco-enzyme reaction networks that produce desired glycomolecules in a predictable and controllable manner. We also highlight novel cell-free methods for shedding light on poorly understood aspects of diverse glycosylation processes and engineering these processes toward desired outcomes. Taken together, cell-free synthetic glycobiology represents a promising set of tools and techniques for accelerating basic glycoscience research (e.g., deciphering the "glycan code") and its application (e.g., biomanufacturing high-value glycomolecules on demand).

Protein Glycoengineering: An Approach for Improving Protein Properties

Natural proteins are an important source of therapeutic agents and industrial enzymes. While many of them have the potential to be used as highly effective medical treatments for a wide range of diseases or as catalysts for conversion of a range of molecules into important product types required by modern society, problems associated with poor biophysical and biological properties have limited their applications. Engineering proteins with reduced side-effects and/or improved biophysical and biological properties is therefore of great importance. As a common protein modification, glycosylation has the capacity to greatly influence these properties. Over the past three decades, research from many disciplines has established the importance of glycoengineering in overcoming the limitations of proteins. In this review, we will summarize the methods that have been used to glycoengineer proteins and briefly discuss some representative examples of these methods, with the goal of providing a general overview of this research area.

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

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 Systematic Review on the Implications of O-linked Glycan Branching and Truncating Enzymes on Cancer Progression and Metastasis

Glycosylation is the most commonly occurring post-translational modifications, and is believed to modify over 50% of all proteins. The process of glycan modification is directed by different glycosyltransferases, depending on the cell in which it is expressed. These small carbohydrate molecules consist of multiple glycan families that facilitate cell-cell interactions, protein interactions, and downstream signaling. An alteration of several types of O-glycan core structures have been implicated in multiple cancers, largely due to differential glycosyltransferase expression or activity. Consequently, aberrant O-linked glycosylation has been extensively demonstrated to affect biological function and protein integrity that directly result in cancer growth and progression of several diseases. Herein, we provide a comprehensive review of several initiating enzymes involved in the synthesis of O-linked glycosylation that significantly contribute to a number of different cancers.

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.

Glycosylation Flux Analysis of Immunoglobulin G in Chinese Hamster Ovary Perfusion Cell Culture

The terminal sugar molecules of the N-linked glycan attached to the fragment crystalizable (Fc) region is a critical quality attribute of therapeutic monoclonal antibodies (mAbs) such as immunoglobulin G (IgG). There exists naturally-occurring heterogeneity in the N-linked glycan structure of mAbs, and such heterogeneity has a significant influence on the clinical safety and efficacy of mAb drugs. We previously proposed a constraint-based modeling method called glycosylation flux analysis (GFA) to characterize the rates (fluxes) of intracellular glycosylation reactions. One contribution of this work is a significant improvement in the computational efficiency of the GFA, which is beneficial for analyzing large datasets. Another contribution of our study is the analysis of IgG glycosylation in continuous perfusion Chinese Hamster Ovary (CHO) cell cultures. The GFA of the perfusion cell culture data indicated that the dynamical changes of IgG glycan heterogeneity are mostly attributed to alterations in the galactosylation flux activity. By using a random forest regression analysis of the IgG galactosylation flux activity, we were further able to link the dynamics of galactosylation with two process parameters: cell-specific productivity of IgG and extracellular ammonia concentration. The characteristics of IgG galactosylation dynamics agree well with what we previously reported for fed-batch cultivations of the same CHO cell strain.

A mechanism for bistability in glycosylation

Glycosyltransferases are a class of enzymes that catalyse the posttranslational modification of proteins to produce a large number of glycoconjugate acceptors from a limited number of nucleotide-sugar donors. The products of one glycosyltransferase can be the substrates of several other enzymes, causing a combinatorial explosion in the number of possible glycan products. The kinetic behaviour of systems where multiple acceptor substrates compete for a single enzyme is presented, and the case in which high concentrations of an acceptor substrate are inhibitory as a result of abortive complex formation, is shown to result in non-Michaelian kinetics that can lead to bistability in an open system. A kinetic mechanism is proposed that is consistent with the available experimental evidence and provides a possible explanation for conflicting observations on the β-1,4-galactosyltransferases. Abrupt switching between steady states in networks of glycosyltransferase-catalysed reactions may account for the observed changes in glycosyl-epitopes in cancer cells.

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

Glyco-Mapper is a novel systems biology product quality prediction tool created using a new framework termed: Discretized Reaction Network Modeling using Fuzzy Parameters (DReaM-zyP). Within Glyco-Mapper, users fix the nutrient feed composition and the glycosylation reaction fluxes to fit the model glycoform to the reference experimental glycoform, enabling cell-line specific glycoform predictions as a result of cell engineering strategies. Glyco-Mapper accurately predicts glycoforms associated with genetic alterations that result in the appearance or disappearance of one or more glycans with an accuracy, sensitivity, and specificity of 96%, 85%, and 97%, respectively, for publications between 1999 and 2014. The modeled glycoforms span a large range of glycoform engineering strategies, including the altered expression of glycosylation, nucleotide sugar transport, and metabolism genes, as well as an altered nutrient feeding strategy. A glycoprotein-producing CHO cell line reference glycoform was modeled and a novel Glyco-Mapper prediction was experimentally confirmed with an accuracy and specificity of 95% and 98%, respectively. Glyco-Mapper is a product quality prediction tool that provides a streamlined way to design host cell line genomes to achieve specific product quality attributes.

Glycoengineering in CHO Cells: Advances in Systems Biology

For several decades, glycoprotein biologics have been successfully produced from Chinese hamster ovary (CHO) cells. The therapeutic efficacy and potency of glycoprotein biologics are often dictated by their post-translational modifications, particularly glycosylation, which unlike protein synthesis, is a non-templated process. Consequently, both native and recombinant glycoprotein production generate heterogeneous mixtures containing variable amounts of different glycoforms. Stability, potency, plasma half-life, and immunogenicity of the glycoprotein biologic are directly influenced by the glycoforms. Recently, CHO cells have also been explored for production of therapeutic glycosaminoglycans (e.g., heparin), which presents similar challenges as producing glycoproteins biologics. Approaches to controlling heterogeneity in CHO cells and directing the biosynthetic process toward desired glycoforms are not well understood. A systems biology approach combining different technologies is needed for complete understanding of the molecular processes accounting for this variability and to open up new venues in cell line development. In this review, we describe several advances in genetic manipulation, modeling, and glycan and glycoprotein analysis that together will provide new strategies for glycoengineering of CHO cells with desired or enhanced glycosylation capabilities.

Glycosylation flux analysis reveals dynamic changes of intracellular glycosylation flux distribution in Chinese hamster ovary fed-batch cultures

N-linked glycosylation of proteins has both functional and structural significance. Importantly, the glycan structure of a therapeutic protein influences its efficacy, pharmacokinetics, pharmacodynamics and immunogenicity. In this work, we developed glycosylation flux analysis (GFA) for predicting intracellular production and consumption rates (fluxes) of glycoforms, and applied this analysis to CHO fed-batch immunoglobulin G (IgG) production using two different media compositions, with and without additional manganese feeding. The GFA is based on a constraint-based modeling of the glycosylation network, employing a pseudo steady state assumption. While the glycosylation fluxes in the network are balanced at each time point, the GFA allows the fluxes to vary with time by way of two scaling factors: (1) an enzyme-specific factor that captures the temporal changes among glycosylation reactions catalysed by the same enzyme, and (2) the cell specific productivity factor that accounts for the dynamic changes in the IgG production rate. The GFA of the CHO fed-batch cultivations showed that regardless of the media composition, galactosylation fluxes decreased with the cultivation time more significantly than the other glycosylation reactions. Furthermore, the GFA showed that the addition of Mn, a cofactor of galactosyltransferase, has the effect of increasing the galactosylation fluxes but only during the beginning of the cultivation period. The results thus demonstrated the power of the GFA in delineating the dynamic alterations of the glycosylation fluxes by local (enzyme-specific) and global (cell specific productivity) factors.

Model-based analysis of N-glycosylation in Chinese hamster ovary cells

The Chinese hamster ovary (CHO) cell is the gold standard for manufacturing of glycosylated recombinant proteins for production of biotherapeutics. The similarity of its glycosylation patterns to the human versions enable the products of this cell line favorable pharmacokinetic properties and lower likelihood of causing immunogenic responses. Because glycan structures are the product of the concerted action of intracellular enzymes, it is difficult to predict a priori how the effects of genetic manipulations alter glycan structures of cells and therapeutic properties. For that reason, quantitative models able to predict glycosylation have emerged as promising tools to deal with the complexity of glycosylation processing. For example, an earlier version of the same model used in this study was used by others to successfully predict changes in enzyme activities that could produce a desired change in glycan structure. In this study we utilize an updated version of this model to provide a comprehensive analysis of N-glycosylation in ten Chinese hamster ovary (CHO) cell lines that include a wild type parent and nine mutants of CHO, through interpretation of previously published mass spectrometry data. The updated N-glycosylation mathematical model contains up to 50,605 glycan structures. Adjusting the enzyme activities in this model to match N-glycan mass spectra produces detailed predictions of the glycosylation process, enzyme activity profiles and complete glycosylation profiles of each of the cell lines. These profiles are consistent with biochemical and genetic data reported previously. The model-based results also predict glycosylation features of the cell lines not previously published, indicating more complex changes in glycosylation enzyme activities than just those resulting directly from gene mutations. The model predicts that the CHO cell lines possess regulatory mechanisms that allow them to adjust glycosylation enzyme activities to mitigate side effects of the primary loss or gain of glycosylation function known to exist in these mutant cell lines. Quantitative models of CHO cell glycosylation have the potential for predicting how glycoengineering manipulations might affect glycoform distributions to improve the therapeutic performance of glycoprotein products.

What can mathematical modelling say about CHO metabolism and protein glycosylation?

Chinese hamster ovary cells have been in the spotlight for process optimization in recent years, due to being the major, long established cell factory for the production of recombinant proteins. A deep, quantitative understanding of CHO metabolism and mechanisms involved in protein glycosylation has proven to be attainable through the development of high throughput technologies. Here we review the most notable accomplishments in the field of modelling CHO metabolism and protein glycosylation.

Fc gamma receptors: glycobiology and therapeutic prospects

Therapeutic antibodies hold great promise for the treatment of cancer and autoimmune diseases, and developments in antibody–drug conjugates and bispecific antibodies continue to enhance treatment options for patients. Immunoglobulin (Ig) G antibodies are proteins with complex modifications, which have a significant impact on their function. The most important of these modifications is glycosylation, the addition of conserved glycans to the antibody Fc region, which is critical for its interaction with the immune system and induction of effector activities such as antibody-dependent cell cytotoxicity, complement activation and phagocytosis. Communication of IgG antibodies with the immune system is controlled and mediated by Fc gamma receptors (FcγRs), membrane-bound proteins, which relay the information sensed and gathered by antibodies to the immune system. These receptors are also glycoproteins and provide a link between the innate and adaptive immune systems. Recent information suggests that this receptor glycan modification is also important for the interaction with antibodies and downstream immune response. In this study, the current knowledge on FcγR glycosylation is discussed, and some insight into its role and influence on the interaction properties with IgG, particularly in the context of biotherapeutics, is provided. For the purpose of this study, other Fc receptors such as FcαR, FcεR or FcRn are not discussed extensively, as IgG-based antibodies are currently the only therapeutic antibody-based products on the market. In addition, FcγRs as therapeutics and therapeutic targets are discussed, and insight into and comment on the therapeutic aspects of receptor glycosylation are provided.

A systematic framework to derive N-glycan biosynthesis process and the automated construction of glycosylation networks

Background

Abnormalities in glycan biosynthesis have been conclusively related to various diseases, whereas the complexity of the glycosylation process has impeded the quantitative analysis of biochemical experimental data for the identification of glycoforms contributing to disease. To overcome this limitation, the automatic construction of glycosylation reaction networks in silico is a critical step.

Results

In this paper, a framework K2014 is developed to automatically construct N-glycosylation networks in MATLAB with the involvement of the 27 most-known enzyme reaction rules of 22 enzymes, as an extension of previous model KB2005. A toolbox named Glycosylation Network Analysis Toolbox (GNAT) is applied to define network properties systematically, including linkages, stereochemical specificity and reaction conditions of enzymes. Our network shows a strong ability to predict a wider range of glycans produced by the enzymes encountered in the Golgi Apparatus in human cell expression systems.

Conclusions

Our results demonstrate a better understanding of the underlying glycosylation process and the potential of systems glycobiology tools for analyzing conventional biochemical or mass spectrometry-based experimental data quantitatively in a more realistic and practical way.

Controlling the time evolution of mAb N-linked glycosylation - Part II: Model-based predictions

N-linked glycosylation is known to be a crucial factor for the therapeutic efficacy and safety of monoclonal antibodies (mAbs) and many other glycoproteins. The nontemplate process of glycosylation is influenced by external factors which have to be tightly controlled during the manufacturing process. In order to describe and predict mAb N-linked glycosylation patterns in a CHO-S cell fed-batch process, an existing dynamic mathematical model has been refined and coupled to an unstructured metabolic model. High-throughput cell culture experiments carried out in miniaturized bioreactors in combination with intracellular measurements of nucleotide sugars were used to tune the parameter configuration of the coupled models as a function of extracellular pH, manganese and galactose addition. The proposed modeling framework is able to predict the time evolution of N-linked glycosylation patterns during a fed-batch process as a function of time as well as the manipulated variables. A constant and varying mAb N-linked glycosylation pattern throughout the culture were chosen to demonstrate the predictive capability of the modeling framework, which is able to quantify the interconnected influence of media components and cell culture conditions. Such a model-based evaluation of feeding regimes using high-throughput tools and mathematical models gives rise to a more rational way to control and design cell culture processes with defined glycosylation patterns. © 2016 American Institute of Chemical Engineers Biotechnol. Prog., 32:1135-1148, 2016.

Comprehensive analysis of the N-glycan biosynthetic pathway using bioinformatics to generate UniCorn: A theoretical N-glycan structure database

Glycan structures attached to proteins are comprised of diverse monosaccharide sequences and linkages that are produced from precursor nucleotide-sugars by a series of glycosyltransferases. Databases of these structures are an essential resource for the interpretation of analytical data and the development of bioinformatics tools. However, with no template to predict what structures are possible the human glycan structure databases are incomplete and rely heavily on the curation of published, experimentally determined, glycan structure data. In this work, a library of 45 human glycosyltransferases was used to generate a theoretical database of N-glycan structures comprised of 15 or less monosaccharide residues. Enzyme specificities were sourced from major online databases including Kyoto Encyclopedia of Genes and Genomes (KEGG) Glycan, Consortium for Functional Glycomics (CFG), Carbohydrate-Active enZymes (CAZy), GlycoGene DataBase (GGDB) and BRENDA. Based on the known activities, more than 1.1 million theoretical structures and 4.7 million synthetic reactions were generated and stored in our database called UniCorn. Furthermore, we analyzed the differences between the predicted glycan structures in UniCorn and those contained in UniCarbKB (www.unicarbkb.org), a database which stores experimentally described glycan structures reported in the literature, and demonstrate that UniCorn can be used to aid in the assignment of ambiguous structures whilst also serving as a discovery database.

N-Glycosylation Design and Control of Therapeutic Monoclonal Antibodies

The N-linked glycan profiles on recombinant monoclonal antibody therapeutics significantly affect antibody biological functions and are largely determined by host cell genotypes and culture conditions. A key step in bioprocess development for monoclonal antibodies (mAbs) involves optimization and control of N-glycan profiles. With pressure from pricing and biosimilars looming, more efficient and effective approaches are sought in the field of glycoengineering. Metabolic studies and mathematical modeling are two such approaches that optimize bioprocesses by better understanding and predicting glycosylation. In this review, we summarize a group of strategies currently used for glycan profile modulation and control. Metabolic analysis and mathematical modeling are then explored with an emphasis on how these two techniques can be utilized to advance glycoengineering.

A control strategy to investigate the relationship between specific productivity and high-mannose glycoforms in CHO cells

The integration of physiological knowledge into process control strategies is a cornerstone for the improvement of biopharmaceutical cell culture technologies. The present contribution investigates the applicability of specific productivity as a physiological control parameter in a cell culture process producing a monoclonal antibody (mAb) in CHO cells. In order to characterize cell physiology, the on-line oxygen uptake rate (OUR) was monitored and the time-resolved specific productivity was calculated as physiological parameters. This characterization enabled to identify the tight link between the deprivation of tyrosine and the decrease in cell respiration and in specific productivity. Subsequently, this link was used to control specific productivity by applying different feeding profiles. The maintenance of specific productivity at various levels enabled to identify a correlation between the rate of product formation and the relative abundance of high-mannose glycoforms. An increase in high mannose content was assumed to be the result of high specific productivity. Furthermore, the high mannose content as a function of cultivation pH and specific productivity was investigated in a design of experiment approach. This study demonstrated how physiological parameters could be used to understand interactions between process parameters, physiological parameters, and product quality attributes.

Dynamics of immature mAb glycoform secretion during CHO cell culture: An integrated modelling framework

Ensuring consistent glycosylation-associated quality of therapeutic monoclonal antibodies (mAbs) has become a priority in pharmaceutical bioprocessing given that the distribution and composition of the carbohydrates (glycans) bound to these molecules determines their therapeutic efficacy and immunogenicity. However, the interaction between bioprocess conditions, cellular metabolism and the intracellular process of glycosylation remains to be fully understood. To gain further insight into these interactions, we present a novel integrated modelling platform that links dynamic variations in mAb glycosylation with cellular secretory capacity. Two alternative mechanistic representations of how mAb specific productivity (qp ) influences glycosylation are compared. In the first, mAb glycosylation is modulated by the linear velocity with which secretory cargo traverses the Golgi apparatus. In the second, glycosylation is influenced by variations in Golgi volume. Within our modelling framework, both mechanisms accurately reproduce experimentally-observed dynamic changes in mAb glycosylation. In addition, an optimisation-based strategy has been developed to estimate the concentration of glycosylation enzymes required to minimise mAb glycoform variability. Our results suggest that the availability of glycosylation machinery relative to cellular secretory capacity may play a crucial role in mAb glycosylation. In the future, the modelling framework presented here may aid in selecting and engineering cell lines that ensure consistent mAb glycosylatio.

Macroscopic modeling of mammalian cell growth and metabolism

We review major modeling strategies and methods to understand and simulate the macroscopic behavior of mammalian cells. These strategies comprise two important steps: the first step is to identify stoichiometric relationships for the cultured cells connecting the extracellular inputs and outputs. In a second step, macroscopic kinetic models are introduced. These relationships together with bioreactor and metabolite balances provide a complete description of a system in the form of a set of differential equations. These can be used for the simulation of cell culture performance and further for optimization of production.

Integration of systems glycobiology with bioinformatics toolboxes, glycoinformatics resources, and glycoproteomics data

The glycome constitutes the entire complement of free carbohydrates and glycoconjugates expressed on whole cells or tissues. 'Systems Glycobiology' is an emerging discipline that aims to quantitatively describe and analyse the glycome. Here, instead of developing a detailed understanding of single biochemical processes, a combination of computational and experimental tools are used to seek an integrated or 'systems-level' view. This can explain how multiple biochemical reactions and transport processes interact with each other to control glycome biosynthesis and function. Computational methods in this field commonly build in silico reaction network models to describe experimental data derived from structural studies that measure cell-surface glycan distribution. While considerable progress has been made, several challenges remain due to the complex and heterogeneous nature of this post-translational modification. First, for the in silico models to be standardized and shared among laboratories, it is necessary to integrate glycan structure information and glycosylation-related enzyme definitions into the mathematical models. Second, as glycoinformatics resources grow, it would be attractive to utilize 'Big Data' stored in these repositories for model construction and validation. Third, while the technology for profiling the glycome at the whole-cell level has been standardized, there is a need to integrate mass spectrometry derived site-specific glycosylation data into the models. The current review discusses progress that is being made to resolve the above bottlenecks. The focus is on how computational models can bridge the gap between 'data' generated in wet-laboratory studies with 'knowledge' that can enhance our understanding of the glycome.

Galactosyltransferase 4 is a major control point for glycan branching in N-linked glycosylation

Protein N-glycosylation is a common post-translational modification that produces a complex array of branched glycan structures. The levels of branching, or antennarity, give rise to differential biological activities for single glycoproteins. However, the precise mechanism controlling the glycan branching and glycosylation network is unknown. Here, we constructed quantitative mathematical models of N-linked glycosylation that predicted new control points for glycan branching. Galactosyltransferase, which acts on N-acetylglucosamine residues, was unexpectedly found to control metabolic flux through the glycosylation pathway and the level of final antennarity of nascent protein produced in the Golgi network. To further investigate the biological consequences of glycan branching in nascent proteins, we glycoengineered a series of mammalian cells overexpressing human chorionic gonadotropin (hCG). We identified a mechanism in which galactosyltransferase 4 isoform regulated N-glycan branching on the nascent protein, subsequently controlling biological activity in an in vivo model of hCG activity. We found that galactosyltransferase 4 is a major control point for glycan branching decisions taken in the Golgi of the cell, which might ultimately control the biological activity of nascent glycoprotein.

Systems glycobiology for glycoengineering

Glycosylation serves essential functions on many proteins produced in biopharmaceutical manufacturing, making it mandatory to thoroughly consider its biogenesis during the production process. Glycoengineering efforts involve the rational design of glycosylation through adjustments in culturing conditions or genetic modifications. Computational models have been developed to aid this process, aiming to offer cheaper and faster alternatives to costly screening strategies. Recently, these models have been successfully utilized to predict glycosylation of products of industrial relevance. Furthermore, systems-level analyses of glycan diversity are elucidating deeper insights into the mechanisms underlying glycosylation. As computational models of glycosylation continue to be expanded, refined, and leveraged for detailed analysis of glycomics data, they will become invaluable resources for cell line development and glycoengineering.

A computational framework for the automated construction of glycosylation reaction networks

Glycosylation is among the most common and complex post-translational modifications identified to date. It proceeds through the catalytic action of multiple enzyme families that include the glycosyltransferases that add monosaccharides to growing glycans, and glycosidases which remove sugar residues to trim glycans. The expression level and specificity of these enzymes, in part, regulate the glycan distribution or glycome of specific cell/tissue systems. Currently, there is no systematic method to describe the enzymes and cellular reaction networks that catalyze glycosylation. To address this limitation, we present a streamlined machine-readable definition for the glycosylating enzymes and additional methodologies to construct and analyze glycosylation reaction networks. In this computational framework, the enzyme class is systematically designed to store detailed specificity data such as enzymatic functional group, linkage and substrate specificity. The new classes and their associated functions enable both single-reaction inference and automated full network reconstruction, when given a list of reactants and/or products along with the enzymes present in the system. In addition, graph theory is used to support functions that map the connectivity between two or more species in a network, and that generate subset models to identify rate-limiting steps regulating glycan biosynthesis. Finally, this framework allows the synthesis of biochemical reaction networks using mass spectrometry (MS) data. The features described above are illustrated using three case studies that examine: i) O-linked glycan biosynthesis during the construction of functional selectin-ligands; ii) automated N-linked glycosylation pathway construction; and iii) the handling and analysis of glycomics based MS data. Overall, the new computational framework enables automated glycosylation network model construction and analysis by integrating knowledge of glycan structure and enzyme biochemistry. All the implemented features are provided as part of the Glycosylation Network Analysis Toolbox (GNAT), an open-source, platform-independent, MATLAB based toolbox for studies of Systems Glycobiology.

Kinetic models in industrial biotechnology - Improving cell factory performance

An increasing number of industrial bioprocesses capitalize on living cells by using them as cell factories that convert sugars into chemicals. These processes range from the production of bulk chemicals in yeasts and bacteria to the synthesis of therapeutic proteins in mammalian cell lines. One of the tools in the continuous search for improved performance of such production systems is the development and application of mathematical models. To be of value for industrial biotechnology, mathematical models should be able to assist in the rational design of cell factory properties or in the production processes in which they are utilized. Kinetic models are particularly suitable towards this end because they are capable of representing the complex biochemistry of cells in a more complete way compared to most other types of models. They can, at least in principle, be used to in detail understand, predict, and evaluate the effects of adding, removing, or modifying molecular components of a cell factory and for supporting the design of the bioreactor or fermentation process. However, several challenges still remain before kinetic modeling will reach the degree of maturity required for routine application in industry. Here we review the current status of kinetic cell factory modeling. Emphasis is on modeling methodology concepts, including model network structure, kinetic rate expressions, parameter estimation, optimization methods, identifiability analysis, model reduction, and model validation, but several applications of kinetic models for the improvement of cell factories are also discussed.

Controllability analysis of protein glycosylation in CHO cells

To function as intended in vivo, a majority of biopharmaceuticals require specific glycan distributions. However, achieving a precise glycan distribution during manufacturing can be challenging because glycosylation is a non-template driven cellular process, with the potential for significant uncontrolled variability in glycan distributions. As important as the glycan distribution is to the end-use performance of biopharmaceuticals, to date, no strategy exists for controlling glycosylation on-line. However, before expending the significant amount of effort and expense required to develop and implement on-line control strategies to address the problem of glycosylation heterogeneity, it is imperative to assess first the extent to which the very complex process of glycosylation is controllable, thereby establishing what is theoretically achievable prior to any experimental attempts. In this work, we present a novel methodology for assessing the output controllability of glycosylation, a prototypical example of an extremely high-dimensional and very non-linear system. We first discuss a method for obtaining the process gain matrix for glycosylation that involves performing model simulations and data analysis systematically and judiciously according to a statistical design of experiments (DOE) scheme and then employing Analysis of Variance (ANOVA) to determine the elements of process gain matrix from the resulting simulation data. We then discuss how to use the resulting high-dimensional gain matrix to assess controllability. The utility of this method is demonstrated with a practical example where we assess the controllability of various classes of glycans and of specific glycoforms that are typically found in recombinant biologics produced with Chinese Hamster Ovary (CHO) cells. In addition to providing useful insight into the extent to which on-line glycosylation control is achievable in actual manufacturing processes, the results also have important implications for genetically engineering cell lines design for enhanced glycosylation controllability.