Computer programs for modeling mammalian cell batch and fed-batch cultures using logistic equations

A MATLAB^® toolbox was developed for applying the logistic modeling approach to mammalian cell batch and fed-batch cultures. The programs in the toolbox encompass sensitivity analyses and simulations of the logistic equations in addition to cell specific rate estimation. The toolbox was first used to generate time courses of the sensitivity equations for characterizing the relationship between the logistic variable and the model parameters. Subsequently, the toolbox was used to describe CHO cell data from batch and fed-batch mammalian cell cultures. Cell density, product, glucose, lactate, glutamine, and ammonia data were analyzed for the batch culture while fed-batch analysis included cell density and product concentration. In all instances, experimental data were well described by the logistic equations and the resulting specific rate profiles were representative of the underlying cell physiology. The 6-variable batch culture data set was also used to compare the logistic specific rates with those from polynomial fitting and discrete derivative methods. The polynomial specific rates grossly misrepresented cell behavior in the initial and final stages of culture while those based on discrete derivatives had high variability due to computational artifacts. The utility of logistic specific rates to guide process development activities was demonstrated using specific protein productivity versus growth rate trajectories for the 3 cultures examined in this study. Overall, the computer programs developed in this study enable rapid and robust analysis of data from mammalian cell batch and fed-batch cultures which can help process development and metabolic flux estimation.     

Re-interpretation of the logistic equation for batch microbial growth in relation to Monod kinetics

Aims:  To determine the underlying substrate utilization mechanism in the logistic equation for batch microbial growth by revealing the relationship between the logistic and Monod kinetics. Also, to determine the logistic rate constant in terms of Monod kinetic constants.
Methods and Results:  The logistic equation used to describe batch microbial growth was related to the Monod kinetics and found to be first-order in terms of the substrate and biomass concentrations. The logistic equation constant was also related to the Monod kinetic constants. Similarly, the substrate utilization kinetic equations were derived by using the logistic growth equation and related to the Monod kinetics.
Conclusion:  It is revaled that the logistic growth equation is a special form of the Monod growth kinetics when substrate limitation is first-order with respect to the substrate concentration. The logistic rate constant ( k ) is directly proportional to the maximum specific growth rate constant ( μ _m) and initial substrate concentration ( S ₀) and also inversely related to the saturation constant ( K _s).
Significance and Impact of the Study:  The semi-empirical logistic equation can be used instead of Monod kinetics at low substrate concentrations to describe batch microbial growth using the relationship between the logistic rate constant and the Monod kinetic constants.     

Transcriptional profiling of apoptotic pathways in batch and fed-batch CHO cell cultures

Chinese Hamster ovary (CHO) cells are regarded as one of the "work-horses" for complex biotherapeutics production. In these processes, loss in culture viability occurs primarily via apoptosis, a genetically controlled form of cellular suicide. Using our "in-house" developed CHO cDNA array and a mouse oligonucleotide array for time profile expression analysis of batch and fed-batch CHO cell cultures, the genetic circuitry that regulates and executes apoptosis induction were examined. During periods of high viability, most pro-apoptotic genes were down-regulated but upon loss in viability, several early pro-apoptotic signaling genes were up-regulated. At later stages of viability loss, we detected late pro-apoptotic effector genes such as caspases and DNases being up-regulated. This sequential regulation of apoptotic genes showed that DNA microarrays could be used as a tool to study apoptosis. We found that in batch and fed-batch cultures, apoptosis signaling occurred primarily via death receptor- and mitochondria-mediated signaling pathways rather than endoplasmic reticulum-mediated signaling. These insights provide a greater understanding of the regulatory circuitry of apoptosis during cell culture and allow for subsequent targeting of relevant apoptosis signaling genes to prolong cell culture.     

Effects of dissolved oxygen levels and the role of extra- and intracellular amino acid concentrations upon the metabolism of mammalian cell lines during batch and continuous cultures

The effects of dissolved oxygen and the concentration of essential amino acids upon the metabolism of two mammalian cell lines (rCHO producing human active (t-PA) and a mouse-mouse hybridoma) were investigated in batch, chemostat, and perfusion cultures. Intracellular amino acid concentrations were measured for both cell lines during repeated batch cultures and the K_S-values for the essential amino acids were calculated using Monod equations via computer simulation. The K_S-values were in the range of 10 mmol L⁻¹ and the pool of most intracellular amino acids remained constant at about 10–100 fold higher in concentration than in the medium. No significant differences were observed between the hybridoma and CHO cell. The specific nutrient uptake rates corresponded with the cell specific growth rate and the effects of reduced dissolved oxygen concentrations only became evident when the DO dropped below 5% of air saturation (critical concentration below 1%). Nevertheless, a correlation between nutrient concentration and specific oxygen uptake was detected. This revised version was published online in June 2006 with corrections to the Cover Date.     

Towards a metabolic and isotopic steady state in CHO batch cultures for reliable isotope-based metabolic profiling

Attaining metabolic and isotopic balanced growth is one critical condition for physiological studies using isotope-labeled tracers, but is very difficult to obtain in batch culture due to the extensive metabolite exchange with the surrounding medium and related physiological changes. In the present study, we investigated metabolic and isotopic behavior of CHO cells in differently designed media. We observed that the assumption of balanced cell growth cannot be justified in batch culture of CHO cells directly using conventional, commercially available media. By systematically redesigning media composition and characterizing metabolic steady state based on mass balances and measurement of labeling dynamics, we achieved balanced cell growth for the main cellular substrates in CHO cells. This was done in a step-by-step analysis of growth and primary metabolism of CHO cells with the use of [U-¹³C]glucose feeding and adjusting concentrations of amino acids in the growth medium. The optimized media obtained at the end of the study provide balanced growth and isotopic steady state or at least asymptotic steady state. As a result, we established a platform to conduct isotope-based physiological studies of mammalian systems more reliably and therefore well suited for later use in metabolic profiling of mammalian systems such as ¹³C-labeled metabolic flux analysis.     

批次和流加培养中国仓鼠卵巢工程细胞的细胞周期相关基因转录谱分析

以表达人重组尿激酶原中国仓鼠卵巢 (CHO) 工程细胞系11G-S为研究对象,运用基因芯片技术比较了CHO工程细胞在批次及流加培养不同生长阶段基因表达水平的差异,在此基础上采用Genmapp软件,同时结合已知的细胞周期信号通路图,着重分析了批次及流加培养CHO工程细胞的细胞周期调控基因转录谱差异。在基因芯片涉及的19 191个目标基因中,批次和流加培养不同生长阶段CHO工程细胞的下调表达的基因数量多于上调表达基因数目;两种培养模式下的基因差异表达有着明显的不同,尤其是在细胞生长的衰退期,流加培养CHO工程细胞中下调表达的基因数量明显多于批次培养。有关调控细胞周期关键基因的转录谱分析表明,CHO工程细胞主要是通过下调表达CDKs、Cyclin及CKI家族中的Cdk6、Cdk2、Cdc2a、Ccne1、Ccne2基因及上调表达Smad4基因,来达到调控细胞增殖及维持自身活力的目的。     

Specific inhibition of caspase-8 and -9 in CHO cells enhances cell viability in batch and fed-batch cultures

In an attempt to investigate the molecular mechanism that leads to apoptotic death in Chinese hamster ovary (CHO) cells in batch and fed-batch cultures, we cloned caspase-2, -8 and -9 from a CHO cDNA library. Recombinant Chinese hamster caspase-2 and -9 expressed in Escherichia coli show highest activities towards commercial peptide substrates Ac-VDVAD-pNA and Ac-LEHD-pNA, the designated commercial substrates for human caspase-2 and -9, respectively. However, Chinese hamster caspase-8 shows a broad specificity profile and it cleaves the caspase-9 substrate more efficiently than it cleaves the caspase-8 substrate. The commercially available fluoromethyl ketone type of caspase inhibitors, such as Z-LEHD-fmk, Z-IETD-fmk, Z-VDVAD-fmk and Z-DEVD-fmk, were shown to completely lack specificity in inhibiting these caspases. The reversible aldehyde form of inhibitors for human caspase-8 and -9, Ac-LEHD-CHO and Ac-IETD-CHO, are equally efficient in inhibiting Chinese hamster caspase-8. Therefore, the wildly used method of utilizing the "caspase-specific" inhibitors to track the role of individual caspases in dying cells can be inaccurate and thus misleading. As an alternative, we stably expressed dominant negative (DN) mutants of Chinese hamster caspase-2, -8 and -9 to specifically inhibit these enzymes in CHO cells. Our results showed that inhibition of either endogenous caspase-8 or caspase-9 enhanced the viability of the CHO cells in both batch and fed-batch suspension cultures, but the inhibition of caspase-2 had minimal effects. These results suggest that caspase-8 and -9 are possibly involved in the apoptotic cell death in batch and fed-batch cultures of CHO cells, whereas caspase-2 is not. These findings can be valuable in the development of strategies for genetically engineering CHO cells to counter apoptotic death in batch and fed-batch cultures.     

Influence of intracellular nucleotide and nucleotide sugar contents on recombinant interferon-gamma glycosylation during batch and fed-batch cultures of CHO cells

Both the macroheterogeneity of recombinant human IFN-gamma produced by CHO cells and intracellular levels of nucleotides and sugar nucleotides, have been characterized during batch and fed-batch cultures carried out in different media. Whereas PF-BDM medium was capable to maintain a high percentage of the doubly- glycosylated glycoforms all over the process, mono-glycosylated and non-glycosylated forms increased during the batch culture using SF-RPMI medium. Intracellular level of UTP was higher in PF-BDM all over the batch culture compared to the SF-RPMI process. UDP-Gal accumulated only during the culture performed in PF-BDM medium, probably as a consequence of the reduced UDP-Glc synthesis flux in SF-RPMI medium. When the recombinant CHO cells were cultivated in fed-batch mode, the UTP level remained at a relatively high value in serum-containing RPMI and its titer increased during the fed-phase indicating an excess of biosynthesis. Besides, an accumulation of UDP-Gal occurred as well. Those results all together indicate that UTP and UDP-Glc syntheses in CHO cells cultivated in SF-RPMI medium in batch process, could be limiting during the glycosylation processes of the recombinant IFN-gamma. At last, the determination of the energetic status of the cells over the three studied processes suggested that a relationship between the adenylate energy charge and the glycosylation macroheterogeneity of the recombinant IFN-gamma may exist.     

Fermentation kinetics of recombinant yeast in batch and fed-batch cultures

Fed-batch cultures of recombinant microorganisms have attracted attention as they can separate cell growth stage from cloned-gene expression phase during fermentations. In this work, the effect of different glucose feeding strategies on cell growth and cloned gene expression was studied during aerobic fed-batch fermentations of recombinant yeast, containing the plasmid pRB58. The plasmid contains the yeast SUC2 gene, which codes for the enzyme invertase. Some feeding policies resulted in a constant glucose concentration inside the fermentor, while others deliberately introduced a cyclic variation. The cell mass yield was found to be higher at low glucose concentrations, thus indicating a shift to the more energy-efficient respiratory pathway. The SUC2 gene expression was derepressed at glucose levels below 2 g/L. The response of specific invertase activity to changes in the medium glucose concentration was found to be almost immediate.     

Controlling glycation of recombinant antibody in fed-batch cell cultures

Protein glycation is a non-enzymatic glycosylation that can occur to proteins in the human body, and it is implicated in the pathogenesis of multiple chronic diseases. Glycation can also occur to recombinant antibodies during cell culture, which generates structural heterogeneity in the product. In a previous study, we discovered unusually high levels of glycation (>50%) in a recombinant monoclonal antibody (rhuMAb) produced by CHO cells. Prior to that discovery, we had not encountered such high levels of glycation in other in-house therapeutic antibodies. Our goal here is to develop cell culture strategies to decrease rhuMAb glycation in a reliable, reproducible, and scalable manner. Because glycation is a post-translational chemical reaction between a reducing sugar and a protein amine group, we hypothesized that lowering the concentration of glucose--the only source of reducing sugar in our fed-batch cultures--would lower the extent of rhuMAb glycation. When we decreased the supply of glucose to bioreactors from bolus nutrient and glucose feeds, rhuMAb glycation decreased to below 20% at both 2-L and 400-L scales. When we maintained glucose concentrations at lower levels in bioreactors with continuous feeds, we could further decrease rhuMAb glycation levels to below 10%. These results show that we can control glycation of secreted proteins by controlling the glucose concentration in the cell culture. In addition, our data suggest that rhuMAb glycation occurring during the cell culture process may be approximated as a second-order chemical reaction that is first order with respect to both glucose and non-glycated rhuMAb. The basic principles of this glycation model should apply to other recombinant proteins secreted during cell culture.     

Apoptosis in CHO cell batch cultures: examination by flow cytometry

Chinese hamster ovary cells grown under conditions which are optimal for the production of a genetically engineered protein in batch culture, lose significant viability shortly after entering the stationary phase. This cell death was investigated morphologically and was found to be almost exclusively via apoptosi. Furthermore, cells were analyzed by flow cytometry using a fluorescent DNA end-labeling assay to label apoptotic cells, in conjunction with cell cycle analysis using propidium iodide. Apoptotic cells could be detected by this method, and by the radioactive end-labeling of extracted DNA, on all days of culture from day 1 to day 7; however, the degree of apoptotic cell death increased dramatically when the cells entered the stationary phase, rising to 50–60% of the total cell number at the termination of the culture. Flow cytometric analysis showed that the majority of cells underwent apoptosis whilst in G₁/G₀ and formed an apoptotic population with high DNA FITC end-labeling and hypodiploid propidium iodide binding. Additionally, the ability or inability to secrete specific protein products did not appear to interfere with the development of the apoptotic population with time.     

Caspase activation, sialidase release and changes in sialylation pattern of recombinant human erythropoietin produced by CHO cells in batch and fed-batch cultures

The activation of caspases represents a crucial turning point during a batch or a fed-batch culture of mammalian cells. It not only affects the quantity but also the quality of the recombinant glycoprotein produced. In this study, the activation of various caspases, the release of intracellular sialidase and the changes in sialylation pattern of a recombinant product, erythropoietin (EPO), in the culture medium were analyzed in both batch and fed-batch cultures. In both setups, all caspase activities peaked at the culture time point at which decline of cell viability was most pronounced. In addition, the release of intracellular lactate dehydrogenase (LDH) was also tracked during these cultures. The increase in LDH activity in the medium coincided with the increase of intracellular caspase activities, the release of sialidase and the observed decline in cell viability, suggesting that the LDH activity in the medium can be used as an indirect indicator of apoptotic cell death in bioreactors. Isoelectric focusing (IEF) coupled with double blotting was employed to analyze the changes in sialylation pattern of the recombinant EPO. This assay resulted in a prompt resolution of secreted EPO isoforms in a time course format. IEF profile of batch culture showed relatively consistent product sialylation compared to fed-batch culture, which showed gradual band shifts towards the isoforms with fewer sialic acid as the culture progressed. These data provided a guideline for the optimal time point to terminate the culture and collect products in batch and fed-batch cultures.Kok Hwee Chuan and Sing Fee Lim contributed equally to this work.     

Improvement of rifemycins production by Amycolatopsis mediterranei in batch and fed-batch cultures

The production of rifamycins B and SV using glucose as main C-source by Amycolatopsis mediterranei in batch and fed-batch culture was investigated. Fed-batch culture using glucose as mono feeding substrate either in the form of pulse addition, in case of shake flask, or with constant feeding rate, in bioreactor level, proved to be an alternative production system with a significant increase in both volumetric and specific antibiotic production. The maximal concentrations of about 1146 mg/l and 2500 mg/l of rifamycins B and SV, respectively, was obtained in fed-batch culture in bioreactor level under non-oxygen limitation. On the other hand, the rate of rifamycins production was increased from 6.58 to 12.13 mg/l x h for rifamycin B and from 9.47 to 31.83 mg/l x h for rifamycin SV on the bioprocess transfer and improvement from the conventional batch cultivation in shake flask to fed-batch cultivation in stirred tank bioreactor.     

Metabolic flux-based modeling of mAb production during batch and fed-batch operations

This paper proposes mathematical models that predict the physiology, growth behavior and productivity of hybridoma cells in both batch and fed-batch systems. Murine hybridoma 130-8F producing anti-F-glycoprotein monoclonal antibody was employed as a model system. A systematic approach based on metabolic flux analysis (MFA) was utilized to yield a dynamic model for predicting the concentration of significant metabolites over time. Correlation analysis was performed to formulate a Biomass Model for predicting cell concentration and viability as a function of the extracellular metabolite concentrations. The coefficients of the model equation were estimated by employing the Metropolis–Hastings algorithm. The Metabolites Model was combined with the Biomass Model to get an Integrated Model capable of predicting concentration values for substrates, extracellular metabolites, and viable and dead cell concentration by utilizing only starting concentrations as input. The prediction accuracy of the model was tested by using experimental data.     

Patterns of RNA synthesis in batch and cyclic fed-batch cultures of tylosin-producingStreptomyces fradiae

Ribonucleic acid (RNA) synthesis was studied in batch and cyclic fed-batch cultures of tylosin-producingStreptomyces fradiae. In batch culture, net RNA synthesis increased with increasing growth rate then decreased as the fermentation progressed. A cyclic response pattern of net RNA synthesis was observed when glucose and glutamate solutions were cyclically fed to the batch culture; RNA synthesis rate increased sharply to a maximum level following the feeding cycles then decreased to lower values between the maxima. Data obtained from an extended cyclic fed-batch culture showed that cyclic feeding of the key nutrients allowed RNA and tylosin syntheses to continue during a period when such activities had ceased in the control batch culture.     

基质穿膜传递过程和Monod方程的耦合模型

在分析基质进入细胞穿膜传质机理的基础上,提出了相应的简单传质模型。以此讨论了传递过程对Monod方程的影响,得出了传递过程不影响Monod方程的形式,但影响其动力学参数的结论。这和文献结果和实验数据一致。     

Robust and efficient design of experiments for the Monod model

In this paper the problem of designing experiments for the Monod model, which is frequently used in microbiology, is studied. The model is defined implicitly by a differential equation and has numerous applications in microbial growth kinetics, environmental research, pharmacokinetics, and plant physiology. The designs presented so far in the literature are local optimal designs, which depend sensitively on a preliminary guess of the unknown parameters, and are for this reason in many cases not robust with respect to their misspecification. Uniform designs and maximin optimal designs are considered as a strategy to obtain robust and efficient designs for parameter estimation. In particular, standardized maximin D- and E-optimal designs are determined and compared with uniform designs, which are usually applied in these microbiological models. It is demonstrated that maximin optimal designs are substantially more efficient than uniform designs. Parameter variances can be decreased by a factor of two by simply sampling at optimal times during the experiment. Moreover, the maximin optimal designs usually provide the possibility for the experimenter to check the model assumptions, because they have more support points than parameters in the Monod model.