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.     

Concomitant reduction of lactate and ammonia accumulation in fed‐batch cultures: Impact on glycoprotein production and quality

Lactate and ammonia accumulation is a major factor limiting the performance of fed‐batch strategies for mammalian cell culture processes. In addition to the detrimental effects of these by‐products on production yield, ammonia also contributes to recombinant glycoprotein quality deterioration. In this study, we tackled the accumulation of these two inhibiting metabolic wastes by culturing in glutamine‐free fed‐batch cultures an engineered HEK293 cell line displaying an improved central carbon metabolism. Batch cultures highlighted the ability of PYC2‐overexpressing HEK293 cells to grow and sustain a relatively high viability in absence of glutamine without prior adaptation to the culture medium. In fed‐batch cultures designed to maintain glucose at high concentration by daily feeding a glutamine‐free concentrated nutrient feed, the maximum lactate and ammonia concentrations did not exceed 5 and 1 mM, respectively. In flask, this resulted in more than a 2.5‐fold increase in IFNα2b titer in comparison to the control glutamine‐supplied fed‐batch. In bioreactor, this strategy led to similar reductions in lactate and ammonia accumulation and an increase in IFNα2b production. Of utmost importance, this strategy did not affect IFNα2b quality with respect to sialylation and glycoform distribution as confirmed by surface plasmon resonance biosensing and LC‐MS, respectively. Our strategy thus offers an attractive and simple approach for the development of efficient cell culture processes for the mass production of high‐quality therapeutic glycoproteins. © 2018 American Institute of Chemical Engineers Biotechnol. Prog., 34:494–504, 2018     

Automated dynamic fed‐batch process and media optimization for high productivity cell culture process development

Current industry practices for large‐scale mammalian cell cultures typically employ a standard platform fed‐batch process with fixed volume bolus feeding. Although widely used, these processes are unable to respond to actual nutrient consumption demands from the culture, which can result in accumulation of by‐products and depletion of certain nutrients. This work demonstrates the application of a fully automated cell culture control, monitoring, and data processing system to achieve significant productivity improvement via dynamic feeding and media optimization. Two distinct feeding algorithms were used to dynamically alter feed rates. The first method is based upon on‐line capacitance measurements where cultures were fed based on growth and nutrient consumption rates estimated from integrated capacitance. The second method is based upon automated glucose measurements obtained from the Nova Bioprofile FLEX® autosampler where cultures were fed to maintain a target glucose level which in turn maintained other nutrients based on a stoichiometric ratio. All of the calculations were done automatically through in‐house integration with a Delta V process control system. Through both media and feed strategy optimization, a titer increase from the original platform titer of 5 to 6.3 g/L was achieved for cell line A, and a substantial titer increase of 4 to over 9 g/L was achieved for cell line B with comparable product quality. Glucose was found to be the best feed indicator, but not all cell lines benefited from dynamic feeding and optimized feed media was critical to process improvement. Our work demonstrated that dynamic feeding has the ability to automatically adjust feed rates according to culture behavior, and that the advantage can be best realized during early and rapid process development stages where different cell lines or large changes in culture conditions might lead to dramatically different nutrient demands. Biotechnol. Bioeng. 2013; 110: 191–205. © 2012 Wiley Periodicals, Inc.     

Equalizing growth in high‐throughput small scale cultivations via precultures operated in fed‐batch mode

An often underestimated problem when working with different clones in microtiter plates and shake flask screenings is the non‐parallel and non‐equal growth of batch cultures. These growth differences are caused by variances of individual clones regarding initial biomass concentration, lag‐phase or specific growth rate. Problems arising from unequal growth kinetics are different induction points in expression studies or uneven cultivation periods at the time of harvest. Screening for the best producing clones of a library under comparable conditions is thus often impractical or even impossible. A new approach to circumvent the problem of unequal growth kinetics of main cultures is the application of fed‐batch mode in precultures in microtiter plates and shake flasks. Fed‐batch operation in precultures is realized through a slow‐release system for glucose. After differently growing cultures turn to glucose‐limited growth, they all consume the same amount of glucose due to the fixed feed profile of glucose provided by the slow‐release system. This leads to equalized growth. Inherent advantages of this method are that it is easy to use and requires no additional equipment like pumps. This new technique for growth equalization in high‐throughput cultivations is simulated and verified experimentally. The growth of distinctly inoculated precultures in microtiter plates and shake flasks could be equalized for different microorganisms such as Escherichia coli and Hansenula polymorpha. Biotechnol. Bioeng. 2009;103: 1095–1102. © 2009 Wiley Periodicals, Inc.     

Modeling kinetics of a large‐scale fed‐batch CHO cell culture by Markov chain Monte Carlo method

Markov chain Monte Carlo (MCMC) method was applied to model kinetics of a fed‐batch Chinese hamster ovary cell culture process in 5,000‐L bioreactors. The kinetic model consists of six differential equations, which describe dynamics of viable cell density and concentrations of glucose, glutamine, ammonia, lactate, and the antibody fusion protein B1 (B1). The kinetic model has 18 parameters, six of which were calculated from the cell culture data, whereas the other 12 were estimated from a training data set that comprised of seven cell culture runs using a MCMC method. The model was confirmed in two validation data sets that represented a perturbation of the cell culture condition. The agreement between the predicted and measured values of both validation data sets may indicate high reliability of the model estimates. The kinetic model uniquely incorporated the ammonia removal and the exponential function of B1 protein concentration. The model indicated that ammonia and lactate play critical roles in cell growth and that low concentrations of glucose (0.17 mM) and glutamine (0.09 mM) in the cell culture medium may help reduce ammonia and lactate production. The model demonstrated that 83% of the glucose consumed was used for cell maintenance during the late phase of the cell cultures, whereas the maintenance coefficient for glutamine was negligible. Finally, the kinetic model suggests that it is critical for B1 production to sustain a high number of viable cells. The MCMC methodology may be a useful tool for modeling kinetics of a fed‐batch mammalian cell culture process. © 2009 American Institute of Chemical Engineers Biotechnol. Prog., 2010     

Perfusion seed cultures improve biopharmaceutical fed‐batch production capacity and product quality

Volumetric productivity and product quality are two key performance indicators for any biopharmaceutical cell culture process. In this work, we showed proof‐of‐concept for improving both through the use of alternating tangential flow perfusion seed cultures coupled with high‐seed fed‐batch production cultures. First, we optimized the perfusion N‐1 stage, the seed train bioreactor stage immediately prior to the production bioreactor stage, to minimize the consumption of perfusion media for one CHO cell line and then successfully applied the optimized perfusion process to a different CHO cell line. Exponential growth was observed throughout the N‐1 duration, reaching >40 × 10⁶ vc/mL at the end of the perfusion N‐1 stage. The cultures were subsequently split into high‐seed (10 × 10⁶ vc/mL) fed‐batch production cultures. This strategy significantly shortened the culture duration. The high‐seed fed‐batch production processes for cell lines A and B reached 5 g/L titer in 12 days, while their respective low‐seed processes reached the same titer in 17 days. The shortened production culture duration potentially generates a 30% increase in manufacturing capacity while yielding comparable product quality. When perfusion N‐1 and high‐seed fed‐batch production were applied to cell line C, higher levels of the active protein were obtained, compared to the low‐seed process. This, combined with correspondingly lower levels of the inactive species, can enhance the overall process yield for the active species. Using three different CHO cell lines, we showed that perfusion seed cultures can optimize capacity utilization and improve process efficiency by increasing volumetric productivity while maintaining or improving product quality. © 2014 American Institute of Chemical Engineers Biotechnol. Prog., 30:616–625, 2014     

Feed development for fed‐batch CHO production process by semisteady state analysis

Semisteady state cultures are useful for studying cell physiology and facilitating media development. Two semisteady states with a viable cell density of 5.5 million cells/mL were obtained in CHO cell cultures and compared with a fed‐batch mode control. In the first semisteady state, the culture was maintained at 5 mM glucose and 0.5 mM glutamine. The second condition had threefold higher concentrations of both nutrients, which led to a 10% increase in lactate production, a 78% increase in ammonia production, and a 30% reduction in cell growth rate. The differences between the two semisteady states indicate that maintaining relatively low levels of glucose and glutamine can reduce the production of lactate and ammonia. Specific amino acid production and consumption indicated further metabolic differences between the two semisteady states and fed‐batch mode. The results from this experiment shed light in the feeding strategy for a fed‐batch process and feed medium enhancement. The fed‐batch process utilizes a feeding strategy whereby the feed added was based on glucose levels in the bioreactor. To evaluate if a fixed feed strategy would improve robustness and process consistency, two alternative feeding strategies were implemented. A constant volume feed of 30% or 40% of the initial culture volume fed over the course of cell culture was evaluated. The results indicate that a constant volumetric‐based feed can be more beneficial than a glucose‐based feeding strategy. This study demonstrated the applicability of analyzing CHO cultures in semisteady state for feed enhancement and continuous process improvement. © 2009 American Institute of Chemical Engineers Biotechnol. Prog., 2010     

Catabolic control of hybridoma cells by glucose and glutamine limited fed batch cultures

Substrate limited fed batch cultures were used to study growth and overflow metabolism in hybridoma cells. A glucose limited fed batch, a glutamine limited fed batch, and a combined glucose and glutamine limited red batch culture were compared with batch cultures. In all cultures mu reaches its maximum early during growth and decreases thereafter so that no exponential growth and decreases thereafter so that no exponential growth rate limiting, although the glutamine concentration (>0.085mM) was lower than reported K(s) vales and glucose was below 0.9mM; but some other nutrients (s) was the cause as verified by simulations. Slightly more cells and antibodies were produced in the combined fed batch compared with the batch culture. The specific rates for consumption of glucose and glutamine were dramatically influenced in fed batch cultures resulting in major metabolic changes. Glucose limitation decreased lactate formation, but increased glutamine consumption and ammonium formation. Glutamine limitation decreased ammonium and alanine formation of lactate, alanine, and ammonium was negligible in the dual-substrate limited fed batch culture. The efficiency of the energy metabolism increased, as judged by the increase in the cellular yield coefficient for glucose by 100% and for glutamine by 150% and by the change in the metabolic ratios lac/glc, ala/ln, and NH(x)/ln, in the combined fed culture. The data indicate that a larger proportion of consumed glutamine enters the TCA cycle through the glutamate dehydrogenase pathway, which releases more energy from glutamine than the transamination pathway. We suggest that the main reasons for these changes are decreased uptake rates of glucose and glutamine, which in turn lead to a reduction of the pyruvate pool and a restriction of the flux through glutaminase and lactate dehydrogenase. There appears to be potential for further cell growth in the dual-substrate-limited fed batch culture as judged by a comparison of mu in the different cultures. (c) 1994 John Wiley & Sons, Inc.     

Fed‐batch production of tetanus toxin by Clostridium tetani

This study deals with the effects of the initial nitrogen source (NZ Case TT) level and the protocol of glucose addition during the fed‐batch production of tetanus toxin by Clostridium tetani. An increase in the initial concentration of NZ Case TT (NZ₀) accelerated cell growth, increased the consumption of the nitrogen source as well as the final yield of tetanus toxin, which achieved the highest values (50–60 L_f/mL) for NZ₀ ≥ 50 g/L. The addition of glucose at fixed times (16, 56, and 88 h) ensured a toxin yield (～60 L_f/mL) about 33% higher than those of fed‐batch runs with addition at fixed concentration (～45 L_f/mL) and about 300% higher than those obtained in reference batch runs nowadays used at industrial scale. The results of this work promise to substantially improve the present production of tetanus toxin and may be adopted for human vaccine production after detoxification and purification. © 2009 American Institute of Chemical Engineers Biotechnol. Prog., 2010     

Model‐based optimization and pO2 control of fed‐batch Escherichia coli and Saccharomyces cerevisiae cultivation processes

A model‐based approach for optimization and cascade control of dissolved oxygen partial pressure (pO₂) and maximization of biomass in fed‐batch cultivations is presented. The procedure is based on the off‐line model‐based optimization of the optimal feeding rate profiles and the subsequent automatic pO₂ control using a proposed cascade control technique. During the model‐based optimization of the process, feeding rate profiles are optimized with respect to the imposed technological constraints (initial and maximal cultivation volume, cultivation time, feeding rate range, maximal oxygen transfer rate and pO₂ level). The cascade pO₂ control is implemented using activation of cascades for agitation, oxygen enrichment, and correction of the preoptimized feeding rate profiles. The proposed approach is investigated in two typical fed‐batch processes with Escherichia coli and Saccharomyces cerevisiae. The obtained results show that it was possible to achieve sufficiently high biomass levels with respect to the given technological constraints and to improve controllability of the investigated processes.     

法夫酵母分批发酵虾青素动力学研究

通过三联30L全自动发酵罐对虾青素产生菌法夫酵母的分批发酵动力学进行了研究,结果表明,法夫酵母的生长与限制性基质葡萄糖浓度之间符合Logistic方程,建立了细胞生长、产物合成和基质消耗随时间变化的数学模型。应用MATLAB软件对发酵动力学模型进行最优参数估计和非线性拟和,获得最大比生长速率（umax）和产物得率（Yp/x）分别为0.1829/h、0.1524g/g,虾青素分批发酵中细胞生长与产物合成属于偶联型,模型模拟计算结果和实验值能较好地吻合,动力学研究结果表明该模型能较好地反映细胞的生长、底物消耗和产物合成过程机制。     

Fed‐batch CHO cell t‐PA production and feed glutamine replacement to reduce ammonia production

Industrial therapeutic protein production has been greatly improved through fed‐batch development. In this study, improvement to the productivity of a tissue‐plasminogen activator (t‐PA) expressing Chinese hamster ovary (CHO) cell line was investigated in shake flask culture through the optimization of the fed‐batch feed and the reduction of ammonia generation by glutamine replacement. The t‐PA titer was increased from 33 mg/L under batch conditions to 250 mg/L with daily feeding starting after three days of culture. A commercially available fed‐batch feed was supplemented with cotton seed hydrolysate and the four depleted amino acids, aspartic acid, asparagine, cysteine, and tyrosine. The fed‐batch operation increased the generation of by‐products such as lactate and ammonia that can adversely affect the fed‐batch performance. To reduce the ammonia production, a glutamine‐containing dipeptide, pyruvate, glutamate, and wheat gluten hydrolysate, were investigated as glutamine substitutes. To minimize the lag phase as the cells adjusted to the new energy source, a feed glutamine replacement process was developed where the cells were initially cultured with a glutamine containing basal medium to establish cell growth followed by feeding with a feed containing the glutamine substitutes. This two‐step feed glutamine replacement process not only reduced the ammonia levels by over 45% but, in the case of using wheat gluten hydrolysate, almost doubled the t‐PA titer to over 420 mg/L without compromising the t‐PA product quality or glycosylation pattern. The feed glutamine replacement process combined with optimizing other feed medium components provided a simple, practical, and effective fed‐batch strategy that could be applied to the production of other recombinant therapeutic proteins. © 2012 American Institute of Chemical Engineers Biotechnol. Prog., 2013     

Multi-objective optimization of glycopeptide antibiotic production in batch and fed batch processes

Fermentation optimization involves potentially conflicting multiple objectives such as product concentration and production media cost. Simultaneous optimization of these objectives would result in a multiobjective optimization problem, which is characterized by a set of multiple solutions, knows as pareto optimal solutions. These solutions gives flexibility in evaluating the trade-offs and selecting the most suitable operating policy. Here, ε-constraint approach was used to generate the pareto solutions for two objectives: product concentration and product per unit cost of media, for batch and fed batch operations using process model for Amycolatopsis balhimycina, a glycopeptide antibiotic producer. This resulted in a set of several pareto optimal solutions with the two objectives ranging from (0.75 g l⁻¹, 3.97 g $^-1) to (0.44 g l⁻¹, 5.19 g $^-1) for batch and from (1.5 g l⁻¹, 5.46 g $^-1) to (1.1 g l⁻¹, 6.34 g $^-1) for fed batch operations. One pareto solution each for batch and for fed batch mode was experimentally validated.     

Elucidation of metabolism in hybridoma cells grown in fed‐batch culture by genome‐scale modeling

Genome‐scale modeling of mouse hybridoma cells producing monoclonal antibodies (mAb) was performed to elucidate their physiological and metabolic states during fed‐batch cell culture. Initially, feed media nutrients were monitored to identify key components among carbon sources and amino acids with significant impact on the desired outcome, for example, cell growth and antibody production. The monitored profiles indicated rapid assimilation of glucose and glutamine during the exponential growth phase. Significant increase in mAb concentration was also observed when glutamine concentration was controlled at 0.5 mM as a feeding strategy. Based on the reconstructed genome‐scale metabolic network of mouse hybridoma cells and fed‐batch profiles, flux analysis was then implemented to investigate the cellular behavior and changes in internal fluxes during the cell culture. The simulated profile of the cell growth was consistent with experimentally measured specific growth rate. The in silico simulation results indicated (i) predominant utilization of glycolytic pathway for ATP production, (ii) importance of pyruvate node in metabolic shifting, and (iii) characteristic pattern in lactate to glucose ratio during the exponential phase. In future, experimental and in silico analyses can serve as a promising approach to identifying optimal feeding strategies and potential cell engineering targets as well as facilitate media optimization for the enhanced production of mAb or recombinant proteins in mammalian cells. Biotechnol. Bioeng. 2009;102: 1494–1504. © 2008 Wiley Periodicals, Inc.     

A robust feeding strategy to maintain set‐point glucose in mammalian fed‐batch cultures when input parameters have a large error

Industrial CHO cell cultures run under fed‐batch conditions are required to be controlled in particular ranges of glucose, while glucose is constantly consumed and must be replenished by a feed. The most appropriate feeding rate is ideally stoichiometric and adaptive in nature to balance the dynamically changing rate of glucose consumption. However, high errors in biomass and glucose estimation as well as limited knowledge of the true metabolic state challenge the control strategy. In this contribution, we take these errors into account and simulate the output with uncertainty trajectories in silico in order to control glucose concentration. Other than many control strategies, which require parameter estimation, our assumptions are founded on two pillars: (i) first principles and (ii) prior knowledge about the variability of fed‐batch CHO cell culture. The algorithm was exposed to an in‐silico Design of Experiments (DoE), in which variations of parameters were changed simultaneously, such as clone‐specific behavior, precision of equipment and desired control range used. The results demonstrate that our method achieved the target of holding the glucose concentration within an acceptable range. A robust and sufficient level of control could be demonstrated even with high errors for biomass or metabolic state estimation. In a time where blockbuster drugs are queuing up for time slots of their production, this transferable control strategy that is independent of tedious establishment runs may be a decisive advantage for rapid implementation during technology transfer and scale up and decrease in campaign change over time. © 2017 American Institute of Chemical Engineers Biotechnol. Prog., 33:317–336, 2017     

Rapid media transition: An experimental approach for steady state analysis of metabolic pathways

Commonly steady state analysis of microbial metabolism is performed under well defined physiological conditions in continuous cultures with fixed external rates. However, most industrial bioprocesses are operated in fed‐batch mode under non‐stationary conditions, which cannot be realized in chemostat cultures. A novel experimental setup—rapid media transition—enables steady state perturbation of metabolism on a time scale of several minutes in parallel to operating bioprocesses. For this purpose, cells are separated from the production process and transferred into a lab‐scale stirred‐tank reactor with modified environmental conditions. This new approach was evaluated experimentally in four rapid media transition experiments with Escherichia coli from a fed‐batch process. We tested the reaction to different carbon sources entering at various points of central metabolism. In all cases, the applied substrates (glucose, succinate, acetate, and pyruvate) were immediately utilized by the cells. Extracellular rates and metabolome data indicate a metabolic steady state during the short‐term cultivation. Stoichiometric analysis revealed distribution of intracellular fluxes, which differs drastically subject to the applied carbon source. For some reactions, the variation of flux could be correlated to changes of metabolite concentrations. © 2009 American Institute of Chemical Engineers Biotechnol. Prog., 2010     

Nisin production of Lactococcus lactis N8 with hemin‐stimulated cell respiration in fed‐batch fermentation system

In this study, nisin production of Lactococcus lactis N8 was optimized by independent variables of glucose, hemin and oxygen concentrations in fed‐batch fermentation in which respiration of cells was stimulated with hemin. Response surface model was able to explain the changes of the nisin production of L. lactis N8 in fed‐batch fermentation system with high fidelity (R² 98%) and insignificant lack of fit. Accordingly, the equation developed indicated the optimum parameters for glucose, hemin, and dissolved oxygen were 8 g L^?1 h^?1, 3 μg mL^?1 and 40%, respectively. While 1711 IU mL^?1 nisin was produced by L. lactis N8 in control fed‐batch fermentation, 5410 IU mL^?1 nisin production was achieved within the relevant optimum parameters where the respiration of cell was stimulated with hemin. Accordingly, nisin production was enhanced 3.1 fold in fed‐batch fermentation using hemin. In conclusion the nisin production of L. lactis N8 was enhanced extensively as a result of increasing the biomass by stimulating the cell respiration with adding the hemin in the fed‐batch fermentation. © 2015 American Institute of Chemical Engineers Biotechnol. Prog., 31:678–685, 2015     

Software sensor design considering oscillating conditions as present in industrial scale fed‐batch cultivations

Investigations of inhomogeneous dynamics in industrial‐scale bioreactors can be realized in laboratory simulators. Such studies will be improved by on line observation of the growth of microorganisms and their product synthesis at oscillating substrate availability which represents the conditions in industrial‐scale fed‐batch cultivations. A method for the kinetic monitoring of such processes, supported by on line measurements accessible in industrial practice, is proposed. It consists of a software sensor (SS) system composed of a cascade structure. Process kinetics are simulated in models with a structure including time‐varying yield coefficients. SSs for measured variable kinetics have classical structures. The SS design of unmeasured variables is realized after a linear transformation using a logarithmic function. For these software sensors, a tuning procedure is proposed, at which an arbitrary choice of one tuning parameter value that guarantees stability of the monitoring system allows the calculation of the optimal values of six parameters. The effectiveness of the proposed monitoring approach is demonstrated with experimental data from a glucose‐limited fed‐batch process of Bacillus subtilis in a laboratory two‐compartment scale down reactor which tries to mimic the conditions present in industrial scale nutrient‐limited fed‐batch cultivations. Biotechnol. Bioeng. 2013; 110: 1945–1955. © 2013 Wiley Periodicals, Inc.