Application of a curated genome-scale metabolic model of CHO DG44 to an industrial fed-batch process

Driving towards digital biomanufacturing by CHO genome-scale models

Genome-scale metabolic models (GEMs) of Chinese hamster ovary (CHO) cells are valuable for gaining mechanistic understanding of mammalian cell metabolism and cultures. We provide a comprehensive overview of past and present developments of CHO-GEMs and in silico methods within the flux balance analysis (FBA) framework, focusing on their practical utility in rational cell line development and bioprocess improvements. There are many opportunities for further augmenting the model coverage and establishing integrative models that account for different cellular processes and data for future applications. With supportive collaborative efforts by the research community, we envisage that CHO-GEMs will be crucial for the increasingly digitized and dynamically controlled bioprocessing pipelines, especially because they can be successfully deployed in conjunction with artificial intelligence (AI) and systems engineering algorithms.

Building blocks needed for mechanistic modeling of bioprocesses: A critical review based on protein production by CHO cells

This paper reviews the key building blocks needed to develop a mechanistic model for use as an operational production tool. The Chinese Hamster Ovary (CHO) cell, one of the most widely used hosts for antibody production in the pharmaceutical industry, is considered as a case study. CHO cell metabolism is characterized by two main phases, exponential growth followed by a stationary phase with strong protein production. This process presents an appropriate degree of complexity to outline the modeling strategy. The paper is organized into four main steps: (1) CHO systems and data collection; (2) metabolic analysis; (3) formulation of the mathematical model; and finally, (4) numerical solution, calibration, and validation. The overall approach can build a predictive model of target variables. According to the literature, one of the main current modeling challenges lies in understanding and predicting the spontaneous metabolic shift. Possible candidates for the trigger of the metabolic shift include the concentration of lactate and carbon dioxide. In our opinion, ammonium, which is also an inhibiting product, should be further investigated. Finally, the expected progress in the emerging field of hybrid modeling, which combines the best of mechanistic modeling and machine learning, is presented as a fascinating breakthrough. Note that the modeling strategy discussed here is a general framework that can be applied to any bioprocess.

Development of a novel tyrosine-based selection system for generation of recombinant Chinese hamster ovary cells

Efficiently expanding Chinese hamster ovary (CHO) cells, which serve as the primary host cells for recombinant protein production, have gained increasing industrial significance. A significant hurdle in stable cell line development is the low efficiency of the target gene integrated into the host genome, implying the necessity for an effective screening and selection procedure to separate these stable cells. In this study, the genes of phenylalanine hydroxylase (PAH) and pterin 4 alpha carbinolamine dehydratase 1 (PCBD1), which are key enzymes in the tyrosine synthesis pathway, were utilized as selection markers and transduced into host cells together with the target genes. This research investigated the enrichment effect of this system and advanced further in understanding its benefits for cell line development and rCHO cell culture. A novel tyrosine-based selection system that only used PCBD1 as a selection marker was designed to promote the enrichment effect. Post 9 days of starvation, positive transductants in the cell pool approached 100%. Applied the novel tyrosine-based selection system, rCHO cells expressing E2 protein were generated and named CHO TS cells. It could continue to grow, and the yield of E2 achieved 95.95 mg/L in a tyrosine-free and chemically-defined (CD) medium. Herein, we introduced an alternative to antibiotic-based selections for the establishment of CHO cell lines and provided useful insights for the design and development of CD medium.

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.

How reliable are Chinese hamster ovary (CHO) cell genome-scale metabolic models?

Genome-scale metabolic models (GEMs) possess the power to revolutionize bioprocess and cell line engineering workflows thanks to their ability to predict and understand whole-cell metabolism in silico. Despite this potential, it is currently unclear how accurately GEMs can capture both intracellular metabolic states and extracellular phenotypes. Here, we investigate this knowledge gap to determine the reliability of current Chinese hamster ovary (CHO) cell metabolic models. We introduce a new GEM, iCHO2441, and create CHO-S and CHO-K1 specific GEMs. These are compared against iCHO1766, iCHO2048, and iCHO2291. Model predictions are assessed via comparison with experimentally measured growth rates, gene essentialities, amino acid auxotrophies, and 13 C intracellular reaction rates. Our results highlight that all CHO cell models are able to capture extracellular phenotypes and intracellular fluxes, with the updated GEM outperforming the original CHO cell GEM. Cell line-specific models were able to better capture extracellular phenotypes but failed to improve intracellular reaction rate predictions in this case. Ultimately, this work provides an updated CHO cell GEM to the community and lays a foundation for the development and assessment of next-generation flux analysis techniques, highlighting areas for model improvements.

A genome-scale dynamic constraint-based modelling (gDCBM) framework predicts growth dynamics, medium composition and intracellular flux distributions in CHO clonal variations

Optimizing mammalian cell growth and bioproduction is a tedious task. However, due to the inherent complexity of eukaryotic cells, heuristic experimental approaches such as, metabolic engineering and bioprocess design, are frequently integrated with mathematical models of cell culture to improve biological process efficiency and find paths for improvement. Constraint-based metabolic models have evolved over the last two decades to be used for dynamic modelling in addition to providing a linear description of steady-state metabolic systems. Formulation and implementation of the underlying optimization problems require special attention to the model's performance and feasibility, lack of defects in the definition of system components, and consideration of optimal alternate solutions, in addition to processing power limitations. Here, the time-resolved dynamics of a genome-scale metabolic network of Chinese hamster ovary (CHO) cell metabolism are shown using a genome-scale dynamic constraint-based modelling framework (gDCBM). The metabolic network was adapted from a reference model of CHO genome-scale metabolic model (GSMM), iCHO_DG44_v1, and dynamic restrictions were imposed to its exchange fluxes based on experimental results. We used this framework for predicting physiological changes in CHO clonal variants. Because of the methodical creation of the components for the flux balance analysis optimization problem and the integration of a switch time, this model can generate sequential predictions of intracellular fluxes during growth and non-growth phases (per hour of culture time) and transparently reveal the shortcomings in such practice. As a result of the differences exploited by various clones, we can understand the relevance of changes in intracellular flux distribution and exometabolomics. The integration of various omics data into the given gDCBM framework, as well as the reductionist analysis of the model, can further help bioprocess optimization.

Systems biology and metabolic modeling for cultivated meat: A promising approach for cell culture media optimization and cost reduction

The cultivated meat industry, also known as cell-based meat, cultured meat, lab-grown meat, or meat alternatives, is a growing field that aims to generate animal tissues ex-vivo in a cost-effective manner that achieves price parity with traditional agricultural products. However, cell culture media costs account for 55%-90% of production costs. To address this issue, efforts are aimed at optimizing media composition. Systems biology-driven approaches have been successfully used to improve the biomass and productivity of multiple bioproduction platforms, like Chinese hamster ovary cells, by accelerating the development of cell line-specific media and reducing research and development and production costs related to cell media and its optimization. In this review, we summarize systems biology modeling approaches, methods for cell culture media and bioprocess optimization, and metabolic studies done in animals of interest to the cultivated meat industry. More importantly, we identify current gaps in knowledge that prevent the identification of metabolic bottlenecks. These include the lack of genome-scale metabolic models for some species (pigs and ducks), a lack of accurate biomass composition studies for different growth conditions, and 13 C-metabolic flux analysis (MFA) studies for many of the species of interest for the cultivated meat industry (only shrimp and duck cells have been subjected to 13 C-MFA). We also highlight the importance of characterizing the metabolic requirements of cells at the organism, breed, and cell line-specific levels, and we outline future steps that this nascent field needs to take to achieve price parity and production efficiency similar to those of other bioproduction platforms. Practical Application: Our article summarizes systems biology techniques for cell culture media design and bioprocess optimization, which may be used to significantly reduce cell-based meat production costs. We also present the results of experimental studies done on some of the species of interest to the cultivated meat industry and highlight why modeling approaches are required for multiple species, cell-types, and cell lines.

Machine learning and circular bioeconomy: Building new resource efficiency from diverse waste streams

Biorefinery systems are playing pivotal roles in the technological support of resource efficiency for circular bioeconomy. Meanwhile, artificial intelligence presents great potential in handling scientific tasks of high-dimensional complexity. This review article scrutinizes the status of machine learning (ML) applications in four critical biorefinery systems (i.e. composting, fermentation, anaerobic digestion, and thermochemical conversions) as well as their advancements against traditional modeling techniques of mechanistic approach. The contents cover their algorithm selections, modeling challenges, and prospective improvements. Perspectives are sketched to further inform collective efforts on crucial aspects. The multidisciplinary interchange of modeling knowledge will enable a more progressive digital transformation of sustainability efforts in supporting sustainable development goals.

Genome-scale modeling of Chinese hamster ovary cells by hybrid semi-parametric flux balance analysis

Flux balance analysis (FBA) is currently the standard method to compute metabolic fluxes in genome-scale networks. Several FBA extensions employing diverse objective functions and/or constraints have been published. Here we propose a hybrid semi-parametric FBA extension that combines mechanistic-level constraints (parametric) with empirical constraints (non-parametric) in the same linear program. A CHO dataset with 27 measured exchange fluxes obtained from 21 reactor experiments served to evaluate the method. The mechanistic constraints were deduced from a reduced CHO-K1 genome-scale network with 686 metabolites, 788 reactions and 210 degrees of freedom. The non-parametric constraints were obtained by principal component analysis of the flux dataset. The two types of constraints were integrated in the same linear program showing comparable computational cost to standard FBA. The hybrid FBA is shown to significantly improve the specific growth rate prediction under different constraints scenarios. A metabolically efficient cell growth feed targeting minimal byproducts accumulation was designed by hybrid FBA. It is concluded that integrating parametric and nonparametric constraints in the same linear program may be an efficient approach to reduce the solution space and to improve the predictive power of FBA methods when critical mechanistic information is missing.

Reconstruction of a generic genome-scale metabolic network for chicken: Investigating network connectivity and finding potential biomarkers

Chicken is the first sequenced avian that has a crucial role in human life for its meat and egg production. Because of various metabolic disorders, study the metabolism of chicken cell is important. Herein, the first genome-scale metabolic model of a chicken cell named iES1300, consists of 2427 reactions, 2569 metabolites, and 1300 genes, was reconstructed manually based on KEGG, BiGG, CHEBI, UNIPROT, REACTOME, and MetaNetX databases. Interactions of metabolic genes for growth were examined for E. coli, S. cerevisiae, human, and chicken metabolic models. The results indicated robustness to genetic manipulation for iES1300 similar to the results for human. iES1300 was integrated with transcriptomics data using algorithms and Principal Component Analysis was applied to compare context-specific models of the normal, tumor, lean and fat cell lines. It was found that the normal model has notable metabolic flexibility in the utilization of various metabolic pathways, especially in metabolic pathways of the carbohydrate metabolism, compared to the others. It was also concluded that the fat and tumor models have similar growth metabolisms and the lean chicken model has a more active lipid and carbohydrate metabolism.

Advances of Glycometabolism Engineering in Chinese Hamster Ovary Cells

As the most widely used mammalian cell line, Chinese hamster ovary (CHO) cells can express various recombinant proteins with a post translational modification pattern similar to that of the proteins from human cells. During industrial production, cells need large amounts of ATP to support growth and protein expression, and since glycometabolism is the main source of ATP for cells, protein production partly depends on the efficiency of glycometabolism. And efficient glycometabolism allows less glucose uptake by cells, reducing production costs, and providing a better mammalian production platform for recombinant protein expression. In the present study, a series of progresses on the comprehensive optimization in CHO cells by glycometabolism strategy were reviewed, including carbohydrate intake, pyruvate metabolism and mitochondrial metabolism. We analyzed the effects of gene regulation in the upstream and downstream of the glucose metabolism pathway on cell’s growth and protein expression. And we also pointed out the latest metabolic studies that are potentially applicable on CHO cells. In the end, we elaborated the application of metabolic models in the study of CHO cell metabolism.

Perfusion culture of Chinese Hamster Ovary cells for bioprocessing applications

Much of the biopharmaceutical industry's success over the past 30 years has relied on products derived from Chinese Hamster Ovary (CHO) cell lines. During this time, improvements in mammalian cell cultures have come from cell line development and process optimization suited for large-scale fed-batch processes. Originally developed for high cell densities and sensitive products, perfusion processes have a long history. Driven by high volumetric titers and a small footprint, perfusion-based bioprocess research has regained an interest from academia and industry. The recent pandemic has further highlighted the need for such intensified biomanufacturing options. In this review, we outline the technical history of research in this field as it applies to biologics production in CHO cells. We demonstrate a number of emerging trends in the literature and corroborate these with underlying drivers in the commercial space. From these trends, we speculate that the future of perfusion bioprocesses is bright and that the fields of media optimization, continuous processing, and cell line engineering hold the greatest potential. Aligning in its continuous setup with the demands for Industry 4.0, perfusion biomanufacturing is likely to be a hot topic in the years to come.

Synergisms of machine learning and constraint-based modeling of metabolism for analysis and optimization of fermentation parameters

Recent noteworthy advances in the development of high-performing microbial and mammalian strains have enabled the sustainable production of bio-economically valuable substances such as bio-compounds, biofuels, and biopharmaceuticals. However, to obtain an industrially viable mass-production scheme, much time and effort are required. The robust and rational design of fermentation processes requires analysis and optimization of different extracellular conditions and medium components, which have a massive effect on growth and productivity. In this regard, knowledge- and data-driven modeling methods have received much attention. Constraint-based modeling (CBM) is a knowledge-driven mathematical approach that has been widely used in fermentation analysis and optimization due to its capabilities of predicting the cellular phenotype from genotype through high-throughput means. On the other hand, machine learning (ML) is a data-driven statistical method that identifies the data patterns within sophisticated biological systems and processes, where there is inadequate knowledge to represent underlying mechanisms. Furthermore, ML models are becoming a viable complement to constraint-based models in a reciprocal manner when one is used as a pre-step of another. As a result, more predictable model is produced. This review highlights the applications of CBM and ML independently and the combination of these two approaches for analyzing and optimizing fermentation parameters. This article is protected by copyright. All rights reserved.

From omics to Cellular mechanisms in mammalian cell factory development

Mammalian cells have been used widely as biopharmaceutical cell factories due to their ability to make complex biotherapeutic proteins with human-compatible modifications. However, their application for some products has been hampered by low protein yields. Numerous studies have aimed to characterize cellular bottlenecks in the hope of boosting protein productivity, but the complexity of the underlying pathways and the diversity of the modifications have complicated cell engineering when relying solely on traditional methodologies. Incorporating omics-based and systems approaches into cell engineering can provide valuable insights into desirable phenotypes of cell factories. Here, we discuss cell engineering strategies for enhancing protein productivity in mammalian cell factories, particularly CHO and HEK293, and the opportunities and limitations of the genome-wide screening and multi-omics approaches for guiding cell engineering. Systems biology strategies will also be discussed to show how they refine our understanding of the cellular mechanisms which will aid in effective engineering strategies.

Strategies for Enhancing in vitro Degradation of Linuron by Variovorax sp. Strain SRS 16 Under the Guidance of Metabolic Modeling

Phenyl urea herbicides are being extensively used for weed control in both agricultural and non-agricultural applications. Linuron is one of the key herbicides in this family and is in wide use. Like other phenyl urea herbicides, it is known to have toxic effects as a result of its persistence in the environment. The natural removal of linuron from the environment is mainly carried through microbial biodegradation. Some microorganisms have been reported to mineralize linuron completely and utilize it as a carbon and nitrogen source. Variovorax sp. strain SRS 16 is one of the known efficient degraders with a recently sequenced genome. The genomic data provide an opportunity to use a genome-scale model for improving biodegradation. The aim of our study is the construction of a genome-scale metabolic model following automatic and manual protocols and its application for improving its metabolic potential through iterative simulations. Applying flux balance analysis (FBA), growth and degradation performances of SRS 16 in different media considering the influence of selected supplements (potential carbon and nitrogen sources) were simulated. Outcomes are predictions for the suitable media modification, allowing faster degradation of linuron by SRS 16. Seven metabolites were selected for in vitro validation of the predictions through laboratory experiments confirming the degradation-promoting effect of specific amino acids (glutamine and asparagine) on linuron degradation and SRS 16 growth. Overall, simulations are shown to be efficient in predicting the degradation potential of SRS 16 in the presence of specific supplements. The generated information contributes to the understanding of the biochemistry of linuron degradation and can be further utilized for the development of new cleanup solutions without any genetic manipulation.

Error propagation in constraint-based modeling of Chinese hamster ovary cells

Chinese hamster ovary (CHO) cells are the most popular mammalian cell factories for the production of glycosylated biopharmaceuticals. To further increase titer and productivity and ensure product quality, rational system-level engineering strategies based on constraint-based metabolic modeling, such as flux balance analysis (FBA), have gained strong interest. However, the quality of FBA predictions depends on the accuracy of the experimental input data, especially on the exchange rates of extracellular metabolites. Yet, it is not standard practice to devote sufficient attention to the accurate determination of these rates. In this work, we investigated to what degree the sampling frequency during a batch culture and the measurement errors of metabolite concentrations influence the accuracy of the calculated exchange rates and further, how this error then propagates into FBA predictions of growth rates. We determined that accurate measurements of essential amino acids with low uptake rates are crucial for the accuracy of FBA predictions, followed by a sufficient number of analyzed time points. We observed that the measured difference in growth rates of two cell lines can only be reliably predicted when both high measurement accuracy and sampling frequency are ensured.

Systematic evaluation of parameters for genome-scale metabolic models of cultured mammalian cells

Genome-scale metabolic models describe cellular metabolism with mechanistic detail. Given their high complexity, such models need to be parameterized correctly to yield accurate predictions and avoid overfitting. Effective parameterization has been well-studied for microbial models, but it remains unclear for higher eukaryotes, including mammalian cells. To address this, we enumerated model parameters that describe key features of cultured mammalian cells - including cellular composition, bioprocess performance metrics, mammalian-specific pathways, and biological assumptions behind model formulation approaches. We tested these parameters by building thousands of metabolic models and evaluating their ability to predict the growth rates of a panel of phenotypically diverse Chinese Hamster Ovary cell clones. We found the following considerations to be most critical for accurate parameterization: (1) cells limit metabolic activity to maintain homeostasis, (2) cell morphology and viability change dynamically during a growth curve, and (3) cellular biomass has a particular macromolecular composition. Depending on parameterization, models predicted different metabolic phenotypes, including contrasting mechanisms of nutrient utilization and energy generation, leading to varying accuracies of growth rate predictions. Notably, accurate parameter values broadly agreed with experimental measurements. These insights will guide future investigations of mammalian metabolism.

A genome-scale metabolic network model and machine learning predict amino acid concentrations in Chinese Hamster Ovary cell cultures

The control of nutrient availability is critical to large-scale manufacturing of biotherapeutics. However, the quantification of proteinogenic amino acids is time-consuming and thus is difficult to implement for real-time in situ bioprocess control. Genome-scale metabolic models describe the metabolic conversion from media nutrients to proliferation and recombinant protein production, and therefore are a promising platform for in silico monitoring and prediction of amino acid concentrations. This potential has not been realized due to unresolved challenges: (1) the models assume an optimal and highly efficient metabolism, and therefore tend to underestimate amino acid consumption, and (2) the models assume a steady state, and therefore have a short forecast range. We address these challenges by integrating machine learning with the metabolic models. Through this we demonstrate accurate and time-course dependent prediction of individual amino acid concentration in culture medium throughout the production process. Thus, these models can be deployed to control nutrient feeding to avoid premature nutrient depletion or provide early predictions of failed bioreactor runs.

Modelling Cell Metabolism: A Review on Constraint-Based Steady-State and Kinetic Approaches

Studying cell metabolism serves a plethora of objectives such as the enhancement of bioprocess performance, and advancement in the understanding of cell biology, of drug target discovery, and in metabolic therapy. Remarkable successes in these fields emerged from heuristics approaches, for instance, with the introduction of effective strategies for genetic modifications, drug developments and optimization of bioprocess management. However, heuristics approaches have showed significant shortcomings, such as to describe regulation of metabolic pathways and to extrapolate experimental conditions. In the specific case of bioprocess management, such shortcomings limit their capacity to increase product quality, while maintaining desirable productivity and reproducibility levels. For instance, since heuristics approaches are not capable of prediction of the cellular functions under varying experimental conditions, they may lead to sub-optimal processes. Also, such approaches used for bioprocess control often fail in regulating a process under unexpected variations of external conditions. Therefore, methodologies inspired by the systematic mathematical formulation of cell metabolism have been used to address such drawbacks and achieve robust reproducible results. Mathematical modelling approaches are effective for both the characterization of the cell physiology, and the estimation of metabolic pathways utilization, thus allowing to characterize a cell population metabolic behavior. In this article, we present a review on methodology used and promising mathematical modelling approaches, focusing primarily to investigate metabolic events and regulation. Proceeding from a topological representation of the metabolic networks, we first present the metabolic modelling approaches that investigate cell metabolism at steady state, complying to the constraints imposed by mass conservation law and thermodynamics of reactions reversibility. Constraint-based models (CBMs) are reviewed highlighting the set of assumed optimality functions for reaction pathways. We explore models simulating cell growth dynamics, by expanding flux balance models developed at steady state. Then, discussing a change of metabolic modelling paradigm, we describe dynamic kinetic models that are based on the mathematical representation of the mechanistic description of nonlinear enzyme activities. In such approaches metabolic pathway regulations are considered explicitly as a function of the activity of other components of metabolic networks and possibly far from the metabolic steady state. We have also assessed the significance of metabolic model parameterization in kinetic models, summarizing a standard parameter estimation procedure frequently employed in kinetic metabolic modelling literature. Finally, some optimization practices used for the parameter estimation are reviewed.

Incorporating RNA-Seq transcriptomics into glycosylation-integrating metabolic network modelling kinetics: Multiomic Chinese hamster ovary (CHO) cell bioreactors

In this work, the kinetic model based on the previously developed metabolic and glycan reaction networks of the ovarian cells of the Chinese hamster ovary (CHO) cell line was improved by the inclusion of transcriptomic data that took into account the values of the RPKM gene (Reads per Kilobase of Exon per Million Reads Mapped). The transcriptomic (RNASeq) data were obtained together with metabolic and glycan data from the literature, and the concentrations with RPKM values were collected at several points in time from two fed-batch processes. First, the fluxes were determined by regression analysis of the metabolic data, then these fluxes were corrected by using the fold change in gene expression as a measure of enzyme concentrations. Next, the corrected fluxes in the kinetic model were used to calculate the concentration profiles of the metabolites, and literature data were used to evaluate the predicted results of the model. Compared to other studies where the concentration profiles of CHO cell metabolites were described using a kinetic model without consideration of RNA-Seq data to correct the fluxes, this model is unique. The additional integration of transcriptomic data led to better predictions of metabolic concentrations in the fed-batch process, which is a significant improvement of the modelling technique used.

The era of big data: Genome-scale modelling meets machine learning

With omics data being generated at an unprecedented rate, genome-scale modelling has become pivotal in its organisation and analysis. However, machine learning methods have been gaining ground in cases where knowledge is insufficient to represent the mechanisms underlying such data or as a means for data curation prior to attempting mechanistic modelling. We discuss the latest advances in genome-scale modelling and the development of optimisation algorithms for network and error reduction, intracellular constraining and applications to strain design. We further review applications of supervised and unsupervised machine learning methods to omics datasets from microbial and mammalian cell systems and present efforts to harness the potential of both modelling approaches through hybrid modelling.

Systematically gap-filling the genome-scale metabolic model of CHO cells

OBJECTIVE

Chinese hamster ovary (CHO) cells are the leading cell factories for producing recombinant proteins in the biopharmaceutical industry. In this regard, constraint-based metabolic models are useful platforms to perform computational analysis of cell metabolism. These models need to be regularly updated in order to include the latest biochemical data of the cells, and to increase their predictive power. Here, we provide an update to iCHO1766, the metabolic model of CHO cells.

RESULTS

We expanded the existing model of Chinese hamster metabolism with the help of four gap-filling approaches, leading to the addition of 773 new reactions and 335 new genes. We incorporated these into an updated genome-scale metabolic network model of CHO cells, named iCHO2101. In this updated model, the number of reactions and pathways capable of carrying flux is substantially increased.

CONCLUSIONS

The present CHO model is an important step towards more complete metabolic models of CHO cells.

BiGG Models 2020: multi-strain genome-scale models and expansion across the phylogenetic tree

The BiGG Models knowledge base (http://bigg.ucsd.edu) is a centralized repository for high-quality genome-scale metabolic models. For the past 12 years, the website has allowed users to browse and search metabolic models. Within this update, we detail new content and features in the repository, continuing the original effort to connect each model to genome annotations and external databases as well as standardization of reactions and metabolites. We describe the addition of 31 new models that expand the portion of the phylogenetic tree covered by BiGG Models. We also describe new functionality for hosting multi-strain models, which have proven to be insightful in a variety of studies centered on comparisons of related strains. Finally, the models in the knowledge base have been benchmarked using Memote, a new community-developed validator for genome-scale models to demonstrate the improving quality and transparency of model content in BiGG Models.

Key Challenges in Designing CHO Chassis Platforms

Following the success of and the high demand for recombinant protein-based therapeutics during the last 25 years, the pharmaceutical industry has invested significantly in the development of novel treatments based on biologics. Mammalian cells are the major production systems for these complex biopharmaceuticals, with Chinese hamster ovary (CHO) cell lines as the most important players. Over the years, various engineering strategies and modeling approaches have been used to improve microbial production platforms, such as bacteria and yeasts, as well as to create pre-optimized chassis host strains. However, the complexity of mammalian cells curtailed the optimization of these host cells by metabolic engineering. Most of the improvements of titer and productivity were achieved by media optimization and large-scale screening of producer clones. The advances made in recent years now open the door to again consider the potential application of systems biology approaches and metabolic engineering also to CHO. The availability of a reference genome sequence, genome-scale metabolic models and the growing number of various “omics” datasets can help overcome the complexity of CHO cells and support design strategies to boost their production performance. Modular design approaches applied to engineer industrially relevant cell lines have evolved to reduce the time and effort needed for the generation of new producer cells and to allow the achievement of desired product titers and quality. Nevertheless, important steps to enable the design of a chassis platform similar to those in use in the microbial world are still missing. In this review, we highlight the importance of mammalian cellular platforms for the production of biopharmaceuticals and compare them to microbial platforms, with an emphasis on describing novel approaches and discussing still open questions that need to be resolved to reach the objective of designing enhanced modular chassis CHO cell lines.

Using Metabolomics to Identify Cell Line-Independent Indicators of Growth Inhibition for Chinese Hamster Ovary Cell-based Bioprocesses

Chinese hamster ovary (CHO) cells are widely used for the production of biopharmaceuticals. Efforts to improve productivity through medium design and feeding strategy optimization have focused on preventing the depletion of essential nutrients and managing the accumulation of lactate and ammonia. In addition to ammonia and lactate, many other metabolites accumulate in CHO cell cultures, although their effects remain largely unknown. Elucidating these effects has the potential to further improve the productivity of CHO cell-based bioprocesses. This study used untargeted metabolomics to identify metabolites that accumulate in fed-batch cultures of monoclonal antibody (mAb) producing CHO cells. The metabolomics experiments profiled six cell lines that are derived from two different hosts, produce different mAbs, and exhibit different growth profiles. Comparing the cell lines’ metabolite profiles at different growth stages, we found a strong negative correlation between peak viable cell density (VCD) and a tryptophan metabolite, putatively identified as 5-hydroxyindoleacetaldehyde (5-HIAAld). Amino acid supplementation experiments showed strong growth inhibition of all cell lines by excess tryptophan, which correlated with the accumulation of 5-HIAAld in the culture medium. Prospectively, the approach presented in this study could be used to identify cell line- and host-independent metabolite markers for clone selection and bioprocess development.

Enzyme capacity-based genome scale modelling of CHO cells

Chinese hamster ovary (CHO) cells are most prevalently used for producing recombinant therapeutics in biomanufacturing. Recently, more rational and systems approaches have been increasingly exploited to identify key metabolic bottlenecks and engineering targets for cell line engineering and process development based on the CHO genome-scale metabolic model which mechanistically characterizes cell culture behaviours. However, it is still challenging to quantify plausible intracellular fluxes and discern metabolic pathway usages considering various clonal traits and bioprocessing conditions. Thus, we newly incorporated enzyme kinetic information into the updated CHO genome-scale model (iCHO2291) and added enzyme capacity constraints within the flux balance analysis framework (ecFBA) to significantly reduce the flux variability in biologically meaningful manner, as such improving the accuracy of intracellular flux prediction. Interestingly, ecFBA could capture the overflow metabolism under the glucose excess condition where the usage of oxidative phosphorylation is limited by the enzyme capacity. In addition, its applicability was successfully demonstrated via a case study where the clone- and media-specific lactate metabolism was deciphered, suggesting that the lactate-pyruvate cycling could be beneficial for CHO cells to efficiently utilize the mitochondrial redox capacity. In summary, iCHO2296 with ecFBA can be used to confidently elucidate cell cultures and effectively identify key engineering targets, thus guiding bioprocess optimization and cell engineering efforts as a part of digital twin model for advanced biomanufacturing in future.

Integration of Time-Series Transcriptomic Data with Genome-Scale CHO Metabolic Models for mAb Engineering

Chinese hamster ovary (CHO) cells are the most commonly used cell lines in biopharmaceutical manufacturing. Genome-scale metabolic models have become a valuable tool to study cellular metabolism. Despite the presence of reference global genome-scale CHO model, context-specific metabolic models may still be required for specific cell lines (for example, CHO-K1, CHO-S, and CHO-DG44), and for specific process conditions. Many integration algorithms have been available to reconstruct specific genome-scale models. These methods are mainly based on integrating omics data (i.e., transcriptomics, proteomics, and metabolomics) into reference genome-scale models. In the present study, we aimed to investigate the impact of time points of transcriptomics integration on the genome-scale CHO model by assessing the prediction of growth rates with each reconstructed model. We also evaluated the feasibility of applying extracted models to different cell lines (generated from the same parental cell line). Our findings illustrate that gene expression at various stages of culture slightly impacts the reconstructed models. However, the prediction capability is robust enough on cell growth prediction not only across different growth phases but also in expansion to other cell lines.

Robust Analytical Methods for the Accurate Quantification of the Total Biomass Composition of Mammalian Cells

Biomass composition is an important input for genome-scale metabolic models and has a big impact on their predictive capabilities. However, researchers often rely on generic data for biomass composition, e.g. collected from similar organisms. This leads to inaccurate predictions, because biomass composition varies between different cell lines, conditions, and growth phases. In this chapter we present protocols for the determination of the biomass composition of Chinese Hamster Ovary (CHO) cells. These methods can easily be adapted to other types of mammalian cells. The protocols include the quantification of cell dry mass and of the main biomass components, namely protein, lipid, DNA, RNA, and carbohydrates. Cell dry mass is determined gravimetrically by weighing a defined number of cells. Amino acid composition and protein content are measured by gas chromatography mass spectrometry. Lipids are quantified by shotgun mass spectrometry, which provides quantities for the different lipid classes and also the distribution of fatty acids. RNA is purified and then quantified spectrophotometrically. The methods for DNA and carbohydrates are simple fluorometric and colorimetric assays adapted to a 96-well plate format. To ensure quantitative results, internal standards or spike-in controls are used in all methods, e.g. to account for possible matrix effects or loss of material. Finally, the last section provides a guide on how to convert the measured data into biomass equations, which can then be integrated into a metabolic model.

Computer-Aided Strategies for Determining the Amino Acid Composition of Medium for Chinese Hamster Ovary Cell-Based Biomanufacturing Platforms

Chinese hamster ovary (CHO) cells are used for the production of the majority of biopharmaceutical drugs, and thus have remained the standard industry host for the past three decades. The amino acid composition of the medium plays a key role in commercial scale biologics manufacturing, as amino acids constitute the building blocks of both endogenous and heterologous proteins, are involved in metabolic and non-metabolic pathways, and can act as main sources of nitrogen and carbon under certain conditions. As biomanufactured proteins become increasingly complex, the adoption of model-based approaches become ever more popular in complementing the challenging task of medium development. The extensively studied amino acid metabolism is exceptionally suitable for such model-driven analyses, and although still limited in practice, the development of these strategies is gaining attention, particularly in this domain. This paper provides a review of recent efforts. We first provide an overview of the widely adopted practice, and move on to describe the model-driven approaches employed for the improvement and optimization of the external amino acid supply in light of cellular amino acid demand. We conclude by proposing the likely prevalent direction the field is heading towards, providing a critical evaluation of the current state and the future challenges and considerations.

Systems biology approach in the formulation of chemically defined media for recombinant protein overproduction

The cell culture medium is an intricate mixture of components which has a tremendous effect on cell growth and recombinant protein production. Regular cell culture medium includes various components, and the decision about which component should be included in the formulation and its optimum amount is an underlying issue in biotechnology industries. Applying conventional techniques to design an optimal medium for the production of a recombinant protein requires meticulous and immense research. Moreover, since the medium formulation for the production of one protein could not be the best choice for another protein, hence, the most suitable media should be determined for each recombinant cell line. Accordingly, medium formulation becomes a laborious, time-consuming, and costly process in biomanufacturing of recombinant protein, and finding alternative strategies for medium development seems to be crucial. In silico modeling is an attractive concept to be adapted for medium formulation due to its high potential to supersede laboratory examinations. By emerging the high-throughput datasets, scientists can disclose the knowledge about the effect of medium components on cell growth and metabolism, and via applying this information through systems biology approach, medium formulation optimization could be accomplished in silico with no need of significant amount of experimentation. This review demonstrates some of the applications of systems biology as a powerful tool for medium development and illustrates the effect of medium optimization with system-level analysis on the production of recombinant proteins in different host cells.

Genetic engineering approaches to improve posttranslational modification of biopharmaceuticals in different production platforms

The number of approved biopharmaceuticals, where product quality attributes remain of major importance, is increasing steadily. Within the available variety of expression hosts, the production of biopharmaceuticals faces diverse limitations with respect to posttranslational modifications (PTM), while different biopharmaceuticals demand different forms and specifications of PTMs for proper functionality. With the growing toolbox of genetic engineering technologies, it is now possible to address general as well as host- or biopharmaceutical-specific product quality obstacles. In this review, we present diverse expression systems derived from mammalians, bacteria, yeast, plants, and insects as well as available genetic engineering tools. We focus on genes for knockout/knockdown and overexpression for meaningful approaches to improve biopharmaceutical PTMs and discuss their applicability as well as future trends in the field.

Strategic Framework for Parameterization of Cell Culture Models

Global Sensitivity Analysis (GSA) is a technique that numerically evaluates the significance of model parameters with the aim of reducing the number of parameters that need to be estimated accurately from experimental data. In the work presented herein, we explore different methods and criteria in the sensitivity analysis of a recently developed mathematical model to describe Chinese hamster ovary (CHO) cell metabolism in order to establish a strategic, transferable framework for parameterizing mechanistic cell culture models. For that reason, several types of GSA employing different sampling methods (Sobol’, Pseudo-random and Scrambled-Sobol’), parameter deviations (10%, 30% and 50%) and sensitivity index significance thresholds (0.05, 0.1 and 0.2) were examined. The results were evaluated according to the goodness of fit between the simulation results and experimental data from fed-batch CHO cell cultures. Then, the predictive capability of the model was tested against four different feeding experiments. Parameter value deviation levels proved not to have a significant effect on the results of the sensitivity analysis, while the Sobol’ and Scrambled-Sobol’ sampling methods and a 0.1 significance threshold were found to be the optimum settings. The resulting framework was finally used to calibrate the model for another CHO cell line, resulting in a good overall fit. The results of this work set the basis for the use of a single mechanistic metabolic model that can be easily adapted through the proposed sensitivity analysis method to the behavior of different cell lines and therefore minimize the experimental cost of model development.