Hypothesis and Theory ARTICLE
Metabolic Reprogramming During Multidrug Resistance in Leukemias
- 1Instituto de Bioquímica Médica Leopoldo de Meis, Universidade Federal do Rio de Janeiro, Rio de Janeiro, Brazil
- 2Instituto de Nutrição Josué de Castro, Universidade Federal do Rio de Janeiro, Rio de Janeiro, Brazil
Cancer outcome has improved since introduction of target therapy. However, treatment success is still impaired by the same drug resistance mechanism of classical chemotherapy, known as multidrug resistance (MDR) phenotype. This phenotype promotes resistance to drugs with different structures and mechanism of action. Recent reports have shown that resistance acquisition is coupled to metabolic reprogramming. High-gene expression, increase of active transport, and conservation of redox status are one of the few examples that increase energy and substrate demands. It is not clear if the role of this metabolic shift in the MDR phenotype is related to its maintenance or to its induction. Apart from the nature of this relation, the metabolism may represent a new target to avoid or to block the mechanism that has been impairing treatment success. In this mini-review, we discuss the relation between metabolism and MDR resistance focusing on the multiple non-metabolic functions that enzymes of the glycolytic pathway are known to display, with emphasis with the diverse activities of glyceraldehyde-3-phosphate dehydrogenase.
Multidrug resistance (MDR) in cancer is the major factor impairing the success of conventional chemotherapy (1). Originally, many hopes were placed on the possibility that by inhibiting the activity of ABC transporters a reversal of the resistance would be attained. This was based on the knowledge that the transporters were shown to be capable of mediating the efflux of many chemotherapeutic drugs. Among the efflux transporters of the ABC superfamily, ABCB1 (P-gp), ABCC1 (MRP1), and ABCG2 (BCRP or MXR) have been described as major players in the development of MDR (2–4), with particular emphasis being given to ABCB1. However, it soon became apparent that inhibition of these transporters was not an effective approach since normal cells may also express these transporters and, therefore, inhibitors could sometimes generate an unacceptable toxicity. Furthermore, in experimental situations, where the down regulation of ABCB1 was achieved, a number of other effects were still observed. Now, it is quite clear that MDR is a multifactorial phenomenon that is involved in the regulation of survival and apoptosis, as well as a number of other cellular pathways. ABCB1 transporter expression, in cells with the MDR phenotype, is but one factor linked to pharmacological evasion of chemotherapeutic drugs (1, 5).
A case in point is the observation that hypoxia participates in the regulation of drug resistance. For instance, ABCB1 gene expression and synthesis of functional proteins are induced by hypoxic environments (6). Furthermore, ABC transporters are expressed not only in MDR cancer cells, but also in a number of stem and progenitor cells. Additionally, it has been reported that hypoxia promotes an undifferentiated cell state in various stem and precursor cell populations, as well as in cancer stem cells (7–9). In this respect, it has also been suggested that NOTCH signaling is involved. However, it must be recalled that when cells are under hypoxic conditions, there is a metabolic shift from oxidative phosphorylation to glycolysis (10). This situation contrasts with cells under normoxia, in which glucose is first anaerobically catabolized to pyruvate which is then further catabolized along the Krebs cycle where NADH and FADH2 are reoxidized by the respiratory chain associated to the electron transport system. Incidentally, glycolysis is a hallmark in many types of tumor cells (11). This phenotype is in fact the basis of the so-called Warburg effect, also known as aerobic glycolysis. The Warburg effect describes a situation in which the glycolytic pathway is fully activated even in the presence of adequate oxygen supply (12). Although Warburg originally proposed that cancer was due to an impairment of mitochondrial function, it is accepted today that these organelles retain full oxidative capacities.
It must be mentioned, however, that apart from red blood cells, aerobic glycolysis is prevalent in highly proliferative cells, whether tumoral or not. Stem cells are a case in point (13). The common belief that cells undergoing glycolysis selected an inefficient form of energy production is misguided. Barring the comparative stoichiometry of ATP formation between glycolysis and OXPHOS, aerobic glycolysis is in fact an efficient form of ATP production due to the kinetic properties of the enzymes participating in the pathway which afford very fast fluxes compatible with the ATP demand of the rapidly growing cells. Beyond its role in bioenergetics, glycolysis constitutes a branch of the pentose phosphate pathway (PPP), since glucose-6-phosphate is also the substrate for glucose-6-phosphate dehydrogenase, the first enzyme of that pathway. Thus, glycolysis also contributes to the production of precursors for the biosynthesis of nucleotides (via generation of ribulose-5-phosphate). In addition, the PPP pathway promotes the formation of NADPH, an essential coenzyme for reductive biosynthetic processes such as that of fatty acids. NADPH also has an important role in maintaining the redox equilibrium. Similarly, glycolysis can be considered as an anaplerotic pathway by way of its participation in amino acids synthesis (via 3-phosphoglycerate or via pyruvate). Thus, from an energetic stand point glycolysis more than compensates the relatively small amounts of ATP produced when compared with oxidative phosphorylation.
However, it must be emphasized that tumors are in fact constituted by a mosaic of different cellular subpopulations. As such, from the biochemical point of view tumors can also be envisaged as being functionally heterogeneous. Accordingly, within the context of types of metabolism, tumors can be perceived as composed of subsets of resistant quiescent/slow-cycling cells that occasionally rely more on mitochondrial respiration and less on glycolysis. Likewise the same tumor could also harbor cells that are exclusively glycolytic (14, 15). Interestingly, the possibility of a switch that regulates mitochondrial function in the case of metastasis has been proposed. The results of Porporato et al. showed that overburdening the electron transport system may be an essential step in enhancing migration of cells in vitro and in vivo (16). The authors concluded that in order to achieve metastasis, mitochondria must be active, although not necessarily functional. By extension such findings suggest that in tumor cells there may be switches that constantly activate/inactivate mitochondrial function depending on changes dictated by the microenvironment that affect, for example, the availability of metabolites. The intermittent switching between anaerobic and oxidative metabolism seems to be a feature of metastasis. According to this scheme, accumulating data show that there is a tradeoff involving growth versus migration, i.e., cells which are proliferating prevalently exhibit a glycolytic type of metabolism, whereas migrating cells which proliferate less, rely more on mitochondria (16). Within this framework, it is known that the switch between the two main types of energy metabolism may be regulated by ATP demand. For example, cells expressing ABC transporters on their surfaces require a considerable amount of ATP in order to sustain the drug efflux activity. Additionally, it has been suggested that the transporter activity might suffer the impact of alterations in pH gradient due to the glycolytic phenotype (17). Normal cells maintain a gradient between acidic vesicular compartments and an alkaline cytoplasm a situation not observed in tumor cells that have an acidic cytoplasm. Nevertheless, it has been observed that MCF-7adr, that has an MDR phenotype, presents a similar pH gradient to that of normal cells (18, 19). This gradient contributes to the sequestration of drugs in acidic organelles and subsequent extrusion from the cell (19).
The mitochondrial electron transport chain generates reactive oxygen species (ROS) (20). In some instances, such as partial disintegration of complex I (21–23), oxidative stress may result from excessive production of ROS which is buffered by redox homeostasis. In turn, homeostasis is achieved by the participation of a number of enzymes such as catalase, superoxide dismutase, gluthatione peroxidase, and cofactors such as reduced gluthatione and NADPH. Control of redox equilibrium is important since an imbalance between the amount of ROS produced and of antioxidant systems may lead to DNA damage, particularly mitochondrial DNA, and other cell lesions.
Notwithstanding, the effect of ROS as agents of oxidative stress affecting the expression of ABCB1 is controversial. Both, downregulation and upregulation of the transporter have been reported (24, 25). Antioxidants (26), as well as products resulting from glycolysis can act as scavengers of free radicals. In this way, the energy metabolism could also play an indirect role as modulators of ABCB1 expression. For instance, it is known that under conditions leading to glycolysis inhibition, ABCB1 expression was observed to be downregulated. Conversely, when exogenous pyruvate was added to the tumor cells there was increased drug resistance and transporter expression (25).
Therefore, there is consensus that under normal conditions, glycolysis, and mitochondrial oxidative phosphorilation operate in concert in so far as energy production in the form of ATP is concerned.
Originally, glyceraldehyde-3-phosphate dehydrogenase (GAPDH) was thought to function exclusively as part of the glycolytic process in the cytoplasm, where its role is well established. However, under conditions of oxidative stress GAPDH may also redirect glyceraldehyde-3-P to glucose 6 phosphate (G6P) as a result from the reversal of part of the glycolytic reactions. The formed G6P then becomes the initial substrate of the PPP leading to the increase of NADPH production (27). GAPDH is particularly sensitive to H2O2-induced oxidation and it has been suggested that cytosolic GAPDH might function as a sensor for redox signals and an information hub to transduce these signals (28). Oxidative stress may also promote GAPDH aggregation leading to mitochondrial dysfunction and necrotic cell death via the permeability transition pore (29). Interestingly, under some circumstances GAPDH was shown to bind to band 3 protein, an anion transporter located on the inner side of the red blood cell membrane. In this membrane bound state, GAPDH along with other glycolytic enzymes contributed toward ATP channeling thus allowing its direct consumption by ion pumps without release into the cytoplasm (30). Such a role for GAPDH would be in keeping with an accessory function within the context of drug efflux. GAPDH, however, plays a number of other roles in different cell compartments, where several pools of GAPDH sense cellular stresses and activate cognate pathways to maintain homeostasis or activate cell death (31). It has also been reported that following exposure to stressors GAPDH translocate into the nucleus (32) where it may suffer ADP-ribosylation by NO (33). Apart from GAPDH, other enzymes of the glycolytic pathway are known to display multiple non-metabolic functions. In fact, of the 10 enzymes that constitute glycolysis, at least 7 have been shown to display extra-glycolytic activities that may bear on the MDR phenotype (34). For example, hexokinase II (HKII) binds to the mitochondrial voltage-gated anion channels and thus relieves the negative feedback effect of G6P (35). HKII overexpression is regulated by HIF-1α and also by c-Myc oncogene. HKII is also under the control of many miRNAs. When phosphohexoseisomerase, or phosphoglucoisomerase (PGI) is secreted by cells it acquires the status of a cytokine and is renamed as the autocrine motility factor (AMF). As such, AMF stimulates cell motility and so PGI is thought to be one of the factors driving metastasis (36). Besides, PGI is involved in many other activities such as apoptosis and EMT. Other glycolytic enzymes such as phosphofructokinase 1, aldolase and triose phosphate isomerase, phosphoglycerate mutase, and pyruvate kinase have all been reported to take part in several cellular functions that have a direct relation to tumorigenesis. Hence, it is entirely plausible that individually or collectively, the glycolytic enzymes and in particular GAPDH may constitute integral parts of the MDR phenotype by acting in a non-canonical fashion (37).
Alternatively, the metabolic rewiring of tumor cells may also result from gene rearrangement. Using yeast as a model system, genomic instability and the reprogramming of central metabolism have been approached (38). Regarding gene rearrangement, factors that normally contribute toward the integrity of the replicative process may be compromised. Drug resistance has been observed in yeast in which genome translocants were investigated (39). A parallel between yeast and tumor cells could thus be established even though the detailed mechanisms are still not understood. Indeed in MDR cells the over expression of mini chromosome maintenance 7 (MCM7) was detected (40).
Oxidative stress induces severe damage to proteins, lipids, and DNA. Normally, the generation of oxyradicals is prevented by the mitochondrial antioxidant system; however, the degree of damage will also depend on the repair capacity of the cell. DNA damage leads to the activation of Poly(ADP-ribose) polymerases or PARPs that are nuclear enzymes responsible for catalyzing the attachment of ADP-ribose units from nicotinamide adenine dinucleotide (NAD+) to acceptor proteins involved in the recognition and repair of DNA strands breaks. After repair is completed the Poly(ADP-ribose) (PAR) chains are degraded (41). During caspase mediated apoptosis PARP is cleaved into fragments that in turn inactivate the enzyme inhibiting repair to proceed (42, 43). It has been suggested that PARP inhibitors (PARPi) could be used in combination with chemotherapeutic drugs, but PARPi tends to be extruded from the cell by the ABCB1 transporter in MDR tumors (44).
Cell survival and energy metabolism in MDR tumor cells will be discussed in the next section.
MDR, Chronic Myeloid Leukemia, and Energy Metabolism
Multidrug resistant cells may be very heterogeneous in their characteristics and survival pathways. In the present review, we will highlight, including our own data, some aspects related to the energy metabolism of two MDR chronic leukemia cell lines derived from the CML cell line, K562. The two cell lines were selected after exposure to different chemotherapeutic drugs. One cell line was exposed to vincristine, originating Lucena-1 (45), and the other exposed to daunorubicin, originating FEPS (46). Despite arising from the same parental cell, they were quite distinct. Comparative microarray analysis identified 130 differentially expressed genes between K562 versus Lucena-1, 1,932 between K562 versus FEPS, and 1,211 between Lucena-1 versus FEPS. ABCB1 was overexpressed in both MDR cell lines, but highly overexpressed in FEPS, which is the most resistant line (47). Similarly, comparative proteomics of the parental cell line and the MDR counterparts indicated that K562 presented 560 unique proteins, Lucena-1 had 38 and FEPS 63 unique proteins. Lucena-1 and FEPS shared 929 proteins. From the latter, 112 were common only to Lucena-1 and FEPS (48). Results from another survey investigating the proteomic profiles of K562 and Lucena-1, identified 36 differentially expressed proteins between these two cell lines (40). From those, the leucine-rich PPR motif-containing protein and MCM7, as well as the expression of ABCB1 could be used as markers to identify patients that would respond or fail to therapy with the tirosine kinase inhibitor, Imatinib (40). The importance of ABCB1 in the resistance to Imatinib in the clinical setting has been demonstrated (49, 50), but it is not the only transporter involved. Tumor cells expressing ABCG2 are capable of Imatinib extrusion with great affinity (51, 52). In a different CML model, comparing the proteomic analysis of another MDR cell line derived from K562 selected with doxorubicin, Qinghong et al. highlighted 44 differentially expressed proteins. Some of the differentially expressed proteins were common to those observed in FEPS (53).
The possibility that overexpression of ABCB1 in these cell lines might reflect stem cell characteristics was analyzed by looking at the NOTCH pathway. Notch was overexpressed in FEPS, the most resistant cell line displaying the highest expression of ABCB1, and slow replication time (47). Another pathway present in both, normal stem cells and CML stem cells is the canonical Wnt pathway. This signaling pathway was shown to be more strongly activated to positively regulate ABCB1 in Lucena-1 cells when compared with the non-MDR K562. However, FEPS was not studied on this occasion (40). Conversely, both MDR cell lines had an increased expression of carbonic anhydrase and hemoglobin (47). This was confirmed by proteomics (48) suggesting that these cells are more differentiated compared with the parental erythroleukemic cell line K562 (54). Furthermore, they have the potential to maintain the intracellular pH (17). It has been reported that Imatinib induces erythrocytic differentiation in K562 cells and this is independent of blockade of apoptosis being also observed in resistant clones (55).
One characteristic of the two MDR cell lines described above (Lucena-1 and FEPS) is their antioxidant capacity. The elevated catalase activity observed in Lucena-1 provides these cells with protection against cytotoxic chemicals as well as UV radiation (56). Catalase activity is also elevated in FEPS and these cells also present increased Glucose 6 phosphate dehydrogenase (G6PD) activity (Vidal RS, Faria G, Maia RC, and Rumjanek VM, unpublished data). Despite the well-known role in redox homeostasis, G6PD is also involved in cell growth and signaling and this might be an equally important role in resistant cells (57). Many redox changes are now perceived to allow, in a localized compartment, a rapid and physiological signaling event that may regulate the activity of certain proteins (26, 58). Important players in regulating these intracellular effects are members of the thioredoxin family, including thioredoxin that has increased levels in Lucena-1 and FEPS (48), glutaredoxins, and peroxiredoxins. The third most upregulated gene in Lucena-1 when compared with K562 is SESN3 that catalyzes peroxiredoxins leading to ROS detoxification (47). Furthermore, several members of peroxiredoxins, a family of antioxidant enzymes, are increased in FEPS (PRDX1, PRDX2, PRDX3, and PRDX6) or Lucena-1 (PDRX1) compared with K562 (47, 48). Therefore, ROS generated by the MDR tumor cell lines Lucena-1 and FEPS, are rapidly reduced. In this way, many anticancer drugs that act via generation of oxidative stress become ineffective.
When the MDR cell lines and their parental counterpart, K562, were tested for oxygen consumption measured by high-resolution respirometry, the most resistant cell line, FEPS, reproducibly displayed comparatively lower values in all parameters measured (Figure 1). This result suggested that in FEPS, energy was probably being obtained via glycolysis. Gene expression (47) and protein expression (48) indicate an increase in pyruvate kinase levels in FEPS, whereas in Lucena-1 only increased protein expression was observed (48). Differences in the glycolytic pathway, with higher expression of hexoquinase 2, GAPDH and LDH, were observed using SKOV3 TR, an ovarian cell line transfected with the ABCB1 gene (59).
Figure 1. High-resolution respirometry of K562, Lucena-1, and FEPS cells. Respiratory parameters of intact K562, Lucena-1, and FEPS cells. Routine respiration—basal respiration of intact cells; leak respiration—rate of oxygen consumption after the addition of oligomycin, that is, uncoupled respiration; coupled respiration (routine–leak); ETS—maximum respiratory capacity (induced by the addition of FCCP); ROX—rate of oxygen consumption in the presence of rotenone and antimycin A, that is, respiration not associated with electron transport system. All the parameters were corrected by ROX values. The bars marked with asterisks denote the values that are significantly different with reference to K562 cells.
Proteomic analysis of another MDR line obtained by doxorubicin selection using K562 as the parental cell, described upregulation of fructose-biphosphate aldolase A, fructose-biphosphate aldolase C, transaldolase, and alpha-enolase suggesting that the cells need more energy to survive chemical stress (53). FEPS cells, selected with daunorubicin also presented upregulated fructose-biphosphate aldolase A (48).
GAPDH, MDR, and Cell Death
Among the differences observed in the present study was the variation found in the levels of GAPDH between K562 and its MDR counterparts. Using a proteomic approach (48), it was possible to verify that in relation to K562, the levels of GAPDH were 2.4 times higher in FEPS and 1.2 times in Lucena-1. These results were validated by western blot (Figure 2). Overexpression of GAPDH in Imatinib-resistant cells has been observed by other authors (60, 61). However, these results differ from those found by Cerezo et al. where no difference in GAPDH was observed in daunorubicin resistant cell lines overexpressing ABCB1 (62). Presumably the inconsistency reflects the well-known intrinsic differences encountered in tumor cells of different origin.
Figure 2. Expression of proteins related to cell death-resistance on K562, Lucena-1, and FEPS cells. Equal amounts of total cellular proteins (100 µg) were loaded in each lane for the detection of glyceraldehyde-3-phosphate dehydrogenase (GAPDH), pro-caspase-7, cleaved PARP, and Hsp70 (loading control) by western blot. Densitometric analysis of each lane was calculated using ImageJ Software. The data are expressed as arbitrary units and represent the mean of two independent experiments.
Interestingly, a recent report describes a strategy to circumvent the ABCB1 transporter activity by transferring constructs that specifically inhibited GAPDH into target tumor cells by using liposomes. Such treatment was effective in vitro and in vivo (63).
Therefore, the possibility exists that GAPDH plays a role in the resistance observed in some MDR tumors. Various non-glycolytic roles have been assigned to GAPDH (64). Post-translational modifications of GAPDH may dictate subcellular localization and different functions (31, 65). It has been proposed that Sirtuin1 expression retains GAPDH in the cytosol (66). However, no differences in gene or protein expression of Sirtuin 1, were observed in Lucena-1 or FEPS when compared with their parental cell line K562. As mentioned before, GAPDH translocates to the nucleus where it binds to DNA and participates in a number of DNA-dependent processes (67). However, the lack of differential expression of Sirtuin 1 in Lucena-1 and FEPS does not invalidate the hypothesis of GAPDH mediated action in the nucleus. GAPDH trafficking may occur by a number of mechanisms. For example, it is known that GAPDH can bind to microtubules and may thus get access to several intracellular organelles, including the nucleus (68).
It has been suggested that GAPDH is an inhibitor of caspase independent cell death (CID) (61, 69, 70). Death evasion to a number of chemotherapeutic drugs is the characteristic of cells displaying the MDR phenotype. Defects in the apoptotic process involving caspase activation to induce cell death have been observed in a number of multidrug resistant cell lines (62, 71).
Apoptosis is a result of caspase activation as a consequence of mitochondrial permeabilization and cytochrome c release. However, cells might also be killed following mitochondrial permeabilization even when caspase activation is inhibited, in the process known as CID. It has been suggested that, in such cases, cell death might occur as a result of a collapse of mitochondrial function, or the release of other proteins that could mediate the death process such as apoptosis-inducing factor (AIF), Smac/Diablo, HtrA2/Omi, Endonuclease G (72).
Mitochondrial permeabilization leading to cytochrome c release, and subsequent caspase activation, involves Bcl-2 family members having pro and anti-apoptotic properties. No difference in Bcl-2 levels could be detected when Lucena-1 and FEPS were compared with their parental non-MDR cell line K562 (46). On the other hand, the inhibitor of apoptosis survivin is increased in these cells, similarly to what has been described in K562/ADR (73) and other MDR cell lines (74, 75). Survivin suppresses cell death via caspase inhibition. Therefore, in a situation where caspases are inhibited, cell death could be a result of CID.
In a study, using CML cells, Lavallard et al. described that the tyrosine kinase inhibitor, Imatinib, was able to induce cell death in Bcr-Abl-positive cells by both caspase-dependent and independent manner (61). To induce CID, Imatinib was added to CML cells treated with caspase-inhibitors. In such situation cells transfected with GAPDH were protected from CID, however, these same transfected cells were not protected from Imatinib-induced apoptosis when no caspase-inhibitors were used. Furthermore, Imatinib-resistant K562 cells spontaneously overexpressed GAPDH compared with parental K562 and were protected from CID (61). This finding is in agreement with the results obtained with both Lucena-1 and FEPS MDR cell lines that also showed resistance to Imatinib and displayed higher levels of GAPDH compared with the parental K562 (Figure 2).
In experiments where caspase activity was inhibited but cytochrome c release occurred, it has been suggested that GAPDH acted by increasing the glycolytic metabolism and generating ATP as well as translocating to the nucleus where it is involved in the expression of Atg12 (70). Using GAPDH mutants, that either supported ATP production, but did not translocate to the nucleus or presented nuclear function albeit unable to produce ATP, these authors verified the requirement for the dual role played by GAPDH (70). The induction of Atg12 expression following mitochondrial permeabilization leads to autophagy with the removal of damaged mitochondria and subsequent cell survival (69).
Another death process involves the release of AIF from the mitochondrial intermembrane followed by its translocation to the nucleus. This step is a result of overactivation of the nuclear enzyme poly (ADP-ribose) synthetase 1 (PARP-1) (76). In the mitochondrial intermembrane space, AIF co-localizes with Hsp60. Due to its oxyredutase activity, AIF might act as a scanveger in the mitochondrial electron transport. However, AIF nuclear activity is independent of oxyredutase activity as mutations in the oxyredutase domain do not inhibit death-induction when it translocates to the nucleus (77). Once present in the nucleus AIF induces chromatin condensation and large-scale DNA fragmentation. The use of AIF mutants lacking the DNA-binding property abrogated cell death in spite of the preservation of AIF’s nuclear translocation (78).
In response to oxidative/nitrosative stress GAPDH binds to Siah and translocates to the nucleus where it activates PARP-1 (79). In the nucleus PARP-1 functions as a sensor to regulate celular DNA repair. PARP-1 overactivation, leading to AIF translocation to the nucleus, is a result of an attempt to restore damaged DNA. Usually, after repair is completed, the PAR chains are degraded. However, in cells with severely damaged DNA, nuclear PARP-1 is extensively activated and promotes the synthesis of an excess of PAR polymer. This reaches toxic levels and in the cytosol constitutes a death signal inducing AIF nuclear translocation (76, 80). The activation of PARP-1 results in the depletion of the cellular NAD+ and ATP pools. PARP-1 may be cleaved in vitro and in vivo by caspase-3 and caspase-7 originating two fragments an 89-kD catalytic fragment and a 24-kD DNA binding domain, capable of arresting the process (81). Similar fragments of 89 and 24 kD were detected in the MDR cell line FEPS, whereas only the 89 was observed in Lucena-1 cells (Figure 2). The relationship between PARP-1 activity and drug resistance is unclear but ABCB1 activity occurs at the expense of ATP. It has been described that cells from animals knockout for PARP-1 showed and increased ABCB1 expression and activity (82).
Intracellular Calcium and MDR
The homeostatic control of cytosolic calcium concentration is of fundamental importance and changes in intracellular microenvironmental calcium levels can impact on cell survival, diverse cell functions, and cell death.
Proteins of the S100 family are calcium sensing proteins important in maintaining the homeostasis of the cell (83). Reports show that the calcium binding protein S100A6, also known as calcyclin or CACYBP, interacts with GAPDH (84). S100A6 is present in tumor cells as well as in normal fibroblasts, smooth, and heart muscle cells (83). This protein, as well as its gene expression, is increased twofold in FEPS compared with K562 (48). However, it is not clear how the interaction of S100A6 and GAPDH might play a role inhibiting CID.
Sorcin is another calcium–protein associated with MDR (85–87). Using K562 cells overexpressing sorcin, it was observed that ABCB1 was upregulated. The opposite was also true. In MDR-K562 cells selected by exposure to doxorubicin it was found that sorcin was upregulated, suggesting that ABCB1 and sorcin may regulate the expression of each other (85). When sorcin was analyzed using the MDR cell line Lucena-1 selected by vincristine, an increased gene and protein expression was also observed compared with the parental K562 (47, 48). This result is in agreement with the proteomics of another doxorubicin-induced K562 where the protein sorcin was also increased (53).
Despite the fact that the relationship between sorcin and ABCB1 has been known for a number of years (88), the meaning of such relation is still unknown. Differences in intracellular calcium levels have been reported in cells presenting the MDR phenotype related to ABCB1 overexpression (89). Some of these differences could be attributed to increased sorcin levels but in most cases, a causal relationship was not established.
Cells transfected with cDNA for ABCB1 did not mobilize calcium when the SERCA inhibitor thapsigargin was used (90). The explanation proposed by Gutheil et al. was that thapsigargin was being extruded. However, no calcium mobilization could be induced in Lucena-1 cells using thapsigargin in a condition in which the inhibitor could not be extruded from the cell, suggesting that drug efflux could not fully explain the lack of mobilization (91). Considering that the main intracellular calcium store is the endoplasmic reticulum (ER) and that thapsigargin elicits ER-stress inhibiting the ER calcium pump, this is another death-inducing mechanism that is altered in MDR cells. Similarly, other workers using two types of resistant cells, one selected with vincristine and one with a stable transfection with a gene encoding ABCB1, also described the reduced sensitivity of MDR cells to thapsigargin (92, 93). Moreover, using immunofluorescence, they observed differences in the localization of the ER proteins ryanodine receptor (RyR), inositol 1,4,5-trisphosphate receptor and calnexin, between the MDR cells and their parental (92). Despite the fact that no statistical difference was observed in the amount of GAPDH comparing the two types of resistant cells with the parental one, treatment with thapsigargin decreased the protein content of GAPDH as well as ABCB1, suggesting a relationship between the two (93).
The MDR phenotype is a very complex phenomenon. There is now a growing awareness that drug resistance transcends the ABCB1 transporters and involves other elements associated to metabolic reprogramming. Accordingly, the enzymes of the glycolytic pathway have been shown to exert several regulatory roles that bear on drug resistance. Therefore, pharmacological interference studies on the MDR phenotype may have a better chance to succeed if they are expanded by include as potential targets GAPDH and other enzymes of the glycolytic pathway. In this review, many ancillary roles of such enzymes have been commented and underlined those that affect, for instance, the supply of ATP for drug extrusion purposes. Furthermore, attention was called to the fact that cells displaying the MDR phenotype also displayed enhanced expression/activity of enzymes involved in the maintenance of redox homeostasis of tumor cells. Thus, chemotherapy aiming at the combined abrogation of drug resistance and careful modulation of the redox equilibrium may open up alternative avenues for the control of tumor growth and metastasis.
VR wrote and reviewed the article. RV and JQ contributed to the writing and data gathering. JQ contributed to the editing and submission of this manuscript. FR contributed to the writing and reviewed the article. MR contributed with the experiments involving the high resolution respirometry (oxygen consumption by the cells) which was included in the paper as Figure 1.
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
The reviewer CVB and the handling Editor declared their past shared affiliation.
This work received financial support from the following Brazilian agencies: Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq), Fundação Carlos Chagas Filho de Amparo à Pesquisa do Estado do Rio de Janeiro (FAPERJ), and Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES-PROCAD). RV and JQ were recipients of PhD fellowships from CNPq.
3. Litman T, Druley TE, Stein WD, Bates SE. From MDR to MXR: new understanding of multidrug resistance systems, their properties and clinical significance. Cell Mol Life Sci (2001) 58:931–59. doi:10.1007/PL00000912
5. Tsuruo T, Naito M, Tomida A, Fujita N, Mashima T, Sakamoto H, et al. Molecular targeting therapy of cancer: drug resistance, apoptosis and survival signal. Cancer Sci (2003) 94:15–21. doi:10.1111/j.1349-7006.2003.tb01345.x
7. Schöning JP, Monteiro M, Gu W. Drug resistance and cancer stem cells: the shared but distinct roles of hypoxia-inducible factors HIF1α and HIF2α. Clin Exp Pharmacol Physiol (2017) 44:153–61. doi:10.1111/1440-1681.12693
9. Gustafsson MV, Zheng X, Pereira T, Gradin K, Jin S, Lundkvist J, et al. Hypoxia requires notch signaling to maintain the undifferentiated cell state. Dev Cell (2005) 9:617–28. doi:10.1016/j.devcel.2005.09.010
10. Milane L, Duan Z, Amiji M. Role of hypoxia and glycolysis in the development of multi-drug resistance in human tumor cells and the establishment of an orthotopic multi-drug resistant tumor model in nude mice using hypoxic pre-conditioning. Cancer Cell Int (2011) 11:3. doi:10.1186/1475-2867-11-3
15. Moreno-Sánchez R, Marín-Hernández A, Saavedra E, Pardo JP, Ralph SJ, Rodríguez-Enríquez S. Who controls the ATP supply in cancer cells? Biochemistry lessons to understand cancer energy metabolism. Int J Biochem Cell Biol (2014) 50:10–23. doi:10.1016/j.biocel.2014.01.025
16. Porporato PE, Payen VL, Pérez-Escuredo J, De Saedeleer CJ, Danhier P, Copetti T, et al. A mitochondrial switch promotes tumor metastasis. Cell Rep (2014) 8:754–66. doi:10.1016/j.celrep.2014.06.043
17. Omran Z, Scaife P, Stewart S, Rauch C. Physical and biological characteristics of multi drug resistance (MDR): an integral approach considering pH and drug resistance in cancer. Semin Cancer Biol (2017) 43:42–8. doi:10.1016/j.semcancer.2017.01.002
18. Schindler M, Grabski S, Hoff E, Simon SM. Defective pH regulation of acidic compartments in human breast cancer cells (MCF-7) is normalized in adriamycin-resistant cells (MCF-7adr). Biochemistry (1996) 35:2811–7. doi:10.1021/bi952234e
21. Ricci JE, Muñoz-Pinedo C, Fitzgerald P, Bailly-Maitre B, Perkins GA, Yadava N, et al. Disruption of mitochondrial function during apoptosis is mediated by caspase cleavage of the p75 subunit of complex I of the electron transport chain. Cell (2004) 117:773–86. doi:10.1016/j.cell.2004.05.008
22. Martinvalet D, Dykxhoorn DM, Ferrini R, Lieberman J. Granzyme A cleaves a mitochondrial complex I protein to initiate caspase-independent cell death. Cell (2008) 133:681–92. doi:10.1016/j.cell.2008.03.032
23. Jacquemin G, Margiotta D, Kasahara A, Bassoy EY, Walch M, Thiery J, et al. Granzyme B-induced mitochondrial ROS are required for apoptosis. Cell Death Differ (2015) 22:862–74. doi:10.1038/cdd.2014.180
24. Cen J, Zhang L, Liu F, Zhang F, Ji BS. Long-term alteration of reactive oxygen species led to multidrug resistance in MCF-7 cells. Oxid Med Cell Longev (2016) 2016:7053451. doi:10.1155/2016/7053451
25. Wartenberg M, Richter M, Datchev A, Günther S, Milosevic N, Bekhite MM, et al. Glycolytic pyruvate regulates P-Glycoprotein expression in multicellular tumor spheroids via modulation of the intracellular redox state. J Cell Biochem (2010) 109:434–46. doi:10.1002/jcb.22422
26. Cort A, Ozben T, Saso L, De Luca C, Korkina L. Redox control of multidrug resistance and its possible modulation by antioxidants. Oxid Med Cell Longev (2016) 2016:4251912. doi:10.1155/2016/4251912
27. Ralser M, Wamelink MM, Kowald A, Gerisch B, Heeren G, Struys EA, et al. Dynamic rerouting of the carbohydrate flux is key to counteracting oxidative stress. J Biol (2007) 6:10. doi:10.1186/jbiol61
28. Hildebrandt T, Knuesting J, Berndt C, Morgan B, Scheibe R. Cytosolic thiol switches regulating basic cellular functions: GAPDH as an information hub? Biol Chem (2015) 396:523–37. doi:10.1515/hsz-2014-0295
29. Nakajima H, Itakura M, Kubo T, Kaneshige A, Harada N, Izawa T, et al. Glyceraldehyde-3-phosphate dehydrogenase (GAPDH) aggregation causes mitochondrial dysfunction during oxidative stress-induced cell death. J Biol Chem (2017) 292:4727–42. doi:10.1074/jbc.M116.759084
32. Sawa A, Khan AA, Hester LD, Snyder SH. Glyceraldehyde-3-phosphate dehydrogenase: nuclear translocation participates in neuronal and nonneuronal cell death. Proc Natl Acad Sci U S A (1997) 94:11669–74. doi:10.1073/pnas.94.21.11669
38. Tosato V, Grüning NM, Breitenbach M, Arnak R, Ralser M, Bruschi CV. Warburg effect and translocation-induced genomic instability: two yeast models for cancer cells. Front Oncol (2013) 2:212. doi:10.3389/fonc.2012.00212
39. Nikitin DV, Bruschi CV, Sims J, Breitenbach M, Rinnerthaler M, Tosato V. Chromosome translocation may lead to PRK1-dependent anticancer drug resistance in yeast via endocytic actin network deregulation. Eur J Cell Biol (2014) 93:145–56. doi:10.1016/j.ejcb.2014.03.003
40. Corrêa S, Pizzatti L, Du Rocher B, Mencalha A, Pinto D, Abdelhay E. A comparative proteomic study identified LRPPRC and MCM7 as putative actors in imatinib mesylate cross-resistance in Lucena cell line. Proteome Sci (2012) 10:23. doi:10.1186/1477-5956-10-23
41. De Vos M, Schreiber V, Dantzer F. The diverse roles and clinical relevance of PARPs in DNA damage repair: current state of the art. Biochem Pharmacol (2012) 84:137–46. doi:10.1016/j.bcp.2012.03.018
42. Nicholson DW, Ali A, Thornberry NA, Vaillancourt JP, Ding CK, Gallant M, et al. Identification and inhibition of the ICE/CED-3 protease necessary for mammalian apoptosis. Nature (1995) 376:37–43. doi:10.1038/376037a0
43. Boulares AH, Yakovlev AG, Ivanova V, Stoica BA, Wang G, Iyer S, et al. Role of poly(ADP-ribose) polymerase (PARP) cleavage in apoptosis. J Biol Chem (1999) 274:22932–40. doi:10.1074/jbc.274.33.22932
45. Rumjanek VM, Trindade GS, Wagner-Souza K, de-Oliveira MC, Marques-Santos LF, Maia RC, et al. Multidrug resistance in tumour cells: characterization of the multidrug resistant cell line K562-Lucena 1. An Acad Bras Cienc (2001) 73:57–69. doi:10.1590/S0001-37652001000100007
46. Daflon-Yunes N, Pinto-Silva FE, Vidal RS, Novis BF, Berguetti T, Lopes RR, et al. Characterization of a multidrug-resistant chronic myeloid leukemia cell line presenting multiple resistance mechanisms. Mol Cell Biochem (2013) 383:123–35. doi:10.1007/s11010-013-1761-0
47. Moreira MA, Bagni C, de Pinho MB, Mac-Cormick TM, dos Santos Mota M, Pinto-Silva FE, et al. Changes in gene expression profile in two multidrug resistant cell lines derived from a same drug sensitive cell line. Leuk Res (2014) 38:983–7. doi:10.1016/j.leukres.2014.06.001
49. da Cunha Vasconcelos F, Mauricio Scheiner MA, Moellman-Coelho A, Mencalha AL, Renault IZ, Rumjanek VM, et al. Low ABCB1 and high OCT1 levels play a favorable role in the molecular response to imatinib in CML patients in the community clinical practice. Leuk Res (2016) 51:3–10. doi:10.1016/j.leukres.2016.10.005
51. Burger H, van Tol H, Boersma AW, Brok M, Wiemer EA, Stoter G, et al. Imatinib mesylate (STI571) is a substrate for the breast cancer resistance protein (BCRP)/ABCG2 drug pump. Blood (2004) 104:2940–2. doi:10.1182/blood-2004-04-1398
52. Ozvegy-Laczka C, Hegedus T, Várady G, Ujhelly O, Schuetz JD, Váradi A, et al. High-affinity interaction of tyrosine kinase inhibitors with the ABCG2 multidrug transporter. Mol Pharmacol (2004) 65:1485–95. doi:10.1124/mol.65.6.1485
53. Qinghong S, Shen G, Lina S, Yueming Z, Xiaoou L, Jianlin W, et al. Comparative proteomics analysis of differential proteins in respond to doxorubicin resistance in myelogenous leukemia cell lines. Proteome Sci (2015) 13:1. doi:10.1186/s12953-014-0057-y
54. Rutherford T, Clegg JB, Higgs DR, Jones RW, Thompson J, Weatherall DJ. Embryonic erythroid differentiation in the human leukemic cell line K562. Proc Natl Acad Sci U S A (1981) 78:348–52. doi:10.1073/pnas.78.1.348
55. Jacquel A, Colosetti P, Grosso S, Belhacene N, Puissant A, Marchetti S, et al. Apoptosis and erythroid differentiation triggered by Bcr-Abl inhibitors in CML cell lines are fully distinguishable processes that exhibit different sensitivity to caspase inhibition. Oncogene (2007) 26:2445–58. doi:10.1038/sj.onc.1210034
56. Trindade GS, Capella MA, Capella LS, Affonso-Mitidieri OR, Rumjanek VM. Differences in sensitivity to UVC, UVB and UVA radiation of a multidrug-resistant cell line overexpressing P-glycoprotein. Photochem Photobiol (1999) 69:694–9. doi:10.1111/j.1751-1097.1999.tb03348.x
58. Hanschmann EM, Godoy JR, Berndt C, Hudemann C, Lillig CH. Thioredoxins, glutaredoxins, and peroxiredoxins—molecular mechanisms and health significance: from cofactors to antioxidants to redox signaling. Antioxid Redox Signal (2013) 19:1539–605. doi:10.1089/ars.2012.4599
59. Milane L, Ganesh S, Shah S, Duan ZF, Amiji M. Multi-modal strategies for overcoming tumor drug resistance: hypoxia, the Warburg effect, stem cells, and multifunctional nanotechnology. J Control Release (2011) 155:237–47. doi:10.1016/j.jconrel.2011.03.032
60. Ohmine K, Nagai T, Tarumoto T, Miyoshi T, Muroi K, Mano H, et al. Analysis of gene expression profiles in an imatinib-resistant cell line, KCL22/SR. Stem Cells (2003) 21:315–21. doi:10.1634/stemcells.21-3-315
61. Lavallard VJ, Pradelli LA, Paul A, Bénéteau M, Jacquel A, Auberger P, et al. Modulation of caspase-independent cell death leads to resensitization of imatinib mesylate-resistant cells. Cancer Res (2009) 69:3013–20. doi:10.1158/0008-5472.CAN-08-2731
62. Cerezo D, Cánovas M, García-Peñarrubia P, Martín-Orozco E. Collateral sensitivity to cold stress and differential Bcl-2 family expression in new daunomycin-resistant lymphoblastoid cell lines. Exp Cell Res (2015) 331:11–20. doi:10.1016/j.yexcr.2014.11.017
63. Guan J, Sun J, Sun F, Lou B, Zhang D, Mashayekhi V, et al. Hypoxia-induced tumor cell resistance is overcome by synergistic GAPDH-siRNA and chemotherapy co-delivered by long-circulating and cationic-interior liposomes. Nanoscale (2017) 9:9190–201. doi:10.1039/c7nr02663c
66. Joo HY, Woo SR, Shen YN, Yun MY, Shin HJ, Park ER, et al. SIRT1 interacts with and protects glyceraldehyde-3-phosphate dehydrogenase (GAPDH) from nuclear translocation: implications for cell survival after irradiation. Biochem Biophys Res Communn (2012) 424:681–6. doi:10.1016/j.bbrc.2012.07.006
70. Colell A, Ricci JE, Tait S, Milasta S, Maurer U, Bouchier-Hayes L, et al. GAPDH and autophagy preserve survival after apoptotic cytochrome c release in the absence of caspase activation. Cell (2007) 129:983–97. doi:10.1016/j.cell.2007.03.045
71. Cerezo D, Lencina M, Ruiz-Alcaraz AJ, Ferragut JA, Saceda M, Sanchez M, et al. Acquisition of MDR phenotype by leukemic cells is associated with increased caspase-3 activity and a collateral sensitivity to cold stress. J Cell Biochem (2012) 113:1416–25. doi:10.1002/jcb.24016
75. Liu F, Xie ZH, Cai GP, Jiang YY. The effect of survivin on multidrug resistance mediated by P-glycoprotein in MCF-7 and its adriamycin resistant cells. Biol Pharm Bull (2007) 30:2279–83. doi:10.1248/bpb.30.2279
77. Daugas E, Susin SA, Zamzami N, Ferri KF, Irinopoulou T, Larochette N, et al. Mitochondrio-nuclear translocation of AIF in apoptosis and necrosis. FASEB J (2000) 14:729–39. doi:10.1096/fasebj.14.5.729
78. Ye H, Cande C, Stephanou NC, Jiang S, Gurbuxani S, Larochette N, et al. DNA binding is required for the apoptogenic action of apoptosis inducing factor. Nat Struct Biol (2002) 9:680–4. doi:10.1038/nsb836
79. Nakajima H, Kubo T, Ihara H, Hikida T, Danjo T, Nakatsuji M, et al. Nuclear-translocated glyceraldehyde-3-phosphate dehydrogenase promotes poly(ADP-ribose) polymerase-1 activation during oxidative/nitrosative stress in stroke. J Biol Chem (2015) 290:14493–503. doi:10.1074/jbc.M114.635607
82. Wurzer G, Herceg Z, Wesierska-Gadek J. Increased resistance to anticancer therapy of mouse cells lacking the poly(ADP-ribose) polymerase attributable to up-regulation of the multidrug resistance gene product P-glycoprotein. Cancer Res (2000) 60:4238–44.
84. Filipek A, Wojda U, Leśniak W. Interaction of calcyclin and its cyanogen bromide fragments with annexin II and glyceraldehyde 3-phosphate dehydrogenase. Int J Biochem Cell Biol (1995) 27:1123–31. doi:10.1016/1357-2725(95)00096-8
85. Qi J, Liu N, Zhou Y, Tan Y, Cheng Y, Yang C, et al. Overexpression of sorcin in multidrug resistant human leukemia cells and its role in regulating cell apoptosis. Biochem Biophys Res Commun (2006) 349:303–9. doi:10.1016/j.bbrc.2006.08.042
86. Zhou Y, Xu Y, Tan Y, Qi J, Xiao Y, Yang C, et al. Sorcin, an important gene associated with multidrug-resistance in human leukemia cells. Leuk Res (2006) 30:469–76. doi:10.1016/j.leukres.2005.08.024
87. Colotti G, Poser E, Fiorillo A, Genovese I, Chiarini V, Ilari A. Sorcin, a calcium binding protein involved in the multidrug resistance mechanisms in cancer cells. Molecules (2014) 19:13976–89. doi:10.3390/molecules190913976
88. Sugawara I, Mizumoto K, Ohkochi E, Hamada H, Tsuruo T, Mori S. Immunocytochemical identification and localization of the Mr 22,000 calcium-binding protein (sorcin) in an adriamycin-resistant myelogenous leukemia cell line. Jpn J Cancer Res (1989) 80:469–74. doi:10.1111/j.1349-7006.1989.tb02338.x
89. Sulová Z, Seres M, Barancík M, Gibalová L, Uhrík B, Poleková L, et al. Does any relationship exist between P-glycoprotein-mediated multidrug resistance and intracellular calcium homeostasis. Gen Physiol Biophys (2009) 28 Spec No Focus:F89–95.
91. Wagner-Souza K, Echevarria-Lima J, Rodrigues LA, Reis M, Rumjanek VM. Resistance to thapsigargin-induced intracellular calcium mobilization in a multidrug resistant tumour cell line. Mol Cell Biochem (2003) 252:109–16. doi:10.1023/A:1025586225941
92. Sereš M, Poláková E, Krizanová O, Hudecová S, Klymenko SV, Breier A, et al. Overexpression of P-glycoprotein in L1210/VCR cells is associated with changes in several endoplasmic reticulum proteins that may be partially responsible for the lack of thapsigargin sensitivity. Gen Physiol Biophys (2008) 27:211–21.
Keywords: multidrug resistance, glycolysis, glyceraldehyde-3-phosphate dehydrogenase, leukemia, reactive oxygen species
Citation: Vidal RS, Quarti J, Rodrigues MF, Rumjanek FD and Rumjanek VM (2018) Metabolic Reprogramming During Multidrug Resistance in Leukemias. Front. Oncol. 8:90. doi: 10.3389/fonc.2018.00090
Received: 10 January 2018; Accepted: 15 March 2018;
Published: 04 April 2018
Edited by:Michael Breitenbach, University of Salzburg, Austria
Reviewed by:Johannes A. Mayr, Paracelsus Medical University Salzburg, Austria
Hamid Morjani, Université de Reims Champagne Ardenne, France
Carlo Vito Bruschi, CVBiotechnologies, E.U., Austria
Copyright: © 2018 Vidal, Quarti, Rodrigues, Rumjanek and Rumjanek. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Vivian M. Rumjanek, firstname.lastname@example.org