Skip to main page content
U.S. flag

An official website of the United States government

Dot gov

The .gov means it’s official.
Federal government websites often end in .gov or .mil. Before sharing sensitive information, make sure you’re on a federal government site.

Https

The site is secure.
The https:// ensures that you are connecting to the official website and that any information you provide is encrypted and transmitted securely.

Access keys NCBI Homepage MyNCBI Homepage Main Content Main Navigation
. 2017 Sep 18;12(9):e0185085.
doi: 10.1371/journal.pone.0185085. eCollection 2017.

The Warburg effect as an adaptation of cancer cells to rapid fluctuations in energy demand

Affiliations

The Warburg effect as an adaptation of cancer cells to rapid fluctuations in energy demand

Tamir Epstein et al. PLoS One. .

Abstract

To maintain optimal fitness, a cell must balance the risk of inadequate energy reserve for response to a potentially fatal perturbation against the long-term cost of maintaining high concentrations of ATP to meet occasional spikes in demand. Here we apply a game theoretic approach to address the dynamics of energy production and expenditure in eukaryotic cells. Conventionally, glucose metabolism is viewed as a function of oxygen concentrations in which the more efficient oxidation of glucose to CO2 and H2O produces all or nearly all ATP except under hypoxic conditions when less efficient (2 ATP/ glucose vs. about 36ATP/glucose) anaerobic metabolism of glucose to lactic acid provides an emergency backup. We propose an alternative in which energy production is governed by the complex temporal and spatial dynamics of intracellular ATP demand. In the short term, a cell must provide energy for constant baseline needs but also maintain capacity to rapidly respond to fluxes in demand particularly due to external perturbations on the cell membrane. Similarly, longer-term dynamics require a trade-off between the cost of maintaining high metabolic capacity to meet uncommon spikes in demand versus the risk of unsuccessfully responding to threats or opportunities. Here we develop a model and computationally explore the cell's optimal mix of glycolytic and oxidative capacity. We find the Warburg effect, high glycolytic metabolism even under normoxic conditions, is represents a metabolic strategy that allow cancer cells to optimally meet energy demands posed by stochastic or fluctuating tumor environments.

PubMed Disclaimer

Conflict of interest statement

Competing Interests: The authors have declared that no competing interests exist.

Figures

Fig 1
Fig 1. Schematic representation of the demand-driven metabolic model.
ATP demand is composed of two types of demand; (dark) slow changing base-load demand, primarily for macromolecules synthesis and (light) peak demand, that rapidly changes primarily to support membrane transporters. In this model oxidative phosphorylation (high efficiency, slow response time) supplies base-load demand, and glycolytic metabolism (less efficient, fast response time) supplies peak ATP demand. The later is less efficient but is more highly responsive.
Fig 2
Fig 2. Typical dynamics of the system for different values of glycolytic capacity, cg.
The optimal level of glycolytic capacity occurs where the curves in Panels E, F and G become flat. The left panels show the temporal dynamics of (A) ATP demand, d(t), (B) glycolytic ATP production, cgPg, (C) oxidative ATP production, copo, and (D) ATP level UATP. Dash line denotes the minimum ATP level that is allowed in the system, U*ATP. Higher glycolytic capacity allows more glycolytic ATP to supply peak demand and reduces the maximum ATP level that is needed to ensure that UATPU*ATP. The right panels show (E) total glycolytic ATP production over one period, Pgtotal, (F) total oxidative ATP production over one period, Pototal, and (G) maximum ATP level over the period max(UATP). Parameters used in this simulation: Demand parameters: f = 10Hz, period = 4, BL = 5, km = 50, kg = 10000, U*ATP = 0.745, Pgmax = 5 and cm = 1.
Fig 3
Fig 3. Normalized cost terms.
(A) Normalized excessive ATP, Anorm, and (B) normalized glucose consumption, Fnorm, as function of normalized glycolytic capacity, cgnorm, for efficiency ration of fg/o = 2. Dots are overlay of all values between PL = 1, 1¼, 1½, …, 5. Red lines are quadratic polynomial fit, which is used for calculating production cost for any value of peak-demand and glycolytic capacity. Whether presented as panel (A) or (B), the optimal level of glycolytic capacity that minimizes costs occurs at the point where normalized excess ATP or normalized glucose consumption reaches 0.
Fig 4
Fig 4. Optimal ATP production mode for peak-demand.
The effects of the costs of excess ATP, a, glycolytic capacity, c, glucose, f, and the production efficiency ratio, fg/o, on whether it is optimal for the cell to use oxidative phosphorylation (o) or glycolysis (g) to meet peak ATP demand. As cost parameters change, there are abrupt shifts from one metabolic strategy to the other. The lines separating regions of oxidative versus glycolytic metabolism are called “isolegs” and their slopes and intercepts indicate how the cost parameters interact. Increasing the production efficiency ratio, fg/o, decreases the region over which glycolysis is optimal. Decreasing a or c generally favors glycolysis. The effect of f on the optimal metabolism for meeting peak demand depends on the relative glucose efficiency of oxidation versus glycolysis, fg/o. The values for all other parameters are the same as those in Fig 2.
Fig 5
Fig 5. Influence of peak-demand parameters on the steady-state population distribution.
Mean of population distribution as function power of the amplitude distribution, p, (bottom) and the mean time between events of demand, λ (left). The graph demonstrates that short time between bursts of peak-demand and high probability of high-amplitude demand results with a selection of population with high glycolytic capacity. Simulation parameters: kdrift = 0.05, Αhaz = 1, Αopp = 1, Αfix = 1 and m = 0.01. Detail description of the simulation process can be found in S3 Appendix.
Fig 6
Fig 6. Influence of opportunistic and hazardous coefficient on the steady-state population distribution.
Mean capacity that is selected as function of the opportunistic coefficient, αopp, and the hazardous coefficient, αhaz for fixed-cost coefficient, αfix = 1 (a) and αfix = 0.5 (b). The graphs demonstrate that hazard is more influential in selection of high glycolytic capacity than opportunity

References

    1. Nelson D, Lehninger A, Cox M. Lehninger principles of biochemistry. 2008;
    1. Hippocrates., Pare A, Harvey W, Jenner E, Holmes OW, Lister J, et al. Scientific papers: physiology, medicine, surgery, geology New York: P.F. Collier & Son; 1910.
    1. Warburg O, Posener K, Negelein E. Ueber den stoffwechsel der tumoren. Biochem Z. 1924;
    1. Hanahan D, Weinberg RA. Hallmarks of cancer: The next generation. Cell. 2011. pp. 646–674. doi: 10.1016/j.cell.2011.02.013 - DOI - PubMed
    1. Moreno-Sánchez R, Marín-Hernández A, Saavedra E, Pardo JP, Ralph SJ, Rodríguez-Enríquez S. Who controls the ATP supply in cancer cells? Biochemistry lessons to understand cancer energy metabolism. Int J Biochem Cell Biol. 2014;50: 10–23. doi: 10.1016/j.biocel.2014.01.025 - DOI - PubMed