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

Mimicking CHO large-scale effects in the single multicompartment bioreactor: A new approach to access scale-up behavior

During the scale-up of biopharmaceutical production processes, insufficiently predictable performance losses may occur alongside gradients and heterogeneities. To overcome such performance losses, tools are required to explain, predict, and ultimately prohibit inconsistencies between laboratory and commercial scale. In this work, we performed CHO fed-batch cultivations in the single multicompartment bioreactor (SMCB), a new scale-down reactor system that offers new access to study large-scale heterogeneities in mammalian cell cultures. At volumetric power inputs of 20.4-1.5 W m-3, large-scale characteristics like long mixing times and dissolved oxygen (DO) heterogeneities were mimicked in the SMCB. Compared to a reference bioreactor (REFB) set-up, the conditions in the SMCB provoked an increase in lactate accumulation of up to 87%, an increased glucose uptake, and reduced viable cell concentrations in the stationary phase. All are characteristic for large-scale performance. The unique possibility to distinguish between the effects of changing power inputs and observed heterogeneities provided new insights into the potential reasons for altered product quality attributes. Apparently, the degree of galactosylation in the evaluated glycan patterns changed primarily due to the different power inputs rather than the provoked heterogeneities. The SMCB system could serve as a potent tool to provide new insights into scale-up behavior and to predict cell line-specific drawbacks at an early stage of process development.

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.

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.

Progress toward rapid, at-line N-glycosylation detection and control for recombinant protein expression

Proteins continue to represent a large fraction of the therapeutics market, reaching over a hundred billion dollars in market size globally. One key feature of protein modification that can affect both structure and function is the addition of glycosylation following protein folding, leading to regulatory requirements for the accurate assessment of protein attributes, including glycan structures. The non-template-driven, innately heterogeneous N-glycosylation process thus requires accurate detection to robustly generate protein therapies. A challenge exists in the timely detection of protein glycosylation without labor-intensive manipulation. In this article, we discuss progress toward N-glycoprotein control, focusing on novel control strategies and the advancement of rapid, high-throughput analysis methods.

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.

Harnessing the potential of machine learning for advancing "Quality by Design" in biomanufacturing

Ensuring consistent high yields and product quality are key challenges in biomanufacturing. Even minor deviations in critical process parameters (CPPs) such as media and feed compositions can significantly affect product critical quality attributes (CQAs). To identify CPPs and their interdependencies with product yield and CQAs, design of experiments, and multivariate statistical approaches are typically used in industry. Although these models can predict the effect of CPPs on product yield, there is room to improve CQA prediction performance by capturing the complex relationships in high-dimensional data. In this regard, machine learning (ML) approaches offer immense potential in handling non-linear datasets and thus are able to identify new CPPs that could effectively predict the CQAs. ML techniques can also be synergized with mechanistic models as a ‘hybrid ML’ or ‘white box ML’ to identify how CPPs affect the product yield and quality mechanistically, thus enabling rational design and control of the bioprocess. In this review, we describe the role of statistical modeling in Quality by Design (QbD) for biomanufacturing, and provide a generic outline on how relevant ML can be used to meaningfully analyze bioprocessing datasets. We then offer our perspectives on how relevant use of ML can accelerate the implementation of systematic QbD within the biopharma 4.0 paradigm.

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.

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.

Probabilistic model by Bayesian network for the prediction of antibody glycosylation in perfusion and fed-batch cell cultures

Glycosylation is a critical quality attribute of therapeutic monoclonal antibodies (mAbs). The glycan pattern can have a large impact on the immunological functions, serum half-life and stability. The medium components and cultivation parameters are known to potentially influence the glycosylation profile. Mathematical modelling provides a strategy for rational design and control of the upstream bioprocess. However, the kinetic models usually contain a very large number of unknown parameters, which limit their practical applications. In this article, we consider the metabolic network of N-linked glycosylation as a Bayesian network (BN) and calculate the fluxes of the glycosylation process as joint probability using the culture parameters as inputs. The modelling approach is validated with data of different Chinese hamster ovary cell cultures in pseudo perfusion, perfusion, and fed batch cultures, all showing very good predictive capacities. In cases where a large number of cultivation parameters is available, it is shown here that principal components analysis can efficiently be employed for a dimension reduction of the inputs compared to Pearson correlation analysis and feature importance by decision tree. The present study demonstrates that BN model can be a powerful tool in upstream process and medium development for glycoprotein productions.

The emerging role of cellular post-translational modifications in modulating growth and productivity of recombinant Chinese hamster ovary cells

Chinese hamster ovary (CHO) cells are one of the most commonly used host cell lines used for the production human therapeutic proteins. Much research over the past two decades has focussed on improving the growth, titre and cell specific productivity of CHO cells and in turn lowering the costs associated with production of recombinant proteins. CHO cell engineering has become of particular interest in recent years following the publication of the CHO cell genome and the availability of data relating to the proteome, transcriptome and metabolome of CHO cells. However, data relating to the cellular post-translational modification (PTMs) which can affect the functionality of CHO cellular proteins has only begun to be presented in recent years. PTMs are important to many cellular processes and can further alter proteins by increasing the complexity of proteins and their interactions. In this review, we describe the research presented from CHO cells to date related on three of the most important PTMs; glycosylation, phosphorylation and ubiquitination.

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.

Control of IgG glycosylation in CHO cell perfusion cultures by GReBA mathematical model supported by a novel targeted feed, TAFE

The N-linked glycosylation pattern is an important quality attribute of therapeutic glycoproteins. It has been reported by our group and by others that different carbon sources, such as glucose, mannose and galactose, can differently impact the glycosylation profile of glycoproteins in mammalian cell culture. Acting on the sugar feeding is thus an attractive strategy to tune the glycan pattern. However, in case of feeding of more than one carbon source simultaneously, the cells give priority to the one with the highest uptake rate, which limits the usage of this tuning, e.g. the cells favor consuming glucose in comparison to galactose. We present here a new feeding strategy (named 'TAFE' for targeted feeding) for perfusion culture to adjust the concentrations of fed sugars influencing the glycosylation. The strategy consists in setting the sugar feeding such that the cells are forced to consume these substrates at a target cell specific consumption rate decided by the operator and taking into account the cell specific perfusion rate (CSPR). This strategy is applied in perfusion cultures of Chinese hamster ovary (CHO) cells, illustrated by ten different regimes of sugar feeding, including glucose, galactose and mannose. Applying the TAFE strategy, different glycan profiles were obtained using the different feeding regimes. Furthermore, we successfully forced the cells to consume higher proportions of non-glucose sugars, which have lower transport rates than glucose in presence of this latter, in a controlled way. In previous work, a mathematical model named Glycan Residues Balance Analysis (GReBA) was developed to model the glycosylation profile based on the fed carbon sources. The present data were applied to the GReBA to design a feeding regime targeting a given glycosylation profile. The ability of the model to achieve this objective was confirmed by a multi-round of leave-one-out cross-validation (LOOCV), leading to the conclusion that the GReBA model can be used to design the feeding regime of a perfusion cell culture to obtain a desired glycosylation profile.

Digital Twins in Pharmaceutical and Biopharmaceutical Manufacturing: A Literature Review

The development and application of emerging technologies of Industry 4.0 enable the realization of digital twins (DT), which facilitates the transformation of the manufacturing sector to a more agile and intelligent one. DTs are virtual constructs of physical systems that mirror the behavior and dynamics of such physical systems. A fully developed DT consists of physical components, virtual components, and information communications between the two. Integrated DTs are being applied in various processes and product industries. Although the pharmaceutical industry has evolved recently to adopt Quality-by-Design (QbD) initiatives and is undergoing a paradigm shift of digitalization to embrace Industry 4.0, there has not been a full DT application in pharmaceutical manufacturing. Therefore, there is a critical need to examine the progress of the pharmaceutical industry towards implementing DT solutions. The aim of this narrative literature review is to give an overview of the current status of DT development and its application in pharmaceutical and biopharmaceutical manufacturing. State-of-the-art Process Analytical Technology (PAT) developments, process modeling approaches, and data integration studies are reviewed. Challenges and opportunities for future research in this field are also discussed.

Glycan Residues Balance Analysis - GReBA: A novel model for the N-linked glycosylation of IgG produced by CHO cells

The structure of N-linked glycosylation is a very important quality attribute for therapeutic monoclonal antibodies. Different carbon sources in cell culture media, such as mannose and galactose, have been reported to have different influences on the glycosylation patterns. Accurate prediction and control of the glycosylation profile are important for the process development of mammalian cell cultures. In this study, a mathematical model, that we named Glycan Residues Balance Analysis (GReBA), was developed based on the concept of Elementary Flux Mode (EFM), and used to predict the glycosylation profile for steady state cell cultures. Experiments were carried out in pseudo-perfusion cultivation of antibody producing Chinese Hamster Ovary (CHO) cells with various concentrations and combinations of glucose, mannose and galactose. Cultivation of CHO cells with mannose or the combinations of mannose and galactose resulted in decreased lactate and ammonium production, and more matured glycosylation patterns compared to the cultures with glucose. Furthermore, the growth rate and IgG productivity were similar in all the conditions. When the cells were cultured with galactose alone, lactate was fed as well to be used as complementary carbon source, leading to cell growth rate and IgG productivity comparable to feeding the other sugars. The data of the glycoprofiles were used for training the model, and then to simulate the glycosylation changes with varying the concentrations of mannose and galactose. In this study we showed that the GReBA model had a good predictive capacity of the N-linked glycosylation. The GReBA can be used as a guidance for development of glycoprotein cultivation processes.

Consequences of trace metal variability and supplementation on Chinese hamster ovary (CHO) cell culture performance: A review of key mechanisms and considerations

Trace metals are supplied to chemically-defined media (CDM) for optimal Chinese hamster ovary (CHO) cell culture performance during the production of monoclonal antibodies and other therapeutic proteins. However, lot-to-lot and vendor-to-vendor variability in raw materials consequently leads to an imbalance of trace metals that are supplied to CDM. This imbalance can yield detrimental effects rooted in several primary mechanisms and pathways including oxidative stress, apoptosis, lactate accumulation, and unfavorable glycan synthesis. Recent research endeavors involve supplying zinc, copper, and manganese to CDM in excess to further maximize culture productivity and product quality. These treatments significantly impact critical quality attributes and furthermore highlight the degree to which trace metal availability can affect CHO cell culture performance. This review highlights the role of trace metal variability, supplementation, and interplay on key cellular mechanisms responsible for overall culture performance and the production and quality of therapeutic proteins.

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.

Dissecting N-Glycosylation Dynamics in Chinese Hamster Ovary Cells Fed-batch Cultures using Time Course Omics Analyses

N-linked glycosylation affects the potency, safety, immunogenicity, and pharmacokinetic clearance of several therapeutic proteins including monoclonal antibodies. A robust control strategy is needed to dial in appropriate glycosylation profile during the course of cell culture processes accurately. However, N-glycosylation dynamics remains insufficiently understood owing to the lack of integrative analyses of factors that influence the dynamics, including sugar nucleotide donors, glycosyltransferases, and glycosidases. Here, an integrative approach involving multi-dimensional omics analyses was employed to dissect the temporal dynamics of glycoforms produced during fed-batch cultures of CHO cells. Several pathways including glycolysis, tricarboxylic citric acid cycle, and nucleotide biosynthesis exhibited temporal dynamics over the cell culture period. The steps involving galactose and sialic acid addition were determined as temporal bottlenecks. Our results show that galactose, and not manganese, is able to mitigate the temporal bottleneck, despite both being known effectors of galactosylation. Furthermore, sialylation is limited by the galactosylated precursors and autoregulation of cytidine monophosphate-sialic acid biosynthesis.

•

Major glycosylated species exhibit temporal dynamics during fed-batch processes

•

Key metabolic pathways linked to N-glycosylation exhibit significant temporal dynamics

•

Dynamics in nucleotide sugar donors (NSDs) directly influences glycoform heterogeneity

•

Glycoform heterogeneity can be mitigated by supplementing NSD biosynthetic precursors

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.