Dynamic multiscale metabolic network modeling of Chinese hamster ovary cell metabolism integrating N-linked glycosylation in industrial biopharmaceutical manufacturing.

Experimental and modeling work, described in this article, is focused on the metabolic pathway of Chinese hamster ovary (CHO) cells, which are the preferred expression system for monoclonal antibody protein production. CHO cells are one of the primary hosts for monoclonal antibodies production, which have extensive applications in multiple fields like biochemistry, biology and medicine. Here, an approach to explain cellular metabolism with in silico modeling of a microkinetic reaction network is presented and validated with unique experimental results. Experimental data of 25 different fed-batch bioprocesses included the variation of multiple process parameters, such as pH, agitation speed, oxygen and CO2 content, and dissolved oxygen. A total of 151 metabolites were involved in our proposed metabolic network, which consisted of 132 chemical reactions that describe the reaction pathways, and include 25 reactions describing N-glycosylation and additional reactions for the accumulation of the produced glycoforms. Additional eight reactions are considered for accumulation of the N-glycosylation products in the extracellular environment and one reaction to correlate cell degradation. The following pathways were considered: glycolysis, pentose phosphate pathway, nucleotide synthesis, tricarboxylic acid cycle, lipid synthesis, protein synthesis, biomass production, anaplerotic reactions, and membrane transport. With the applied modeling procedure, different operational scenarios and fed-batch techniques can be tested.

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.