Alteration of growth yield by overexpression of phosphoenolpyruvate carboxylase and phosphoenolpyruvate carboxykinase in Escherichia coli.

Phosphoenolpyruvate and oxaloacetate are key intermediates at the junction between catabolism and biosynthesis. Alteration of carbon flow at these branch points will affect the growth yield and the formation of products. We attempted to modulate the metabolic flow between phosphoenolpyruvate and oxaloacetate by overexpressing phosphoenolpyruvate carboxylase and phosphoenolpyruvate carboxykinase from a multicopy plasmid under the control of the tac promoter. It was found that overexpression of phosphoenolpyruvate carboxylase decreased the rates of glucose consumption and organic acid excretion, but the growth and respiration rates remained unchanged. Consequently, the growth yield on glucose was improved. This result indicates that the wild-type level of phosphoenolpyruvate carboxylase is not optimal for the most efficient glucose utilization in batch cultures. On the other hand, overexpression of phosphoenolpyruvate carboxykinase increased glucose consumption and decreased oxygen consumption relative to those levels required for growth. Therefore, the growth yield on glucose was reduced because of a higher rate of fermentation product excretion. These data provide useful insights into the regulation of central metabolism and facilitate further manipulation of pathways for metabolite production.     

Fundamental Escherichia coli biochemical pathways for biomass and energy production: identification of reactions

Cells grow by oxidizing nutrients using a complex network of biochemical reactions. During this process new biological material is produced along with energy used for maintaining cellular organization. Because the metabolic network is highly branched, these tasks can be accomplished using a wide variety of unique reaction sequences. However, evolutionary pressures under carbon-limited growth conditions likely select organisms that utilize highly efficient pathways. Using elementary-mode analysis, we demonstrate that the metabolism of the bacterium Escherichia coli contains four unique pathways that most efficiently convert glucose and oxygen into new cells and maintenance energy under any level of oxygen limitation. Observed regulatory patterns and experimental findings suggest growing cells use these highly efficient pathways. It is predicted that five knockout mutations generate a strain that supports growth using only the most efficient reaction sequence. The analysis approach should be generally useful for predicting metabolic capabilities and efficient network designs based on only genomic information.     

Metabolic analysis of the synthesis of high levels of intracellular human SOD in Saccharomyces cerevisiae rhSOD 2060 411 SGA122

The synthesis of human superoxide dismutase (SOD) in batch cultures of a Saccharomyces cerevisiae strain using a glucose-limited minimal medium was studied through metabolic flux analysis. A stoichiometric model was built, which included 78 reactions, according to metabolic pathways operative in these strains during respirofermentative and oxidative metabolism. It allowed calculation of the distribution of metabolic fluxes during diauxic growth on glucose and ethanol. Fermentation profiles and metabolic fluxes were analyzed at different phases of diauxic growth for the recombinant strain (P+) and for its wild type (P-). The synthesis of SOD by the strain P+ resulted in a decrease in specific growth rate of 34 and 54% (growth on glucose and ethanol respectively) in comparison to the wild type. Both strains exhibited similar flux of glucose consumption and ethanol synthesis but important differences in carbon distribution with biomass/substrate yields and ATP production 50% higher in P-. A higher contribution of fermentative metabolism, with 64% of the energy produced at the phosphorylation level, was observed during SOD production. The flux of precursors to amino acids and nucleotides was higher in the recombinant strain, in agreement with the higher total RNA and protein levels. Lower specific growth rates in strain P+ appear to be related to the decrease in the rate of synthesis of nonrecombinant protein, as well as a decrease in the activities of the pentose phosphate (PP) pathway and TCA cycle. A very different way of entry into the stationary phase was observed for each strain: in the wild-type strain most metabolic fluxes decreased and fluxes related to energy reserve synthesis increased, while in the P+ strain the flux of 22 reactions (including PP pathway and amino acids biosynthesis) related to SOD production increased their fluxes. Changes in SOD production rates at different physiological states appear to be related to the differences in building blocks availability between respirofermentative and oxidative metabolism. Using the present expression system, ideal conditions for SOD synthesis are represented by either active growth during respirofermentative metabolism or transition from a growing to a nongrowing state. An increase in SOD flux could be achieved using an expression system nonassociated to growth and potentially eliminating part of the metabolic burden.     

Physiological response of central metabolism in Escherichia coli to deletion of pyruvate oxidase and introduction of heterologous pyruvate carboxylase

We studied the physiological response of Escherichia coli central metabolism to the expression of heterologous pyruvate carboxylase (PYC) in the presence and absence of pyruvate oxidase (POX). These studies were complemented with expression analysis of central and intermediary metabolic genes and conventional in vitro enzyme assays to evaluate glucose metabolism at steady-state growth conditions (chemostats). The absence of POX activity reduced nongrowth-related energy metabolism (maintenance coefficient) and increased the maximum specific rate of oxygen consumption. The presence of PYC activity (i.e., with POX activity) increased the biomass yield coefficient and reduced the maximum specific oxygen consumption rate compared to the wildtype. The presence of PYC in a poxB mutant resulted in a 42% lower maintenance coefficient and a 42% greater biomass yield compared to the wildtype. Providing E. coli with PYC or removing POX increased the threshold specific growth rate at which acetate accumulation began, with an 80% reduction in acetate accumulation observed at a specific growth rate of 0.4 h-1 in the poxB-pyc+ strain. Gene expression analysis suggests utilization of energetically less favorable glucose metabolism via glucokinase and the Entner-Doudoroff pathway in the absence of functional POX, while the upregulation of the phosphotransferase glucose uptake system and several amino acid biosynthetic pathways occurs in the presence of PYC. The physiological and expression changes resulting from these genetic perturbations demonstrate the importance of the pyruvate node in respiration and its impact on acetate overflow during aerobic growth.     

Stoichiometric flux balance models quantitatively predict growth and metabolic by-product secretion in wild-type Escherichia coli W3110.

Flux balance models of metabolism use stoichiometry of metabolic pathways, metabolic demands of growth, and optimality principles to predict metabolic flux distribution and cellular growth under specified environmental conditions. These models have provided a mechanistic interpretation of systemic metabolic physiology, and they are also useful as a quantitative tool for metabolic pathway design. Quantitative predictions of cell growth and metabolic by-product secretion that are experimentally testable can be obtained from these models. In the present report, we used independent measurements to determine the model parameters for the wild-type Escherichia coli strain W3110. We experimentally determined the maximum oxygen utilization rate (15 mmol of O2 per g [dry weight] per h), the maximum aerobic glucose utilization rate (10.5 mmol of Glc per g [dry weight] per h), the maximum anaerobic glucose utilization rate (18.5 mmol of Glc per g [dry weight] per h), the non-growth-associated maintenance requirements (7.6 mmol of ATP per g [dry weight] per h), and the growth-associated maintenance requirements (13 mmol of ATP per g of biomass). The flux balance model specified by these parameters was found to quantitatively predict glucose and oxygen uptake rates as well as acetate secretion rates observed in chemostat experiments. We have formulated a predictive algorithm in order to apply the flux balance model to describe unsteady-state growth and by-product secretion in aerobic batch, fed-batch, and anaerobic batch cultures. In aerobic experiments we observed acetate secretion, accumulation in the culture medium, and reutilization from the culture medium. In fed-batch cultures acetate is cometabolized with glucose during the later part of the culture period.(ABSTRACT TRUNCATED AT 250 WORDS)     

A metabolic model of cellular energetics and carbon flux during aerobic Escherichia coli fermentation

An integrated metabolic model for the production of acetate by Escherichia coli growing on glucose under aerobic conditions was presented previously (Ko et al., 1993). The resulting model equations can be used to explain phenomena often observed with industrial fermentations, i.e., increased acetate production which follows from high glucose uptake rate, a low dissolved oxygen concentration, a high specific growth rate, or a combination of these conditions. However, several questions still need to be addressed. First, cell composition is growth rate and media dependent. Second, the macromolecular composition varied between E. coli strains. And finally, a model that represents the carbon fluxes between the Embden-Meyerhof-Parnas (EMP) and the hexose monophosphate (HMP) pathways when cells are subject to internal and/or external stresses is still not well defined. In the present work, we have made an effort to account for these effects, and the resulting model equations show good agreement for wild-type and recombinant E. coli experimental data for the acetate concentration, the onset of acetate secretion, and cell yield based on glucose. These results are useful for optimizing aerobic E. coli fermentation processes. More specifically, we have determined the EMP pathway carbon flux profiles required by the integrated metabolic model for an accurate fit of the acetic acid profile data from a wild-type E. coli strain ML308. These EMP carbon flux profiles were correlated with a dimensionless measurement of biomass and then used to predict the acetic acid profiles for E. coli strain F-122 expressing human immunodeficiency virus-(HIV(528)) beta-galactosidase fusion protein. The effect of different macromolecular compositions and growth rates between these two E. coli strains required a constant scaling factor for improved quantitative predictions.     

Fluxes of carbon, phosphorylation, and redox intermediates during growth of saccharomyces cerevisiae on different carbon sources

In the present work we develop a method for estimating anabolic fluxes when yeast are growing on various carbon substrates (glucose, glycerol, lactate, pyruvate, acetate, or ethanol) in minimal medium. Fluxes through the central amphibolic pathways were calculated from the product of the total required amount of a specified carbon intermediate times the growth rate. The required amount of each carbon intermediate was estimated from the experimentally determined macromolecular composition of cells grown in each carbon source and the monomer composition of macromolecules.Substrates sharing most metabolic pathways such as ethanol and acetate, despite changes in the macromolecular composition, namely carbohydrate content (34% +/- 1 and 21% +/- 3, respectively), did not show large variations in the overall fluxes through the main amphibolic pathways. For instance, in order to supply anabolic precursors to sustain growth rates in the range of 0.16/h to 0.205/h, similar large fluxes through Acetyl CoA synthase were required by acetate (4.2 mmol/hr g dw) or ethanol (5.2 mmol/h g dw).The V(max) activities of key enzymes of the main amphibolic pathways measured in permeabilized yeast cells allowed to confirm, qualitatively, the operation of those pathways for all substrates and were consistent on most substrates with the estimated fluxes required to sustain growth.When ATP produced from oxidation of the NADH synthesized along with the key intermediary metabolites was taken into account, higher Y(ATP) (max) values (36 with respect to 24 g dw/mol ATP) were obtained for glucose. The same result was obtained for glycerol, ethanol, and acetate. A yield index (YI) was defined as the ratio of the theoretically estimated substrate flux required to sustain a given growth rate over the experimentally measured flux of substrate consumption. Comparison of Yl between growth on various carbon sources led us to conclude that ethanol (Yl = 0.84), acetate (Yl = 0.77), and lactate (Yl = 0.77) displayed the most efficient use of substrate for biomass production. For the other substrates, the Yl decayed in the following order: pyruvate > glycerol > glucose.An improvement of the quantitative understanding of yeast metabolism, energetics, and physiology is provided by the present analysis. The methodology proposed can be applied to other eukaryotic organisms of known chemical composition. (c) 1995 John Wiley & Sons, Inc.     

基于代谢流量分布信息理解大肠杆菌中异柠檬酸裂解酶调节因子的代谢调控作用

基因的表达受不同的转录调节因子调节。大肠杆菌中的异柠檬酸裂解酶调节因子(IclR)能够抑制编码乙醛酸支路酶的aceBAK操纵子的表达。本研究基于代谢物的13C同位体物质分布来定量解析代谢反应,主要研究了iclR基因在大肠杆菌生理和代谢中的作用。大肠杆菌iclR基因缺失突变株的生长速率、糖耗速率和乙酸的产量相对于原始菌株都有所降低,但菌体得率略有增加。通过代谢途径的流量比率分析发现基因缺失株的乙醛酸支路得到了激活,33%的异柠檬酸流经了乙醛酸支路;戊糖磷酸途径的流量变小,使得CO2的生成量减少。同时,乙醛酸支路激活,但草酰乙酸形成磷酸烯醇式丙酮酸的流量基本不变,说明磷酸烯醇式丙酮酸-乙醛酸循环没有激活,没有过多的碳原子在磷酸烯醇式丙酮酸羧化激酶反应中以CO2形式排出,从而确保了菌体得率。葡萄糖利用速率的降低、乙酰辅酶A的代谢效率提高等使得iclR基因敲除菌的乙酸分泌较原始菌株有所降低。     

Icilin induces G1 arrest through activating JNK and p38 kinase in a TRPM8-independent manner

Mechanistic models of glucose stimulated insulin secretion (GSIS) established in minimal media in vitro, may not accurately describe the complexity of coupling metabolism with insulin secretion that occurs in vivo. As a first approximation, we have evaluated metabolic pathways in a typical growth media, DMEM as a surrogate in vivo medium, for comparison to metabolic fluxes observed under the typical experimental conditions using the simple salt-buffer of KRB. Changes in metabolism in response to glucose and amino acids and coupling to insulin secretion were measured in INS-1 832/13 cells. Media effects on mitochondrial function and the coupling efficiency of oxidative phosphorylation were determined by fluorometrically measured oxygen consumption rates (OCRs) combined with ³¹P NMR measured rates of ATP synthesis. Substrate preferences and pathways into the TCA cycle, and the synthesis of mitochondrial 2nd messengers by anaplerosis were determined by ¹³C NMR isotopomer analysis of the fate of [U-¹³C] glucose metabolism.Despite similar incremental increases in insulin secretion, the changes of OCR in response to increasing glucose from 2.5 to 15 mM were blunted in DMEM relative to KRB. Basal and stimulated rates of insulin secretion rates were consistently higher in DMEM, while ATP synthesis rates were identical in both DMEM and KRB, suggesting greater mitochondrial uncoupling in DMEM. The relative rates of anaplerosis, and hence synthesis and export of 2nd messengers from the mitochondria were found to be similar in DMEM to those in KRB. And, the correlation of total PC flux with insulin secretion rates in DMEM was found to be congruous with the correlation in KRB. Together, these results suggest that signaling mechanisms associated with both TCA cycle flux and with anaplerotic flux, but not ATP production, may be responsible for the enhanced rates of insulin secretion in more complex, and physiologically-relevant media.     

Effect of oxygen on steady-state product distribution in Bacillus polymyxa fermentations

Bacillus polymyxa ferments glucose to 1-2,3 butanediol, acetoin, ethanol, acetic acid, lactic acid, and formic acid. This research investigates product formation as a function of oxygen availability. A predictive model that simulates product distribution at known oxygen transfer rates is developed on the hypothesis that, in an energy-limited environment, B. polymyxa utilizes glucose and oxygen in the most efficient manner. The efficiency of utilization of glucose and oxygen is measured in terms of the ATP yields of each oxidative pathway. The identity of the products constituting the profile at the given oxygen transfer rate is determined by comparing the ATP production and consumption rates. While the ATP generated is calculated from a knowledge of the oxygen transfer rate and ATP yields of the oxidative pathways, the ATP consumption is estimated by the Pirt expression in terms of growth- and nongrowth-associated components. The product formation rates are obtained by solving ATP and NAD balance equations. They equate the production and consumption rates of these intermediates and are derived from the pseudo-steady-state hypothesis. The model is applied to continuous culture systems that are both open and closed with respect to biomass. At a given oxygen transfer rate, dilution rate, and inlet glucose concentration, the model predicts steady-state concentrations of two dominant fermentation endproducts with the help of four parameters that can be determined from independent experiments. In contrast with earlier approaches, the experimental studies are carried out in continuous culture. Product profiles are obtained at various oxygen transfer rates, fer rates, inlet glucose concentrations, and dilution rates. The effect of pH on the relative distribution of products is also demonstrated. Results indicate that the model is fairly successful in predicting product profiles as a function of oxygen availability. (c) 1992 John Wiley & Sons, Inc.     

Enhanced antibody production following intermediate addition based on flux analysis in mammalian cell continuous culture

Generally, mammalian cells utilize glucose and glutamine as primary energy sources. To investigate the effect of energy sources on metabolic fluxes and antibody production, glucose- or glutamine-limited serum-free continuous culture of hybridoma 3A21 cells, which produce anti-ribonuclease A antibody, was carried out. The cell volume and dry cell weight were evaluated under various steady-state conditions. The specific consumption and production rates were evaluated on the basis of dry cell weight. On the basis of these results, the fluxes of the metabolic pathway were calculated. It was found that increasing the specific growth rate causes the specific ATP and antibody production rates to decrease. The fluxes between malate and pyruvate also decreased with the increase in specific growth rate. To increase the ATP production rate under steady-state conditions by the enhancement of fluxes between malate and pyruvate, the reduced metabolic fluxes were increased by an intermediate (pyruvate, malate, and citrate) addition. As a result, higher specific ATP and antibody production rates were achieved following the intermediate addition at a constant dilution rate.     

Response of fluxome and metabolome to temperature-induced recombinant protein synthesis in Escherichia coli

The response of the central carbon metabolism of Escherichia coli to temperature-induced recombinant production of human fibroblast growth factor was studied on the level of metabolic fluxes and intracellular metabolite levels. During production, E. coli TG1:plambdaFGFB, carrying a plasmid encoded gene for the recombinant product, revealed stress related characteristics such as decreased growth rate and biomass yield and enhanced by-product excretion (acetate, pyruvate, lactate). With the onset of production, the adenylate energy charge dropped from 0.85 to 0.60, indicating the occurrence of a severe energy limitation. This triggered an increase of the glycolytic flux which, however, was not sufficient to compensate for the increased ATP demand. The activation of the glycolytic flux was also indicated by the readjustment of glycolytic pool sizes leading to an increased driving force for the reaction catalyzed by phosphofructokinase. Moreover, fluxes through the TCA cycle, into the pentose phosphate pathway and into anabolic pathways decreased significantly. The strong increase of flux into overflow pathways, especially towards acetate was most likely caused by a flux redirection from pyruvate dehydrogenase to pyruvate oxidase. The glyoxylate shunt, not active during growth, was the dominating anaplerotic pathway during production. Together with pyruvate oxidase and acetyl CoA synthase this pathway could function as a metabolic by-pass to overcome the limitation in the junction between glycolysis and TCA cycle and partly recycle the acetate formed back into the metabolism.     

Metabolic flux analysis gives an insight on verapamil induced changes in central metabolism of HL-1 cells

Verapamil has been shown to inhibit glucose transport in several cell types. However, the consequences of this inhibition on central metabolism are not well known. In this study we focused on verapamil induced changes in metabolic fluxes in a murine atrial cell line (HL-1 cells). These cells were adapted to serum free conditions and incubated with 4 μM verapamil and [U-¹³C₅] glutamine. Specific extracellular metabolite uptake/production rates together with mass isotopomer fractions in alanine and glutamate were implemented into a metabolic network model to calculate metabolic flux distributions in the central metabolism. Verapamil decreased specific glucose consumption rate and glycolytic activity by 60%. Although the HL-1 cells show Warburg effect with high lactate production, verapamil treated cells completely stopped lactate production after 24 h while maintaining growth comparable to the untreated cells. Calculated fluxes in TCA cycle reactions as well as NADH/FADH₂ production rates were similar in both treated and untreated cells. This was confirmed by measurement of cell respiration. Reduction of lactate production seems to be the consequence of decreased glucose uptake due to verapamil. In case of tumors, this may have two fold effects; firstly depriving cancer cells of substrate for anaerobic glycolysis on which their growth is dependent; secondly changing pH of the tumor environment, as lactate secretion keeps the pH acidic and facilitates tumor growth. The results shown in this study may partly explain recent observations in which verapamil has been proposed to be a potential anticancer agent. Moreover, in biotechnological production using cell lines, verapamil may be used to reduce glucose uptake and lactate secretion thereby increasing protein production without introduction of genetic modifications and application of more complicated fed-batch processes.     

Pyruvate and oxygen consumption throughout the growth and development of murine oocytes

Growing oocytes in vitro from the most immature stages until they are developmentally competent is a major goal of reproductive technology, requiring fundamental knowledge of metabolic processes. Carbohydrate metabolism and oxygen consumption have been analysed in a series of experiments designed to investigate important energy substrates for mouse oocytes and to reveal any qualitative or quantitative changes between the primordial and ovulatory follicle stages. Primordial follicles were incubated in groups in modified-KSOM medium, whereas growing or ovulated oocytes were studied singly and, in both cases, the depletion or accumulation of metabolites in spent medium were analysed using ultramicrofluorometric assays. The rates of glucose (0.014 +/- 0.006 pmol/hr) and pyruvate (0.028 +/- 0.009 pmol/hr) consumption and l-lactate (0.058 +/- 0.023 pmol/hr) production by primordial follicles suggested that energy production was supported by a combination of metabolic pathways, including glycolysis. Pyruvate and oxygen consumption per oocyte increased two- and ninefold, respectively, between the primary and pre-ovulatory stages (0.82 +/- 0.1 and 1.67 +/- 0.1 pmol pyruvate/hr, respectively and 1.4 +/- 0.3 and 7 +/- 0.6 pmol oxygen/hr) after which oxygen (12.7 +/- 1.1 pmol/hr) utilisation nearly doubled. Oxygen consumption by fully grown oocytes was in excess of oxidation requirements for pyruvate. When pyruvate and oxygen consumption rates were normalised for oocyte cellular volume, which increased over 130-fold during growth, oocyte metabolism was higher in primary follicles than at any subsequent stage, indicating that energy needs are greater during a developmental transition. To conclude, pyruvate and oxygen were consumed throughout oocyte development at increasing rates. When oocyte cellular volume was accounted for, oocytes from primary follicles displayed greatest metabolic rates.     

Design, construction and performance of the most efficient biomass producing E. coli bacterium

Inverse metabolic engineering based on elementary mode analysis was applied to maximize the biomass yield of Escherchia coli MG1655. Elementary mode analysis was previously employed to identify among 1691 possible pathways for cell growth the most efficient pathway with maximum biomass yield. The metabolic network analysis predicted that deletion of only 6 genes reduces the number of possible elementary modes to the most efficient pathway. We have constructed a strain containing these gene deletions and we evaluated its properties in batch and in chemostat growth experiments. The results show that the theoretical predictions are closely matched by the properties of the designed strain.     

Growth and metabolism of marine fish Chinook salmon embryo cells: response to lack of glucose and glutamine

A peculiar phenomenon, differing from the response of mammalian cells, occurred when Chinook salmon embryo (CHSE) cells were passaged in the medium lacking of both glucose and glutamine. To elucidate metabolic mechanism of CHSE cells, the metabolism parameters, key metabolic enzymes, and ATP levels were measured at different glucose and glutamine concentrations. In the glutamine-free culture, hexokinase activity kept constant, and lactate dehydrogenase (LDH) activity decreased. This indicated that lack of glutamine did not expedite glucose consumption but made it shift to lower lactate production and more efficient energy metabolism. The results coincided with the experimental results of unaltered specific glucose consumption rate and decreased yield coefficients of lactate to glucose. In the glucose-free culture, simultaneous increase of glutaminase activity and of specific ammonia production rate suggested an increased flux into the glutaminolysis pathway, and increases of both glutamate dehydrogenase activity and yield coefficient of ammonia to glutamine showed an increased flux into deamination pathway. However, when glucose and glutamine were both lacking, the specific consumption rates of most of amino acids increased markedly, together with decrease of LDH activity, indicating that pyruvate derived from amino acids, away from lactate production, remedied energy deficiency. When both glucose and glutamine were absent, intracellular ATP contents and the energy charge remained virtually unaltered.Revisions requested 16 December 2004; Revisions received 24 January 2005     

Effects of dissolved oxygen on hybridoma cell growth, metabolism, and antibody production kinetics in continuous culture.

The effects of dissolved oxygen concentration (DO) on hybridoma cell physiology were examined in a continuous stirred tank bioreactor with a murine hybridoma cell line (167.4G5.3). Dissolved oxygen concentration was varied between 0% and 100% air saturation. Cell growth and viability, carbohydrate, amino acid, and energy metabolism, oxygen uptake, and antibody production rates were investigated. Cell growth was inhibited at both high and low DO. Cells could grow at 0% DO and maintain viability under a nitrogen atmosphere. Cell viability was higher at low DO. Glucose, glutamine, and oxygen consumption rates changed little at DO above 1% air saturation. However, the metabolic uptake rates changed below 1% DO, where growth became oxygen limited, and a Km value of 0.6% DO was obtained for the specific oxygen uptake rate. The metabolic rates of glucose, glutamine, lactate, and ammonia increased 2-3-fold as the DO dropped from 1% to 0%. Amino acid metabolism followed the same general pattern as that of glutamine and glucose. Alanine was the only amino acid produced. The consumption rates of amino acids changed little above 1% DO, but under anaerobic conditions the consumption rates of all amino acids increased severalfold. Cells obtained most of their metabolic energy from glutamine oxidation except under oxygen limitation, when glucose provided most of the energy. The calculated ATP production rate was only slightly influenced by DO and rose at 0% DO. Antibody concentration was highest at 35% DO, while the specific antibody production rate was insensitive to DO.     

Continuous modeling of metabolic networks with gene regulation in yeast and in vivo determination of rate parameters

A continuous model of a metabolic network including gene regulation to simulate metabolic fluxes during batch cultivation of yeast Saccharomyces cerevisiae was developed. The metabolic network includes reactions of glycolysis, gluconeogenesis, glycerol and ethanol synthesis and consumption, the tricarboxylic acid cycle, and protein synthesis. Carbon sources considered were glucose and then ethanol synthesized during growth on glucose. The metabolic network has 39 fluxes, which represent the action of 50 enzymes and 64 genes and it is coupled with a gene regulation network which defines enzyme synthesis (activities) and incorporates regulation by glucose (enzyme induction and repression), modeled using ordinary differential equations. The model includes enzyme kinetics, equations that follow both mass-action law and transport as well as inducible, repressible, and constitutive enzymes of metabolism. The model was able to simulate a fermentation of S. cerevisiae during the exponential growth phase on glucose and the exponential growth phase on ethanol using only one set of kinetic parameters. All fluxes in the continuous model followed the behavior shown by the metabolic flux analysis (MFA) obtained from experimental results. The differences obtained between the fluxes given by the model and the fluxes determined by the MFA do not exceed 25% in 75% of the cases during exponential growth on glucose, and 20% in 90% of the cases during exponential growth on ethanol. Furthermore, the adjustment of the fermentation profiles of biomass, glucose, and ethanol were 95%, 95%, and 79%, respectively. With these results the simulation was considered successful. A comparison between the simulation of the continuous model and the experimental data of the diauxic yeast fermentation for glucose, biomass, and ethanol, shows an extremely good match using the parameters found. The small discrepancies between the fluxes obtained through MFA and those predicted by the differential equations, as well as the good match between the profiles of glucose, biomass, and ethanol, and our simulation, show that this simple model, that does not rely on complex kinetic expressions, is able to capture the global behavior of the experimental data. Also, the determination of parameters using a straightforward minimization technique using data at only two points in time was sufficient to produce a relatively accurate model. Thus, even with a small amount of experimental data (rates and not concentrations) it was possible to estimate the parameters minimizing a simple objective function. The method proposed allows the obtention of reasonable parameters and concentrations in a system with a much larger number of unknowns than equations. Hence a contribution of this study is to present a convenient way to find in vivo rate parameters to model metabolic and genetic networks under different conditions.     

Metabolic and morphological differences between rapidly proliferating cancerous and normal breast epithelial cells

The metabolic and morphological characteristics of two human epithelial breast cell populations--MCF7 cells, a cancerous cell line, and 48R human mammary epithelial cells (48R HMECs), a noncancerous, finite lifespan cell strain--were compared at identical growth rates. Both cell types were induced to grow rapidly in nutrient-rich media containing 13C-labeled glucose, and the isotopic enrichment of cellular metabolites was quantified to calculate metabolic fluxes in key pathways. Despite their similar growth rates, the cells exhibited distinctly different metabolic and morphological profiles. MCF7 cells have an 80% smaller exposed surface area and contain 26% less protein per cell than the 48R cells. Surprisingly, rapidly proliferating 48R cells exhibited a 225% higher per-cell glucose consumption rate, a 250% higher per-cell lactate production rate, and a nearly identical per-cell glutamine consumption rate relative to the cancer cell line. However, when fluxes were considered on the basis of exposed area, the cancer cells were observed to have higher glucose, lactate, and glutamine fluxes, demonstrating superior transport capabilities per unit area of cell membrane. MCF7 cells also consumed amino acids at rates much higher than are generally required for protein synthesis, whereas 48R cells generally did not. Pentose phosphate pathway activity was higher in MCF7 cells, and the flux of glutamine to glutamate was less reversible. Energy efficiency was significantly higher in MCF7 cells, as a result of a combination of their smaller size and greater reliance on the TCA cycle than the 48R cells. These observations support evolutionary models of cancer cell metabolism and suggest targets for metabolic drugs in metastatic breast cancers.     

Bonded Cumomer Analysis of Human Melanoma Metabolism Monitored by 13C NMR Spectroscopy of Perfused Tumor Cells

A network model for the determination of tumor metabolic fluxes from ¹³C NMR kinetic isotopomer data has been developed and validated with perfused human DB-1 melanoma cells carrying the BRAF V600E mutation, which promotes oxidative metabolism. The model generated in the bonded cumomer formalism describes key pathways of tumor intermediary metabolism and yields dynamic curves for positional isotopic enrichment and spin-spin multiplets. Cells attached to microcarrier beads were perfused with 26 mm [1,6-¹³C₂]glucose under normoxic conditions at 37 °C and monitored by ¹³C NMR spectroscopy. Excellent agreement between model-predicted and experimentally measured values of the rates of oxygen and glucose consumption, lactate production, and glutamate pool size validated the model. ATP production by glycolytic and oxidative metabolism were compared under hyperglycemic normoxic conditions; 51% of the energy came from oxidative phosphorylation and 49% came from glycolysis. Even though the rate of glutamine uptake was ∼50% of the tricarboxylic acid cycle flux, the rate of ATP production from glutamine was essentially zero (no glutaminolysis). De novo fatty acid production was ∼6% of the tricarboxylic acid cycle flux. The oxidative pentose phosphate pathway flux was 3.6% of glycolysis, and three non-oxidative pentose phosphate pathway exchange fluxes were calculated. Mass spectrometry was then used to compare fluxes through various pathways under hyperglycemic (26 mm) and euglycemic (5 mm) conditions. Under euglycemic conditions glutamine uptake doubled, but ATP production from glutamine did not significantly change. A new parameter measuring the Warburg effect (the ratio of lactate production flux to pyruvate influx through the mitochondrial pyruvate carrier) was calculated to be 21, close to upper limit of oxidative metabolism.