Macroscopic modeling of mammalian cell growth and metabolism

Dynamic optimization of an integrated cultivation-aggregation model for mAb production

Due to their proteinaceous structure, monoclonal antibodies (mAbs) are susceptible to irreversible aggregation, with harmful consequences on drug efficacy and patient safety. To mitigate this risk in modern biopharmaceutical processes, it is critical to comply with current good manufacturing practices (cGMP) and pursue operating strategies minimizing irreversible aggregation whilst also maximizing mAb throughput. These conflicting objectives are targeted in this study by formulating and analyzing an integrated dynamic model accounting for both cultivation and aggregation of mAbs from a Chinese Hamster Ovary (CHO) cell line. Two manipulated dynamic variables are considered here in simulation studies: firstly temperature manipulation within a batch reactor, and secondly feed flow manipulation within a series of isothermal fed-batch reactors. Following this, dynamic optimization investigations have been conducted, firstly with the single objective of maximizing mAb throughput and secondly with multiple (two) objectives of maximizing mAb throughput while also minimizing irreversible aggregate content, simultaneously. The study provides key insight into tradeoffs of how simultaneous temperature and feed flowrate manipulation affects mAb throughput and aggregation inside bioreactors.

Challenges of Spatially Resolved Metabolism in Cancer Research

Stable isotope-resolved metabolomics comprises a critical set of technologies that can be applied to a wide variety of systems, from isolated cells to whole organisms, to define metabolic pathway usage and responses to perturbations such as drugs or mutations, as well as providing the basis for flux analysis. As the diversity of stable isotope-enriched compounds is very high, and with newer approaches to multiplexing, the coverage of metabolism is now very extensive. However, as the complexity of the model increases, including more kinds of interacting cell types and interorgan communication, the analytical complexity also increases. Further, as studies move further into spatially resolved biology, new technical problems have to be overcome owing to the small number of analytes present in the confines of a single cell or cell compartment. Here, we review the overall goals and solutions made possible by stable isotope tracing and their applications to models of increasing complexity. Finally, we discuss progress and outstanding difficulties in high-resolution spatially resolved tracer-based metabolic studies.

A perspective-driven and technical evaluation of machine learning in bioreactor scale-up: A case-study for potential model developments

Bioreactor scale-up and scale-down have always been a topical issue for the biopharmaceutical industry and despite considerable effort, the identification of a fail-safe strategy for bioprocess development across scales remains a challenge. With the ubiquitous growth of digital transformation technologies, new scaling methods based on computer models may enable more effective scaling. This study aimed to evaluate the potential application of machine learning (ML) algorithms for bioreactor scale-up, with a specific focus on the prediction of scaling parameters. Factors critical to the development of such models were identified and data for bioreactor scale-up studies involving CHO cell-generated mAb products collated from the literature and public sources for the development of unsupervised and supervised ML models. Comparison of bioreactor performance across scales identified similarities between the different processes and primary differences between small- and large-scale bioreactors. A series of three case studies were developed to assess the relationship between cell growth and scale-sensitive bioreactor features. An embedding layer improved the capability of artificial neural network models to predict cell growth at a large-scale, as this approach captured similarities between the processes. Further models constructed to predict scaling parameters demonstrated how ML models may be applied to assist the scaling process. The development of data sets that include more characterization data with greater variability under different gassing and agitation regimes will also assist the future development of ML tools for bioreactor scaling.

Integration Approaches to Model Bioreactor Hydrodynamics and Cellular Kinetics for Advancing Bioprocess Optimisation

Large-scale bioprocesses are increasing globally to cater to the larger market demands for biological products. As fermenter volumes increase, the efficiency of mixing decreases, and environmental gradients become more pronounced compared to smaller scales. Consequently, the cells experience gradients in process parameters, which in turn affects the efficiency and profitability of the process. Computational fluid dynamics (CFD) simulations are being widely embraced for their ability to simulate bioprocess performance, facilitate bioprocess upscaling, downsizing, and process optimisation. Recently, CFD approaches have been integrated with dynamic Cell reaction kinetic (CRK) modelling to generate valuable information about the cellular response to fluctuating hydrodynamic parameters inside large production processes. Such coupled approaches have the potential to facilitate informed decision-making in intelligent biomanufacturing, aligning with the principles of "Industry 4.0" concerning digitalisation and automation. In this review, we discuss the benefits of utilising integrated CFD-CRK models and the different approaches to integrating CFD-based bioreactor hydrodynamic models with cellular kinetic models. We also highlight the suitability of different coupling approaches for bioprocess modelling in the purview of associated computational loads.

Integrated Experimental and Mathematical Exploration of Modular Tissue Cultures for Developmental Engineering

Developmental engineering (DE) involves culturing various cells on modular scaffolds (MSs), yielding modular tissues (MTs) assembled into three-dimensional (3D) tissues, mimicking developmental biology. This study employs an integrated approach, merging experimental and mathematical methods to investigate the biological processes in MT cultivation and assembly. Human dermal fibroblasts (HDFs) were cultured on tissue culture plastics, poly(lactic acid) (PLA) discs with regular open structures, or spherical poly(methyl methacrylate) (PMMA) MSs, respectively. Notably, HDFs exhibited flattened spindle shapes when adhered to solid surfaces, and complex 3D structures when migrating into the structured voids of PLA discs or interstitial spaces between aggregated PMMA MSs, showcasing coordinated colonization of porous scaffolds. Empirical investigations led to power law models simulating density-dependent cell growth on solid surfaces or voids. Concurrently, a modified diffusion model was applied to simulate oxygen diffusion within tissues cultured on solid surfaces or porous structures. These mathematical models were subsequently combined to explore the influences of initial cell seeding density, culture duration, and oxygen diffusion on MT cultivation and assembly. The findings underscored the intricate interplay of factors influencing MT design for tissue assembly. The integrated approach provides insights into mechanistic aspects, informing bioprocess design for manufacturing MTs and 3D tissues in DE.

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.

A genome-scale nutrient minimization forecast algorithm for controlling essential amino acid levels in CHO cell cultures

Mammalian cell culture processes rely heavily on empirical knowledge in which process control remains a challenge due to the limited characterization/understanding of cell metabolism and inability to predict the cell behaviors. This study facilitates control of Chinese hamster ovary (CHO) processes through a forecast-based feeding approach that predicts multiple essential amino acids levels in the culture from easily acquired viable cell density data. Multiple cell growth behavior forecast extrapolation approaches are considered with logistic curve fitting found to be the most effective. Next, the nutrient-minimized CHO genome-scale model is combined with the growth forecast model to generate essential amino acid forecast profiles of multiple CHO batch cultures. Comparison of the forecast with the measurements suggests that this algorithm can accurately predict the concentration of most essential amino acids from cell density measurement with error mitigated by incorporating off-line amino acids concentration measurements. Finally, the forecast algorithm is applied to CHO fed-batch cultures to support amino acid feeding control to control the concentration of essential amino acids below 1-2 mM for lysine, leucine, and valine as a model over a 9-day fed batch culture while maintaining comparable growth behavior to an empirical-based culture. In turn, glycine production was elevated, alanine reduced and lactate production slightly lower in control cultures due to metabolic shifts in branched-chain amino acid degradation. With the advantage of requiring minimal measurement inputs while providing valuable and in-advance information of the system based on growth measurements, this genome model-based amino acid forecast algorithm represent a powerful and cost-effective tool to facilitate enhanced control over CHO and other mammalian cell-based bioprocesses.

White paper on high-throughput process development for integrated continuous biomanufacturing

Continuous manufacturing is an indicator of a maturing industry, as can be seen by the example of the petrochemical industry. Patent expiry promotes a price competition between manufacturing companies, and more efficient and cheaper processes are needed to achieve lower production costs. Over the last decade, continuous biomanufacturing has had significant breakthroughs, with regulatory agencies encouraging the industry to implement this processing mode. Process development is resource and time consuming and, although it is increasingly becoming less expensive and faster through high-throughput process development (HTPD) implementation, reliable HTPD technology for integrated and continuous biomanufacturing is still lacking and is considered to be an emerging field. Therefore, this paper aims to illustrate the major gaps in HTPD and to discuss the major needs and possible solutions to achieve an end-to-end Integrated Continuous Biomanufacturing, as discussed in the context of the 2019 Integrated Continuous Biomanufacturing conference. The current HTPD state-of-the-art for several unit operations is discussed, as well as the emerging technologies which will expedite a shift to continuous biomanufacturing.

Predictive macroscopic modeling of cell growth, metabolism and monoclonal antibody production: Case study of a CHO fed-batch production

We describe a systematic approach to establish predictive models of CHO cell growth, cell metabolism and monoclonal antibody (mAb) formation during biopharmaceutical production. The prediction is based on a combination of an empirical metabolic model connecting extracellular metabolic fluxes with cellular growth and product formation with mixed Monod-inhibition type kinetics that we generalized to every possible external metabolite. We describe the maximum specific growth rate as a function of the integral viable cell density (IVCD). Moreover, we also take into account the accumulation of metabolites in intracellular pools that can influence cell growth. This is possible even without identification and quantification of these metabolites as illustrated with fed-batch cultures of Chinese Hamster Ovary (CHO) cells producing a mAb. The impact of cysteine and tryptophan on cell growth and cell productivity was assessed, and the resulting macroscopic model was successfully used to predict the impact of new, untested feeding strategies on cell growth and mAb production. This model combining piecewise linear relationships between metabolic rates, growth rate and production rate together with Monod-inhibition type models for cell growth did well in predicting cell culture performance in fed-batch cultures even outside the range of experimental data used for establishing the model. It could therefore also successfully be applied for in silico prediction of optimal operating conditions.

In-silico media optimization for continuous cultures using genome scale metabolic networks: The case of CHO-K1

The cell culture is the central piece of a biotechnological industrial process. It includes upstream (e.g. media preparation, fixed costs, etc.) and downstream steps (e.g. product purification, waste disposal, etc.). In the continuous mode of cell culture, a constant flow of fresh media replaces culture fluid until the system reaches a steady state. This steady state is the standard operation mode which, under very general conditions, is a function of the ratio between the cell density and the dilution rate and depends on the media supplied to the culture. To optimize the production process it is widely accepted that the concentration of the metabolites in this media should be carefully tuned. A poor media may not provide enough nutrients to the culture, while a media too rich in nutrients may be a waste of resources because, either the cells do not use all of the available nutrients, or worse, they over-consume them producing toxic byproducts. In this study, we show how an in-silico study of a genome scale metabolic network coupled to the dynamics of a chemostat could guide the strategy to optimize the media to be used in a continuous process. Given a known media we model the concentrations of the cells in a chemostat as a function of the dilution rate. Then, we cast the problem of optimizing the production process within a linear programming framework in which the goal is to minimize the cost of the media keeping fixed the cell concentration for a given dilution rate in the chemostat. We evaluate our results in two metabolic models: first a simplified model of mammalian cell metabolism, and then in a realistic genome-scale metabolic network of mammalian cells, the Chinese hamster ovary cell line. We explore the latter in more detail given specific meaning to the predictions of the concentrations of several metabolites.

Digital Twins in Biomanufacturing

In recent years process modelling has become an established method which generates digital twins of manufacturing plant operation with the aid of numerically solved process models. This article discusses the benefits of establishing process modelling, in-house or by cooperation, in order to support the workflow from process development, piloting and engineering up to manufacturing. The examples are chosen from the variety of botanicals and biologics manufacturing thus proving the broad applicability from variable feedstock of natural plant extracts of secondary metabolites to fermentation of complex molecules like mAbs, fragments, proteins and peptides.Consistent models and methods to simulate whole processes are available. To determine the physical properties used as model parameters, efficient laboratory-scale experiments are implemented. These parameters are case specific since there is no database for complex molecules of biologics and botanicals in pharmaceutical industry, yet.Moreover, Quality-by-Design approaches, demanded by regulatory authorities, are integrated within those predictive modelling procedures. The models could be proven to be valid and predictive under regulatory aspects. Process modelling does earn its money from the first day of application. Process modelling is a key-enabling tool towards cost-efficient digitalization in chemical-pharmaceutical industries.

A unifying modelling formalism for the integration of stoichiometric and kinetic models

Current research on systems and synthetic biology relies heavily on mathematical models of the systems under study. The usefulness of such models depends on the quantity and quality of biological data, and on the availability of appropriate modelling formalisms that can gather and accommodate such data so that they can be exploited properly. Given our incomplete knowledge of biological systems and the fact that they consist of many subsystems, biological data are usually uncertain and heterogeneous. These facts hinder the use of mathematical models and computational methods. In the scope of dynamic biological systems, e.g. metabolic networks, this difficulty can be overcome by the novel modelling formalism of flexible nets (FNs). We show that an FN can combine, in a natural way, a stoichiometric model and a kinetic model. Moreover, the resulting net admits nonlinear dynamics and can be analysed in both transient and steady states.

Modelling of hollow fiber membrane bioreactor for mammalian cell cultivation using computational hydrodynamics

The hollow fiber membrane bioreactor (HFMB) has been investigated for the cultivation of mammalian Chinese hamster ovary cell expansion. The experiments were carried out in Petri's dishes and in the hollow fiber membrane bioreactor having 20 fibers (S2025 from FiberCell Systems). The approach to HFMB modelling which combines the model of cell growth kinetics and hydrodynamics has been proposed. The hydrodynamic model is made using ANSYS Fluent software. The mathematical model of HFMB was developed, allowing the study of the hydrodynamics into the lumen and the extracapillary spaces, the filtration through the membrane fiber with the cell expansion on outer membrane surface. The direct nutrient medium flow variant into the extracapillary space was suggested. Based on the numerical simulations, the optimal parameters were selected for daily changes in the medium flow-rate into the lumen space. The HFMB scaling up was performed for the larger size HFMB (60 fibers).

Early integration of Design of Experiment (DOE) and multivariate statistics identifies feeding regimens suitable for CHO cell line development and screening

In Chinese Hamster Ovary (CHO) cell lines, the establishment of the ideal fed-batch regimen promotes metabolic conditions advantageous for the bioproduction of therapeutic molecules. A tailored, cell line-specific feeding scheme is typically defined during process development (PD) activities, through the incorporation of Design of Experiment (DOE) and late stage cell culture approaches. The feeding during early stage cell line development (CLD) was a simplified "one-fits-all" design, inherited from PD lab, that didn't account for CLD needs of throughput and streamlined workflow. The "one-fits-all" efficiency was not routinely verified when novel technologies were incorporated in CLD and sub-optimal feeding carried the risk of not selecting the most desirable cell lines amenable to late stage PD. In our work we developed the DOE-feed method; a streamlined, three-stages framework for identifying efficient feeding schemes as the CLD technologies evolved. We combined early stage cell culture input data with late-stage techniques, such as statistical modelling, principal component analysis (PCA), DOE and Prediction Profiler. Novel in our DOE-feed work, we deliberately anticipated the application of statistics and approached the method development as an early-stage, continuously updated process, by building iterative datasets and statistically interpreting their responses. We capitalized on the statistical models defined by the DOE-feed methodology to study the influence of feeds on daily productivity and growth and to extrapolate feeding-schemes that improved the cell line screening. The DOE-feed became a methodology suited for CLD needs at AbbVie, and optimized the early stage screening, reduced the operational hands-on time and improved the overall workstream efficiency.

Model-based optimization of temperature and pH shift to increase volumetric productivity of a Chinese hamster ovary fed-batch process

Industrial production processes, which utilize mammalian cells for the production of therapeutic proteins are routinely designed to include temperature and pH shifts. The process conditions are shifted away from growth promoting conditions towards a state of decreased metabolic activity. This results in the extension of the cultivation duration and therefore in a higher final product concentration. Although the correct timing of these shifts is essential for peak process performance, not many studies have been investigating this topic. In this work temperature and pH shift were optimized with a mechanistic model to increase the final product concentration in comparison to an established industrial fed-batch process. The major advantage of the mechanistic in comparison to a data-driven approach lies in the reduced number of experiments, which is needed. Therefore process development is faster, which decreases the time of the product to the market. Based on the optimization, an increased final product concentration of 14% was achieved in comparison to an already established industrial fed-batch process with the same cell line. Furthermore, the space-time-yield of the process did increase in comparison, resulting in a 20% increase of the final volumetric product concentration.

Accelerating Biologics Manufacturing by Upstream Process Modelling

Intensified and accelerated development processes are being demanded by the market, as innovative biopharmaceuticals such as virus-like particles, exosomes, cell and gene therapy, as well as recombinant proteins and peptides will possess no available platform approach. Therefore, methods that are able to accelerate this development are preferred. Especially, physicochemical rigorous process models, based on all relevant effects of fluid dynamics, phase equilibrium, and mass transfer, can be predictive, if the model is verified and distinctly quantitatively validated. In this approach, a macroscopic kinetic model based on Monod kinetics for mammalian cell cultivation is developed and verified according to a general valid model validation workflow. The macroscopic model is verified and validated on the basis of four decision criteria (plausibility, sensitivity, accuracy and precision as well as equality). The process model workflow is subjected to a case study, comprising a Chinese hamster ovary fed-batch cultivation for the production of a monoclonal antibody. By performing the workflow, it was found that, based on design of experiments and Monte Carlo simulation, the maximum growth rate µmax exhibited the greatest influence on model variables such as viable cell concentration XV and product concentration. In addition, partial least squares regressions statistically evaluate the correlations between a higher µmax and a higher cell and product concentration, as well as a higher substrate consumption.

Accelerating Biologics Manufacturing by Modeling or: Is Approval under the QbD and PAT Approaches Demanded by Authorities Acceptable Without a Digital-Twin?

Innovative biologics, including cell therapeutics, virus-like particles, exosomes,recombinant proteins, and peptides, seem likely to substitute monoclonal antibodies as the maintherapeutic entities in manufacturing over the next decades. This molecular variety causes agrowing need for a general change of methods as well as mindset in the process development stage,as there are no platform processes available such as those for monoclonal antibodies. Moreover,market competitiveness demands hyper-intensified processes, including accelerated decisionstoward batch or continuous operation of dedicated modular plant concepts. This indicates gaps inprocess comprehension, when operation windows need to be run at the edges of optimization. Inthis editorial, the authors review and assess potential methods and begin discussing possiblesolutions throughout the workflow, from process development through piloting to manufacturingoperation from their point of view and experience. Especially, the state-of-the-art for modeling inred biotechnology is assessed, clarifying differences and applications of statistical, rigorousphysical-chemical based models as well as cost modeling. “Digital-twins” are described and effortsvs. benefits for new applications exemplified, including the regulation-demanded QbD (quality bydesign) and PAT (process analytical technology) approaches towards digitalization or industry 4.0based on advanced process control strategies. Finally, an analysis of the obstacles and possiblesolutions for any successful and efficient industrialization of innovative methods from processdevelopment, through piloting to manufacturing, results in some recommendations. A centralquestion therefore requires attention: Considering that QbD and PAT have been required byauthorities since 2004, can any biologic manufacturing process be approved by the regulatoryagencies without being modeled by a “digital-twin” as part of the filing documentation?

Novel column generation-based optimization approach for poly-pathway kinetic model applied to CHO cell culture

Mathematical modelling can provide precious tools for bioprocess simulation, prediction, control and optimization of mammalian cell-based cultures. In this paper we present a novel method to generate kinetic models of such cultures, rendering complex metabolic networks in a poly-pathway kinetic model. The model is based on subsets of elementary flux modes (EFMs) to generate macro-reactions. Thanks to our column generation-based optimization algorithm, the experimental data are used to identify the EFMs, which are relevant to the data. Here the systematic enumeration of all the EFMs is eliminated and a network including a large number of reactions can be considered. In particular, the poly-pathway model can simulate multiple metabolic behaviors in response to changes in the culture conditions.

We apply the method to a network of 126 metabolic reactions describing cultures of antibody-producing Chinese hamster ovary cells, and generate a poly-pathway model that simulates multiple experimental conditions obtained in response to variations in amino acid availability. A good fit between simulated and experimental data is obtained, rendering the variations in the growth, product, and metabolite uptake/secretion rates. The intracellular reaction fluxes simulated by the model are explored, linking variations in metabolic behavior to adaptations of the intracellular metabolism.

•

Novel method to model multiple states by a poly-pathway kinetic model.

•

EFMs relevant to data identified by column generation (CG)-based optimization.

•

CG optimization enables use of networks much larger than systematic enumeration.

•

A kinetic model simulates changes in metabolic rates linked to available amino acids.

•

The flux distribution of each metabolic state is visualized in the original network.

Advanced Cell and Tissue Biomanufacturing

This position paper assesses state-of-the-art advanced biomanufacturing and identifies paths forward to advance this emerging field in biotechnology and biomedical engineering, including new research opportunities and translational and corporate activities. The vision for the field is to see advanced biomanufacturing emerge as a discipline in academic and industrial communities as well as a technological opportunity to spur research and industry growth. To navigate this vision, the paths to move forward and to identify major barriers were a focal point of discussions at a National Science Foundation-sponsored workshop focused on the topic. Some of the major needs include but are not limited to the integration of specific scientific and engineering disciplines and guidance from regulatory agencies, infrastructure requirements, and strategies for reliable systems integration. Some of the recommendations, major targets, and opportunities were also outlined, including some "grand challenges" to spur interest and progress in the field based on the participants at the workshop. Many of these recommendations have been expanded, materialized, and adopted by the field. For instance, the formation of an initial collaboration network in the community was established. This report provides suggestions for the opportunities and challenges to help move the field of advanced biomanufacturing forward. The field is in the early stages of effecting science and technology in biomanufacturing with a bright and important future impact evident based on the rapid scientific advances in recent years and industry progress.

Relevant biological processes for tissue development with stem cells and their mechanistic modeling: A review

A potential alternative for tissue transplants is tissue engineering, in which the interaction of cells and biomaterials can be optimized. Tissue development in vitro depends on the complex interaction of several biological processes such as extracellular matrix synthesis, vascularization and cell proliferation, adhesion, migration, death, and differentiation. The complexity of an individual phenomenon or of the combination of these processes can be studied with phenomenological modeling techniques. This work reviews the main biological phenomena in tissue development and their mathematical modeling, focusing on mesenchymal stem cell growth in three-dimensional scaffolds.

Process Analytical Technology for Advanced Process Control in Biologics Manufacturing with the Aid of Macroscopic Kinetic Modeling

Productivity improvements of mammalian cell culture in the production of recombinant proteins have been made by optimizing cell lines, media, and process operation. This led to enhanced titers and process robustness without increasing the cost of the upstream processing (USP); however, a downstream bottleneck remains. In terms of process control improvement, the process analytical technology (PAT) initiative, initiated by the American Food and Drug Administration (FDA), aims to measure, analyze, monitor, and ultimately control all important attributes of a bioprocess. Especially, spectroscopic methods such as Raman or near-infrared spectroscopy enable one to meet these analytical requirements, preferably in-situ. In combination with chemometric techniques like partial least square (PLS) or principal component analysis (PCA), it is possible to generate soft sensors, which estimate process variables based on process and measurement models for the enhanced control of bioprocesses. Macroscopic kinetic models can be used to simulate cell metabolism. These models are able to enhance the process understanding by predicting the dynamic of cells during cultivation. In this article, in-situ turbidity (transmission, 880 nm) and ex-situ Raman spectroscopy (785 nm) measurements are combined with an offline macroscopic Monod kinetic model in order to predict substrate concentrations. Experimental data of Chinese hamster ovary cultivations in bioreactors show a sufficiently linear correlation (R² ≥ 0.97) between turbidity and total cell concentration. PLS regression of Raman spectra generates a prediction model, which was validated via offline viable cell concentration measurement (RMSE ≤ 13.82, R² ≥ 0.92). Based on these measurements, the macroscopic Monod model can be used to determine different process attributes, e.g., glucose concentration. In consequence, it is possible to approximately calculate (R² ≥ 0.96) glucose concentration based on online cell concentration measurements using turbidity or Raman spectroscopy. Future approaches will use these online substrate concentration measurements with turbidity and Raman measurements, in combination with the kinetic model, in order to control the bioprocess in terms of feeding strategies, by employing an open platform communication (OPC) network—either in fed-batch or perfusion mode, integrated into a continuous operation of upstream and downstream.

Characterizing steady states of genome-scale metabolic networks in continuous cell cultures

In the continuous mode of cell culture, a constant flow carrying fresh media replaces culture fluid, cells, nutrients and secreted metabolites. Here we present a model for continuous cell culture coupling intra-cellular metabolism to extracellular variables describing the state of the bioreactor, taking into account the growth capacity of the cell and the impact of toxic byproduct accumulation. We provide a method to determine the steady states of this system that is tractable for metabolic networks of arbitrary complexity. We demonstrate our approach in a toy model first, and then in a genome-scale metabolic network of the Chinese hamster ovary cell line, obtaining results that are in qualitative agreement with experimental observations. We derive a number of consequences from the model that are independent of parameter values. The ratio between cell density and dilution rate is an ideal control parameter to fix a steady state with desired metabolic properties. This conclusion is robust even in the presence of multi-stability, which is explained in our model by a negative feedback loop due to toxic byproduct accumulation. A complex landscape of steady states emerges from our simulations, including multiple metabolic switches, which also explain why cell-line and media benchmarks carried out in batch culture cannot be extrapolated to perfusion. On the other hand, we predict invariance laws between continuous cell cultures with different parameters. A practical consequence is that the chemostat is an ideal experimental model for large-scale high-density perfusion cultures, where the complex landscape of metabolic transitions is faithfully reproduced.

Segmented linear modeling of CHO fed-batch culture and its application to large scale production

We describe a systematic approach to model CHO metabolism during biopharmaceutical production across a wide range of cell culture conditions. To this end, we applied the metabolic steady state concept. We analyzed and modeled the production rates of metabolites as a function of the specific growth rate. First, the total number of metabolic steady state phases and the location of the breakpoints were determined by recursive partitioning. For this, the smoothed derivative of the metabolic rates with respect to the growth rate were used followed by hierarchical clustering of the obtained partition. We then applied a piecewise regression to the metabolic rates with the previously determined number of phases. This allowed identifying the growth rates at which the cells underwent a metabolic shift. The resulting model with piecewise linear relationships between metabolic rates and the growth rate did well describe cellular metabolism in the fed‐batch cultures. Using the model structure and parameter values from a small‐scale cell culture (2 L) training dataset, it was possible to predict metabolic rates of new fed‐batch cultures just using the experimental specific growth rates. Such prediction was successful both at the laboratory scale with 2 L bioreactors but also at the production scale of 2000 L. This type of modeling provides a flexible framework to set a solid foundation for metabolic flux analysis and mechanistic type of modeling. Biotechnol. Bioeng. 2017;114: 785–797. © 2016 The Authors. Biotechnology and Bioengineering Published by Wiley Periodicals, Inc.