The avian cell line AGE1.CR.pIX characterized by metabolic flux analysis
- Verena Lohr†1, 2Email author,
- Oliver Hädicke†1,
- Yvonne Genzel1,
- Ingo Jordan2,
- Heino Büntemeyer3,
- Steffen Klamt1 and
- Udo Reichl1, 4
© Lohr et al.; licensee BioMed Central Ltd. 2014
Received: 31 October 2013
Accepted: 16 July 2014
Published: 30 July 2014
In human vaccine manufacturing some pathogens such as Modified Vaccinia Virus Ankara, measles, mumps virus as well as influenza viruses are still produced on primary material derived from embryonated chicken eggs. Processes depending on primary cell culture, however, are difficult to adapt to modern vaccine production. Therefore, we derived previously a continuous suspension cell line, AGE1.CR.pIX, from muscovy duck and established chemically-defined media for virus propagation.
To better understand vaccine production processes, we developed a stoichiometric model of the central metabolism of AGE1.CR.pIX cells and applied flux variability and metabolic flux analysis. Results were compared to literature dealing with mammalian and insect cell culture metabolism focusing on the question whether cultured avian cells differ in metabolism. Qualitatively, the observed flux distribution of this avian cell line was similar to distributions found for mammalian cell lines (e.g. CHO, MDCK cells). In particular, glucose was catabolized inefficiently and glycolysis and TCA cycle seem to be only weakly connected.
A distinguishing feature of the avian cell line is that glutaminolysis plays only a minor role in energy generation and production of precursors, resulting in low extracellular ammonia concentrations. This metabolic flux study is the first for a continuous avian cell line. It provides a basis for further metabolic analyses to exploit the biotechnological potential of avian and vertebrate cell lines and to develop specific optimized cell culture processes, e.g. vaccine production processes.
KeywordsAvian cell line AGE1.CR.pIX Biomass composition Flux variability analysis Metabolic network modeling Glutaminolysis
Novel therapeutic options have been made possible with production of animal cell culture-derived biopharmaceuticals. An enormous amount of insight has been published on production of recombinant proteins (including antibodies) since the first products have been licensed for application in humans 1986 and 1987. CHO cells are the main substrate for recombinant proteins and therefore most metabolic studies describe these cells [1, 2]. Regarding vaccine manufacturing, literature on properties and metabolism of host cells during growth and virus replication as well as process design and optimization is less abundant and not as focused. One reason is that safety considerations and execution of pre-clinical and clinical trials are paramount for vaccines because this medication usually is administered to healthy recipients. Another reason is the multitude of process options and host cell systems used for propagation of pathogens. For historic reasons, primary cultures (today usually embryonated chicken eggs or chicken embryo fibroblasts), diploid cell strains (human embryonic WI-38 and MRC-5), and few continuous cell lines (i.e. macaque-derived Vero cells or MDCK cells) are associated with the lowest possible risks to the vaccine recipients. Production with such cellular substrates often is performed in the presence of calf serum and adherent culture with little room for optimization. Furthermore, standardization of experiments is difficult with primary chicken material.
Several virus strains and viral vectors, e.g. Modified Vaccinia Virus Ankara (MVA), only replicate efficiently in avian cells . Accordingly, primary chicken fibroblasts are still a commonly accepted production substrate. As vectored vaccines are gaining increasing importance , the potential of continuous suspension cells that proliferate in chemically-defined media are attracting more and more attention. Therefore, continuous avian cell lines are being investigated as a viable option to replace primary material. One possible new cell candidate that meets these criteria is the avian designer cell line AGE1.CR.pIX (in the following: CR.pIX) that proliferates in fully scalable suspension culture and is adapted to growth in a chemically-defined medium [5, 6]. The latter property is an advantage for metabolic flux analysis as no unknown or complex components such as animal sera or hydrolysates that complicate carbon and nitrogen balance closure are present.
Metabolic models have been developed for a variety of cell lines to study physiological states or changes in metabolism as a response to different cultivation strategies, medium composition or stimuli like virus infection or accumulation of toxins. For example, the human cell lines HEK293 [7, 8] and AGE1.HN [9, 10], other mammalian cells like CHO , BHK , MDCK [13, 14] and hybridoma cells [15, 16] or insect cell lines like Sf9 [17, 18] have been studied because of their relevance for the production of biopharmaceuticals. For all these cell lines, the two main pathways for energy generation and precursor supply were found to be glucose catabolism via glycolysis and glutamine catabolism via the TCA cycle, referred to as glutaminolysis [19–21]. Therefore, cell culture media usually contain substantial amounts of both substrates. Associated with an overflow metabolism based on glucose and a high glutaminolysis activity is the accumulation of lactate and ammonia in the cultivation broth with a negative impact on cell growth and product formation [22, 23].
Results and discussion
To allow metabolic flux analysis, some cell characteristics need to be known, i.e. specific dry cell weight and biomass composition. As it was not known whether similar results to other eukaryotic cells could be expected for an avian cell, most characteristics were determined experimentally rather than presumed.
The specific cell dry weight of CR.pIX cells was measured with 314 pg/cell. Amounts of DNA (2.3 ± 0.5% of biomass) and RNA (3.1 ± 0.2%) per CR.pIX cell were comparable to published values obtained for other eukaryotic cells. The protein content (55.2 ± 8.4%) was lower in CR.pIX cells than the typical range of 70–75% determined for other cells [7, 16, 29]. However, there are reports supporting a low protein content such as Zupke et al. who found 60% protein in mouse hybridoma cells and Carnicer et al. who measured a protein content of only 37% for yeast cells [30, 31]. The relative amounts of amino acids of whole cell protein determined from CR.pIX cells are in general similar to published data on yeast , mammalian cells [9, 29, 32] and insect cell lines  (see Additional file 1: Table S1). The remaining fraction of biomass was assigned to lipids and carbohydrates in a ratio of 1:2. The estimated lipid content (13.1%) agrees with previous reports that vary between 9 and 20%, whereas the assumed carbohydrate content (26.3%) is at the upper limit of the wide range of 3.5–25.0% reported in literature [7, 16, 29, 30, 34]. However, a sensitivity analysis showed that absolute fractions of lipids and carbohydrates have only a minor impact on rate values of the calculated flux distributions (data not shown).
Growth phases of CR.pIX cells in STR
The region chosen for MFA and FVA was 24–97 h, where a growth rate of 0.0132 h-1 was observed. In the following, only this time frame is considered when describing and discussing extracellular rates and intracellular flux distributions as it is very likely that the quasi steady state assumption holds. Furthermore, it is also the typical time frame in which cells are grown before infection in virus production processes (MVA or influenza A and B virus) [6, 35]. Studies in this time frame are therefore most relevant in terms of vaccine manufacturing using CR.pIX cells.
Exchange rates from measured extracellular metabolites
Determined extracellular fluxes for cultivation of CR.pIX cells in comparison with literature data
Initial glc conc. [mM]
Initial gln conc. [mM]
-211.2 ± 10.8g
297.6 ± 8.8
-30.1 ± 0.8
-0.1 ± 0.1
12.6 ± 0.9
-1.4 ± 1.1
31.6 ± 4.9
-9.9 ± 2.8
-17.1 ± 8.5
-17.8 ± 3.2
-4.8 ± 4.8
-2.3 ± 0.4
-6.2 ± 4.4
-7.9 ± 3.2
-5.5 ± 0.3
-4.1 ± 2.0
-2.5 ± 1.0
-0.1 ± 0.1
-8.3 ± 8.3
-3.4 ± 1.9
-1.5 ± 0.2
-1.8 ± 1.8
-4.3 ± 2.8
First, glutamine uptake is zero with CR.pIX cells (Figure 2C) as the observed glutamine concentration decrease is solely due to glutamine hydrolysis. Associated with this, only low amounts of ammonia accumulated in the medium (1.4 mM after 97 h: maximal 2.8 mM at the end of the batch cultivation at 210 h). For production of cell culture-derived products this is a beneficial characteristic as higher ammonia concentrations are known to be toxic with negative impact on cell growth and product formation . The observed ammonia concentration is at the lower limit compared to other cells, where concentrations between 3 and 8 mM were reported to lead to impaired growth [37–39].
The second difference is observable for glycine which is taken up by CR.pIX cells, but released by other cell lines. This observation can be explained by the fact that glycine is reported to be an essential amino acid for birds and therefore needs to be provided for CR.pIX cells by the medium . The third finding is that aspartate and asparagine are consumed by CR.pIX cells in higher amounts or at least at higher rates compared to other cell lines. Both amino acids can be converted to oxaloacetate and thus channeled into the TCA cycle. One could hypothesize that CR.pIX cells compensate in that way for the negligible glutamine uptake to fill the TCA cycle intermediate pools. We also observed that the measured uptake rates for the essential amino acids histidine, lysine and tryptophan were slightly below their minimal requirements for biomass formation with the determined growth rate. For the following analyses we therefore adjusted their upper boundaries slightly to fit the minimal stoichiometric demands for the measured specific growth rate.
Intracellular flux distribution
Using a metabolic network model of the central metabolism of avian cells that we developed according to literature data and database entries, the metabolism of CR.pIX cells during cultivation in 1 L STR was studied using FVA (details see section “Flux variability and metabolic flux analysis”). The time span of 24 h and 97 h was chosen for analysis, as within this time interval during exponential growth, the quasi steady state assumption holds. FVA uses as few constraints for the metabolic network as possible in order to get an unbiased look at the metabolism. This generated a distribution of flux ranges that is named scenario 1 in this manuscript. However, flux ranges are sometimes difficult to interpret. Therefore, after having performed FVA, we set few reasonable constraints, e.g. uptake and excretion rates were set as fixed and the pyruvate carboxylase activity was assumed to be inactive. With these constraints a unique flux distribution could be calculated which is discussed as scenario 2 in the following.
All calculated flux ranges (scenario 1; determined by FVA) and fluxes (for scenario 2, determined by MFA) are given in the Additional file (Additional file 1: Table S2). Together with the network, main fluxes are also depicted in Figure 1. Although only flux ranges can be calculated with the constraints of scenario 1 (Figure 1A), some conclusions can already be drawn from them as we will discuss below. However, most of the following discussion will refer to the unique flux distribution calculated for scenario 2 (Figure 1B). With the applied constraints, the linear equation system had two degrees of redundancy. Therefore, the sum of the variance weighted squared residuals (h) need to be below 5.99 to exclude significant measurement errors (with an applied significance level of α = 0.05). Since we obtained an h-value of 0.63 we concluded that the model fitted the data sufficiently well.
Glycolysis and cytosolic pyruvate metabolism
Glucose is taken up with a high rate (211 μmol/gDW/h) and then processed via glycolytic reactions to pyruvate. Only 1% of the glucose is fed into the pentose phosphate pathway for nucleotide synthesis. Such low pentose phosphate influxes were also observed in other studies with transformed cells, e.g. 1.5% and 2.9% in human AGE1.HN cells [9, 36]. The calculated ranges for the rates of glycolysis are relatively narrow as they mainly depend on the glucose uptake rate that was given by the measurement of extracellular glucose. Cytosolic pyruvate is mainly derived from glycolysis (92%); the rest of the cytoplasmic pyruvate is taken up by the cells and in small parts obtained from cysteine degradation and the malic enzyme reaction (together 1%). Similar to other transformed cells, cytosolic pyruvate is mainly converted to lactate (73%) or alanine (10%). Only a small fraction (16%) is transported to the mitochondria where it is converted to acetyl-CoA and further oxidized. However, if a complete pentose phosphate pathway (PPP) is included in the model, the exact ratio between glycolysis and PPP is unresolvable and the corresponding flux ranges would be notably wider. Since CO2 is produced via the PPP, a high PPP flux would reduce the amount of carbon that re-enters the lower part of glycolysis, and therefore even less pyruvate would enter the TCA cycle.
In general, the flux distribution at the cytosolic pyruvate node shows that the majority of the glucose-associated metabolism of CR.pIX cells is used to generate ATP and reducing equivalents via glycolysis. The known Warburg effect of enhanced aerobic glycolysis, leading mostly to lactate, can thus also be observed in transformed cells from avian origin. Furthermore, there seems to be a weak coupling to the mitochondrial TCA cycle as only a minority of the cytosolic pyruvate is channeled into the mitochondria. Other studies have shown that several transformed cell lines show low activities of the three pyruvate dependent enzymes pyruvate dehydrogenase (PDH), pyruvate decarboxylase (PC) and phosphoenolpyruvate carboxykinase (PEPCK) that can serve as connectors between glycolysis and TCA cycle [41, 42]. For our cells, we calculated that the PDH (r33) must carry some flux (since the lower boundary of this reaction is strictly positive), but the non-resolvable PC reaction (r46) is not obligatory as its flux range also includes zero. Measurements on activities of anaplerotic enzymes will be discussed in a later subsection.
Mitochondrial pyruvate metabolism and TCA cycle
The mitochondrial pyruvate pool is fed from the cytosolic pyruvate pool (64%) and the anaplerotic malate conversion (34%). 100% of the efflux is catalyzed by the pyruvate dehydrogenase thereby generating 88% of the mitochondrial acetyl-CoA. Other influxes towards this TCA precursor come from amino acids, namely isoleucine (4%) and leucine (7%), tryptophan, and lysine (together 1%). The subsequent reactions of the TCA cycle are fed from other catabolic reactions: α-ketoglutarate from cytosolic influx and succinyl-CoA from isoleucine, valine, serine and methionine degradation. This shows the dependency of the TCA cycle activity on relatively high amino acid uptake and catabolic rates in CR.pIX cells.
Oxidative phosphorylation & ATP balance
We calculated high fluxes for the oxidation of NADH (177 μmol/gDW/h) and FADH2 (46 μmol/gDW/h). Since a transport reaction for NADH between the cytosol and mitochondria was included in the model, the pools of reducing equivalents for both compartments are balanced simultaneously. A minor flux for the corresponding transport reaction shifting reducing equivalents from the mitochondria to the cytosol was calculated with 30 μmol/gDW/h. By setting this transport reaction to zero, cytosolic and mitochondrial reducing equivalents have to be balanced separately. Besides an increased oxygen uptake and ATP synthesis rate, this leads to no significant differences in the flux distributions (data not shown).
ATP generation via the TCA cycle and the electron transport chain makes up to 74% of the total ATP generated. A smaller fraction (26%) is derived from glycolysis. This is similar to results from studies on mouse hybridoma and CHO cells [43, 44] where values of 60–80% of ATP production via the TCA cycle have been reported. Although considerable ATP amounts are needed for lipid synthesis, amino acid transporters, upper glycolysis and biomass, the high (unspecific) ATP consumption through reaction r77 (966 μmol/gDW/h) indicates that there seems to be an overload with ATP which is either used for certain maintenance processes or passed through futile cycles.
Amino acid and nitrogen metabolism
Degradation rates of some amino acids, e.g. tryptophan and histidine, are low, indicating that these amino acids are taken up to match the anabolic demand for synthesizing biomass rather than to further metabolize them via catabolic routes. The calculated rate ranges for most fluxes of the amino acid catabolism are constrained to relatively small intervals (e.g. for degradation of valine and histidine). There is no uric acid synthesis (r68) as no uric acid could be measured experimentally. Although this reaction is known to take place in avian cells, this is not a fully unexpected finding. Uric acid was detected in livers and kidneys of chicken and turkey, but not in all bird tissues . CR.pIX cells, as other transformed cells, use ammonia as a sink and release the excess cytosolic ammonia (that is produced by deamination reactions of amino acids) into the medium where it accumulates. However, due to low glutamine exploitation, the ammonium concentration stays at a low level so that toxic concentrations are not reached.
Interestingly, the fluxes through the glutamine synthetase (r48) (i.e. in the direction of glutamate synthesis) are negative as glutamine uptake was not sufficient for biomass synthesis and glutamine therefore has to be synthesized from the glutamate pool. Glutamate itself is likely to be mainly derived from degradation reactions of isoleucine, valine and aspartate. Especially aspartate is taken up in large parts (18 μmol/gDW/h) and can be converted to glutamate. Other amino acids that could compensate for glutamine conversion to glutamate are arginine, histidine and proline. However, only arginine is taken up in considerable amounts by the cells so that one could speculate that the path from glutamate to α-ketoglutarate is generally not used extensively by the CR.pIX cells and thus, the classical dependency on glutamine is not given for CR.pIX cells. To elucidate unambiguously the source of carbon for glutamate synthesis, 13C labelling experiments would be required, which are beyond the scope of the current contribution.
In vitroenzyme activities
To validate the two hypotheses from the previous section that i) glycolysis is only weakly connected to the TCA, and that ii) glutaminolysis plays a minor role in energy supply and precursor generation, we measured the maximum in vitro enzyme activities of the related enzymes.
For the first hypothesis, we measured maximum activities of PC and PDH. For both reactions, low enzyme activities of 0.2 nmol/106cells/min (i.e. 39 μmol/gDW/h) were measured as for MDCK cells . These values supported the applied constraint of a missing carboxylating anaplerotic reaction from pyruvate towards oxaloacetate. However, the measured low PDH activity seems to contradict the calculated flux of 102 μmol/gDW/h. One possible explanation is given by the significantly smaller calculated lower boundary of this flux when applying fewer and less stringent assumptions with scenario 1. However, at least a basal activity is required as the lower boundary is strictly positive (53 μmol/gDW/h).
The second hypothesis, stating that glutaminolysis is inversed, was validated via measurements of glutaminase (converting glutamine to glutamate) and glutamine synthetase (converting glutamate to glutamine). Glutaminase showed a maximal activity of 0.1 nmol/106cells/min (i.e. 19 μmol/gDW/h), which is an even lower value than the measured 1.5–7.4 nmol/106cells/min for MDCK cells [45, 46]. The flux from glutamate to glutamine was predicted by flux analysis to be active towards glutamine in order to account for the unmet demand of glutamine for biomass formation. The measured enzyme activity of the catalyzing enzyme glutamine synthetase showed a maximal activity of 5.9 nmol/106cells/min (i.e. 1135 μmol/gDW/h). This value is on the one hand significantly higher than the value for the reverse reaction and on the other hand at the upper limit of reported values for MDCK cells (0.6–6.5 nmol/106cells/min, [45, 46]). Thus, although values from in vitro enzyme activity measurement and flux analysis are not directly convertible, the measured enzyme activities support the calculated flux distribution towards glutamine and the low level of glutamine usage.
Growth in glutamine-free medium
Avian cells are important substrates for producing a number of licensed vaccines and various viral vector-based processes under development. The only avian substrates that are currently approved for production are embryonated chicken eggs or primary chicken embryo fibroblasts. Because processes using primary cultures are difficult to optimize and standardize, access to well characterized continuous avian cell lines, such as the CR.pIX cell line derived from a muscovy duck embryo, would be a significant step forward.
Here, we determined metabolic properties of the CR.pIX cell line to complement previous developments on scalable vaccine production processes in chemically defined media. The current study comprised three parts: i) analysis of biomass composition and measurements of cell and extracellular metabolite concentrations, ii) construction of a stoichiometric network for the central metabolism of avian cells and subsequent analysis with flux variability and metabolic flux analysis, and iii) experiments and analytics to validate hypotheses derived from the stoichiometric analysis.
The first part expands the available data base on the metabolism of immortalized cell lines. Results on biomass composition and growth properties of the avian CR.pIX cell line agree with data determined for other transformed mammalian and insect cell lines used in biopharmaceutical production.
Flux variability and metabolic flux analyses in the second part reveal an overflow metabolism similar to mammalian cells with high uptake of glucose followed by conversion to and release of lactate and alanine. The TCA cycle was only weakly connected to glycolysis and appeared to depend on the influx of amino acids, mainly isoleucine, valine and aspartate. Finally, our measurements together with flux (variability) analysis indicated that glutamine supplementation is not an essential requirement for CR.pIX cell proliferation.
The derived hypothesis on low glutaminolysis activity could be confirmed by enzyme activity measurements. However, the hypothesis of a weak connection between the cytosolic and mitochondrial pyruvate pools could not be confirmed entirely. To elucidate whether cytosolic pyruvate is transported in significant amounts from the cytosol to the mitochondria, e.g. 13C labelling experiments could be helpful, which was beyond the scope of the current contribution. Low dependence on glutamine was confirmed by passaging cells in glutamine-free medium and agrees well with the observation that ammonia is accumulating to very low concentrations in the medium. This is a very beneficial property of CR.pIX cells especially for fed-batch cultivation strategies where ammonia often accumulates to toxic levels.
Overall, the present study provides a basis for a more focused metabolic analysis of avian cell lines and thus helps to further optimize medium and feed requirements. Acute infection and associated take-over of the metabolism by a virus may shift the requirements for certain nutrients. The derived model may therefore also provide a context for future studies on the metabolism of CR.pIX cells during the virus replication phase.
CR.pIX suspension cells were routinely cultivated in 250 mL shaker flasks (working volume 150 mL, with membrane cap for aeration) in a chemically-defined medium (CD-U2, PAA) supplemented with glutamine (2 mM final concentration, Sigma), alanine (2 mM final concentration, Sigma) and a recombinant insulin-like growth factor (Long-R3 IGF, 10 ng/mL final concentration, Sigma). Flasks were incubated at 185 rpm with 5 cm amplitude, 37°C and 5% CO2. Pre-cultures for bioreactor experiments, cultures for determination of biomass composition as well as the comparative cultivations in glutamine-containing and glutamine-free medium were performed in shaker flasks. For inoculation of a 1 L stirred tank bioreactor (STR, DasGip AG, vessel: Spinner Type BS, cellferm-pro® system), 8 × 105 cells/mL in 1 L CD-U2 medium were transferred into the bioreactor. The bioreactor was equipped with a pitched-blade stirrer that was operated at 120 rpm. The temperature was controlled at 37°C; pO2 was set to 50% by pulsed aeration with air enriched with variable contents of O2 and CO2. During the first 60 h, pH was slowly reduced from 7.4 (initial value of the medium) to 7.0 and controlled at this value subsequently by addition of CO2 or 1 M Na2CO3.
Measurement of biomass composition
All assays and analytics were performed with cells in the exponential growth phase. Therefore, samples were taken from shaker flask cultures at 72 h post inoculation. For each analysis, 3 biological replicates were measured with at least 3 technical replicates each.
For determination of dry cell weight, cell suspension containing about 1 × 108 cells/mL was centrifuged (500 × g, 10 min) and re-suspended in 10 mL of fresh medium. Samples were dried at 60°C until weight remained constant. The specific dry cell weight was calculated against weight of empty tubes and against a medium control (10 mL of fresh medium dried in tubes).
The BCA assay (Pierce/ThermoScientific) was used to determine the cellular protein content. Samples containing 4 × 106 cells were centrifuged (500 × g, 10 min) and re-suspended in 1 mL lysis buffer (from CytoTox 96® Non-Radioactive Cytotoxicity Assay, Promega). After 15 min incubation at room temperature, samples were sonicated for 3 min in a water bath. Standard curves were prepared using BSA and analyzed together with samples by measuring the absorption at 562 nm (inifinite200 plate reader, TecanGroup).
For DNA purification, a culture volume containing 5 × 106 cells was taken and purified with the QIAamp DNA Blood Mini kit (Qiagen) according to the procedure “protocol for cultured cells”. After purification, 2 μL of assay buffer (in which the samples were dissolved later) were pipetted onto 16 spots on a NanoQuantPlate (TecanGroup) as blank. After blanking, 2 μL of samples were measured. DNA concentration was then determined by multiplication of the measured absorption coefficient with 50 as a dsDNA specific factor. The same method as for DNA content analysis was employed for RNA content analysis. Here, purification was done with the NucleoSpin RNAII kit (Macherey-Nagel) according to the procedure “total RNA purification from cultured cells and tissue”. RNA concentration was determined from absorption measurements at 260 nm, but by multiplication with a factor of 40. Measurement procedure was analogous to the DNA analysis.
Amino acid contents were determined as follows: about 1 × 108 cells were washed with 0.9% NaCl solution and re-suspended in 2 M HCl. Cells were incubated in acid at 100°C for 24 h. After this hydrolyzation, the solution was neutralized by adding 2 M NaOH and subsequently analysed by HPLC. Aspartate and glutamate percentages were corrected for a hydrolysis rate of 47.5% as described in literature .
A HPLC method was used to measure the amino acid composition of cell protein . Due to harsh hydrolytic conditions, some amino acid concentrations could not be measured, but were estimated from literature values (see Additional file 1: Table S1).
Cell and metabolite concentration measurements
Samples were taken once or twice a day. Cell concentration, viability and average cell diameter were determined using the Vi-CELL™ XR device (Beckman Coulter). Samples for measuring extracellular metabolite concentrations were centrifuged after sampling and the supernatant was stored at -80°C until analysis. Measurement of glucose, lactate and ammonia concentrations was done with a Bioprofile 100 Plus (Nova Biomedical). A HPLC method was used to measure the concentration of pyruvate. Extracellular amino acid concentrations for metabolic flux analysis were determined using a derivatization method and subsequent measurement in a HPLC system . Uric acid concentration was measured with an enzymatic reaction assay kit (Uric Acid Assay Kit, BioVision) including standard curves and controls as described in the manual. Uric acid was only measured at culture start and end. Standard deviations for all concentration measurements (determined by method validations) can be found in Additional file 1: Table S3.
In vitrodetermination of maximum enzyme activities
Maximum activities of pyruvate carboxylase (PC), pyruvate dehydrogenase (PDH), glutamine synthetase (GS) and glutaminase (GLNase) in exponentially growing cells were measured. Therefore, cells were harvested from shaker flasks 72 h after inoculation. After a washing step in 0.9% NaCl, dry cell pellets were frozen at -80°C and stored until analysis. Enzyme extraction and assays were performed as described for MDCK cells . Three assay parameters had to be optimized before measuring CR.pIX cell extracts: pH, amount of cells that are used and – for GLNase and GS activity assay – incubation time. After these pre-tests, CR.pIX cell extracts were measured with the following modifications: a two-fold higher cell concentration as for MDCK cells was necessary for the GLNase and GS assay and the incubation time of the GLNase assay was set to 30 min instead of 60 min. The assay plates were analyzed in a 96 well plate reader (Tecan) and enzyme activities were calculated from absorption values.
Reactions included in the metabolic network model of CR.pIX cells
Glc → Glccyt
Pyr + 0.33 ATPcyt → Pyrcyt
O2 → O2,cyt
Gln + 0.33 ATPcyt → Glncyt
Glu + ATPcyt → Glucyt
Ala + 0.33 ATPcyt → Alacyt
Asp + ATPcyt → Aspcyt
Arg + 0.33 ATPcyt → Argcyt
Asn + 0.33 ATPcyt → Asncyt
Cys + 0.33 ATPcyt → Cyscyt
Gly + 0.33 ATPcyt → Glycyt
His + 0.33 ATPcyt → Hiscyt
Ile + 0.33 ATPcyt → Ilecyt
Leu + 0.33 ATPcyt → Leucyt
Lys + 0.33 ATPcyt → Lyscyt
Val + 0.33 ATPcyt → Valcyt
Met + 0.33 ATPcyt → Metcyt
Phe + 0.33 ATPcyt → Phecyt
Pro + 0.33 ATPcyt → Procyt
Ser + 0.33 ATPcyt → Sercyt
Thr + 0.33 ATPcyt → Thrcyt
Trp + 0.33 ATPcyt → Trpcyt
Tyr + 0.33 ATPcyt → Tyrcyt
Glccyt + ATPcyt → G6Pcyt
G6Pcyt ↔ F6Pcyt
F6Pcyt + ATPcyt ↔ FBPcyt
FBPcyt ↔ GAPcyt + DHAPcyt
DHAPcyt ↔ GAPcyt
GAPcyt ↔ NADHcyt + ATPcyt + PGcyt
PEPcyt → Pyrcyt + ATPcyt
Pyrcyt + NADHcyt ↔ Laccyt
Pyrmit + CoAmit → AcCoAmit + CO2,mit + NADHmit
Pentose phosphate pathway
G6Pcyt → R5Pcyt + CO2,cyt + 2 NADPHcyt
Malcyt ↔ OAAcyt + NADHcyt
AcCoAcyt + OAAcyt → Citcyt + CoAcyt
Fumcyt ↔ Malcyt
AcCoAmit + OAAmit → Citmit + CoAmit
Malmit ↔ OAAmit + NADHmit
Fummit ↔ Malmit
aKGmit + CoAmit → SuccCoAmit + CO2,mit + NADHmit
SuccCoAmit ↔ Fummit + CoAmit + ATPmit + FADH2,mit
Citmit ↔ aKGmit + CO2,mit + NADHmit
Malcyt → Pyrcyt + CO2,cyt + NADPHcyt
Malmit ↔ Pyrmit + CO2,mit + NADPHmit
Pyrmit + CO2,mit + ATPmit → OAAmit
Amino acid catabolism
Glumit ↔ aKGmit + Ammmit + NADPHmit
Glncyt ↔ Glucyt + Ammcyt + ATPcyt
Alacyt + aKGcyt ↔ Pyrcyt + Glucyt
Asncyt ↔ Aspcyt + Ammcyt
Hiscyt → Glucyt + 2 Ammcyt + CO2,cyt
Ilecyt + aKGcyt + 2 CoAmit + ATPmit → SuccCoAmit + AcCoAmit + Glucyt + NADHmit + FADH2,mit
Leucyt + aKGcyt + 3 CoAmit + 2 ATPmit → 3 AcCoAmit + Glucyt + NADHmit + FADH2,mit
Lyscyt + 2 aKGcyt + 2 CoAmit + NADPHcyt → 2 AcCoAmit + 2 Glucyt + 2 CO2,mit + 2 NADHcyt + 2 NADHmit + FADH2,mit
Metcyt + Sercyt + CoAmit + 3 ATPcyt + ATPmit → Cyscyt + SuccCoAmit + Ammcyt + NADHmit + CO2,cyt
Phecyt + O2,cyt + NADPHcyt → Tyrcyt
Procyt ↔ Glucyt + 2 NADHcyt
Thrcyt → Pyrmit + Ammmit + CO2,mit + NADHcyt + NADHmit + FADH2,mit
Trpcyt + 2 CoAmit + 3 O2,cyt + NADPHcyt → Alacyt + 2 AcCoAmit + Ammcyt + 2 CO2,mit + 2 CO2,cyt + NADHcyt + 2 NADHmit + FADH2,mit
Valcyt + aKGcyt + CoAmit + ATPmit → SuccCoAmit + Glucyt + CO2,cyt + 3 NADHmit + FADH2,mit
Tyrcyt + aKGcyt + 2 CoAcyt + 2 O2,cyt → Fumcyt + 2 AcCoAcyt + Glucyt + CO2,cyt
Sercyt ↔ Pyrcyt + Ammcyt
Cyscyt + aKGcyt + O2,cyt → Pyrcyt + Glucyt
Aspcyt + aKGmit ↔ Glucyt + OAAmit
Argcyt + aKGcyt → 2 Glucyt + Ureacyt + NADHcyt
MTHF & uric acid synthesis
Sercyt + THFcyt → Glycyt + MTHFcyt
Glycyt + THFcyt → NADHcyt + CO2,cyt + Ammcyt + MTHFcyt
Aspcyt + 2 Glncyt + Glycyt + 2 MTHFcyt + 7 ATPcyt + CO2,cyt → UricAcidcyt + Fumcyt + 2 Glucyt + 2 THFcyt
18 AcCoAcyt + 18 ATPcyt + 11 O2,cyt + 27 NADPHcyt → CH + 18 CoAcyt + 9 CO2,cyt
GAPcyt + Sercyt + 27.6 ATPcyt + 17.6 AcCoAcyt + 4 MTHFcyt + 2 NADHcyt + 31.2 NADPHcyt → PC + 17.6 CoAcyt + 4 THFcyt
GAPcyt + Sercyt + 18.6 ATPcyt + 17.6 AcCoAcyt + MTHFcyt + 2 NADHcyt + 31.2 NADPHcyt → PE + 17.6 CoAcyt + THFcyt
GAPcyt + Sercyt + 18.6 ATPcyt + 17.6 AcCoAcyt + 2 MTHFcyt + 2 NADHcyt + 31.2 NADPHcyt → PS + 17.6 CoAcyt + 2 THFcyt
2 GAPcyt + 17.6 ATPcyt + 17.6 AcCoAcyt + 4 NADHcyt + 31.2 NADPHcyt → PGL + 17.6 CoAcyt
GAPcyt + G6Pcyt + 17.6 ATPcyt + 17.6 AcCoAcyt + 2 NADHcyt + 31.2 NADPHcyt → PI + 17.6 CoAcyt
2 Sercyt + 27.8 ATPcyt + 16.8 AcCoAcyt + 3 MTHFcyt + 2 NADHcyt + 29.6 NADPHcyt → SM + 16.8 CoAcyt + 3 THFcyt
3 GAPcyt + 35.2 ATPcyt + 35.2 AcCoAcyt + 6 NADHcyt + 62.4 NADPHcyt → DPG + 35.2 CoAcyt
ATPcyt → maintenance
Laccyt → Lac
Alacyt → Ala
Uric acidcyt → Uric acid
Ureacyt → Urea
Ammcyt → Amm
r83, CO2 out
CO2,cyt → CO2
Pyrcyt → Pyr
Transport reactions, oxidative phosphorylation
NADHcyt ↔ NADHmit
ATPcyt ↔ ATPmitt
r87, CO2 trans
CO2,cyt ↔ CO2,mit
Malcyt + Citmit ↔ Malmit + Citcyt
Glucyt ↔ Glumit
Pyrcyt ↔ Pyrmit
aKGcyt ↔ aKGmit
Ammcyt ↔ Ammmit
O2,cyt + 2 FADH2,mit → 3 ATPmit
NADHmit ↔ NADPHmit
O2,cyt + 2 NADHmit → 5 ATPmit
NADHcyt ↔ NADPHcyt
Synthesis of macromolecules and biomass
0.552 proteins + 0.263 carbohydrates + 0.131 lipids + 0.023 DNA + 0.031 RNA → biomass
Proteins [1 g] = 955.79 Aspcyt + 1344.29 Alacyt + 543.16 Glncyt + + 817.99 Glucyt + 1024.72 Argcyt + 362.10 Asncyt + 19.61 Cyscyt + 913.16 Glycyt + 271.58 Hiscyt + 240.54 Ilecyt + 588.41 Leucyt + 724.21 Lyscyt + 191.97 Valcyt + 114.57 Metcyt + 201.91 Phecyt + 205.97 Procyt + 247.44 Sercyt + 202.07 Thrcyt + 21.17 Trpcyt + 61.87 Tyrcyt + 24046.3 ATPcyt
Lipids [1 g] = 181 CH + 661.4 PC + 250.3 PE + 90.9 PI + 24.9 PS + 12.6 PGL + 81.4 SM + 26.8 DPG
DNA [1 g] = 3009 R5Pcyt + 3912 Aspcyt + 5717 Glncyt + -2106 NADHcyt + 1505 Glycyt + 22569 ATPcyt + 903 NADHmit + -2407 Malcyt + 903 NADPHcyt + -5717 Glucyt + 5417 MTHFcyt + -5417 THFcyt
RNA [1 g] = 3020 R5Pcyt + 3606 Aspcyt + 6316 Glncyt + 293 O2,cyt + -2435 NADHcyt + 1477 Glycyt + 22614 ATPcyt + 586 NADHmit + -2069 Malcyt + -2954 NADPHcyt + -6316 Glucyt + 4431 MTHFcyt + -4431 THFcyt
Carbohydrates [1 g] = 6172.8 G6Pcyt + 21605 ATPcyt
Generally, the network includes
Two compartments: cytosol and mitochondria
Main pathways of energy metabolism: glycolysis, TCA cycle and oxidative phosphorylation
Catabolism of amino acids
Transport of metabolites and corresponding energy demands (as described before )
Synthesis of macromolecules (biomass)
Synthesis of urea and uric acid
A simplified/lumped pentose-phosphate pathway (as in )
Calculation of extracellular rates
where we applied the experimentally determined cell dry weight and the calculated μ from Equation (1).
For glutamine, a first order degradation process to ammonia and 5-oxopyrrolidine-2-carboxylic acid was assumed. The rate constant k Gln = 0.0032 h-1 of this process was determined experimentally (n = 3) using medium without cells under standard cultivation conditions in a bioreactor. Based on this degradation rate, glutamine and ammonia rates (Equation (2)) were corrected as described elsewhere . To compute standard deviations of the fluxes, we performed a Monte-Carlo simulation and added a normally distributed error with mean zero and corresponding standard deviation (derived from assay validations) at each time point and for each measurement (metabolites and cell concentration). With this, we generated 100,000 random sample data sets and recalculated for each data set the rates of uptake, release and growth. From these simulated data sets, the empirical variances and standard deviations for each of these rates were computed.
Flux variability and metabolic flux analysis
As a result from the 2∗u optimizations (u = number of unknown fluxes), we obtained the physiologically feasible flux range for the unknown reactions. Moreover, if the computed minimal and maximal rate of a reaction coincide (ru,imin = ru,imax), the reaction rate follows to be uniquely determined. Note that FVA as described above does not make any assumption about biological objectives in contrast to flux balance analysis . The objective function in equation (4) only serves as a tool to identify the feasible flux ranges. This approach has also been used to estimate flux distributions in CHO cells and was introduced under the term flux-spectrum .
We define two scenarios where we consider different sets of constraints. In scenario 1, we applied the measured rate values plus/minus their corresponding standard deviations as upper/lower boundaries for the corresponding rates. In case that the lower boundary of an uptake rate becomes negative, we set it to zero. Rates that were measured not to be active (e.g. uric acid excretion and alanine uptake) were also fixed to zero. No further constraints were used so that the results obtained with this scenario are robust against measurement errors and assumptions on reaction activities. However, since the system is underdetermined with these constraints, only flux ranges can be calculated via FVA.
In scenario 2, we set the measured uptake/excretion rates as fixed constraints. To resolve the anaplerotic fluxes and since the corresponding flux range contains zero (see Results), we assumed the pyruvate-carboxylase (r46) to be inactive. This assumption is biologically reasonable as it is often described and assumed to be negligible in transformed cell cultures (e.g. in [25, 61]) and supported by our enzyme measurements. With these constraints the network is over-determined with two degrees of redundancy. By applying MFA with a variance-weighted least squares estimation , we can calculate a unique flux distribution. To avoid a significant deviation of the estimated against the measured growth rate, we set the measurement variance of the growth rate to 1 × 10-9 in scenario 2.
All computations presented in this study were performed with our software CellNetAnalyzer, a MATLAB toolbox with graphical user interface facilitating metabolic network analysis . It can be downloaded from http://www.mpi-magdeburg.mpg.de/projects/cna/cna.html, and the network project files will be made available on this site within CellNetAnalyzer’s model repository.
The authors thank C. Best, N. Wynserski, I. Behrendt and S. König for their excellent technical assistance and A. Lagoda for her advice in RNA, DNA and protein quantification. Also, advice and assistance by Anne-Kareen Blechert in enzyme activity assays was highly appreciated.
- Butler M: Optimisation of the cellular metabolism of glycosylation for recombinant proteins produced by Mammalian cell systems. Cytotechnology. 2006, 50 (1–3): 57-76.View ArticleGoogle Scholar
- Lim Y, Wong NS, Lee YY, Ku SC, Wong DC, Yap MG: Engineering mammalian cells in bioprocessing - current achievements and future perspectives. Biotechnol Appl Biochem. 2010, 55 (4): 175-189.View ArticleGoogle Scholar
- Blanchard TJ, Alcami A, Andrea P, Smith GL: Modified vaccinia virus Ankara undergoes limited replication in human cells and lacks several immunomodulatory proteins: implications for use as a human vaccine. J Gen Virol. 1998, 79 (Pt 5): 1159-1167.View ArticleGoogle Scholar
- Draper SJ, Heeney JL: Viruses as vaccine vectors for infectious diseases and cancer. Nat Rev Microbiol. 2010, 8 (1): 62-73.View ArticleGoogle Scholar
- Jordan I, Vos A, Beilfuss S, Neubert A, Breul S, Sandig V: An avian cell line designed for production of highly attenuated viruses. Vaccine. 2009, 27 (5): 748-756.View ArticleGoogle Scholar
- Lohr V, Rath A, Genzel Y, Jordan I, Sandig V, Reichl U: New avian suspension cell lines provide production of influenza virus and MVA in serum-free media: studies on growth, metabolism and virus propagation. Vaccine. 2009, 27 (36): 4975-4982.View ArticleGoogle Scholar
- Nadeau I, Sabatie J, Koehl M, Perrier M, Kamen A: Human 293 cell metabolism in low glutamine-supplied culture: interpretation of metabolic changes through metabolic flux analysis. Metab Eng. 2000, 2 (4): 277-292.View ArticleGoogle Scholar
- Nadeau I, Jacob D, Perrier M, Kamen A: 293SF metabolic flux analysis during cell growth and infection with an adenoviral vector. Biotechnol Prog. 2000, 16 (5): 872-884.View ArticleGoogle Scholar
- Niklas J, Schrader E, Sandig V, Noll T, Heinzle E: Quantitative characterization of metabolism and metabolic shifts during growth of the new human cell line AGE1.HN using time resolved metabolic flux analysis. Bioprocess Biosyst Eng. 2011, 34 (5): 533-545.View ArticleGoogle Scholar
- Priesnitz C, Niklas J, Rose T, Sandig V, Heinzle E: Metabolic flux rearrangement in the amino acid metabolism reduces ammonia stress in the alpha1-antitrypsin producing human AGE1.HN cell line. Metab Eng. 2012, 14 (2): 128-137.View ArticleGoogle Scholar
- Altamirano C, Illanes A, Casablancas A, Gamez X, Cairo JJ, Godia C: Analysis of CHO cells metabolic redistribution in a glutamate-based defined medium in continuous culture. Biotechnol Prog. 2001, 17 (6): 1032-1041.View ArticleGoogle Scholar
- Teixeira A, Cunha AE, Clemente JJ, Moreira JL, Cruz HJ, Alves PM, Carrondo MJ, Oliveira R: Modelling and optimization of a recombinant BHK-21 cultivation process using hybrid grey-box systems. J Biotechnol. 2005, 118 (3): 290-303.View ArticleGoogle Scholar
- Sidorenko Y, Wahl A, Dauner M, Genzel Y, Reichl U: Comparison of metabolic flux distributions for MDCK cell growth in glutamine- and pyruvate-containing media. Biotechnol Prog. 2008, 24 (2): 311-320.View ArticleGoogle Scholar
- Wahl A, Sidorenko Y, Dauner M, Genzel Y, Reichl U: Metabolic flux model for an anchorage-dependent MDCK cell line: characteristic growth phases and minimum substrate consumption flux distribution. Biotechnol Bioeng. 2008, 101 (1): 135-152.View ArticleGoogle Scholar
- Bonarius HP, Houtman JH, Schmid G, de Gooijer CD, Tramper J: Metabolic-flux analysis of hybridoma cells under oxidative and reductive stress using mass balances. Cytotechnology. 2000, 32 (2): 97-107.View ArticleGoogle Scholar
- Bonarius HP, Hatzimanikatis V, Meesters KP, de Gooijer CD, Schmid G, Tramper J: Metabolic flux analysis of hybridoma cells in different culture media using mass balances. Biotechnol Bioeng. 1996, 50 (3): 299-318.View ArticleGoogle Scholar
- Bernal V, Carinhas N, Yokomizo AY, Carrondo MJ, Alves PM: Cell density effect in the baculovirus-insect cells system: a quantitative analysis of energetic metabolism. Biotechnol Bioeng. 2009, 104 (1): 162-180.View ArticleGoogle Scholar
- Carinhas N, Bernal V, Yokomizo AY, Carrondo MJ, Oliveira R, Alves PM: Baculovirus production for gene therapy: the role of cell density, multiplicity of infection and medium exchange. Appl Microbiol Biotechnol. 2009, 81 (6): 1041-1049.View ArticleGoogle Scholar
- Kovacevic Z, Brkljac O, Bajin K: Control and function of the transamination pathways of glutamine oxidation in tumour cells. Biochem J . 1991, 273 (Pt 2): 271-275.View ArticleGoogle Scholar
- Moreadith RW, Lehninger AL: The pathways of glutamate and glutamine oxidation by tumor cell mitochondria. Role of mitochondrial NAD(P)+-dependent malic enzyme. J Biol Chem. 1984, 259 (10): 6215-6221.Google Scholar
- Reitzer LJ, Wice BM, Kennell D: Evidence that glutamine, not sugar, is the major energy source for cultured HeLa cells. J Biol Chem. 1979, 254 (8): 2669-2676.Google Scholar
- Glacken MW, Fleischaker RJ, Sinskey AJ: Reduction of waste product excretion via nutrient control: possible strategies for maximizing product and cell yields on serum in cultures of mammalian cells. Biotechnol Bioeng. 1986, 28 (9): 1376-1389.View ArticleGoogle Scholar
- Miller WM, Blanch HW, Wilke CR: A kinetic analysis of hybridoma growth and metabolism in batch and continuous suspension culture: effect of nutrient concentration, dilution rate, and pH. Biotechnol Bioeng. 1988, 32 (8): 947-965.View ArticleGoogle Scholar
- Boghigian BA, Seth G, Kiss R, Pfeifer BA: Metabolic flux analysis and pharmaceutical production. Metab Eng. 2010, 12 (2): 81-95.View ArticleGoogle Scholar
- Quek LE, Dietmair S, Kromer JO, Nielsen LK: Metabolic flux analysis in mammalian cell culture. Metab Eng. 2010, 12 (2): 161-171.View ArticleGoogle Scholar
- Ahn WS, Antoniewicz MR: Metabolic flux analysis of CHO cells at growth and non-growth phases using isotopic tracers and mass spectrometry. Metab Eng. 2011, 13 (5): 598-609.View ArticleGoogle Scholar
- Mahadevan R, Schilling CH: The effects of alternate optimal solutions in constraint-based genome-scale metabolic models. Metab Eng. 2003, 5 (4): 264-276.View ArticleGoogle Scholar
- Klamt S, Saez-Rodriguez J, Gilles ED: Structural and functional analysis of cellular networks with Cell NetAnalyzer. BMC Syst Biol. 2007, 1: 2View ArticleGoogle Scholar
- Xie L, Wang DI: Applications of improved stoichiometric model in medium design and fed-batch cultivation of animal cells in bioreactor. Cytotechnology. 1994, 15 (1–3): 17-29.View ArticleGoogle Scholar
- Zupke C, Stephanopoulos G: Intracellular flux analysis in hybridomas using mass balances and in vitro (13)C nmr. Biotechnol Bioeng. 1995, 45 (4): 292-303.View ArticleGoogle Scholar
- Carnicer M, Baumann K, Toplitz I, Sanchez-Ferrando F, Mattanovich D, Ferrer P, Albiol J: Macromolecular and elemental composition analysis and extracellular metabolite balances of Pichia pastoris growing at different oxygen levels. Microb Cell Fact. 2009, 8: 65View ArticleGoogle Scholar
- Selvarasu S, Ho YS, Chong WP, Wong NS, Yusufi FN, Lee YY, Yap MG, Lee DY: Combined in silico modeling and metabolomics analysis to characterize fed-batch CHO cell culture. Biotechnol Bioeng. 2012, 109 (6): 1415-1429.View ArticleGoogle Scholar
- Ferrance JP, Goel A, Ataai MM: Utilization of glucose and amino acids in insect cell cultures: quantifying the metabolic flows within the primary pathways and medium development. Biotechnol Bioeng. 1993, 42 (6): 697-707.View ArticleGoogle Scholar
- Xie L, Wang DI: Fed-batch cultivation of animal cells using different medium design concepts and feeding strategies. Biotechnol Bioeng. 1994, 43 (11): 1175-1189.View ArticleGoogle Scholar
- Lohr V, Genzel Y, Jordan I, Katinger D, Mahr S, Sandig V, Reichl U: Live attenuated influenza viruses produced in a suspension process with avian AGE1.CR.pIX cells. BMC Biotechnol. 2012, 12: 79View ArticleGoogle Scholar
- Niklas J, Priesnitz C, Rose T, Sandig V, Heinzle E: Primary metabolism in the new human cell line AGE1.HN at various substrate levels: increased metabolic efficiency and alpha1-antitrypsin production at reduced pyruvate load. Appl Microbiol Biotechnol. 2012, 93 (4): 1637-1650.View ArticleGoogle Scholar
- Schneider M, Marison IW, von Stockar U: The importance of ammonia in mammalian cell culture. J Biotechnol. 1996, 46 (3): 161-185.View ArticleGoogle Scholar
- Yang M, Butler M: Effect of ammonia on the glycosylation of human recombinant erythropoietin in culture. Biotechnol Prog. 2000, 16 (5): 751-759.View ArticleGoogle Scholar
- Slivac I, Blajic V, Radosevic K, Kniewald Z, Gaurina Srcek V: Influence of different ammonium, lactate and glutamine concentrations on CCO cell growth. Cytotechnology. 2010, 62 (6): 585-594.View ArticleGoogle Scholar
- Stevens L: Avian Biochemistry and Molecular Biology. 1996, Cambridge: Cambridge University PressView ArticleGoogle Scholar
- Fitzpatrick L, Jenkins HA, Butler M: Glucose and glutamine metabolism of a murine B-lymphocyte hybridoma grown in batch culture. Appl Biochem Biotechnol. 1993, 43 (2): 93-116.View ArticleGoogle Scholar
- Neermann J, Wagner R: Comparative analysis of glucose and glutamine metabolism in transformed mammalian cell lines, insect and primary liver cells. J Cell Physiol. 1996, 166 (1): 152-169.View ArticleGoogle Scholar
- Xie L, Wang DI: Energy metabolism and ATP balance in animal cell cultivation using a stoichiometrically based reaction network. Biotechnol Bioeng. 1996, 52 (5): 591-601.View ArticleGoogle Scholar
- Zamorano F, Wouwer AV, Bastin G: A detailed metabolic flux analysis of an underdetermined network of CHO cells. J Biotechnol. 2010, 150 (4): 497-508.View ArticleGoogle Scholar
- Janke R, Genzel Y, Handel N, Wahl A, Reichl U: Metabolic adaptation of MDCK cells to different growth conditions: effects on catalytic activities of central metabolic enzymes. Biotechnology and bioengineering. 2011, 108 (11): 2691-2704.View ArticleGoogle Scholar
- Janke R, Genzel Y, Wahl A, Reichl U: Measurement of key enzyme activities in mammalian cell lines using rapid and sensitive microplate-based assays. Biotechnol Bioeng. 2010, 107 (3): 566-81.View ArticleGoogle Scholar
- Taschwer M, Hackl M, Hernandez Bort JA, Leitner C, Kumar N, Puc U, Grass J, Papst M, Kunert R, Altmann F, Borth N: Growth, productivity and protein glycosylation in a CHO EpoFc producer cell line adapted to glutamine-free growth. J Biotechnol. 2012, 157 (2): 295-303.View ArticleGoogle Scholar
- Capiaumont J, Legrand C, Carbonell D, Dousset B, Belleville F, Nabet P: Methods for reducing the ammonia in hybridoma cell cultures. J Biotechnol. 1995, 39 (1): 49-58.View ArticleGoogle Scholar
- Hassell T, Gleave S, Butler M: Growth inhibition in animal cell culture. The effect of lactate and ammonia. Appl Biochem Biotechnol. 1991, 30 (1): 29-41.View ArticleGoogle Scholar
- Fuchs BC, Bode BP: Stressing out over survival: glutamine as an apoptotic modulator. J Surg Res. 2006, 131 (1): 26-40.View ArticleGoogle Scholar
- Niederlechner S, Klawitter J, Baird C, Kallweit AR, Christians U, Wischmeyer PE: Fibronectin-integrin signaling is required for L-glutamine’s protection against gut injury. PLoS One. 2012, 7 (11): e50185View ArticleGoogle Scholar
- Chen P, Harcum SW: Effects of amino acid additions on ammonium stressed CHO cells. J Biotechnol. 2005, 117 (3): 277-286.View ArticleGoogle Scholar
- Genzel Y, Ritter JB, König S, Alt R, Reichl U: Substitution of glutamine by pyruvate to reduce ammonia formation and growth inhibition of mammalian cells. Biotechnol Prog. 2005, 21 (1): 58-69.View ArticleGoogle Scholar
- Buentemeyer H: Methods for off-line analysis of nutrients and products in mammalian cell culture. Methods in Biotechnology. Volume 24. Edited by: Poertner R. 2007, New York: Humana Press Inc, 253-268. 2Google Scholar
- Stephanopoulos GN, Aristidou AA, Nielsen J: Metabolic Engineering - Principles and Methodologies. 1998, Orlando: Academic PressGoogle Scholar
- Ozturk SS, Palsson BO: Chemical decomposition of glutamine in cell culture media: effect of media type, pH, and serum concentration. Biotechnol Prog. 1990, 6 (2): 121-128.View ArticleGoogle Scholar
- van der Heijden RT, Heijnen JJ, Hellinga C, Romein B, Luyben KC: Linear constraint relations in biochemical reaction systems: I. Classification of the calculability and the balanceability of conversion rates. Biotechnol Bioeng. 1994, 43 (1): 3-10.View ArticleGoogle Scholar
- Klamt S, Schuster S, Gilles ED: Calculability analysis in underdetermined metabolic networks illustrated by a model of the central metabolism in purple nonsulfur bacteria. Biotechnol Bioeng. 2002, 77 (7): 734-751.View ArticleGoogle Scholar
- Orth JD, Thiele I, Palsson BO: What is flux balance analysis?. Nat Biotechnol. 2010, 28 (3): 245-248.View ArticleGoogle Scholar
- Llaneras F, Pico J: An interval approach for dealing with flux distributions and elementary modes activity patterns. J Theor Biol. 2007, 246 (2): 290-308.View ArticleGoogle Scholar
- Mancuso A, Sharfstein ST, Tucker SN, Clark DS, Blanch HW: Examination of primary metabolic pathways in a murine hybridoma with carbon-13 nuclear magnetic resonance spectroscopy. Biotechnol Bioeng. 1994, 44 (5): 563-585.View ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.