Skip to main content
  • AACR Publications
    • Blood Cancer Discovery
    • Cancer Discovery
    • Cancer Epidemiology, Biomarkers & Prevention
    • Cancer Immunology Research
    • Cancer Prevention Research
    • Cancer Research
    • Clinical Cancer Research
    • Molecular Cancer Research
    • Molecular Cancer Therapeutics

AACR logo

  • Register
  • Log in
  • My Cart
Advertisement

Main menu

  • Home
  • About
    • The Journal
    • AACR Journals
    • Subscriptions
    • Permissions and Reprints
  • Articles
    • OnlineFirst
    • Current Issue
    • Past Issues
    • Meeting Abstracts
    • Collections
      • COVID-19 & Cancer Resource Center
      • Focus on Computer Resources
      • Highly Cited Collection
      • Editors' Picks
      • "Best of" Collection
  • For Authors
    • Information for Authors
    • Author Services
    • Early Career Award
    • Best of: Author Profiles
    • Submit
  • Alerts
    • Table of Contents
    • Editors' Picks
    • OnlineFirst
    • Citations
    • Author/Keyword
    • RSS Feeds
    • My Alert Summary & Preferences
  • News
    • Cancer Discovery News
  • COVID-19
  • Webinars
  • Search More

    Advanced Search

  • AACR Publications
    • Blood Cancer Discovery
    • Cancer Discovery
    • Cancer Epidemiology, Biomarkers & Prevention
    • Cancer Immunology Research
    • Cancer Prevention Research
    • Cancer Research
    • Clinical Cancer Research
    • Molecular Cancer Research
    • Molecular Cancer Therapeutics

User menu

  • Register
  • Log in
  • My Cart

Search

  • Advanced search
Cancer Research
Cancer Research
  • Home
  • About
    • The Journal
    • AACR Journals
    • Subscriptions
    • Permissions and Reprints
  • Articles
    • OnlineFirst
    • Current Issue
    • Past Issues
    • Meeting Abstracts
    • Collections
      • COVID-19 & Cancer Resource Center
      • Focus on Computer Resources
      • Highly Cited Collection
      • Editors' Picks
      • "Best of" Collection
  • For Authors
    • Information for Authors
    • Author Services
    • Early Career Award
    • Best of: Author Profiles
    • Submit
  • Alerts
    • Table of Contents
    • Editors' Picks
    • OnlineFirst
    • Citations
    • Author/Keyword
    • RSS Feeds
    • My Alert Summary & Preferences
  • News
    • Cancer Discovery News
  • COVID-19
  • Webinars
  • Search More

    Advanced Search

Integrated Systems and Technologies

Modeling the Genetic Regulation of Cancer Metabolism: Interplay between Glycolysis and Oxidative Phosphorylation

Linglin Yu, Mingyang Lu, Dongya Jia, Jianpeng Ma, Eshel Ben-Jacob, Herbert Levine, Benny Abraham Kaipparettu and José N. Onuchic
Linglin Yu
1Center for Theoretical Biological Physics, Rice University, Houston, Texas.
2Applied Physics Program, Rice University, Houston, Texas.
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Mingyang Lu
1Center for Theoretical Biological Physics, Rice University, Houston, Texas.
3The Jackson Laboratory, Bar Harbor, Maine.
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • For correspondence: jonuchic@rice.edu kaippare@bcm.edu mingyang.lu@jax.org
Dongya Jia
1Center for Theoretical Biological Physics, Rice University, Houston, Texas.
4Systems, Synthetic and Physical Biology Program, Rice University, Houston, Texas.
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Jianpeng Ma
1Center for Theoretical Biological Physics, Rice University, Houston, Texas.
5Department of Biochemistry and Molecular Biology, Baylor College of Medicine, Houston, Texas.
7Department of Bioengineering, Rice University, Houston, Texas.
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Eshel Ben-Jacob
1Center for Theoretical Biological Physics, Rice University, Houston, Texas.
6School of Physics and Astronomy, Tel-Aviv University, Tel-Aviv, Israel.
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Herbert Levine
1Center for Theoretical Biological Physics, Rice University, Houston, Texas.
7Department of Bioengineering, Rice University, Houston, Texas.
8Department of Biosciences, Rice University, Houston, Texas.
9Department of Physics and Astronomy, Rice University, Houston, Texas.
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Benny Abraham Kaipparettu
10Department of Molecular and Human Genetics, Baylor College of Medicine, Houston, Texas.
11Dan L. Duncan Cancer Center, Baylor College of Medicine, Houston, Texas.
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • For correspondence: jonuchic@rice.edu kaippare@bcm.edu mingyang.lu@jax.org
José N. Onuchic
1Center for Theoretical Biological Physics, Rice University, Houston, Texas.
8Department of Biosciences, Rice University, Houston, Texas.
9Department of Physics and Astronomy, Rice University, Houston, Texas.
12Department of Chemistry, Rice University, Houston, Texas.
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • For correspondence: jonuchic@rice.edu kaippare@bcm.edu mingyang.lu@jax.org
DOI: 10.1158/0008-5472.CAN-16-2074 Published April 2017
  • Article
  • Figures & Data
  • Info & Metrics
  • PDF
Loading

Abstract

Abnormal metabolism is a hallmark of cancer, yet its regulation remains poorly understood. Cancer cells were considered to utilize primarily glycolysis for ATP production, referred to as the Warburg effect. However, recent evidence suggests that oxidative phosphorylation (OXPHOS) plays a crucial role during cancer progression. Here we utilized a systems biology approach to decipher the regulatory principle of glycolysis and OXPHOS. Integrating information from literature, we constructed a regulatory network of genes and metabolites, from which we extracted a core circuit containing HIF-1, AMPK, and ROS. Our circuit analysis showed that while normal cells have an oxidative state and a glycolytic state, cancer cells can access a hybrid state with both metabolic modes coexisting. This was due to higher ROS production and/or oncogene activation, such as RAS, MYC, and c-SRC. Guided by the model, we developed two signatures consisting of AMPK and HIF-1 downstream genes, respectively, to quantify the activity of glycolysis and OXPHOS. By applying the AMPK and HIF-1 signatures to The Cancer Genome Atlas patient transcriptomics data of multiple cancer types and single-cell RNA-seq data of lung adenocarcinoma, we confirmed an anticorrelation between AMPK and HIF-1 activities and the association of metabolic states with oncogenes. We propose that the hybrid phenotype contributes to metabolic plasticity, allowing cancer cells to adapt to various microenvironments. Using model simulations, our theoretical framework of metabolism can serve as a platform to decode cancer metabolic plasticity and design cancer therapies targeting metabolism. Cancer Res; 77(7); 1564–74. ©2017 AACR.

Major Findings

We developed a theoretical framework for modeling gene regulation of cancer metabolism. We found that, in addition to glycolytic and oxidative metabolism in normal cells, cancer cells have a new hybrid phenotype in which both metabolic modes coexist. Cells in the hybrid state have enhanced metabolic plasticity to facilitate tumorigenesis and metastasis. Guided by the modeling, we developed the AMPK and HIF-1 signatures to quantify the activity of OXPHOS and glycolysis, respectively. We propose a new cancer therapeutic strategy by targeting the hybrid state.

Introduction

Cells can utilize multiple metabolic pathways for energy production and biosynthesis depending on the requirements for cellular function and the availability of metabolites. In the presence of glucose, cells typically uptake glucose and convert it to pyruvate inside the cytosol by glycolysis. Under normoxic conditions, pyruvate is further transported into mitochondria where it undergoes oxidative phosphorylation (OXPHOS) to produce ATP through the tricarboxylic acid (TCA) cycle and the electron transport chain (ETC). On the other hand, under hypoxia, cells utilize anaerobic glycolysis instead, which converts pyruvate into lactate, and produces ATP in a much faster but less efficient way. If fatty acid or glutamine is available, cells under normoxia can also undergo fatty acid oxidation (also called β-oxidation) or glutamine oxidation.

Quick Guide to Equations and Assumptions

In this study, we constructed a regulatory network of cell metabolism that consists of regulatory proteins and metabolites. To understand the behavior of the network, we coarse-grained the whole network into a minimalist circuit of three components —AMPK, HIF-1, and ROS. Here we assumed that the essential feature of the network could be captured by this circuit, while the effects of many other relevant genes can be regarded as inputs to this circuit. The dynamic behaviors of the regulatory circuit AMPK:HIF-1:ROS (Fig. 2) were analyzed by solving nonlinear differential rate equations. Typically, the deterministic rate equation for a protein or metabolite has a generic form: Embedded Image Here, Embedded Image represents the level of the protein or the metabolite. Embedded Image and Embedded Image represent the basal production and degradation rates of Embedded Image, respectively. Embedded Image and Embedded Image represent the regulation on the production and degradation of Embedded Image, respectively.

Overall, we assumed that the regulation of Embedded Image by a component Embedded Image can be described by a nonlinear function, namely the shifted Hill function (1, 2), Embedded Image where Embedded Image is the threshold level, Embedded Image is the Hill coefficient, and Embedded Image is the fold change.

Embedded Image is the inhibitory Hill function, defined as Embedded Image Embedded Image is the excitatory Hill function, defined as Embedded Image Embedded Image represents an inhibitory regulation and Embedded Image represents an excitatory regulation.

When the production of Embedded Image is regulated by Embedded Image. When the degradation of Embedded Image is regulated by Embedded Image.

When the production (or degradation) of Embedded Image is regulated by two components Embedded Image and Embedded Image simultaneously, Embedded Image can be expressed as: Embedded Image Details of the functional form of Embedded Imagecan be found in Supplementary Information Section S3.

In our modeling, AMPK and HIF-1 competitively regulate the production of ROS (both mitochondria ROS and cytosol ROS); ROS and oxygen competitively regulate the degradation of HIF-1. All the other combined regulatory processes are assumed to be independent. The parameters were estimated per experimental evidence, as shown in Supplementary Tables S3–S6. Furthermore, the model allowed evaluation of the effects of oncogenes and metabolic therapies. Details can be found in Supplementary Information Section S3–S6, S9.

Abnormal metabolism is a hallmark of cancer (3–5). Unlike normal cells, cancer cells largely depend on glycolysis to produce energy even in the presence of oxygen, which is referred to as the Warburg effect (6, 7) or aerobic glycolysis. Aerobic glycolysis is an aggressive metabolic phenotype in that it has the advantage to produce ATP at a high rate and prepare biomass for amino acids and fatty acid synthesis that are required for rapid cell proliferation (7). Although aerobic glycolysis is generally regarded as a dominant metabolism in cancer cells (8), recent evidence suggests that mitochondrial OXPHOS is also utilized by cancer cells (9–12). Indeed, the influence of mitochondria on cancer cells is well documented (13–19), suggesting that mitochondria are actively functioning in cancer cells. More importantly, it has been revealed that OXPHOS contributes to cancer metastasis (11–15, 19–21). Yet, it is still unclear how these metabolic modes are regulated in cancer cells, particularly the advantages of having each one of them.

To shed light on these puzzles, we utilize a computational systems biology approach to study the genetic regulation of glycolysis and OXPHOS (including glucose oxidation, fatty acid oxidation, and glutamine oxidation). Certainly, cancer metabolism has been studied by several computational modeling approaches. However, each of these earlier works mostly focused on only a part of the interplay between glycolysis and OXPHOS. For example, there are studies of the Warburg effect, including modeling of the regulation by reactive oxygen species (ROS) of hypoxia-inducible factor 1 (HIF-1) in response to hypoxia (22), and modeling of the genome-scale metabolic network based on metabolic flux balance analysis (23). Some other studies focused on mitochondrial functions, for example, in the study of the cooperative effects of TrxR and Prx in defending oxidative stresses (24). The evolutionary advantage of metabolic heterogeneity of cancer has also been investigated (25) with emphasis on glycolysis and acidity. However, we still lack a comprehensive picture of how different metabolic pathways contribute to the survival and oncogenic potential of cancer cells. Therefore, there is an urgent need to explore in detail how genes and metabolites regulate both aerobic glycolysis and OXPHOS.

In this study, we establish a theoretical framework for modeling genetic regulation of the interplay between glycolysis and OXPHOS in cancer metabolism. On the basis of an extensive literature survey, we construct a metabolic network model for both glycolysis and OXPHOS. The network includes both regulatory genes and metabolites that are involved in these metabolic pathways. To capture the two regimes of cancer metabolism, we further coarse-grained the network into a core regulatory circuit, composed of HIF-1, 5′ AMP-activated protein kinase (AMPK) and ROS. The circuit contains HIF-1 and AMPK, as they are the master regulators of glycolysis and OXPHOS, respectively (7); and ROS is included because both cytosolic and mitochondrial ROS (mtROS) mediate the interplay between the regulatory genes and metabolites (20, 26). This reduced representation provides a clear picture of the interplay between glycolysis and OXPHOS. Similar HIF-1:AMPK:ROS representation has also been proposed in the study of longevity in Caenorhabditis elegans (C. elegans; ref. 26), suggesting that the circuit is evolutionarily robust. The core regulatory circuit is modulated by several oncogenic signaling pathways, including MYC, RAS, and c-SRC (27–29). Many metabolic drugs (30) directly target the genes that are associated with specific metabolic reactions, such as GLUT1, HK and FASN, or regulatory genes in the current study, that is, AMPK and HIF-1. Therefore, the regulatory circuit allows us to model the effects of the oncogenes and metabolic drugs on cellular metabolism. As we show later, this model can explain important experimental observations, and is sufficient to capture the main differences between cancer and normal cells.

Computational modeling of our reduced circuit shows that in general cancer cells can have three stable states—a Warburg state (W: high HIF-1, low pAMPK), an oxidative state (O: low HIF-1, high pAMPK), and a hybrid state (W/O: high HIF-1, high pAMPK). Here, the AMPK activity is represented by the level of phosphorylated AMPK (pAMPK) at threonine-172 of the α subunit. Compared with cancer cells, normal cells usually lack the hybrid phenotype due to their lower mtROS production and higher rate of HIF-1 degradation. The activation of oncogenic pathways is an alternative way to induce cells to lie in the hybrid state. The discovery of the hybrid metabolic state for cancer cells from the network analysis has been supported by recent experimental evidence (summarized in Supplementary Table S1; refs. 10, 11, 16, 19, 20). The association of cancer metabolism and oncogenic pathways are also supported by extensive The Cancer Genome Atlas (TCGA) and single-cell transcriptomics data analysis of multiple cancer types under this theoretical framework. Particularly, we develop two signatures to quantify AMPK and HIF-1 activities and find that AMPK activity anticorrelates with HIF-1 activity in several tumors but not in the others. We propose that the W/O hybrid state enhances metabolic plasticity, therefore allowing the cancer cells to adapt to the tumor microenvironment and promote cancer progression and metastasis. Guided by the simulations of targeting AMPK or HIF-1, we make predictions on the effectiveness of various therapeutic strategies, such as hyperbaric oxygen, metformin, and 3-bromopyrvate (3BP), in reducing the metabolic plasticity of cancer cells. All told, our model provides not only a theoretical framework for charactering the regulatory mechanism of glycolysis and OXPHOS and its connection to tumorigenesis, but also new insights into the optimal therapeutic strategies that specifically target cancer metabolism.

Results

The regulatory network of glycolysis and OXPHOS

After an extensive literature analysis, we constructed a comprehensive network (Fig. 1) featuring the regulation of oxidative respiration and glycolysis by both genes and metabolites. The network contains the following four types of regulatory interactions. First, the three major metabolic pathways, aerobic glycolysis, glucose oxidation and fatty acid oxidation (green rectangles), directly inhibit each other because they compete for shared metabolites (see Supplementary Information Section S1 for details). Second, the activity of these metabolic pathways is directly regulated by specific genes (red ovals). Most prominently, the fatty acid oxidation and glucose oxidation are regulated by an energy sensor gene, AMPK and the glycolytic pathway is regulated by the hypoxia-inducible factor, HIF-1. Third, some metabolites produced by the metabolic pathways (yellow rectangles) in turn regulate the activities of some regulatory genes. HIF-1 is stabilized by both NOX-derived ROS (noxROS) from cytosol and OXPHOS-derived ROS from mitochondria (mtROS). Meanwhile, ROS induces phosphorylation and activation of AMPK, yet excessive production of ATP by either metabolic pathway could block the AMPK activity (see Supplementary Table S1 for details). Fourth, the regulatory genes are also coupled to several oncogenic pathways (blue ovals), including RAS, MYC, and c-SRC (Supplementary Table S2).

Figure 1.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 1.

Schematic illustration of the regulatory network of metabolism. The network includes both regulatory genes (ellipses) and metabolites (rectangles). The arrows represent positive regulation and the bars represent negative regulation. AMPK and HIF-1 play a central role in regulating different metabolic pathways, while the metabolic pathways regulate activities of AMPK and HIF-1 partially through ROS. The oncogenic pathways directly regulate the activity of AMPK, HIF-1, and noxROS.

The core AMPK:HIF-1:ROS regulatory circuit

To capture the principles of how genes and metabolites modulate metabolism, we coarse-grained the extensive metabolic network (Fig. 1) to a minimalist regulatory circuit consisting of AMPK, HIF-1, and ROS (including noxROS and mtROS; Fig. 2). These components were chosen, because they play critical roles in regulating the decision-making of both glycolysis and OXPHOS. Moreover, the core regulatory circuit captures the main features of the more complete network. As we show later, the reduced circuit is sufficient to explain important experimental observations of glycolysis and OXPHOS. It is worth noting that the core circuit is directly coupled with both the oncogenic pathways and metabolic pathways. Its dynamic behavior can shed insight into the interplay of the various metabolic modes in cancer.

Figure 2.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 2.

The AMPK:HIF-1:ROS regulatory circuit. AMPK and HIF-1 are the master regulators of OXPHOS and glycolysis, respectively. ROS represents both mtROS and noxROS. RAS, MYC, and c-SRC modulates the balance of glycolysis and OXPHOS (details in Fig. 4).

As shown in Fig. 2, the regulatory links among these four components, AMPK, HIF-1, noxROS, and mtROS, are derived from either direct or indirect regulatory interactions (see Supplementary Information Section S2 for details). Various oncogenic pathways are modeled as input signals to the circuit, and the activities of OXPHOS and glycolysis are modeled as the readout of the circuit. The detailed experimental evidence for the regulatory interactions of the full network and the associated reduced circuit are listed in the Supplementary Table S2. The corresponding chemical rate equations are provided in Supplementary Information Section S3. The parameters and their experimental justifications are provided in Supplementary Tables S3–S6.

The metabolic circuit allows multiple cell phenotypes

To identify the possible cell states that are allowed by the core circuit, we performed an analysis of the circuit equations with two sets of parameter that correspond to normal and cancer cells, respectively. Here we illustrate the results by considering the differences between cancer cells and normal cells in two aspects. First, cancer cells have higher mtROS production due to reprogrammed mitochondria (20). Second, cancer cells have more stable HIF-1 because less oxygen is available to each cell due to abnormally rapid proliferation. The effects of the oncogenic pathways will be evaluated in the next section. Therefore, as compared with normal cells, cancer cells produce more mtROS in response to the AMPK activation (represented by Embedded Image, maximum fold change of mtROS by the AMPK activation), and have a lower HIF-1 degradation rate (denoted by Embedded Image). As mentioned above, the activity of AMPK was quantified by the concentration of its phosphorylated form, that is, pAMPK.

Our simulations show that normal cells have two stable steady states, which correspond to the high HIF-1, low pAMPK state and the low HIF-1, high pAMPK state (Fig. 3A), which are associated with the Warburg effect (W) and oxidative respiration (O) respectively (the “W” state and “O” state in Fig. 3A). This result is consistent with the fact that while cells usually use glucose oxidation to produce energy, they switch to glycolysis during anaerobic exercise. Next, we performed an analogous analysis for the cancer cells, as reflected by larger Embedded Image and lower Embedded Image. Interestingly, we found that cancer cells have a new hybrid state – (high pAMPK, high HIF-1; “W/O” state in Fig. 3B) in addition to the “W” and “O” states. The new hybrid state is found to be a robust feature of the cancer cells by parameter sensitivity analysis (Supplementary Fig. S1). In the hybrid state, the levels of pAMPK and HIF-1 are both high because the cancer cells maintain high levels of ROS that both stabilizes HIF-1 and activates AMPK. The existence of the hybrid metabolic phenotype means that cancer cells have the flexibility to simultaneously utilize both glycolysis and mitochondrial OXPHOS.

Figure 3.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 3.

The nullclines and steady states in the phase space of AMPK and HIF-1. The red line represents the nullcline of Embedded Image, and the blue line represents the nullcline of Embedded Image (see Supplementary Information Eq. S1). The green solid dots denote stable steady states and the green hallow dots denote unstable steady state. Each stable state is associated with a metabolic phenotype. For normal cells (A), the circuit allows a Warburg state, denoted by “W,” and an OXPHOS state, denoted by “O.” For cancer cells (B), the circuit allows an additional hybrid state, denoted as “W/O.”

The role of mtROS production and ontogenetic pathways in cancer metabolism

We further evaluated the effects of mtROS production and several oncogenes on the dynamic behavior of the AMPK:HIF-1:ROS circuit. In Fig. 4A, we show that, when the rest of the parameters remain unchanged, the increase of Embedded Image shifts the cells from bistability (with two stable steady states) to tristability (with three stable steady states) and gives rise to the hybrid metabolic state – “W/O,” in addition to “W” and “O” phenotypes. This is because, once again, higher production of mtROS by mitochondria stimulates both AMPK activation and HIF-1 stabilization, which is consistent with the experimental results that partial inhibition of ETC leads to the coexistence of glycolysis and OXPHOS by producing excessive mtROS (20). In addition to ROS production, stabilization of HIF-1, as reflected by lower kh, can also give rise to the hybrid W/O state (Supplementary Fig. S2).

Figure 4.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 4.

Bifurcation diagrams of the levels of phosphorylated AMPK (pAMPK) and HIF-1 driven by mtROS production (A), MYC (B), c-SRC (C), and RAS (D). Embedded Image represents the fold change of mtROS by AMPK activation. The blue, green, and magenta solid lines represent the phenotypes O, W/O, and W, respectively. The red dashed line represents the unstable states. Different background colors represent different phases.

Similarly, we simulated the response of the core circuit to the changes in the activity of the following oncogenic pathways, that is, MYC (Fig. 4B), c-SRC (Fig. 4C), and RAS (Fig. 4D; See Supplementary Information Section S4–S6 for the modeling details). Interestingly, in all these cases, we found that moderate activation of RAS, MYC, or c-SRC induces the circuit from a two-state system to a three-state system containing the hybrid state “W/O.” Strong activation of RAS can drive the circuit to either the “W” or the “W/O” phenotype, while strong activation of c-SRC can drive the circuit to either the “O” or the “W/O” phenotype. However, strong activation of MYC can drive the circuit to the “W” state. The results are consistent with the experimental evidences on the role of some oncogenic pathways in modulating both glycolysis and OXPHOS. Previously, we show that cells with high c-SRC activity have substantial levels of glycolysis and fatty acid oxidation (31). In addition, MYC overexpression has been shown to promote fatty acid oxidation in triple-negative (TN) breast cancer (32). Moreover, RAS can simultaneously activate glycolysis and OXPHOS (33). It is likely that cancer cells can be in the hybrid phenotype when the effects of the excessive production of ROS or the activation of the oncogenes outweigh the mutual inhibitions between AMPK and HIF-1.

Quantification of metabolic state by the activities of AMPK and HIF-1

To validate the metabolic circuit model, we performed an extensive data analysis on the TCGA data from eight cancer types. Here we quantified the activities of AMPK and HIF-1 by evaluating the expression of downstream targets of both AMPK and HIF-1 (a total of 33 genes for AMPK, and 23 genes for HIF-1, see Supplementary Information Section S7 for details). For each cancer type, we performed principal component analysis (PCA) on the RNA-seq data independently for either the AMPK downstream genes or the HIF-1 downstream genes, from which we assigned the first principal component as the axis to quantify the activities of AMPK or HIF-1. The AMPK or the HIF-1 axes (signatures) derived from different cancer types are surprisingly similar (Supplementary Table S7 and S8; Supplementary Fig. S3), indicating the consistent regulatory functions of AMPK and HIF-1 across these tissues and cancer types. Most of the genes have positive contributions to the respective principal component (Supplementary Fig. S4), indicating that they are positively regulated by AMPK or HIF-1. For a few genes that have substantial negative contributions, we identified experimental evidence for the negative regulation (Supplementary Table S9).

Strikingly, even though the AMPK and HIF-1 gene sets are independent and the corresponding principal axes were obtained independently, we observed strong anticorrelations between the HIF-1 and AMPK activities for several major cancer types—liver hepatocellular carcinoma (HCC), lung adenocarcinoma (LUAD), breast invasive carcinoma (Fig. 5A–C). We obtained similar results for three other types of cancer—stomach adenocarcinoma, acute myeloid leukemia, and pancreatic adenocarcinoma (PAC; Supplementary Fig. S5A–S5C), but not in clear cell renal cell carcinoma (CCRCC), prostate adenocarcinoma, and colorectal adenocarcinoma (Supplementary Fig. S5D–S5F). From the data analysis, we observed a spectrum of cases with “W,” “W/O” and “O” phenotypes. Notably, we did not observe two clearly separated clusters representing glycolysis and OXPHOS respectively; therefore, the AMPK and HIF-1 activities are strongly anticorrelated but not mutually exclusive in cancer. However, the TCGA data analysis was not sufficient to explain the existence of the hybrid “W/O” state due to the mixture of normal and cancer cells in the tissues and possible heterogeneity of metabolic states in the cell populations. Thus, we looked into single-cell RNA-seq data of LUAD (34), which shows a similar anticorrelation between the AMPK and HIF-1 activities and mainly a bimodal distribution with substantial amount of hybrid cells (Fig. 5D; see Supplementary Figs. S8–S11 for more details on the single-cell RNA-seq data analysis). This finding suggests that the hybrid metabolic phenotype can also be observed for individual cancer cell.

Figure 5.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 5.

Evaluation of the AMPK and HIF-1 activities using TCGA patient data and single-cell RNA-seq data of lung adenocarcinoma in HCC (n = 373; A), LUAD (n = 517; B), and breast-invasive carcinoma (n = 971; C). Left, each point represents the AMPK and HIF-1 activities of one sample. For each dataset, the standard k-mean analysis was applied to group the cases into the W, W/O, and O states. Middle, the fraction of the top 20% samples with high oncogene activities in each metabolic group (W, W/O, and O). χ2 test was used to test the significance (See Supplementary Information Section S10). *, P < 0.05; **, P < 0.01; ***, P < 0.001. Right, Kaplan–Meier (KM) overall survival curves of patients that were stratified by their metabolic states. D, Evaluating the AMPK and HIF-1 activities of single LUAD cells (n = 77). Left, each point represents the AMPK and HIF-1 activities of a single cell. Histogram represents the distribution of single cells projected to the first principal component of their AMPK and HIF-1 signatures. Three Gaussians are required to fit to the histogram (Supplementary Fig. S11). Right, box plot for HIF-1 activity. Evaluation of the AMPK and HIF-1 activities of luminal A (n = 531), luminal B (n = 135), HER2+ (n = 36), and TN (n = 141) breast carcinoma (E), and HCC (n = 373), CCRCC (n = 534), and PAC (n = 179) from TCGA (F; see Supplementary Information Section S7 for details of methods).

Next, we applied the AMPK and HIF-1 signatures to evaluate the metabolic activities of different subtypes of breast carcinoma and different types of tumors across TCGA. HER2+ and TN breast tumors clearly show higher HIF-1 activity and lower AMPK activity than luminal tumors (Fig. 5E), consistent with the experimental observation that TN cells are more dependent on glycolysis than luminal cells (35). PAC and CCRCC show higher HIF-1 activity and lower AMPK activity compared with HCC (Fig. 5F; Supplementary Fig. S6), consistent with experimental observations (see Supplementary Fig. S6 for details).

For each of the three cancer types shown in Fig. 5A–C, we further evaluated the activities of the oncogenic pathways in samples with various metabolic phenotypes. First, each cancer sample was assigned a metabolic state (“W,” “W/O” and “O”) by k-mean clustering analysis of the whole dataset (Fig. 5, left). The activities of the oncogenic pathways (MYC, c-SRC, and AKT) were quantified by the oncogenic scores derived experimentally in previous studies [see Supplementary Information Section S8 for details regarding the calculation of the oncogene score (36)]. The classification and scoring methods allow us to evaluate the enrichment of samples with high oncogenic activity for each metabolic phenotype (Fig. 5A–C). Interestingly, high MYC activity was observed to be enriched in the “W” phenotype, whereas high c-SRC and AKT (representing RAS) activities are enriched in both the “W” and the “W/O” phenotypes. Low oncogenic activity is associated with the “O” phenotype and low RS group (better prognosis; ref. 34) of single LUAD cells (Supplementary Figs. S7 and S8). These results are consistent with the modeling results and experimental evidence discussed in the previous section. Moreover, we evaluated the 5-year overall survival rates for different metabolism phenotypes and found that patient samples in “O” phenotype are significantly correlated with better survival results compared with patients' samples in either “W” or “W/O” (Fig. 5A–C). Consistently, single LUAD cells with higher risk scores (worse prognosis; ref. 34) maintain higher HIF-1 activity (Fig. 5D). Thus, we believe that our HIF-1 and AMPK signatures provide an easy and reliable scoring method to quantify the metabolic states of glycolysis and OXPHOS based on RNA-seq data from both bulk tumors and single cells.

The advantages of the hybrid metabolic phenotype and experimental support

We argue that cancer cells in the “W/O” hybrid state can have a significant advantage in supporting their survival, proliferation, and metastasis. First, cells in the hybrid state have higher metabolic plasticity to better adapt to various microenvironments, such as hypoxia and acidic conditions, because these cells are flexible in their use of available nutrients, such as glucose and fatty acids, to produce energy. Second, the hybrid state allows cancer cells to efficiently produce energy by both OXPHOS and glycolysis while at the same time using lactate and pyruvate, the byproducts of glycolysis, to generate biomass for cell proliferation (7). Third, cells in the hybrid state can modulate ROS at a moderate level, so that cells can take advantage of ROS signaling (20, 37) to promote metastasis meanwhile avoid excessive DNA damage (38). Here, ROS scavenging is achieved by pyruvate and NADPH, the byproducts of glycolysis and fatty acid oxidation, respectively. Fourth, the hybrid state might be specifically associated with metastasis, as suggested by some experimental evidence listed below. In summary, cancer cells at the hybrid metabolic phenotype have advantage in multiple aspects over cells with either the “W” or the “O” phenotypes.

Our metabolic circuit model is consistent with many recent observations that, although glycolysis is a typical feature of cancer cells, mitochondrial respiration plays a crucial role in tumor invasion and metastasis (9, 16–21, 39). It has been observed that some aggressive tumor cells, such as SiHa and HeLa, have not only oxidative respiration but also glycolysis due to stronger HIF-1 activation from lactate (21). Moreover, the inhibition of the mitochondrial respiratory chain contributes to the reduction of multidrug resistance of slow cycling melanoma cells (14). It has also been shown that the surviving pancreatic cancer cells after doxycycline withdrawal depend on OXPHOS and are highly sensitive to OXPHOS inhibitors (13). Also, the super-metastatic tumor cells obtained by experimental selection in vitro (SiHa-F3 cells) and in vivo (B16F10 and B16-M1 to M5 tumor cells) have increased OXPHOS with higher ROS production (20). TGFβ1 treatment of non–small lung carcinoma A459 cells induces metastasis, and they were found to have decreased fatty acid synthesis and increased oxygen consumption (11). The inhibition of OXPHOS by graphene for multiple cancer cells, including breast cancer cells and hepatocellular carcinoma cells, can effectively inhibit tumor migration and invasion (39). In the context of 4T1 mammary epithelial cancer cells, PGC-1α promotes metastasis through activating mitochondrial biogenesis and OXPHOS (10). The breast cancer cells can also shift from glycolysis to mitochondrial OXPHOS after radiation exposure to generate more ATP for survival (16). A detailed list of these experimental findings is presented in Supplementary Table S1. These data taken together indicates a critical role of OXPHOS in tumorigenesis, and supports our prediction of the hybrid “W/O” state. Our model and its experimental support suggest that the hybrid metabolic phenotype could be a good target for metabolic drugs.

Modeling therapeutic strategies targeting cancer metabolism

It has been shown that targeting cellular metabolism is a promising strategy for fighting against cancer (details in Supplementary Table S10; ref. 30). Certain metabolic drugs have been shown to be effective in treating cancers in some cases. These drugs, for example, S-trans, trans-farnesylthiosalicylic acid (FTS; ref. 40), 2-deoxy-d-glucose (2DG; ref. 41), 3-bromopyruvate (3BP; ref. 42), metformin (43), and AICAR (44), typically target glycolysis or mitochondrial OXPHOS. It has also been suggested that combinations of different drugs could prove more effective (40, 45, 46). For example, administrating metformin and 2DG together can induce massive ATP depletion in cancer cells and further trigger cellular processes to induce cell death, such as autophagy or p53-dependent apoptosis (45). Often, metabolic drugs work well in some cases but not in others, and the underlying mechanisms are not completely understood. Here we propose that these therapeutic strategies might be effective partly because they target the hybrid metabolic phenotype of cancer cells. An effective treatment outcome could be achieved by shifting the metabolic phenotype from the hybrid state to the other allowed states of the cancer cells. By doing so, the drugs might sensitize the cancer cells, therefore, other therapies, such as chemotherapy, can be more effective in killing the cancer cells.

We evaluated the dynamic response of the AMPK:HIF-1:ROS circuit to various treatment strategies, each of which exerts its own regulatory mechanism on the circuit (See Supplementary Information Section S9 for details). In particular, hyperbaric oxygen therapy can effectively reduce hypoxia and accelerate degradation of HIF-1; 3BP, similar to 2DG, targets the glycolytic enzymes such as glucose transporters (GLUT) and hexokinase (HK), thus effectively reducing cellular glycolysis; metformin activates AMPK, inhibits ETC Complex-1, and inhibits mTOR in an AMPK-independent manner, which further inhibits HIF-1. The effective circuit diagrams for each case are shown in Fig. 6, and details of the modeling procedures are shown in Supplementary Information Section S5 with the relevant parameters given in Supplementary Table S11. The model allows us to calculate the steady states and simulate the time course of the levels of pAMPK, HIF-1, and ROS. We are especially interested in how different treatments drive phenotypic transitions among different metabolic states.

Figure 6.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 6.

The bifurcation diagrams of the levels of pAMPK in response to the oxygen level in the hyperbaric oxygen therapy (A), 3BP level (B), metformin level (C), and combined 3BP and metformin therapy (D). In A, the normal condition of oxygen level is 5%. In B–D, the drug (1/2/3) represents a hypothetic signal with opposite effect of the corresponding drug(s). In A–D, blue, green, and magenta solid lines represent phenotype O, W/O, and W, respectively. The red dashed lines represent the unstable steady states. The dashed lines with arrows indicate transitions between two phenotypes. Different background colors represent different phases.

A comparison of possible metabolic therapies

Starting from the model where cancer cells have the “W,” the “W/O,” and the “O” states, we evaluated the effects of various treatment strategies on allowing the cell to escape from the tristability of the metabolic circuit. Each of the panel of Fig. 6A–D shows the bifurcation diagram of the levels of pAMPK and HIF-1 with respect to the level of the drugs, for each treatment. The results show that these treatment strategies can shift the cancer metabolism from the tristable phase to the monostable phase (“O”) by increasing the dose level.

From the bifurcation analysis, metformin is less effective compared with the hyperbaric oxygen therapy and 3BP. The inefficiency of metformin is caused by the decrease of mitochondrial potential and increase of mtROS production by the drug (46). Various treatment strategies were further evaluated by simulating the time course of the levels of pAMPK, HIF-1, and ROS (Supplementary Figs. S12–S15). We assume that an effective therapy needs to drive cancer cells away from the hybrid states. We found that the administration of metformin alone is less effective in avoiding the hybrid state, while the combined therapies are more effective. The prediction is consistent with the finding that, in human gastric and esophageal cell lines, the administration of both metformin and 2DG (whose effect is similar to 3BP) is more effective than the administration of each drug alone (46).

Discussion

In this study, we established a theoretical framework for modeling genetic regulation of cancer metabolism. By integrating existent data, a network was constructed that features the regulation of mitochondrial OXPHOS and glycolysis by both regulatory protein and metabolites. We further coarse-grained the network into a core regulatory circuit that is composed of HIF-1, AMPK, and ROS. Although we made several simplifications, the reduced circuit was still sufficient to capture regulation of both glycolysis and OXPHOS, as well as the major differences in metabolism between normal and cancer cells. As the regulatory links in the core circuit are supported by multiple experiments (Supplementary Table S1), the circuit model is expected to be robust for studying the behavior of glycolysis and OXPHOS.

Normal cells typically use either OXPHOS or glycolysis at a fixed time depending on the availability of energy sources, because of tight regulation of the metabolic circuit and the competition between different metabolic modes. Cancer cells, however, kidnap the same gene regulatory circuit of metabolism for their own advantage. Because of either high energetic and oxidative stresses or the activation of specific oncogenes, cancer cells can be in a hybrid metabolic phenotype utilizing both glycolysis and OXPHOS. As cells in the hybrid “W/O” state have increased plasticity, they could have an advantage in survival over cells with the other phenotypes. The modeling results are supported by the TCGA patient data from multiple cancer types. The model also explains the phenomenon of oxygen shock (47). Cells that largely rely on glycolysis are exposed to large amount of oxygen when they reach the blood vessel. These cells could switch the metabolic phenotype to the “W/O” state, which further induces metastasis (10). In future work, the landscape approach (48, 49) could be utilized to quantify the transition processes among glycolysis, OXPHOS, and hybrid metabolism phenotypes. Note that, in addition to the use of glycolysis and OXPHOS, cancer cells can also utilize glutamine as a nitrogen source (5), while cancer-associated fibroblasts can perform glycolysis to fuel cancer cells by producing lactate, known as reverse Warburg effects (50). These aspects are out of the scope of this current study, but they are worth further investigation by using modeling approach.

To help test the basic motions of our core circuit model, we developed two metabolic signatures, one for glycolysis based on the expressions of HIF-1 downstream genes and the other for OXPHOS based on the expression of AMPK downstream genes. These metabolic signatures were applied to distinguish patient samples from multiple cancer types or single cells with different metabolic phenotypes, to identify the anticorrelations between the AMPK and HIF-1 activities, and to elucidate the association of oncogenic pathways with the metabolic states. In general, we expect these metabolic signatures to be a powerful tool to predict the metabolic phenotypes of cancer cells directly from gene expression data.

We proposed to design metabolic therapies by considering the hybrid metabolic phenotype (W/O). A putative strategy would be to drive cancer cells away from the hybrid states. Guided by this assumption, we evaluated the effectiveness of several therapeutic strategies by model simulation. We found that therapy using metformin is more likely to drive cells into the hybrid metabolic phenotype during the treatment time-course than the therapies using glycolysis blockers, such as 3BP and 2DG. Our model also provides an explanation of why some combined therapies are more effective (45, 46), because these therapies are efficient in driving cells out of the hybrid state.

Disclosure of Potential Conflicts of Interest

No potential conflicts of interest were disclosed.

Authors' Contributions

Conception and design: L. Yu, M. Lu, E. Ben-Jacob, B.A. Kaipparettu, J.N. Onuchic

Development of methodology: L. Yu, M. Lu, D. Jia, E. Ben-Jacob, H. Levine, B.A. Kaipparettu, J.N. Onuchic

Analysis and interpretation of data (e.g., statistical analysis, biostatistics, computational analysis): L. Yu, M. Lu, D. Jia, E. Ben-Jacob, H. Levine, B.A. Kaipparettu

Writing, review, and/or revision of the manuscript: L. Yu, M. Lu, D. Jia, H. Levine, B.A. Kaipparettu, J.N. Onuchic

Administrative, technical, or material support (i.e., reporting or organizing data, constructing databases): L. Yu, B.A. Kaipparettu

Study supervision: M. Lu, J. Ma, E. Ben-Jacob, B.A. Kaipparettu, J.N. Onuchic

Grant Support

J.N. Onuchic and H. Levine were awarded the Physics Frontiers Center NSF grant PHY-1427654. H. Levine was awarded the NSF grant DMS-1361411. J.N. Onuchic and H. Levine were awarded the Cancer Prevention and Research Institute of Texas (CPRIT) grants R1110 and R1111, respectively. M. Lu had a training fellowship from the Keck Center for Interdisciplinary Bioscience Training of the Gulf Coast Consortia (CPRIT Grant RP140113). J. Ma is thankful for support from the National Institutes of Health (R01-GM067801, R01-GM116280) and the Welch Foundation (Q-1512). E. Ben-Jacob was also supported by the Tauber Family Funds and the Maguy-Glass Chair in Physics of Complex Systems. B.A. Kaipparettu was awarded by the National Cancer Institute (NCI) grants R21CA173150 and R21CA179720.

The costs of publication of this article were defrayed in part by the payment of page charges. This article must therefore be hereby marked advertisement in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.

Footnotes

  • Note: Supplementary data for this article are available at Cancer Research Online (http://cancerres.aacrjournals.org/).

  • This work is a joint collaboration of the two laboratories of B.A. Kaipparettu and J.N. Onuchic.

  • ↵†Deceased

  • Received August 5, 2016.
  • Revision received December 2, 2016.
  • Accepted December 19, 2016.
  • ©2017 American Association for Cancer Research.

References

  1. 1.↵
    1. Lu M,
    2. Jolly MK,
    3. Ben-Jacob E
    . Toward decoding the principles of cancer metastasis circuits. Cancer Res 2014;74:4574–87.
    OpenUrlAbstract/FREE Full Text
  2. 2.↵
    1. Lu M,
    2. Jolly MK,
    3. Gomoto R,
    4. Huang B,
    5. Onuchic J,
    6. Ben-Jacob E
    . Tristability in cancer-associated microRNA-TF chimera toggle switch. J Phys Chem B 2013;117:13164–74.
    OpenUrl
  3. 3.↵
    1. Hanahan D,
    2. Weinberg RA
    . Hallmarks of cancer: the next generation. Cell 2011;144:646–74.
    OpenUrlCrossRefPubMed
  4. 4.↵
    1. Boroughs LK,
    2. DeBerardinis RJ
    . Metabolic pathways promoting cancer cell survival and growth. Nat Cell Biol 2015;17:351–9.
    OpenUrlCrossRefPubMed
  5. 5.↵
    1. Pavlova NN,
    2. Thompson CB
    . The emerging hallmarks of cancer metabolism. Cell Metab 2016;23:27–47.
    OpenUrlCrossRefPubMed
  6. 6.↵
    1. Warburg O
    . On the origin of cancer cells. Science 1956;123:309–14.
    OpenUrlFREE Full Text
  7. 7.↵
    1. Vander Heiden MG,
    2. Cantley LC,
    3. Thompson CB
    . Understanding the Warburg effect: the metabolic requirements of cell proliferation. Science 2009;324:1029–33.
    OpenUrlAbstract/FREE Full Text
  8. 8.↵
    1. Hsu PP,
    2. Sabatini DM
    . Cancer cell metabolism: Warburg and beyond. Cell 2008;134:703–07.
    OpenUrlCrossRefPubMed
  9. 9.↵
    1. Mathupala SP,
    2. Ko YH,
    3. Pedersen PL
    . The pivotal roles of mitochondria in cancer: Warburg and beyond and encouraging prospects for effective therapies. Biochim Biophys Acta 2010;1797:1225–30.
    OpenUrlPubMed
  10. 10.↵
    1. LeBleu VS,
    2. O'Connell JT,
    3. Herrera KNG,
    4. Wikman H,
    5. Pantel K,
    6. Haigis MC,
    7. et al.
    PGC-1α mediates mitochondrial biogenesis and oxidative phosphorylation in cancer cells to promote metastasis. Nat Cell Biol 2014;16:992–1003.
    OpenUrlCrossRefPubMed
  11. 11.↵
    1. Jiang L,
    2. Xiao L,
    3. Sugiura H,
    4. Huang X,
    5. Ali A,
    6. Kuro-o M,
    7. et al.
    Metabolic reprogramming during TGFβ1-induced epithelial-to-mesenchymal transition. Oncogene 2014;34:3908–16.
    OpenUrlPubMed
  12. 12.↵
    1. Viale A,
    2. Corti D,
    3. Draetta GF
    . Tumors and mitochondrial respiration: a neglected connection. Cancer Res 2015;75:3687–91.
    OpenUrlAbstract/FREE Full Text
  13. 13.↵
    1. Viale A,
    2. Pettazzoni P,
    3. Lyssiotis CA,
    4. Ying H,
    5. Sánchez N,
    6. Marchesini M,
    7. et al.
    Oncogene ablation-resistant pancreatic cancer cells depend on mitochondrial function. Nature 2014;514:628–32.
    OpenUrlCrossRefPubMed
  14. 14.↵
    1. Roesch A,
    2. Vultur A,
    3. Bogeski I,
    4. Wang H,
    5. Zimmermann KM,
    6. Speicher D,
    7. et al.
    Overcoming intrinsic multidrug resistance in melanoma by blocking the mitochondrial respiratory chain of slow-cycling JARID1B high cells. Cancer Cell 2013;23:811–25.
    OpenUrlCrossRefPubMed
  15. 15.↵
    1. Strohecker AM,
    2. White E
    . Targeting mitochondrial metabolism by inhibiting autophagy in BRAF-driven cancers. Cancer Discov 2014;4:766–72.
    OpenUrlAbstract/FREE Full Text
  16. 16.↵
    1. Lu C-L,
    2. Qin L,
    3. Liu H-C,
    4. Candas D,
    5. Fan M,
    6. Li JJ
    . Tumor cells switch to mitochondrial oxidative phosphorylation under radiation via mTOR-mediated hexokinase II inhibition-A Warburg-reversing effect. PLoS One 2015;10:e0121046. doi: 10.1371/journal.pone.0121046.
    OpenUrlCrossRefPubMed
  17. 17.↵
    1. Maiuri MC,
    2. Kroemer G
    . Essential Role for Oxidative Phosphorylation in Cancer Progression. Cell Metab 2015;21:11–12.
    OpenUrlCrossRefPubMed
  18. 18.↵
    1. Tan AS,
    2. Baty JW,
    3. Dong L-F,
    4. Bezawork-Geleta A,
    5. Endaya B,
    6. Goodwin J,
    7. et al.
    Mitochondrial genome acquisition restores respiratory function and tumorigenic potential of cancer cells without mitochondrial DNA. Cell Metab 2015;21:81–94.
    OpenUrlCrossRefPubMed
  19. 19.↵
    1. Xu Q,
    2. Biener-Ramanujan E,
    3. Yang W,
    4. Ramanujan VK
    . Targeting metabolic plasticity in breast cancer cells via mitochondrial complex I modulation. Breast Cancer Res Treat 2015;150:43–56.
    OpenUrl
  20. 20.↵
    1. Porporato PE,
    2. Payen VL,
    3. Pérez-Escuredo J,
    4. De Saedeleer CJ,
    5. Danhier P,
    6. Copetti T,
    7. et al.
    A mitochondrial switch promotes tumor metastasis. Cell Rep 2014;8:754–66.
    OpenUrlCrossRefPubMed
  21. 21.↵
    1. De Saedeleer CJ,
    2. Copetti T,
    3. Porporato PE,
    4. Verrax J,
    5. Feron O,
    6. Sonveaux P
    . Lactate activates HIF-1 in oxidative but not in Warburg-phenotype human tumor cells. PLoS One 2012;7:e46571.
    OpenUrlCrossRefPubMed
  22. 22.↵
    1. Qutub AA,
    2. Popel AS
    . Reactive oxygen species regulate hypoxia-inducible factor 1α differentially in cancer and ischemia. Mol Cell Biol 2008;28:5106–19.
    OpenUrlAbstract/FREE Full Text
  23. 23.↵
    1. Yizhak K,
    2. Le Dévédec SE,
    3. Rogkoti VM,
    4. Baenke F,
    5. Boer VC,
    6. Frezza C,
    7. et al.
    A computational study of the Warburg effect identifies metabolic targets inhibiting cancer migration. Mol Syst Biol 2014;10:744.
  24. 24.↵
    1. Pannala VR,
    2. Dash RK
    . Mechanistic characterization of the thioredoxin system in the removal of hydrogen peroxide. Free Radical Biol Med 2015;78:42–55.
    OpenUrlCrossRefPubMed
  25. 25.↵
    1. Robertson-Tessi M,
    2. Gillies RJ,
    3. Gatenby RA,
    4. Anderson AR
    . Impact of metabolic heterogeneity on tumor growth, invasion, and treatment outcomes. Cancer Res 2015;75:1567–79.
    OpenUrlAbstract/FREE Full Text
  26. 26.↵
    1. Hwang AB,
    2. Ryu E-A,
    3. Artan M,
    4. Chang H-W,
    5. Kabir MH,
    6. Nam H-J,
    7. et al.
    Feedback regulation via AMPK and HIF-1 mediates ROS-dependent longevity in Caenorhabditiselegans. Proc Natl Acad Sci U S A 2014;111:E4458–E67.
    OpenUrlAbstract/FREE Full Text
  27. 27.↵
    1. Dang CV
    . Rethinking the Warburg effect with Myc micromanaging glutamine metabolism. Cancer Res 2010;70:859–62.
    OpenUrlAbstract/FREE Full Text
  28. 28.↵
    1. Plas DR,
    2. Thompson CB
    . Cell metabolism in the regulation of programmed cell death. Trends Endocrinol Metab 2002;13:75–78.
    OpenUrlCrossRefPubMed
  29. 29.↵
    1. Caino MC,
    2. Ghosh JC,
    3. Chae YC,
    4. Vaira V,
    5. Rivadeneira DB,
    6. Faversani A,
    7. et al.
    PI3K therapy reprograms mitochondrial trafficking to fuel tumor cell invasion. Proc Natl Acad Sci U S A 2015;112:8638–43.
    OpenUrlAbstract/FREE Full Text
  30. 30.↵
    1. Zhao Y,
    2. Butler E,
    3. Tan M
    . Targeting cellular metabolism to improve cancer therapeutics. Cell Death Dis 2013;4:e532.
    OpenUrlCrossRefPubMed
  31. 31.↵
    1. Park JH,
    2. Vithayathil S,
    3. Kumar S,
    4. Sung P-L,
    5. Dobrolecki LE,
    6. Putluri V,
    7. et al.
    Fatty acid oxidation-driven Src links mitochondrial energy reprogramming and oncogenic properties in triple-negative breast cancer. Cell Rep 2016;14:2154–65.
    OpenUrl
  32. 32.↵
    1. Camarda R,
    2. Zhou AY,
    3. Kohnz RA,
    4. Balakrishnan S,
    5. Mahieu C,
    6. Anderton B,
    7. et al.
    Inhibition of fatty acid oxidation as a therapy for MYC-overexpressing triple-negative breast cancer. Nat Med 2016;22:427–32.
    OpenUrl
  33. 33.↵
    1. Telang S,
    2. Lane AN,
    3. Nelson KK,
    4. Arumugam S,
    5. Chesney J
    . The oncoprotein H-Ras V12 increases mitochondrial metabolism. Mol Cancer 2007;6:77.
    OpenUrlCrossRefPubMed
  34. 34.↵
    1. Kim K-T,
    2. Lee HW,
    3. Lee H-O,
    4. Kim SC,
    5. Seo YJ,
    6. Chung W,
    7. et al.
    Single-cell mRNA sequencing identifies subclonal heterogeneity in anti-cancer drug responses of lung adenocarcinoma cells. Genome Biol 2015;16:127.
    OpenUrlCrossRefPubMed
  35. 35.↵
    1. McCleland ML,
    2. Adler AS,
    3. Shang Y,
    4. Hunsaker T,
    5. Truong T,
    6. Peterson D,
    7. et al.
    An integrated genomic screen identifies LDHB as an essential gene for triple-negative breast cancer. Cancer Res 2012;72:5812–23.
    OpenUrlAbstract/FREE Full Text
  36. 36.↵
    1. Creighton CJ
    . Multiple oncogenic pathway signatures show coordinate expression patterns in human prostate tumors. PLoS One 2008;3:e1816.
    OpenUrlCrossRefPubMed
  37. 37.↵
    1. Ishikawa K,
    2. Takenaga K,
    3. Akimoto M,
    4. Koshikawa N,
    5. Yamaguchi A,
    6. Imanishi H,
    7. et al.
    ROS-generating mitochondrial DNA mutations can regulate tumor cell metastasis. Science 2008;320:661–64.
    OpenUrlAbstract/FREE Full Text
  38. 38.↵
    1. Piskounova E,
    2. Agathocleous M,
    3. Murphy MM,
    4. Hu Z,
    5. Huddlestun SE,
    6. Zhao Z,
    7. et al.
    Oxidative stress inhibits distant metastasis by human melanoma cells. Nature 2015;527:186–91.
    OpenUrlCrossRefPubMed
  39. 39.↵
    1. Zhou H,
    2. Zhang B,
    3. Zheng J,
    4. Yu M,
    5. Zhou T,
    6. Zhao K,
    7. et al.
    The inhibition of migration and invasion of cancer cells by graphene via the impairment of mitochondrial respiration. Biomaterials 2014;35:1597–607.
    OpenUrlCrossRefPubMed
  40. 40.↵
    1. Goldberg L,
    2. Israeli R,
    3. Kloog Y
    . FTS and 2-DG induce pancreatic cancer cell death and tumor shrinkage in mice. Cell Death Dis 2012;3:e284.
    OpenUrlCrossRefPubMed
  41. 41.↵
    1. Wang Q,
    2. Liang B,
    3. Shirwany NA,
    4. Zou M-H
    . 2-Deoxy-D-glucose treatment of endothelial cells induces autophagy by reactive oxygen species-mediated activation of the AMP-activated protein kinase. PLoS ONE 2011;6:e17234.
    OpenUrlCrossRefPubMed
  42. 42.↵
    1. Ganapathy-Kanniappan S,
    2. Vali M,
    3. Kunjithapatham R,
    4. Buijs M,
    5. Syed L,
    6. Rao P,
    7. et al.
    3-bromopyruvate: a new targeted antiglycolytic agent and a promise for cancer therapy. Curr Pharm Biotechnol 2010;11:510–17.
    OpenUrlCrossRefPubMed
  43. 43.↵
    1. Dowling RJ,
    2. Goodwin PJ,
    3. Stambolic V
    . Understanding the benefit of metformin use in cancer treatment. BMC Med 2011;9:33.
    OpenUrlCrossRefPubMed
  44. 44.↵
    1. Jose C,
    2. Hébert-Chatelain E,
    3. Bellance N,
    4. Larendra A,
    5. Su M,
    6. Nouette-Gaulain K,
    7. et al.
    AICAR inhibits cancer cell growth and triggers cell-type distinct effects on OXPHOS biogenesis, oxidative stress and Akt activation. Biochim Biophys Acta 2011;1807:707–18.
    OpenUrlCrossRefPubMed
  45. 45.↵
    1. Sahra IB,
    2. Laurent K,
    3. Giuliano S,
    4. Larbret F,
    5. Ponzio G,
    6. Gounon P,
    7. et al.
    Targeting cancer cell metabolism: the combination of metformin and 2-deoxyglucose induces p53-dependent apoptosis in prostate cancer cells. Cancer Res 2010;70:2465–75.
    OpenUrlAbstract/FREE Full Text
  46. 46.↵
    1. Cheong J-H,
    2. Park ES,
    3. Liang J,
    4. Dennison JB,
    5. Tsavachidou D,
    6. Nguyen-Charles C,
    7. et al.
    Dual inhibition of tumor energy pathway by 2-deoxyglucose and metformin is effective against a broad spectrum of preclinical cancer models. Mol Cancer Ther 2011;10:2350–62.
    OpenUrlAbstract/FREE Full Text
  47. 47.↵
    1. Mantel CR,
    2. O'Leary HA,
    3. Chitteti BR,
    4. Huang X,
    5. Cooper S,
    6. Hangoc G,
    7. et al.
    Enhancing hematopoietic stem cell transplantation efficacy by mitigating oxygen shock. Cell 2015;161:1553–65.
    OpenUrlCrossRefPubMed
  48. 48.↵
    1. Li C,
    2. Wang J
    . Quantifying the landscape for development and cancer from a core cancer stem cell circuit. Cancer Res 2015;75:2607–18.
    OpenUrlAbstract/FREE Full Text
  49. 49.↵
    1. Li C,
    2. Wang J
    . Quantifying the underlying landscape and paths of cancer. J Royal Soc Interface 2014;11:20140774.
    OpenUrl
  50. 50.↵
    1. Pavlides S,
    2. Whitaker-Menezes D,
    3. Castello-Cros R,
    4. Flomenberg N,
    5. Witkiewicz AK,
    6. Frank PG,
    7. et al.
    The reverse Warburg effect: aerobic glycolysis in cancer associated fibroblasts and the tumor stroma. Cell Cycle 2009;8:3984–4001.
    OpenUrlCrossRefPubMed
PreviousNext
Back to top
Cancer Research: 77 (7)
April 2017
Volume 77, Issue 7
  • Table of Contents
  • Table of Contents (PDF)
  • About the Cover
  • Editorial Board (PDF)

Sign up for alerts

View this article with LENS

Open full page PDF
Article Alerts
Sign In to Email Alerts with your Email Address
Email Article

Thank you for sharing this Cancer Research article.

NOTE: We request your email address only to inform the recipient that it was you who recommended this article, and that it is not junk mail. We do not retain these email addresses.

Enter multiple addresses on separate lines or separate them with commas.
Modeling the Genetic Regulation of Cancer Metabolism: Interplay between Glycolysis and Oxidative Phosphorylation
(Your Name) has forwarded a page to you from Cancer Research
(Your Name) thought you would be interested in this article in Cancer Research.
CAPTCHA
This question is for testing whether or not you are a human visitor and to prevent automated spam submissions.
Citation Tools
Modeling the Genetic Regulation of Cancer Metabolism: Interplay between Glycolysis and Oxidative Phosphorylation
Linglin Yu, Mingyang Lu, Dongya Jia, Jianpeng Ma, Eshel Ben-Jacob, Herbert Levine, Benny Abraham Kaipparettu and José N. Onuchic
Cancer Res April 1 2017 (77) (7) 1564-1574; DOI: 10.1158/0008-5472.CAN-16-2074

Citation Manager Formats

  • BibTeX
  • Bookends
  • EasyBib
  • EndNote (tagged)
  • EndNote 8 (xml)
  • Medlars
  • Mendeley
  • Papers
  • RefWorks Tagged
  • Ref Manager
  • RIS
  • Zotero
Share
Modeling the Genetic Regulation of Cancer Metabolism: Interplay between Glycolysis and Oxidative Phosphorylation
Linglin Yu, Mingyang Lu, Dongya Jia, Jianpeng Ma, Eshel Ben-Jacob, Herbert Levine, Benny Abraham Kaipparettu and José N. Onuchic
Cancer Res April 1 2017 (77) (7) 1564-1574; DOI: 10.1158/0008-5472.CAN-16-2074
del.icio.us logo Digg logo Reddit logo Twitter logo CiteULike logo Facebook logo Google logo Mendeley logo
  • Tweet Widget
  • Facebook Like
  • Google Plus One

Jump to section

  • Article
    • Abstract
    • Introduction
    • Results
    • Discussion
    • Disclosure of Potential Conflicts of Interest
    • Authors' Contributions
    • Grant Support
    • Footnotes
    • References
  • Figures & Data
  • Info & Metrics
  • PDF
Advertisement

Related Articles

Cited By...

More in this TOC Section

Integrated Systems and Technologies

  • FLT PET Imaging of Capecitabine and Trifluridine/Tipiracil
  • In Vivo Vascular Imaging for Head and Neck Cancer Diagnosis
  • Decoding Tumor Growth with Different Treatment Schemes
Show more Integrated Systems and Technologies

Mathematical Oncology

  • Chemotherapeutic Dose Scheduling Based on Tumor Growth Rates
  • Immune Cell–Stroma Interactions in Colorectal Cancer
  • Multiscale Modeling of Inflammation-Induced Tumorigenesis
Show more Mathematical Oncology
  • Home
  • Alerts
  • Feedback
  • Privacy Policy
Facebook  Twitter  LinkedIn  YouTube  RSS

Articles

  • Online First
  • Current Issue
  • Past Issues
  • Meeting Abstracts

Info for

  • Authors
  • Subscribers
  • Advertisers
  • Librarians

About Cancer Research

  • About the Journal
  • Editorial Board
  • Permissions
  • Submit a Manuscript
AACR logo

Copyright © 2021 by the American Association for Cancer Research.

Cancer Research Online ISSN: 1538-7445
Cancer Research Print ISSN: 0008-5472
Journal of Cancer Research ISSN: 0099-7013
American Journal of Cancer ISSN: 0099-7374

Advertisement