Inclusion of maintenance energy improves the intracellular flux predictions of CHO

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.

Characterising the metabolic rewiring of extremely slow growing Komagataella phaffii

Retentostat cultivations have enabled investigations into substrate-limited near-zero growth for a number of microbes. Quantitative physiology at these near-zero growth conditions has been widely discussed, yet characterisation of the fluxome is relatively under-reported. We investigated the rewiring of metabolism in the transition of a recombinant protein-producing strain of Komagataella phaffii to glucose-limited near-zero growth rates. We used cultivation data from a 200-fold range of growth rates and comprehensive biomass composition data to integrate growth rate dependent biomass equations, generated using a number of different approaches, into a K. phaffii genome-scale metabolic model. Here, we show that a non-growth-associated maintenance value of 0.65 mmol ATP g CDW - 1 h - 1 and a growth-associated maintenance value of 108 mmol ATP g CDW - 1 lead to accurate growth rate predictions. In line with its role as energy source, metabolism is rewired to increase the yield of ATP per glucose. This includes a reduction of flux through the pentose phosphate pathway, and a greater utilisation of glycolysis and the TCA cycle. Interestingly, we observed activity of an external, non-proton translocating NADH dehydrogenase in addition to the malate-aspartate shuttle. Regardless of the method used for the generation of biomass equations, a similar, yet different, growth rate dependent rewiring was predicted. As expected, these differences between the different methods were clearer at higher growth rates, where the biomass equation provides a much greater constraint than at slower growth rates. When placed on an increasingly limited glucose diet, the metabolism of K. phaffii adapts, enabling it to continue to drive critical processes sustaining its high viability at near-zero growth rates.

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 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.

Genome-Scale Metabolic Model's multi-objective solving algorithm based on the inflexion point of Pareto front including maximum energy utilization and its application in A.niger DS03043

The solution of genome-scale metabolic model (GSMM) directly affects the simulation accuracy of the metabolic process in digital cells. Single-objective optimization methods, such as Flux Balance Analysis (FBA) which is widely used in solving GSMM, have limitations when simulating actual biological processes, which leads to unrealistic results due to other biological constraints being ignored. A novel multi-objective Differential Evolution algorithm based on general FBA (i.e., DEFBA) is hence proposed to solve GSMM. First, in accordance with to the assumption that cells minimize resource consumption and maximize resource utilization, the maximum specific growth rate and the minimum cellular production rate of ATP, NADPH, and NADH are defined as the multi-objective functions of DEFBA. Second, FBA is used to produce the initial individuals of DEFBA by changing the upper bound of biomass reaction in GSMM. Third, mutation and selection operations help in generating new individuals in the solution space to search the Pareto front. Finally, the optimal solution is selected by analyzing the inflexion point of the Pareto front. In DEFBA, multi-objective technology and optimal solution judging technology can introduce the biological constraints into the GSMM solving method, such that the solution can be more consistent with the essential biological mechanism. DEFBA is applied to solve Aspergillus niger's GSMM. The improved results show that DEFBA can be an effective general solving algorithm for GSMM. This article is protected by copyright. All rights reserved.