Cancer proliferation and therapy: the Warburg effect and quantum metabolism
© Demetrius et al; licensee BioMed Central Ltd. 2010
Received: 3 November 2009
Accepted: 19 January 2010
Published: 19 January 2010
Most cancer cells, in contrast to normal differentiated cells, rely on aerobic glycolysis instead of oxidative phosphorylation to generate metabolic energy, a phenomenon called the Warburg effect.
Quantum metabolism is an analytic theory of metabolic regulation which exploits the methodology of quantum mechanics to derive allometric rules relating cellular metabolic rate and cell size. This theory explains differences in the metabolic rates of cells utilizing OxPhos and cells utilizing glycolysis. This article appeals to an analytic relation between metabolic rate and evolutionary entropy - a demographic measure of Darwinian fitness - to: (a) provide an evolutionary rationale for the Warburg effect, and (b) propose methods based on entropic principles of natural selection for regulating the incidence of OxPhos and glycolysis in cancer cells.
The regulatory interventions proposed on the basis of quantum metabolism have applications in therapeutic strategies to combat cancer. These procedures, based on metabolic regulation, are non-invasive, and complement the standard therapeutic methods involving radiation and chemotherapy
Cancer is an age-dependent disease characterized by five key hallmarks in cell physiology that drive the progressive change of normal differentiated cells into diverse states of malignancy : autonomous growth-replication in the absence of growth signals; insensitivity to anti-growth signals; apoptosis-evasion of programmed cell death, angiogenesis-the induction of the growth of new blood vessels; invasion and metastasis
The age-dependency of cancer  and the relatively rare incidence of the disease during an average human life time suggest that adaptive mechanisms exist in cells and tissues to prevent this multi-step transition from a normal differentiated cell into malignancy. Consequently, each of these physiological changes constitutes the rupture of anti-tumor defences developed during the evolutionary history of the organism. It may be worth observing that a graph of the logarithm of the total cancer incidence against age approximates to a straight line with a gradient of 6-7 (the value of the power-law exponent) suggesting that 6-7 separate events are required for neoplastic transformation of a 'typical' human cell .
Cancer cells may be considered as autonomous units which have an impaired capacity to maintain the metabolic stability of the organism in which they reside. Anti-cancer therapies are corrective measures designed to remedy this impairment by eliminating the errant cells.
The nature of these corrective measures has undergone a significant development, starting with surgical resection of solid tumors followed by radiation and then chemotherapy, as the understanding of the biology of cancer has increased. The non-surgical therapies which came to dominate the treatment of the disease were based on the proposition that the disease is primarily the result of dynamic changes in the genome. This gene oriented perspective led to the notion that decoding the genetic instruction that determines the cancer phenotype would elucidate its origin and thus provide an effective biological basis for therapy. Genes represent the blueprint for phenotypic expression. Accordingly, the genetic model entailed therapies based on the complete elimination of cancer cells. Radiation and chemotherapy were the first class of anti-cancer strategies which this model invoked. These two therapeutic methods were designed to eliminate cancer cells from tissues. However, due to the low selectivity of this approach, non-cancer cells are also killed or damaged leading to severe side effects.
The next generation of cancer drugs which developed from this genomic viewpoint explicitly recognized the multi-step progression towards malignancy, and that in most instances death only occurs when the metastatic state is attained. Metastasis is often triggered by angiogenesis, the proliferation of a network of blood vessels that penetrates into cancerous tissue, supplying nutrients and oxygen .
Drugs that impede the formation of tumor blood vessels were therefore proposed as therapeutic agents in the combat against malignancy . There exist, however, some disadvantages to this mode of therapy as angiogenic inhibitors sometimes trigger side effects and induce a more invasive type of tumor .
Studies in recent years have led to a re-evaluation of the genomic model of cancer and the development of a model based on cell metabolism . The research which triggered this shift from genes to metabolic reactions was done in 1924 by Warburg  who recognized certain critical differences between energy regulation in normal differentiated cells and cancer cells. Warburg analyzed the ratio of oxidative phosphorylation (OxPhos) to glycolysis in different tissues of cancer cells and normal cells. Glycolysis under aerobic conditions was found to be particularly high in aggressive tumors when compared with benign tumors and normal tissues. These observations led Warburg to propose deficiency in OxPhos and elevated glycolysis as the primary cause of cancer.
The discovery of the double helix by Watson and Crick in 1953 and its implications for the understanding of molecular processes in biology diverted interest away from research into the significance of Warburg's metabolic hypothesis. However, the failure of the genomic approach to provide effective therapies for certain types of aggressive cancer, and recent studies  creating a rapprochement of genetic and metabolic views have revived interest in the Warburg hypothesis.
The hypothesis, in its simplest form, asserts that cancer is primarily a disease of metabolic dysregulation: a switch, inducible by various agents-genetic, nutritional and environmental, from an OxPhos pathway to a glycolytic mode of energy processing. This focus on metabolism as the primary cause for the progressive transition from normalcy to malignancy suggests a radically new approach to cancer therapy. The focus is to influence metabolic regulation in cancer cells so that the autonomy, proliferative capacity, invasiveness and metastasis which define the aggressive cancer phenotype, is never attained.
A therapeutic strategy which emphasizes containment rather than annihilation is a radical departure from methods designed to eradicate cancer cells completely from tissues. Therapies based on metabolic interventions involve two complementary programs: the down-regulation of glycolysis and the up-regulation of OxPhos.
According to Warburg the aggressiveness of a tumor derives from the elevation of the glycolytic mode of energy processing. This elevation is considered to be the result of competition between cells utilizing the glycolytic mode, and cells adopting OxPhos. Hence, the principle that underlies these complementary programs of disease control is Darwinian: the modulation of the selective advantage of cells using glycolysis and OxPhos, respectively.
It is well known that ATP generation through glycolysis is less efficient than through mitochondrial respiration. Hence a long-standing paradox is how cancer cells with their metabolic disadvantage can survive the competition with normal cells. In terms of biochemical reactions , mitochondrial respiration defects lead to activation of the Akt survival pathway through a mechanism mediated by NADH. Respiration-deficient cells harboring mitochondrial defects exhibit dependency on glycolysis, increased NADH, and activation of Akt, leading to survival advantage and also drug resistance in hypoxia .
The efforts to implement these therapeutic programs have generated certain significant questions regarding the analytical characterization of Darwinian fitness, the capacity of a cell type to displace related types in competition for resources. The problem can be formulated as follows: (a) What class of physiological, biochemical and biophysical properties of cells confers a selective advantage during evolution of the cancer phenotype? (b) To what extent can these properties be analytically described in terms of bio-energetic and kinetic variables associated with the regulatory circuits that describe the metabolic networks?
These questions are consistent with the view that cancer development proceeds according to an evolutionary process in which a succession of genetic and epigenetic changes, due to the selective advantage conferred, leads to the progressive transformation of normal cells to cancer cells [10–12]. The resolution of the problems addressed in (a) and (b) would evidently provide an analytic framework for cancer therapy based on containment rather than eradication. The analysis would also yield a rationale based on natural selection for Warburg's hypothesis, and consequently, an evolutionary understanding of the origin of cancer.
The analytical framework for a theory of metabolic regulation in cells capable of addressing problems of cellular adaptation and somatic evolution within living organism was proposed in a series of articles [13, 14], and called Quantum Metabolism in view of the quantum mechanics methodology the theory invoked. Quantum Metabolism gives a molecular level explanation of certain empirically derived allometric laws relating metabolic rate with cell size.
Here, P is the metabolic rate, the rate of ATP production, W the cell size. The dimensionality parameter d in the scaling exponent, β = d/(d+1), describes the number of degrees of freedom of the enzymes which catalyze the redox reactions within the energy transducing organelles: mitochondria (in the case of OxPhos), metabolosomes (in the case of glycolysis). The proportionality constant, α, depends on the mode of coupling between the electron transport chain and ADP phosphorylation. This mode of coupling is electrical in the case of OxPhos and chemical in the case of glycolysis.
A mathematical theory of evolution by natural selection which, for the first time, considered the effect of resource constraints and finite population size on the outcome of competition between related types, was described in . The cornerstone of this model was the statistical parameter, evolutionary entropy, a measure of the stability of population numbers, and an index of Darwinian Fitness.
In this article we will integrate Quantum Metabolism with certain analytic relations between evolutionary entropy and metabolic rate to show that selective advantage in cellular evolution is contingent on the resource constraints - its abundance and distribution and predicted by the metabolic rate, the rate at which cells transform resources into metabolic work. Quantum Metabolism predicts that the metabolic rate of cells utilizing OxPhos and cells utilizing glycolysis will have the same scaling exponents but will differ in terms of the proportionality constants. We will exploit this prediction, and the characterization of selective advantage in terms of resource constraints and metabolic rate, to provide an evolutionary rationale for the cancer phenotype. The evolutionary argument rests on differences in the metabolic rate of cells utilizing OxPhos and glycolytic pathways, respectively. Cellular metabolic rate can be influenced by perturbing the geometry of the metabolic network or the mode of coupling hence the incidence of OxPhos and glycloysis in normal and cancer cells can be metabolically regulated. We will appeal to these notions of metabolic intervention to propose anti-cancer strategies based on arresting the transformation from a benign tumor, to a malignant tissue characterized by enhanced glycolysis.
This article will give a brief account of Quantum Metabolism and the scaling laws for metabolic rate which the new theory derives. We discuss the evolutionary perspective this theory entails, and then apply the new class of models to propose therapies based on altering the selective advantage of normal and cancer cells during the transition to the cancer phenotype.
Cellular metabolism is the totality of all chemical reactions in cells carried out by an organism. The characteristics of living organisms, such as their growth, the maintenance of their structure and mass transport, depend on the input of energy from the environment. Metabolism designates the series of chemical reactions that transform substrates such as glucose into cellular building blocks, and energy in the form of ATP. A quantitative understanding of the rules which regulate this energy transformation is critical for a quantitative characterization of the selective advantage associated with different modes of energy processing.
Quantum Metabolism exploits the methodology of the quantum theory of solids, as developed by Einstein and Debye, to derive a class of analytic rules relating metabolic rate, with cell size. The variables which define these rules are bio-energetic parameters, whose values depend on the phospholipid composition of the bio-membranes, and enzymatic reaction rates, which depend on the concentration of substrates in metabolic reactions.
The allometric laws of metabolism
Metabolic rate, the rate at which an organism transforms nutrients into thermal energy and biological work required for sustaining life, is highly dependent on the organism's size. This observation draws in large part from the experimental studies of Lavoisier and Laplace who were the first to demonstrate the relationship between combustion and respiration. The systematic empirical study of the relation between metabolic rate, P, and body size, W, which began with Rubner , and was extended later by Kleiber  for non-domesticated mammals, and by Hemmigsen  to uni-cells, led to a set of allometric laws relating metabolic rate with body size.
Quantum Metabolism exploits the formalism of quantum mechanics to study the dynamics of energy flow in the electron transfer chain in cells, and provides a molecular level explanation of the empirical rules documented in Kleiber  and Hemmingsen .
(1) The redox potential difference, that is the actual redox potential between the donor and acceptor couples in the electron transport chain.
(2) The phosphorylation potential for ATP synthesis.
The parameter d in the scaling exponent β = d/(d+1), characterizes the number of degrees of freedom of the enzymes which are embedded in the energy transducing organelles. Enzymatic reactions have an intrinsic direction and the enzymes localized in the organelles have a given orientation. The enzymes are subject to oscillations, due to the redox potential. We will assume that these vibrations can be approximated by harmonic oscillators. We will also assume that the enzymatic vibrations are coupled and inherited by the energy transducing organelles in which the enzymes reside.
There exists a diverse body of empirical support for these assumptions. Some of this support can be annotated as follows. Experimental studies of the organization of mitochondrial networks show that these systems can be regarded as coupled oscillators . Synchronisation of metabolic cycles through gene and enzyme regulation within and between cells has been shown to involve co-ordinated transcriptional cycles not only in cultured yeast but also importantly in mammalian cells [22–24]. Tsong and collaborators have demonstrated that dynamical processes govern the function of metabolic enzymes which can capture and transmit energy from oscillating electric fields  involving electro-conformational coupling  and electric modulation of membrane proteins .
Our model rests on the hypothesis that the metabolic energy of the cell is characterized by the coupled oscillations of the energy transducing organelles. Consequently, there are three levels of metabolic organization to be considered: (a) the energy contained in the vibrations at the level of individual enzymes, (b) the coupling of the enzymes within the organelles, and (c) the coupling of the energy transducing organelles within the cell.
Because the energy which drives the process of metabolic regulation depends on coupling at two distinct levels, the enzymatic level and the level of the energy transducing organelles, we can assume that the dimensionality parameter d will depend on physico-chemical properties of the cellular matrix at: (a) the level of the enzymes, (b) the level of the mitochondria and metabolosomes. Hence, we can assume that d, an index of the number of degrees of freedom of the enzymes aligned in the organelles, may not satisfy the condition 1<d<3, as in the Debye model, but may vary in principle between 1 and infinity.
Cycle time and metabolic rate
Quantum Metabolism establishes that the mode of energy transfer may either be quantized or classical, contingent on: (a) the cycle time, τ, of the metabolic processes transforming substrates into products within the network of chemical reactions, and (b) the relaxation time, τ *, which describes the return time of enzyme concentrations to their steady state condition after a random perturbation.
The mean cycle time, τ, is highly dependent on the external resources, their concentration and diversity, whereas the relaxation time, τ *, depends on the density of the metabolic enzymes and the physical properties of cellular medium .
The mean cycle time is the mean turnover time of the enzymes in the reaction process. Based on typical data for ion pumps in biological membranes, the range of values expected for τ is between 10-6s and 10-3s [27, 28]. Recently, the proton-driven ATP synthase rotor has been reported to have Ohmic conductance on the order of 10 fS . Hence, we can estimate for this value the number of protons involved in each cycle to range from 10 at the high frequency limit to 10, 000 at the low-frequency limit. More recent papers provided a novel description of mitochondria as individual oscillators whose dynamics may obey collective, network properties in terms of high-amplitude, self-sustained and synchronous oscillations of bio-energetic parameters under both physiological and patho-physiological conditions . This was demonstrated through the analysis of their power spectrum that exhibits an exponent vastly different from random behavior. Therefore, a proposed description of the metabolic activity involving mitochondrial proteins as coupled quantum oscillators of the Debye type appears to be supported by recent observations.
According to Quantum Metabolism, the quantized or classical modes of energy transfer are determined by the extreme values assumed by the ratio τ/τ*. We have the following two scenarios :
(a) When τ << τ*, the mode of energy transfer is discrete and the dynamic of energy flow is described by quantum laws. In this case d is finite and the scaling exponent, β, satisfies 1/2 < β < 1.
(b) When τ >> τ*, energy transduction is continuous and the dynamic of energy transfer is described by classical laws. Here, the parameter d is infinite and the scaling exponent β = 1.
Proportionality constant and metabolic rate
The proportionality constant, α, is determined by the mechanism - electrical or chemical - which describes the coupling between electron transfer and ATP assembly. The coupling between the electron transport chain and ADP phosphorylation is electrical, in the case of OxPhos, and chemical in the case of glycolytic processes. The metabolic rate in cells utilizing each of these processes will differ, the difference being due to the proportionality constant.
(I) Oxidative phosphorylation. OxPhos is a mode of energy processing whereby the cell carries out phosphorylation of ADP during the oxygen-dependent transmission of electrons down the respiratory chain in the relevant membrane. This reaction produces NADH which then fuels OxPhos to maximize ATP production with minimal production of lactate. The coupling between the electron transport chain and ADP phosphorylation is generated by the flow of protons across the bio-membrane. The proportionality constant α is a function of the bio-energetic parameters, proton conductance C, and the proton-motive potential Δp. These parameters will depend on the phospholipid composition of the membrane .
(IIa) Anaerobic glycolysis (glycolysis which is suppressed by oxygen). Anaerobic glycolysis describes a mode of energy generation when oxygen supply is limited. In this process, cells direct the pyruvate generated by glycolysis away from the mitochondrial OxPhos pathway. This is achieved by generating lactate. This process does not occur within the mitochondria but in the liquid protoplasm of the cell. This generation of lactate allows glycolysis to continue but results in a lowered ATP production rate. The rate-limiting factors for net ATP production are the efficiencies of glucose delivery and lactate removal.
(IIb) Aerobic glycolysis (glycolysis, which is not suppressed by oxygen). Aerobic glycolysis refers to the conversion of glucose to pyruvate and then to lactate regardless of whether oxygen is present or not. This process is also localized in the liquid cytoplasm. This property is shared by normal proliferative tissue.
The proportionality constant in both aerobic and anerobic glycolysis is a function of the activity of the glycolytic enzymes.
We wish to mention that the terms aerobic and anaerobic glycolysis are imprecise, although they have been used extensively in the literature (see for example ). The term aerobic glycolysis is somehow misleading since both types of glycolysis are in fact anaerobic. Warburg used the term aerobic glycolysis to emphasize the important difference from normal glycolysis that it is not suppressed by oxygen. Any type of glycolysis is always substantially the same pathway. It runs much faster in the absence than in the presence of oxygen because of feedback regulation. For example, when oxygen is available, the highly efficient ADP-phosphorylating system in the energy-transducing bio-membrane ensures a high ATP/AMP ratio in the cell, which down-regulates the glycolytic enzymes. The distinction between aerobic and anaerobic glycolysis is therefore only one of net flux rate. Consequently, some researchers consider the distinction redundant. Glycolysis is not oxygen-dependent; but it is always (partly) suppressed by oxygen.
An empirical observation
Quantum Metabolism predicts that the range of values assumed by the scaling exponent, β, β < 1 (an allometric relation), β = 1 (an isometry), is highly dependent on the cycle time, τ, a quantity which varies with the resource flux - an environmentally regulated property. This observation indicates that changes in environmental factors can induce significant changes in the scaling exponent with concomitant changes in metabolic rate. Empirical support for this prediction is given in Nakaya et. al . The experiments in , reported a shift from a scaling law with β = 3/4. to a scaling law described by β = 1, contingent on changes in the environmental conditions. The existence of such environmental switches is particularly pertinent in applications of quantum metabolism to somatic evolution. Such changes in the scaling exponent indicate a mechanism for regulating the metabolic profile of cells by imposing various environmental constraints.
Oxidative phosphorylation and glycolysis: a bio-energetic comparison
Quantum Metabolism predicts that the scaling laws for cells utilizing OxPhos and glycolysis will be described by similar scaling exponents but different proportionality constants. The proportionality constant is determined by the mode of coupling between the electron transport chain and ADP phosphorylation. In OxPhos coupling is achieved by a single common intermediate between the oxidation of a variety of substrates and ATP formation. This intermediate is the trans-membrane proton gradient. For glycolysis, there is a single set of enzymes for every coupled reaction.
The differences in the mode of coupling entail significant differences in the metabolic efficiency of cells utilizing OxPhos and cells using glycolysis. With glucose as substrate, OxPhos generates about 17 times more ATP than glycolysis. Consequently, the metabolic rate of OxPhos (the respiration rate) will also be greater than the metabolic rate of glycolysis (the fermentation rate). The evolutionary history of the different modes of energy production will provide some perspective on the selective advantage which respiration and fermentation conferred as the environmental conditions and resource constraints changed during the history of life on Earth.
The fermentative way of energy generation is now accepted as the primordial mode of energy processing . The invention of photosynthesis changed the atmosphere, since bacteria and plants used light energy to split water into hydrogen and oxygen. After accumulation of oxygen in the atmosphere, OxPhos as a new way of energy generation was established by bacteria. Such free living bacteria have been integrated into eukaryotic cells and represent the symbiosis of two different cell types. By acquiring bacteria capable of OxPhos a new organism emerged having the choice between energy release based on substrate phosphorylation or OxPhos. During evolution of higher vertebrates duplication and modification of genes  led to a fermentative energy generation which is not suppressed by oxygen . Since this type of glycolysis is performed even in the presence of oxygen, Warburg called it aerobic glycolysis. The term aerobic glycolysis is somehow misleading since both types of glycolysis are anaerobic. Warburg's term aerobic glycolysis tries to emphasize the fact that the important difference to normal glycolysis is that it is not suppressed by oxygen.
These observations suggest that in the course of evolution, anaerobic glycolysis arose first. The primitive nature of glycolysis, in contrast to OxPhos, is indicated by the fact that the glycolytic enzymes exist free in solution in the soluble portion of the cytoplasm. This is in sharp contrast to the enzyme systems responsible for respiration and photosynthesis. These enzymes are grouped and arranged in an intracellular structure organized in the mitochondria and chloroplasts .
Modes of energy processing: their incidence
Predominant energy metabolism in different types of tumor cells
Tissue of Tumor
Predominant Energy Metabolism
Glycolytic ATP contribution in selected normal cell types.
Percentage of Glycolytic ATP contribution
Rat coronary endothelial cells
In humans, aerobic glycolysis is extremely important in cells like neurons, retinal cells, stem cells and germ cells. Neurons, for example, depend absolutely on glucose as an energy source and cannot function anaerobically. The utility of anaerobic glycolysis, to a muscle cell for example, when it needs to utilize large amounts of energy in a short period of time, stems from the fact that the rate of ATP production from glycolysis is approximately 100 times faster than from oxidative phosphorylation. During exertion muscle cells do not need to activate anabolic reaction pathways. The requirement is to generate the maximum amount of ATP, for muscle contraction, in the shortest time frame. This is the reason why muscle cells derive almost all of the ATP consumed during exertion from anaerobic glycolysis. As these types of cells do not display the cancer phenotype, it is evident that elevated glycolysis is not a sufficient condition for tumorigenesis.
Elevated glycolysis is also not always observed at all stages of tumorigenesis. The studies reported in de Groof et. al.  based on transformed fibroblasts in mice, indicate the multi-step nature of the progression from normalcy to the malignant state. These studies document a high rate of OxPhos in newly transformed cells with a subsequent high glycolytic rate in the malignant state. The experimental system shows that the proliferation of malignant cells hinges on the evolution of different metabolic profiles as an adaptive response to the changes in the resource conditions during the transition from normal to malignant cells.
We will now show how these evolutionary changes in metabolic profile can be understood in terms of a measure of selective advantage - a prescription derived from Quantum Metabolism and evolutionary dynamics.
Darwinian fitness in somatic evolution
The concept Darwinian fitness describes the capacity of a variant type to increase in frequency in competition with an incumbent population for the available resources. This notion is of fundamental importance in all analytical studies of the Darwinian process at molecular, cellular and organismic levels.
Here, Δr = r* - r, where r* is growth rate of the variant and incumbent type, respectively. These models have provided qualitative insight into several studies of the evolutionary process, and since 1930 have been the dominant approach in evolutionary genetics - in spite of several inconsistencies in their predictive and explanatory power.
Measure of selective advantage: evolutionary entropy
The analytical and conceptual basis for the limitations of the classical models proposed by Fisher, was only recently recognized. The studies reported in  showed that the population growth rate as a measure of fitness is only valid when population size and resources are infinite. Consequently, the Malthusian parameter will only be a meaningful selective index when population sizes are large.
The studies initiated by Demetrius  and later developed as a general model of the evolutionary process showed that the outcome of competition between an incumbent and a variant type is conditional on the magnitude and the variation in resource constraints, and is predicted by the robustness or the demographic stability of the population . Robustness describes the rate at which the population returns to its steady state condition after a random perturbation in the individual birth and death rates. This property can be analytically described in terms of the quantity evolutionary entropy, an information theoretic measure which describes the uncertainty in the state - age, size, metabolic energy - of the immediate ancestor of a randomly chosen newborn. Evolutionary entropy in cellular populations describes the variability in the rate at which individual cells pass through the various stages of the cell cycle. A synchronous population has small entropy, an asynchronous population has large entropy.
Φ < 0 corresponds to limited resource
Φ > 0 corresponds to unlimited resource
γ < 0 corresponds to a variable resource distribution
γ > 0 characterizes a constant resource distribution
The measure of selective advantage given by Eq. (3) is a far-reaching generalization of the measure given by Eq. (2). Indeed, Eq. (3) reduces to Eq. (2) when γ = 0 or when the population size N tends to infinity. The condition γ = 0 corresponds to a complete correlation between the resource variability and the demographic variability.
In view of the characterization of the statistical parameters Φ and γ as resource constraints, the measure of selective advantage given by Eq. (3) can be qualitatively described in terms of the following two tenets.
(Ia) When resources are constant and limited, variants with increased entropy will have a selective advantage and increase in frequency in competition with the resident type.
An increased entropy describes an increased capacity to maintain a constant population size when subject to random perturbations in the individual birth and death rates. This stability condition confers a selective advantage when resources are limited.
(Ib) When resources fluctuate in abundance, variants with decreased entropy will have a selective advantage and increase in frequency in competition with the resident type.
A decreased entropy represents an increased sensitivity of population numbers to random perturbations in the individual birth and death rates. This increased sensitivity - which entails a decreased resilience or stability - now confers a selective advantage when resources fluctuate in abundance, as it enables the population to track the changes in the resource abundance.
The importance of this measure of selective advantage for the study of cellular evolution and the evolution of the cancer phenotype rests on the fact that evolutionary entropy, a population parameter, is analytically related to physiological and morphometric properties of individual cells . Consequently, the selective advantage can also be described in terms of changes in individual properties such as metabolic rate and cell size.
Metabolic rate as fitness
where the taxon-specific constants a2 and b2 are positive.
This result asserts that the outcome of competition between a resident and a variant cell type will be dependent on the resource abundance, as described by Φ, the resource variability, as expressed by γ, the population size, N, and the metabolic rate P.
The significance of metabolic rate in the description of selective advantage, as given by Eq. (7), revolves around the observation that cellular metabolic rate is constrained by two main factors: (a) the rate at which the cells acquire resources, and (b) the rate at which resources acquired are transformed into cellular work.
In this context, when population size, N, is large, the rules describing the outcome of competition between resident and variant types which Eq. (7) prescribes, can be qualitatively annotated in the following terms.
(IIa) When resources are constant and limited, variants with increased metabolic rate have a selective advantage and increase in frequency in competition with the resident type.
An increased metabolic rate entails an increased efficiency in the acquisition of resources. This condition will confer a selective advantage when resources are limited. Hence, when limited resource conditions prevail, types with higher metabolic rates will be favoured.
(IIb) When resources vary in abundance, and the mean abundance is large, resource density is no longer a limiting constraint and types with decreased metabolic rate will have a selective advantage and increase in frequency.
When the mean abundance is large, selective outcome will no longer be contingent on the efficiency in resource acquisition. The dynamics of selection will now be determined by the capacity of cells to transform resources into cellular work. In view of the allometric relations between metabolic rate and cell size, cells with decreased metabolic rates will typically have smaller sizes; consequently their capacity to transform resources into cellular activity will be enhanced. Therefore, under abundant resource conditions cells with decreased metabolic rate will prevail.
These general rules provide a basis for predicting changes in metabolic profile-metabolic rate and metabolic efficiency - as one cell type replaces a related type by natural selection.
Metabolic rate is the rate of ATP production per unit time. Metabolic efficiency is the amount of ATP produced per unit of substrate material. The metabolic rate is the minimum rate at which the cell uses energy to stay alive. All living systems exist in a steady state relatively far from thermodynamic equilibrium and this state is maintained by sustaining non-equilibrium concentration gradients across membranes. This activity requires a constant influx of metabolic energy, and the consequent generation of entropy. The ratio of the energy stored in the various bio-molecules to the free energy released in redox reactions is the efficiency of the metabolic process, which can be as high as about 95% in the case of the Krebs cycle. This high efficiency is characteristic of cyclic processes.
Our analysis rests on certain critical differences in the metabolic rate of cells utilizing fermentation and respiration. First, we observe that there are significant differences in the metabolic efficiency between glycolysis and OxPhos. The metabolism of glucose to lactate generates only 2 ATP (net) per molecule of glucose, whereas OxPhos generates up to 38 ATP molecules (net) . Now, metabolic efficiency and metabolic rate in cells are determined by the same processes and the same reactions. Consequently, metabolic rate under glycolysis will be less than that under OxPhos. These observations, together with the expression for selective advantage given by Eq. (7), imply that the metabolic profile - fermentation or respiration which defines a cell type will be determined by the resource constraints the populations experience during their evolutionary history.
We can, therefore, infer from the principles described in (IIa) and (IIb) the following rules relating resource constraints and the selective advantage of the two classes of metabolic pathways.
(IIIa) When resources are constant and limited, cells utilizing oxidative phosphorylation will have a selective advantage.
(IIIb) When resources vary in abundance and the mean abundance is large, cells utilizing the glycolytic pathway will have a selective advantage.
These principles, which derive from the higher metabolic efficiency and concomitantly the higher metabolic rate of OxPhos, provide an evolutionary explanation of the mode of energy processing observed in microbes and higher organisms. The principles given by (IIIa) and (IIIb) entail that microbes and higher organisms will be characterized by the OxPhos mode of energy processing when evolving under limited resource conditions, and a glycolytic mode when evolving under exponential growth conditions.
Empirical observations on the metabolic phenotypes of microbes and cells from multi-cellular organisms are consistent with these predictions. Unicellular organisms starved of nutrients respire aerobically, as is the case of cells in multi-cellular organisms that are maintained in a non-proliferative state. Unicellular organisms undergoing exponential growth metabolize glucose by fermentation into a small organic molecule such as ethanol. Proliferating cells in multi-cellular organisms exhibit a similar metabolic profile: these cells utilize glycolysis and excrete large amounts of lactate.
The selection principle based on metabolic rate can also be cited to explain the use of aerobic glycolysis in certain classes of cells, namely, neurons, retinal cells, stem cells and germ cells. Aerobic glycolysis confers a selective advantage in cells evolving under conditions where the resources are abundant. Neurons, retinal cells and stem cells require for their metabolic functioning, a large resource energy. The presence of glycolysis as the mode of energy processing in these cells is the selective outcome of the environmental conditions under which these types of cells evolve. These cells do not develop the cancer phenotype since the energy which the cells absorb is constantly being used to maintain the homeostatic integrity of the organism.
The Warburg Effect: Competition between cancer cells and normal cells
Quantum Metabolism predicts that differences in the metabolic rate of cells using OxPhos and cells using glycolysis derive primarily from the differences in the mode of coupling between the energy generated by the flow of electrons in the electron transport chain and ADP phosphorylation.
In OxPhos, the coupling is electrical and the metabolic rate is determined by the flow of protons across the mitochondrial inner membrane. In glycolysis, the coupling is chemical and the metabolic rate is regulated by the activity of the glycolytic enzymes.
Carcinogenesis is a multi-step process that implicates the activation of oncogenes, the inactivation of tumor suppressor genes and the ultimate emergence and proliferation of the cancer phenotype. The transition to the glycolytic phenotype of malignant cells which define the Warburg effect is a consequence of an evolutionary process. Our synthesis of Quantum Metabolism  and evolutionary dynamics  establishes that the index of selective advantage which describes this process is a function of the resource constraints and the metabolic rate of the competing cell types.
Outcomes of competition between resident population using OxPhos and mutant cells using glycolysis.
(I) Limited resources
Ox Phos [(almost always) = (a.s.)]
(II) Abundant resources
(III) Limited Resources
Large population size
Ox Phos (a.s.)
Small population size
Ox Phos (with probability that increases with population size)
(IV) Abundant resources
Large population size
Small population size
Glycolysis (with probability that increases with population size)
Table 3 underscores the stochastic nature of the outcome of competition between the two types of cells and consequently the stochastic nature of the evolution of the glycolytic phenotype. The selective outcome depends not only on the resource constraints but also on the population size of the cells from which the tumor originated. Items (III) and (IV) in Table 3 indicate that when the population size is small, the selective outcome is highly stochastic - a glycolytic phenotype may outcompete the OxPhos pathway, even though the resources are limited; and an OxPhos phenotype may prevail in competition even though the resources are fluctuating. The stochasticity which small population size generates - a property emphasized in the general theoretical study in  - indicates that when resources are limited, tumor cells may even invade a population of normal differentiated cells.
The stochasticity which Table 3 delineates partly explains the non-universality of the Warburg effect and the caution which must be exercised in using the effect as a basis for detecting cancer. Most but not all cancer cells display the glycolytic phenotype.
The Warburg hypothesis asserts that cancer is primarily a disease of metabolic dysregulation. The activation of oncogenes and the inactivation of tumor suppressor genes have engineered alterations in metabolism and a shift from an oxidative mode of energy processing to a glycolytic mode .
This phenomenon suggests that therapeutic strategies based on arresting the transition from a normal differentiated cell defined by oxidative mode of energy processing, to a malignant cell defined by up-regulation of glycolysis, may be effective in complementing traditional methods based on radiation and chemotherapy.
The density and variability of the resources the cells utilize.
The population size of the tumor cells.
The relative metabolic rate of the normal and tumor cells.
Our evolutionary study has discovered two critical features of the competitive dynamics of tumor cells and their normal ancestors.
(1) When resources are constant and limited, the OxPhos phenotype, in view of its higher metabolic rate, will have a selective advantage. Consequently, cells displaying the glycolytic phenotype will have a relatively low probability of becoming established and proliferating. In essence, limited resource condition will inhibit the proliferation of the cancer phenotype.
(2) When resources vary in abundance, the glycolytic phenotype, in view of its lower metabolic rate, will have a selective advantage. Accordingly, cells displaying the glycolytic phenotype will have an enhanced probability of becoming established. Hence resources which show large variation in abundance will facilitate the development of the cancer phenotype, its invasiveness and metastases.
These observations suggest at least four distinct classes of therapeutic strategies to contain the shift of tumors towards the malignant state. These strategies are based on various types of metabolic interventions, some of which have been extensively investigated in the past oblivious, however, to the theoretical rationale which quantum metabolism confers.
Modifying the metabolic rate of cancer cells
Cancer cells exploit the glycolytic pathway as the mode of energy processing. The metabolic rate of cells using this mode of energy processing is determined primarily by the glycolytic enzymes : hexokinase HXK, phosphofructokinase, (PFK), lactate dehydrogenase, (LDH), transketolase (TKTL1). These enzymes are generally up regulated in various cancer cells. This class of enzymes represent potential therapeutic targets as the modification of their activity will have significant effects on the metabolic rate of cancer cells and hence their selective advantage.
where ΔΨ is the membrane potential, F denotes Faraday's constant and ΔpH is the change in pH across the membrane. The values of Δp, ΔΨ and ΔpH depend on the geometry and chemical composition of the energy transducing membrane. Consequently, the metabolic rates can be controlled by not only the membrane composition and its geometry but also the variations in temperature and pH. One of the well-known characteristics of cancer cells is their difference in pH (acidic) compared to normal cells. Interestingly, cancer cells appear to exhibit greater heat sensitivity than normal cells which may also lead to therapeutic applications. Recent biochemical and clinical studies have revealed a profound and selective toxic effect of elevated temperatures on tumor cells .
Modifying the metabolic rate of normal cells
The metabolic rate of normal cells is primarily determined by the proton conductance and the proton potential of the metabolic systems. These bio-energetic parameters are regulated by the phospholipid composition of the mitochondrial membrane. This composition can be modified by exercise and diet . Such modification can increase the metabolic rate of normal cells, thus enhancing its selective advantage, under limited resource conditions, with cells utilizing the glycolytic pathway.
It is also tantalizing to speculate about the high therapeutic value of metabolism-oriented drugs, such as the recently headlined molecule dichloroacetic acid (DCA) which attempts to restore OxPhos and hence mitochondrial function, ultimately restoring apoptosis which has been demonstrated to kill cancer cells in vitro, and shrink tumors in rats . DCA is currently in clinical trials and it is hoped that it may become the first effective drug in dealing with the glycolytic pathway activated in cancers. A recent report  has shown that inhibition of glycolysis severely depletes ATP in cancer cells, especially in clones of cancer cells with mitochondrial respiration defects, and leads to rapid dephosphorylation of the glycolysis-apoptosis integrating molecule BAD, and cell death. Importantly, inhibition of glycolysis effectively kills colon cancer cells and lymphoma cells in a hypoxic environment in which the cancer cells exhibit high glycolytic activity and decreased sensitivity to common anticancer agents.
Changing the cycle time using uncoupling agents and inhibitors
Cycle time can be experimentally manipulated, for instance, using uncoupling agents and inhibitors of electron transport. There are six distinct types of poisons which may affect mitochondrial function : (i) respiratory chain inhibitors (e.g. cyanide, antimycin, rotenone, TTFA, malonate) which block respiration in the presence of either ADP or uncoupling agents; (ii) phosphorylation inhibitors (e.g. oligomycin) which abolish the burst of oxygen consumption; (iii) uncoupling agents (e.g. DNP, CCCP, FCCP, oligomycin) which disrupt the linkage between the respiratory chain and the phosphorylation system; (iv) transport inhibitors (e.g. atractyloside, bongkrekic acid, NEM) which either prevent the export of ATP, or the import of raw materials across the mitochondrial inner membrane; (v) ionophores (e.g. valinomycin, nigericin) which make the inner membrane permeable to compounds which are ordinarily unable to cross it; and (vi) Krebs cycle inhibitors (e.g. arsenite, aminooxyacetate) which block one or more of the TCA cycle enzymes, or an ancillary reaction.
According to Quantum Metabolism, changing the cycle time can induce a change of the scaling exponent β  which, in turn, will affect the metabolic rate of the cell type. This can further alter the selective advantage of tumor cells and thereby impede the transition to a malignant state. However, the toxicity of these compounds makes clinical applications to cancer therapy very challenging, except perhaps for the arsenic compounds. Hence novel uncoupling agents are still much in demand.
Discussion and Conclusions
These differences in the factors that determine metabolic rate in the two modes of energy processing have been exploited to explain the Warburg effect and the selective advantage of the glycolytic phenotype. The understanding of the metabolic basis of carcinogenesis which this evolutionary study generates suggests new methods for regulating cellular metabolic rate and thereby impeding the transition from a benign tumor to a malignant state. In this connection, it is interesting to note that Western diet, which is very calorie-rich, confers evolutionary advantage to cancer cells, hence the increased incidence of cancer in North America and Europe compared to developing countries. Furthermore, increased longevity contributes to this effect as well, since it is well-known that metabolic rate decreases with age of individuals, and in the West the caloric intake does not fall with age at the same rate, thus providing an additional factor for the growth of tumors.
Quantum metabolism, a new bio-energetic theory of metabolic regulation in cells, shows that the proportionality constants in the scaling laws for metabolic rate of cells utilizing OxPhos and glycolysis pathways are contingent on the different modes of coupling between the electron transfer chain and ADP phosphorylation.
(1) In OxPhos, the coupling between the redox reactions (in the mitochondria) and ADP phosphorylation is electrical. Hence the metabolic rate is determined by bio-energetic parameters such as the proton conductance and proton potential of the metabolic network.
(2) In glycolysis, ADP phosphorylation is linked chemically to individual enzyme-catalyzed reactions, so the rate of ATP production is a function of the flux through the pathway as a whole.
Detailed analysis of internal regulation of ATP turnover through both glycolysis and oxidative phosphorylation in terms of specific biochemical mechanisms has been presented in [43, 44] in the context of normal cells. Metabolic interventions which focus on modulating the efficiency of these modes of coupling should provide a general framework for regulating the incidence of these two modes of energy processing - and thereby yield a new paradigm for cancer therapy . A recent article  has advocated a combination therapy utilizing vitamins as a strategy aimed at energy modulation in cancer leading to the prevention of life-limiting cachexia.
J.A.T. acknowledges generous support for this research provided by NSERC, Alberta Cancer Foundation, the Allard Foundation, the Canadian Breast Cancer Foundation and the Alberta Advanced Education and Technology. Insights gained through fruitful discussions with Dr. J.A. Turner are gratefully acknowledged. L.A.D. acknowledges the support of the Max Planck Institute for Molecular Genetics, Berlin.
- Hanahan D, Weinberg R: The hallmarks of cancer. Cell. 2000, 100: 57-70. 10.1016/S0092-8674(00)81683-9.View ArticlePubMedGoogle Scholar
- DePinho RA: The age of cancer. Nature. 2000, 408: 248-254. 10.1038/35041694.View ArticlePubMedGoogle Scholar
- Armitage P, Doll R: The age distribution of cancer and a multi-stage theory of carcinogenesis, International Journal of. Epidemiology. 2004, 33: 1174-1179. 10.1093/ije/dyh216.Google Scholar
- Zetter BR: Angiogenesis and Tumor Metastasis. Ann Rev Med. 1998, 49: 407-424. 10.1146/annurev.med.49.1.407.View ArticlePubMedGoogle Scholar
- Siemann DW, Chaplin DJ, Horsman MR: Vascular-targeting therapies for treatment of malignant disease. Cancer. 2004, 100: 241-249. 10.1002/cncr.20299.View ArticleGoogle Scholar
- Ebos JML, Lee CR, Cruz-Munoz Bjarnason WGA, Christensen JG, Kerbel RS: Accelerated Metastasis after Short-Term Treatment with a Potent Inhibitor of Tumor Angiogenesis. Cancer Cell. 2009, 15: 232-239. 10.1016/j.ccr.2009.01.021.PubMed CentralView ArticlePubMedGoogle Scholar
- Ramanathan A, Wang G, Schreiber SL: Perturbational profiling of a cell-line model of tumorigenesis by using metabolic measurements. Proc Natl Acad Sci USA. 2005, 102: 5992-5997. 10.1073/pnas.0502267102.PubMed CentralView ArticlePubMedGoogle Scholar
- Warburg O, Posener K, Negelein E: Uber den Stoffwechsel der Carcinomzelle. Biochem Z. 1924, 152: 309-344.Google Scholar
- Pelicano H, Xu R-H, Du M, Feng L, Sasaki R, Carew JS, Hu Y, Ramdas L, Hu L, Keating MJ, Zhang W, Plunkett W, Huang P: Mitochondrial respiration defects in cancer cells cause activation of Akt survival pathway through a redox-mediated mechanism. J Cell Biol. 2006, 175: 913-923. 10.1083/jcb.200512100.PubMed CentralView ArticlePubMedGoogle Scholar
- Fearon ER, Vogelstein B: A genetic model for colorectal tumorigenesis. Cell. 1991, 61: 759-767. 10.1016/0092-8674(90)90186-I.View ArticleGoogle Scholar
- Nunney L: Lineage selection and the evolution of multistage carcinogenesis. Proc R Soc Lond. 1999, 226: 493-498. 10.1098/rspb.1999.0664.View ArticleGoogle Scholar
- Gatenby RA, Gillies RJ: Why do cancers have high aerobic glycolysis?. Nature Reviews Cancer. 2004, 4: 891-899. 10.1038/nrc1478.View ArticlePubMedGoogle Scholar
- Demetrius L: Quantum statistics and allometric scaling relations. Physica A. 2003, 322: 477-490. 10.1016/S0378-4371(03)00013-X.View ArticleGoogle Scholar
- Demetrius L, Tuszynski J: Quantum metabolism explains the allometric scaling of metabolic rates. Jour Royal Soc Interface. 2009, 7: 507-514. 10.1098/rsif.2009.0310. March 6, 2010View ArticleGoogle Scholar
- Demetrius L: Directionality principles in thermodynamics and evolution. Proc Natl Acad Sci USA. 1997, 94: 3491-3498. 10.1073/pnas.94.8.3491.PubMed CentralView ArticlePubMedGoogle Scholar
- Rubner M: The problem of life and its relationship to growth and nutrition. 1908, Berlin: OldenbergGoogle Scholar
- Kleiber M: The fire of life. An introduction to animal energetics. 1961, Wiley, New YorkGoogle Scholar
- Hemmingsen AM: Energy metabolism as related to body size and respiratory surfaces, and its evolution. Rep Steno Mem Hosp Copenhagen. 1960, 9: 1-110.Google Scholar
- Lehninger AL: Bioenergetics. 1965, W. A. Benjamin, Inc., New YorkGoogle Scholar
- Harold FM: The way of the cell. 2001, Oxford University Press, OxfordGoogle Scholar
- Aon MA, Cortassa S, O'Rourke B: The fundamental organization of cardiac mitochondria as a network of coupled oscillators. Biophys J. 2006, 91: 4317-4327. 10.1529/biophysj.106.087817.PubMed CentralView ArticlePubMedGoogle Scholar
- Bianchi M: Collective behavior in gene regulation: Metabolic clocks and cross talking. FEBS Journal. 2008, 275: 2356-2363. 10.1111/j.1742-4658.2008.06397.x.View ArticlePubMedGoogle Scholar
- Klevecz RR, Li CM, Marcus I, Frankel PH: Collective behavior in gene regulation: The cell is an oscillator, the cell cycle a developmental process. FEBS Journal. 2008, 275: 2372-2384. 10.1111/j.1742-4658.2008.06399.x.PubMed CentralView ArticlePubMedGoogle Scholar
- Palumbo MC, Farina L, De Santis A, Giuliani A, Colosimo A, Morelli G, Ruberti I: Collective behavior in gene regulation: Post-transcriptional regulation and the temporal compartmentalization of cellular cycles. FEBS Journal. 2008, 275: 2364-2371. 10.1111/j.1742-4658.2008.06398.x.View ArticlePubMedGoogle Scholar
- Westerhoff H, Tsong TY, Chock PB, Chen YD, Astumian RD: How Enzyme Can Capture and Transmit Free Energy from an Oscillating Electric Field. Proc Natl Acad Sci. 1986, 83: 4734-4738. 10.1073/pnas.83.13.4734.PubMed CentralView ArticlePubMedGoogle Scholar
- Tsong TY: Electroconformational Coupling and Membrane Protein Function. Prog Biophys Mole Biol. 1987, 50: 1-45. 10.1016/0079-6107(87)90002-2.View ArticleGoogle Scholar
- Tsong TY: Electrical Modulation of Membrane Proteins: Enforced Conformational Oscillations and Biological Energy and Signal Transductions. Annual Review of Biophysics and Biophysical Chemistry. 1990, 19: 83-106. 10.1146/annurev.bb.19.060190.000503.View ArticlePubMedGoogle Scholar
- Tsong TY, Astumian RD: Electroconformational Coupling: How Membrane-Bound ATPase Transduces Energy from Dynamic Electric Fields. Ann Rev Physiol. 1988, 50: 273-290. 10.1146/annurev.ph.50.030188.001421.View ArticleGoogle Scholar
- Liu DS, Astumian RD, Tsong TY: Activation of Na+ and K+ pumping modes of (Na,K)-ATPase by an oscillating electric field. J Biol Chem. 1990, 265: 7260-7267.PubMedGoogle Scholar
- Feniouk BA, Kozlova MA, Knorre DA, Cherepanov DA, Armen Y, Mulkidjanian AY, Junge W: The Proton-Driven Rotor ATP Synthase: Ohmic Conductance (10 fS), and Absence of Voltage Gating. Biophys J. 2004, 86: 4094-4109. 10.1529/biophysj.103.036962.PubMed CentralView ArticlePubMedGoogle Scholar
- Hochachka PW, Somero GW: Biochemical Adaptation. 2000, Oxford University Press, OxfordGoogle Scholar
- Nakaya F, Saito Y, Motokawa T: Experimental allometry: effect of size manipulation on metabolic rate of colonial ascidians. Proc Roy Soc B: Biol Sci. 2005, 272: 1963-1969. 10.1098/rspb.2005.3143.View ArticleGoogle Scholar
- Coy J, Wittig R: The Role of Glucose Metabolism and Glucose-Associated Signalling in Cancer. Persp Medicinal Chem. 2007, 1: 64-84.Google Scholar
- de Groof AJC, te Lindert MM, van Dommelen MM, Wu M, Willemse M, Smift AL, Winer M, Oerlemans F, Pluk H, Fransen JA, Wieringa B: Increased OXPHOS activity precedes rise in glycolytic rate in H-RasV12/E1A transformed fibroblasts that develop a Warburg phenotype. Mol Cancer. 2009, 8: 54-10.1186/1476-4598-8-54.PubMed CentralView ArticlePubMedGoogle Scholar
- Fisher RA: The Genetical Theory of Natural Selection. 1930, Clarendon Press, OxfordView ArticleGoogle Scholar
- Demetrius L, Gundlach M, Ochs G: Invasion exponents in biological networks. Physica A. 2009, 388: 651-672. 10.1016/j.physa.2008.10.048.View ArticleGoogle Scholar
- Demetrius L, Legendre S, Harremöes P: Evolutionary entropy: A predictor of body mass, metabolic rate and maximal life span. Bulletin of Mathematical Biology. 2009, 9: 800-818. 10.1007/s11538-008-9382-6.View ArticleGoogle Scholar
- Heiden Vander MG, Cantley LC, Thomson CB: Understanding the Warburg effect: the metabolic requirements of cell proliferation. Science. 2009, 324: 1029-1033. 10.1126/science.1160809.View ArticleGoogle Scholar
- Roskoski R: Biochemistry. 1996, Elsevier, AmsterdamGoogle Scholar
- Cavaliere R, Ciocatto EC, Giovanella BC, Heidelberger C, Johnson RO, Margottini M, Mondovi B, Moricca G, Rossi-Fanelli A: Selective heat sensitivity of cancer cells. Biochemical and clinical studies. Cancer. 209: 1351-1381.
- Quiles JL, Huertas JR, Manãas Battino MM, Mataix J: Physical exercise affects the lipid profile of mitochondrial membranes in rats fed with virgin olive oil or sunflower oil. British Journal of Nutrition. 1999, 81: 21-24.PubMedGoogle Scholar
- Michelakis ED, Webster L, Mackey JR: Dichloroacetate (DCA) as a potential metabolic-targeting therapy for cancer. Br J Cancer. 2008, 99: 989-994. 10.1038/sj.bjc.6604554.PubMed CentralView ArticlePubMedGoogle Scholar
- Brand MD: Regulation analysis of energy metabolism. Journal of Experimental Biology. 1997, 200: 193-202.PubMedGoogle Scholar
- Ainscow EK, Brand MD: Internal regulation of ATP turnover, glycolysis and oxidative phosphorylation in rat hepatocytes. Eur J Biochem. 1999, 266: 737-749. 10.1046/j.1432-1327.1999.00856.x.View ArticlePubMedGoogle Scholar
- Moreno-Sanchez R, Rodriguez-Enriquez S, Marin-Hernandez A, Saavedra E: Energy metabolism in tumor cells. FEBS J. 2007, 274: 1393-1418. 10.1111/j.1742-4658.2007.05686.x.View ArticlePubMedGoogle Scholar
- Perumal SS, Shanthi P, Sachdanandam P: Energy-modulating vitamins - a new combinatorial therapy prevents cancer cachexia in rat mammary carcinoma. British Journal of Nutrition. 2005, 93: 901-909. 10.1079/BJN20051439.View ArticlePubMedGoogle 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 cited.