Cecilie Abildgaard1 , Christina Dahl1 , Astrid L Basse2 , Tao Ma2 and Per Guldberg1*
1 Danish Cancer Society Research Center, Copenhagen, Denmark
2 Department of Biology, University of Copenhagen, Copenhagen, Denmark.
Correspondence: [email protected]
Received: 14 September 2020
Accepted: 4 December 2020
Published: 9 December 2020
Background: Advances in melanoma treatment through targeted inhibition of oncogenic BRAF are limited owing to the development of acquired resistance. The involvement of BRAFV600E in metabolic reprogramming of melanoma cells provides a rationale for co-targeting metabolism as a therapeutic approach.
Methods: We examined the effects of dichloroacetate (DCA), an inhibitor of pyruvate dehydrogenase kinase, on the growth and metabolic activity of human melanoma cell lines. The combined effect of DCA and the BRAF inhibitor vemurafenib was investigated in BRAFV600E -mutated melanoma cell lines. Vemurafenib-resistant cell lines were established in vitro and their sensitivity to DCA was tested.
Results: DCA induced a reduction in glycolytic activity and intracellular ATP levels, and inhibited cellular growth. Co-treatment of BRAFV600E-mutant melanoma cells with DCA and vemurafenib induced a greater reduction in intracellular ATP levels and cellular growth than either compound alone. In addition, melanoma cells with in vitro acquired resistance to vemurafenib retained their sensitivity to DCA.
Conclusions: These results suggest that DCA potentiates the effect of vemurafenib through a cooperative attenuation of energy production. Furthermore, the demonstration of retained sensitivity to DCA in melanoma cells with acquired resistance to vemurafenib could have implications for melanoma treatment.
Keywords: Dichloroacetate, Melanoma, BRAF, Bioenergetics, Metabolism, ATP
Abbreviations: (Acetyl-CoA): Acetyl coenzyme A, (AMPK): AMP activated protein kinase, (DCA): Dichloroacetate, (ECAR): Extracellular acidification rate, (HEMn-LP): Human epidermal melanocytes(IC50): half maximal inhibitory concentration, (LKB1): Liver kinase B1, (MITF): Microphthalmia-associated transcription factor, (OCR): Oxygen consumption rate, (PDH): Pyruvate dehydrogenase, (PDK): Pyruvate dehydrogenase kinase
© 2014 Abildgaard et al.; licensee BioMed Central Ltd.
A hallmark of cancer is the reprogramming of cellular metabolism towards aerobic glycolysis. This metabolic pattern is characterized by increased glucose uptake and highly up-regulated glycolytic activity with fermentation of glucose into lactic acid instead of complete aerobic decomposition in the mitochondria. Aerobic glycolysis, also referred to as the Warburg effect, resembles the anaerobic metabolism of normal cells, but occurs in the context of an adequate oxygen supply . The reprogramming of metabolism in cancer cells is a highly complex and heterogeneous process, which is driven by a wide variety of genetic and non-genetic strategies to overcome energy restriction –.
The BRAF V600E oncogene, present in more than 50% of melanomas , has been directly implicated in the reprogramming of cellular metabolism. The constitutive activity of mutant BRAF reduces the expression of oxidative enzymes and the number of mitochondria, while increasing the expression of glycolytic enzymes and lactic acid production ,. Furthermore, a molecular link was recognized between the RAS-RAF-MEK-ERK-MAPK pathway and the energetic-stress check-point mediated by the liver kinase B1 (LKB1)-AMP activated protein kinase (AMPK) pathway, suggesting a role of BRAFV600E in mediating resistance to energetic stress ,. BRAF affects oxidative metabolism through microphthalmia-associated transcription factor (MITF)-dependent control of the mitochondrial master regulator PGC1α . Previous studies have shown that melanomas expressing PGC1α have a more oxidative phenotype than PGC1α-negative melanomas ,. In addition, BRAFV600E was shown to mediate oncogene-induced senescence through metabolic regulation. This mechanism involves an increase in pyruvate dehydrogenase (PDH) activity through the suppression of pyruvate dehydrogenase kinase (PDK) . PDH controls the coupling between glycolysis and mitochondrial respiration by facilitating the influx of pyruvate into the mitochondria, promoting complete utilization of glucose. The PDK-PDH axis is often dysregulated in cancer, where PDK over-expression reduces the coupling between the two energy systems and thereby contributes to the Warburg effect ,. On the basis of these findings, targeted inhibition of PDK was proposed as a therapeutic option for melanoma, with a possible synergistic effect of chemical BRAFV600E inhibitors, such as vemurafenib ,.
Dichloroacetate (DCA) is an inhibitor of the four isoforms of PDK and was previously used for treatment of lactic acidosis ,, with low toxicity at effective dose levels ,. Several studies have demonstrated that DCA reverses the Warburg effect in cancer cells and negatively affects their growth and survival ,–. This effect was attributed to a normalization of the mitochondrial membrane potential from the hyperpolarized state that characterizes cancer cells. The changes in membrane potential result in the reopening of voltage-gated anion channels and were shown to introduce a re-sensitization to apoptosis, due to a regained ability to release pro-apoptotic mediators . Here we have investigated the effect of DCA on melanoma cells. Specifically, we analyzed cellular responses with regards to metabolism, bioenergetics, growth, proliferation and cell death in melanoma cell lines, primary human melanocytes, and BRAFV600E-mutant melanoma cells with acquired resistance to vemurafenib.
DCA (sodium dichloroacetate) and 2-Deoxy-D-glucose (2-DG) were purchased from Sigma-Aldrich and dissolved in dH2O to working stock concentrations of 1 M. Vemurafenib (PLX4032) was purchased from Selleck Chemicals and dissolved in DMSO to a working stock concentration of 0.05 M.
The melanoma cell lines ED-007, ED-013, ED-024, ED-027, ED-029, ED-034, ED-050, ED-070, ED-071, ED-117, ED-140, ED-179 and ED-196 were obtained from the European Searchable Tumour line Database (ESTDAB, ED) . The melanoma cell line SK-MEL-28 was purchased from ATCC. Primary human epidermal melanocytes (neonatal) from lightly pigmented tissue (HEMn-LP) were purchased from Invitrogen. The melanoma cell lines were cultured at 37°C under 5% CO2 in RPMI-1640 medium supplemented with 10% fetal bovine serum and 1% penicillin/streptomycin. HEMn-LP cells were cultured under the same conditions in 254CF medium supplemented with 1% human melanocyte growth supplement (HMGS-2) and 12-O-tetradecanoyl-phorbol-13-acetate (TPA; 10 ng/ml). All media and supplements were purchased from Invitrogen.
Metabolic characterization was performed on melanoma cell lines and primary human melanocytes using a Seahorse XF96 extracellular flux analyzer (Seahorse Bioscience, Billerica, MA), which performs real-time measurements of extracellular acidification rate (ECAR) and oxygen consumption rate (OCR). An assay was designed to study the capacity of the mitochondrial and glycolytic energy systems. The ECAR and OCR were measured under basal conditions and during successive addition of five metabolic modulators: The ATP synthase inhibitor, oligomycin (1 μM); the mitochondrial membrane permeabilizer, carbonyl cyanide-4-(trifluoromethoxy)phenylhydrazone (FCCP) (1 μM); the inhibitors of mitochondrial respiration, rotenone (1 μM) and antimycin A (1 μM); and the glycolytic inhibitor, 2-DG (100 mM). The XF Cell Mito Stress Kit, containing oligomycin, FCCP, rotenone and antimycin A, was purchased from Seahorse Bioscience.
Intracellular ATP levels were measured using the ATPlite, 1 step Luminescence Assay System (Perkin Elmer), a method based on the reaction of ATP with luciferase and D-luciferin. The cells were seeded in triplicates with 10,000 cells per well and treated with the indicated compounds and vehicle control for 2 or 24 hours. Luminescence was measured with Spectra Max Gemini EM luminescence microplate reader (Molecular Devices) and normalized to background levels.
Crystal violet assay
A crystal violet assay was applied to evaluate the effect of the studied compounds on cell growth. Cells were seeded in duplicates at a suitable density and then treated with DCA, vemurafenib, the two compounds combined and vehicle control. Medium and the treatment compounds were replaced every 48 hours. The experiment was repeated three times independently. To terminate the experiment, medium and unattached cells were removed, and the remaining cells were washed in PBS and fixed with glutaraldehyde for 15 minutes. The fixed cells were incubated with crystal violet solution (0.1% crystal violet, 20% CH3OH) for 1 hour. The amount of dye taken up by the monolayer, proportional to the number of viable cells attached to the well bottom, was quantified by extracting the color with 10% acetic acid and measuring the absorbance at a wavelength of 595 nm. The linear correlation between the absorbance and the number of cells was verified by performing a standard curve. Relative cell growth was determined by normalizing to the untreated controls after background (without cells) subtraction.
Cell proliferation assay
Melanoma cells were seeded in triplicates with 500–1,000 cells per well and treated with DCA at the given concentrations and vehicle control for 96 hours. Proliferation was then measured by detecting BrdU after 12 hours of incorporation into cellular DNA. The procedure was conducted according to the protocol provided with the BrdU Cell Proliferation Assay Kit (Cell Signaling Technology®).
Annexin V-FITC apoptosis detection
Apoptosis detection was performed using an Annexin V-FITC apoptosis detection kit (BD Bioscience), according to the provided protocol. Cells were harvested and washed twice in cold PBS. The cells were then transferred to another tube, spun down and resuspended in binding buffer. From the resuspension, 5 × 105 cells were transferred to FACS tubes and stained with Annexin V-FITC and propidium iodide (PI). After 30 minutes incubation, flow cytometry was performed on a Cytomics FC 500 MPL instrument (Beckman Coulter). Unstained cells were included as control.
Induction of in vitro acquired vemurafenib resistance
Acquired resistance to vemurafenib was induced in seven cultures derived from four BRAFV600E-mutant, vemurafenib-sensitive melanoma cell lines (ED-013, ED-071, ED-196 and SK-MEL-28). Cells were cultured in increasing concentrations of vemurafenib until they grew steadily in a concentration above the IC50, and were then maintained in medium containing vemurafenib.
Pyrosequencing of mutation hotspots in BRAF and NRAS was performed on a PyroMark Q24 platform (Qiagen), using PyroMark Gold Q24 Reagents (Qiagen). The primer sequences are listed in Additional file 1: Table S1.
|Target site||Primer name||Primer sequence (5’-3’)|
PGC1α expression analysis
Total RNA was isolated using RNeasy mini kit (Qiagen) and cDNA was synthesized with the SuperScript™ III Reverse Transcriptase kit (Invitrogen). Oligo dT24 and random hexamers were used as primers for cDNA synthesis. Gene expression of PGC1α was determined with quantitative real-time PCR on Roche LightCycler 2.0 using LigthCycler FastStart DNA MasterPLUS SYBR Green I kit (Roche). The primer sequences were: PPARGC1A_2241F: 5′-GCTGTACTTTTGTGGACGCA-3′ and PPARGC1A_2306R: 5′-GGAAGCAGGGTCAAAGTCAT-3′. The expression was normalized to the expression of the housekeeping gene RPLP0. The primer sequences were: RPLP0_433F: 5′-ACTAAAATCTCCAGGGGCACC-3′ and RPLP0_547R: 5′-ATGACCAGCCCAAAGGAGAA-3′. The two melanoma cell lines ED-050 and SK-MEL-28 were included as positive and negative controls, respectively .
Differences between independent data sets were determined with Student’s t-test. One-way matched-samples ANOVA was used for statistical analysis of variance between different treatments (vehicle control, DCA, vemurafenib and the combination of DCA and vemurafenib). Tukey’s honest significance difference (HSD) multi-comparison test was used to determine statistical significance. Pearson’s correlation coefficient was used to determine correlation between DCA sensitivity and metabolic parameters. A value of 1 indicated a positive correlation, 0 no correlation, and −1 a negative correlation.
The experiments performed in this study involved only commercially available cell lines and therefore required no ethics committee approval.
Metabolic characterization of melanoma cell lines and primary melanocytes
Metabolic profiling of 12 melanoma cell lines and primary human melanocytes (HEMn-LP) was performed using the Seahorse FX96 analyzer. This instrument performs real-time measurements of the extracellular acidification rate (ECAR) and the oxygen consumption rate (OCR), which are indirect measures of glycolytic activity and mitochondrial respiration, respectively . The measurements were performed under basal conditions and during the successive addition of five metabolic modulators (Figure 1A). Compared with normal melanocytes, 11 out of 12 cell lines presented with higher glycolytic rates, as indicated by higher basal glycolytic ECARs, showing that the Warburg effect is a general characteristic of melanoma cells. Furthermore, nine of the cell lines exhibited higher maximal glycolytic capacities compared with melanocytes (Figure 1B). With few exceptions, there were no significant differences in basal and maximal mitochondrial respiration between melanocytes and melanoma cells (Additional file 2: Figure S2). According to the OCR-to-ECAR ratios (OCR/ECAR; Figure 1C), the relative contribution from mitochondrial respiration to ATP production was lower in 10 of the melanoma cell lines compared with melanocytes.
Mitochondrial efficiency (ATP coupling) and performance (respiratory control ratio) were estimated from measurements of two parameters for mitochondrial function; proton leak and the maximal mitochondrial OCR . The proton leak was determined after the addition of the ATP synthase inhibitor oligomycin, and the maximal mitochondrial OCR was determined after the addition of the mitochondrial uncoupler FCCP. This analysis showed a significantly lower ATP coupling in melanoma cells compared with melanocytes (p < 0.05; ED-007: p = 0.07), suggesting a higher proton leak and a less efficient production of ATP relative to the level of oxygen consumption (Figure 1D). Furthermore, six of the melanoma cell lines also had significantly lower respiratory control ratios than melanocytes (p?<?0.01, Figure 1E), indicative of a poor mitochondrial performance.
DCA shifts the metabolism towards mitochondrial respiration and reduces ATP levels
To determine the effect of DCA on melanoma cell metabolism, we analyzed the panel of 12 cell lines using the Seahorse XF96 analyzer. After treatment with 10 mM DCA for 2 hours, all cell lines responded with a reduction in ECAR and an increase in OCR (Figure 2A), indicating a shift towards mitochondrial respiration. The ECAR response was similar among the cell lines, whereas there was a large variation in the OCR response (Figure 2A). The relative changes in ECAR, OCR and OCR/ECAR in response to DCA were concentration dependent (Additional file 3: Figure S3).
To determine if the shift in metabolism interfered with energetic homeostasis, we measured ATP levels in melanoma cells after treatment with 0.1, 1 or 10 mM DCA for 2 hours. All three cell lines tested were able to sustain ATP levels after treatment with low concentrations of DCA (0.1-1 mM), whereas a trend towards reduced ATP was observed in cultures treated with 10 mM DCA (Additional file 3: Figure S3). When cells were treated with DCA (10 or 20 mM) for 24 hours, a significant concentration-dependent decrease in ATP was observed (Figure 2B), indicating a gradual exhaustion of the metabolic system.
DCA reduces the growth of melanoma cells independent of genetic-driver and PGC1α-expression status
To assess a possible clinical application of DCA for treatment of melanoma, we studied its effects on different parameters for cellular growth. The entire panel of cell lines and primary melanocytes were treated with a range of DCA concentrations (0.5-100 mM) for 96 hours (see Figure 3A for representative results). All cell lines showed a concentration-dependent reduction in growth, with IC50 values in the range of 9–38 mM (Table 1), compared with an IC50 value of 70 mM for primary melanocytes (p < 0.001). There was no correlation between the response to DCA and either the BRAF/NRAS mutation status or the expression levels of PGC1α (Table 1).
|Cell line||DCA IC50(mM)4||BRAF/ NRASStatus1||PGC1α expression2,4|
|ED–070||8.9 ± 0.6***||NRASQ61L||1.0 ± 0.2|
|ED–007||12.2 ± 2.2***||WT3||0.6 ± 0.0|
|ED–071||12.3 ± 3.6***||BRAFV600E||0.0 ± 0.0|
|ED–034||12.7 ± 1.7***||BRAFL597S||1.0 ± 0.2|
|ED–013||14.4 ± 2.0***||BRAFV600E||1.9 ± 0.3|
|ED–027||17.7 ± 2.1***||BRAFV600E||0.5 ± 0.1|
|SK–MEL–28||20.0 ± 4.5***||BRAFV600E||0.0 ± 0.0|
|ED–179||20.6 ± 1.8***||NRASQ61R||0.3 ± 0.0|
|ED–024||21.9 ± 1.6***||NRASQ61L||0.0 ± 0.0|
|ED–140||23.9 ± 2.0***||WT||0.4 ± 0.1|
|ED–050||24.1 ± 3.1***||WT||2.7 ± 0.5|
|ED–029||29.7 ± 5.1***||BRAFV600K||1.3 ± 0.2|
|ED–196||35.8 ± 3.2***||BRAFV600E||1.5 ± 0.5|
|ED–117||37.6 ± 2.2***||BRAFV600E||1.2 ± 0.1|
|HEMn–LP||69.1 ± 6.4||WT||1|
|ED–013–R1||12.6 ± 3.0***||BRAFV600E|
|ED–013–R2||13.6 ± 2.4***||BRAFV600E|
|ED–071–R1||12.2 ± 0.9***||BRAFV600E|
|ED–071–R2||13.8 ± 3.7***||BRAFV600E|
|SK–MEL–28–R1||23.1 ± 4.6***||BRAFV600E|
|SK–MEL–28–R2||26.2 ± 8.0**||BRAFV600E|
**p < 0.01; ***p < 0.001 when compared with HEMn-LP.
1BRAF/NRAS status are in accordance with published results .
2Relative to PGC1α expression in normal melanocytes.
4DCA IC50 values and PGC1α expression represent the means of three independent measurements ± standard deviation.
To further characterize the effects of DCA on cellular growth, we measured the incorporation of BrdU in cells treated with 1 or 10 mM DCA for 96 hours. As shown in Figure 3B, all four cell lines tested responded with reduced proliferation, in the range of 11-30% at 1 mM and of 42-88% at 10 mM. We also measured the apoptotic response to DCA by flow-cytometric analysis of annexin V levels. At concentrations of DCA below IC50, the number of annexin V-positive cells did not increase after 96 hours and up to 3 weeks. In contrast, treatment with concentrations above the IC50 increased the number of cells positive for both annexin V and PI, indicating induction of cell death already after 96 hours (Figure 3C).
DCA potentiates the effect of vemurafenib on BRAFV600E-mutant melanoma cells
To investigate if DCA can be used to improve the efficiency of chemical BRAF inhibitors for treatment of melanoma, we tested various combinations of DCA and vemurafenib on cellular growth. Treatment of four BRAFV600E-mutant cell lines with vemurafenib (0.05-5 μM) for 96 hours revealed IC50 values from 0.5 to 4.5 μM, consistent with data from previous studies . When exposing primary melanocytes to the same treatment, we did not reach the IC50 value for these cells, even with the highest concentration tested (5 μM), confirming the specificity of the compound.
When cells were treated with 1 mM DCA in combination with low concentrations of vemurafenib (<IC50), the reduction in cellular growth was more pronounced than with DCA or vemurafenib alone (p < 0.05; Figure 4A, B). DCA did not potentiate the effect of vemurafenib in ED-117 cells, which may be attributed to the inherent resistance of these cells to DCA (IC50 38 mM; Table 1). At IC50 concentrations, both DCA and vemurafenib caused a reduction in intracellular ATP levels when used as single agents and a further reduction when used in combination, although not statistically significant for all cell lines (Figure 4C).
Vemurafenib-resistant cell lines have improved oxidative capacity and retain sensitivity to DCA
Seven cultures of vemurafenib-resistant cells were generated from four BRAFV600E-mutant cell lines by exposing the cells to increasing concentrations of vemurafenib. Cells were considered resistant when they could continuously be propagated at a concentration of vemurafenib above the IC50. DNA was isolated from all seven cultures and tested for BRAF copy number gain and secondary NRAS mutations, which are two well-described mechanisms of acquired resistance to vemurafenib ,. Pyrosequencing revealed an increase in the ratio of BRAF V600E to BRAF WT in one of the resistant cell lines (ED-013-R2) compared with the parental cell line (Additional file 4: Figure S4). No BRAF or NRAS alterations were found in the remaining resistant cell lines.
Metabolic characterization of two of the resistant cell lines (ED-013-R1 and ED-196-R) using the Seahorse XF96 analyzer showed that both resistant cell lines had a transformed metabolic profile with a significantly increased maximal respiratory capacity (Figure 5A), but no changes in basal respiratory OCR, ATP coupling or non-mitochondrial OCR. Interestingly, when tested for sensitivity to DCA, the vemurafenib-resistant cell lines all had IC50 values similar to those for the parental cell lines (Table 1).
Figure 5B illustrates the sensitivity of vemurafenib-resistant cell lines to DCA and vemurafenib, as single agents or in combination, compared with the parental cell lines. The growth of the resistant cell lines was slightly decreased, unaffected or even increased in the presence of vemurafenib after 96 hours, whereas the sensitivity to DCA was similar to the parental cells, both in the presence and absence of vemurafenib (Figure 5B). Similar results were obtained when measuring changes in ATP levels after treatment with the same drugs for 24 hours (Figure 5C). All together, these data show that the vemurafenib-resistant melanoma cells retained sensitivity to DCA, despite the change in their metabolic profile.
Metabolic targeted therapy for cancer has been primarily focused on targeting the energy supply through inhibition of glycolysis. However, the recognition that mitochondria may be active contributors to melanoma progression has increased the attention on oxidative metabolism as a potential therapeutic target ,,. DCA promotes PDK-dependent activation of PDH, reversing lactate production in favor of influx of pyruvate into the mitochondria ,. Through this mechanism, DCA improves the coupling between glycolysis and mitochondrial respiration, which will have a greater impact on cells with a deficient coupling, such as cancer cells . All melanoma cell lines examined in our study responded to DCA with reduced lactate production and an increased OCR. This shift towards mitochondrial respiration was expected to optimize substrate utilization and lead to a more efficient energy yield, but instead led to a significant drop in ATP levels despite an unaffected or even increased mitochondrial ATP coupling. The observed reduction in ECAR in response to DCA suggests that inhibition of glycolysis could be a major contributor to energy deprivation. A glycolysis-inhibitory mechanism of DCA has not been previously described. However, it has been demonstrated that pyruvate kinase, the last ATP-producing site in the glycolytic pathway, is negatively regulated by acetyl coenzyme A (acetyl-CoA) . As PDH activation directly increases the formation of acetyl-CoA , this could explain the DCA-mediated inhibition of glycolysis. The structural similarity between DCA and pyruvate  could also imply a direct inhibition of glycolysis by DCA, possibly through an allosteric feedback mechanism.
The metabolic response to DCA was accompanied by reduced proliferation of melanoma cells, independent of the genetic driver status and metabolic profiles of these cells. Several previous studies have demonstrated an apoptotic effect of DCA on cancer cells ,,,–. However, in accordance with our results, the apoptotic response was only triggered at concentrations too high to be clinically relevant . To further explore the clinical relevance of DCA to melanoma treatment, we examined the efficacy of this agent in combination with the BRAF inhibitor vemurafenib. These experiments demonstrated a potentiating effect of DCA on the growth inhibition of BRAFV600E-mutant melanoma cells. At low concentrations of DCA that alone had no effect on cell growth, the combination with low concentrations of vemurafenib had a significantly stronger growth-reducing effect than vemurafenib alone. This potentiating effect of DCA was also reflected in the reduction of ATP levels. Biochemical analysis has demonstrated the ability of BRAFV600E to uncouple the LKB1-AMPK energy sensing pathway, promoting resistance to energy deprivation and preventing an apoptotic response ,. Treatment with BRAF inhibitors restores this pathway  and may, therefore, potentiate the response to compounds that reduce the generation of ATP. Both DCA and vemurafenib suppress glycolytic activity in melanoma cells and thus render them more dependent on mitochondrial respiration . As glycolysis accounts for a large fraction of the total energy production in these cells, inhibition of this process will place a high demand on the oxidative system for ATP production. The lower performance of the mitochondria in melanoma cells could explain the inability of these cells to sustain ATP levels in the presence of DCA and vemurafenib. The cooperative effect of these compounds in lowering ATP levels suggests that the energetic threshold promoting growth arrest or cell death in melanoma cells can be reached with lower concentrations of vemurafenib in the presence of DCA.
Previous studies have investigated the ability of metabolic modulators to improve the therapeutic effect of BRAF inhibitors for treatment of melanoma. The combination of PLX4720 (a vemurafenib analogue) with either of the two anti-diabetic biguanides, metformin and phenformin, showed synergistic inhibition of melanoma cell viability ,. Both agents impair ATP synthesis through inhibition of the mitochondrial complex I activity, leading to a reduction in the ATP to ADP ratio and activation of the LKB1-AMPK pathway to suppress growth ,. Unlike DCA, metformin and phenformin both stimulate glycolysis and lactic acid production ,, which could explain the growth-stimulating effects of metformin on some melanoma cell lines when used as a single agent. In addition, the concentrations at which metformin was effective were above a therapeutically relevant level . Phenformin was significantly more potent than metformin , but has been associated with a high risk of lactic acidosis , and was taken off the market for treatment of type 2 diabetes in many countries. DCA, on the other hand, was here demonstrated to potentiate the effect of vemurafenib at concentrations down to 1 mM, and was previously shown to have few adverse effects when administered to patients ,,. These findings allude to a therapeutic potential of DCA as a co-drug for vemurafenib treatment of BRAFV600E-mutant melanoma. This was reinforced by the demonstration that sensitivity to DCA was retained in melanoma cell lines with acquired resistance to vemurafenib. Although resistant cells showed an altered metabolic profile with significantly increased maximal mitochondrial respiration, as also shown by Corazao-Rozas et al. , they were as sensitive to DCA as the parental vemurafenib-sensitive cells. Therefore, DCA could possibly provide a strategy to prevent the appearance of vemurafenib-tolerant subpopulations during initial treatment and thereby postpone or prevent the development of resistance.
We here provide a more elaborate understanding of the effects of DCA on the metabolism and growth of melanoma cells. The ability of DCA to lower ATP levels and melanoma growth appears to potentiate the effect of vemurafenib, a drug already used in the clinic for treatment of BRAFV600E-mutant metastatic melanomas. Importantly, melanoma cells with acquired resistance to vemurafenib retained their sensitivity to DCA. These findings should encourage further investigation of this drug combination and the in vivo application of DCA.
Additional file 1: Table S1.(33K, doc)
Pyrosequencing primers for amplification and sequencing of BRAF and NRAS mutation hotspots. Forward, reverse and sequencing primers are denoted F, R and S, respectively.
Additional file 2: Figure S2.(8.5M, tiff)
Basal and maximal mitochondrial OCR values for melanoma cell lines and human epidermal melanocytes (HEMn-LP). The OCR measured after addition of rotenone/antimycin A (non-mitochondrial OCR) was subtracted from all values. The dashed line indicates the basal OCR of HEMn-LP. The indicated values are means of three independent measurements ± standard deviation. Students t-test was used to determine differences between HEMn-LP and the melanoma cell lines (*p < 0.05; **p < 0.01; ***p < 0.001). The maximal mitochondrial OCR for ED-034 was not indicated due to a very large variation among four independent experiments, ranging from 6.19 to 58.36 pmoles/min/1,000 cells.
Additional file 3: Figure S3.(12M, tiff)
Metabolic response of melanoma cells to DCA Relative response in ECAR, OCR, OCR/ECAR and ATP levels after treatment with DCA (0.1, 1 and 10 mM) for 2 h. The error bars in the first three panels represent standard deviations of three repeated measurements of six parallel samples. The error bars in the lower panel represent the standard deviation of three independent experiments. One-way matched-samples ANOVA was used for statistical analysis and Tukey’s HSD test was used to determine statistical significance (*p < 0.05; **p < 0.01).
Additional file 4: Figure S4.(5.6M, tiff)
Allele status of BRAF. Pyrosequencing of the BRAF c.1799 T > A mutation site in ED-013 and the vemurafenib resistant derivative ED-013-R2. The increased BRAF V600E-to-BRAF WT ratio in ED-013-R2 indicates a copy number gain, which could explain the resistance to vemurafenib.
The authors declare that they have no competing interests.
PG and CA planned and organized the study. CA performed the majority of the experiments and the processing of data. CD planned and performed the cell proliferation assay and helped interpret the results. AB and TM helped preparing and optimizing the design for the metabolic analysis on the Seahorse XF instrument. CA, AB and TM discussed and interpreted the results from the metabolic analysis. CA and PG wrote the manuscript with contributions and edits from CD, AB and TM. The final manuscript was read and approved by all authors.
We thank Professor Karsten Kristiansen (The Department of Biology, University of Copenhagen) and Associate Professor Jacob B. Hansen (The Department of Biomedical Sciences, University of Copenhagen) for providing the equipment and technical expertise with respect to the Seahorse XF analysis. This study was supported by grants from the Danish Cancer Society and the Danish Cancer Research Foundation.
REFERENCES1 Warburg O, Wind F, Negelein E: The metabolism of tumors in the body. J Gen Physiol. 1927, 8: 519-530. 10.1085/jgp.8.6.519.
2 Gatenby RA, Gillies RJ: Why do cancers have high aerobic glycolysis?. Nat Rev Cancer. 2004, 4: 891-899. 10.1038/nrc1478.
3 Vander Heiden MG, Cantley LC, Thompson CB: Understanding the Warburg effect: the metabolic requirements of cell proliferation. Science. 2009, 324: 1029-1033. 10.1126/science.1160809.
4 Vazquez F, Lim JH, Chim H, Bhalla K, Girnun G, Pierce K, Clish CB, Granter SR, Widlund HR, Spiegelman BM, Puigserver P: PGC1alpha expression defines a subset of human melanoma tumors with increased mitochondrial capacity and resistance to oxidative stress. Cancer Cell. 2013, 23: 287-301. 10.1016/j.ccr.2012.11.020.
5 vDavies H, Bignell GR, Cox C, Stephens P, Edkins S, Clegg S, Teague J, Woffendin H, Garnett MJ, Bottomley W, Davis N, Dicks E, Ewing R, Floyd Y, Gray K, Hall S, Hawes R, Hughes J, Kosmidou V, Menzies A, Mould C, Parker A, Stevens C, Watt S, Hooper S, Wilson R, Jayatilake H, Gusterson BA, Cooper C, Shipley J: Mutations of the BRAF gene in human cancer. Nature. 2002, 417: 949-954. 10.1038/nature00766.
6 Hall A, Meyle KD, Lange MK, Klima M, Sanderhoff M, Dahl C, Abildgaard C, Thorup K, Moghimi SM, Jensen PB, Bartek J, Guldberg P, Christensen C: Dysfunctional oxidative phosphorylation makes malignant melanoma cells addicted to glycolysis driven by the (V600E)BRAF oncogene. Oncotarget. 2013, 4: 584-599.
7 Haq R, Shoag J, Andreu-Perez P, Yokoyama S, Edelman H, Rowe GC, Frederick DT, Hurley AD, Nellore A, Kung AL, Wargo JA, Song JS, Fisher DE, Arany Z, Widlund HR: Oncogenic BRAF regulates oxidative metabolism via PGC1alpha and MITF. Cancer Cell. 2013, 23: 302-315. 10.1016/j.ccr.2013.02.003.
8 Esteve-Puig R, Canals F, Colome N, Merlino G, Recio JA: Uncoupling of the LKB1-AMPKalpha energy sensor pathway by growth factors and oncogenic BRAF.PLoS One 2009, 4:e4771.,
9 Zheng B, Jeong JH, Asara JM, Yuan YY, Granter SR, Chin L, Cantley LC: Oncogenic B-RAF negatively regulates the tumor suppressor LKB1 to promote melanoma cell proliferation. Mol Cell. 2009, 33: 237-247. 10.1016/j.molcel.2008.12.026.
10 Kaplon J, Zheng L, Meissl K, Chaneton B, Selivanov VA, Mackay G, van der Burg SH, Verdegaal EM, Cascante M, Shlomi T, Gottlieb E, Peeper DS: A key role for mitochondrial gatekeeper pyruvate dehydrogenase in oncogene-induced senescence. Nature. 2013, 498: 109-112. 10.1038/nature12154.
11 Baumunk D, Reichelt U, Hildebrandt J, Krause H, Ebbing J, Cash H, Miller K, Schostak M, Weikert S: Expression parameters of the metabolic pathway genes pyruvate dehydrogenase kinase-1 (PDK-1) and DJ-1/PARK7 in renal cell carcinoma (RCC). World J Urol. 2013, 31: 1191-1196. 10.1007/s00345-012-0874-5.
12 Hur H, Xuan Y, Kim YB, Lee G, Shim W, Yun J, Ham IH, Han SU: Expression of pyruvate dehydrogenase kinase-1 in gastric cancer as a potential therapeutic target. Int J Oncol. 2013, 42: 44-54.
13 Kluza J, Corazao-Rozas P, Touil Y, Jendoubi M, Maire C, Guerreschi P, Jonneaux A, Ballot C, Balayssac S, Valable S, Corroyer-Dulmont A, Bernaudin M, Malet-Martino M, de Lassalle EM, Maboudou P, Formstecher P, Polakowska R, Mortier L, Marchetti P: Inactivation of the HIF-1alpha/PDK3 signaling axis drives melanoma toward mitochondrial oxidative metabolism and potentiates the therapeutic activity of pro-oxidants. Cancer Res. 2012, 72: 5035-5047. 10.1158/0008-5472.CAN-12-0979.
14 Stacpoole PW, Lorenz AC, Thomas RG, Harman EM: Dichloroacetate in the treatment of lactic acidosis. Ann Intern Med. 1988, 108: 58-63. 10.7326/0003-4819-108-1-58.
15 Stacpoole PW: The pharmacology of dichloroacetate. Metabolism. 1989, 38: 1124-1144. 10.1016/0026-0495(89)90051-6.
16 Michelakis ED, Webster L, Mackey JR: Dichloroacetate (DCA) as a potential metabolic-targeting therapy for cancer. Br J Cancer. 2008, 99: 989-994. 10.1038/sj.bjc.6604554.
17 Stacpoole PW, Kerr DS, Barnes C, Bunch ST, Carney PR, Fennell EM, Felitsyn NM, Gilmore RL, Greer M, Henderson GN, Hutson AD, Neiberger RE, O’Brien RG, Perkins LA, Quisling RG, Shroads AL, Shuster JJ, Silverstein JH, Theriaque DW, Valenstein E: Controlled clinical trial of dichloroacetate for treatment of congenital lactic acidosis in children. Pediatrics. 2006, 117: 1519-1531. 10.1542/peds.2005-1226.
18 Bonnet S, Archer SL, lalunis-Turner J, Haromy A, Beaulieu C, Thompson R, Lee CT, Lopaschuk GD, Puttagunta L, Bonnet S, Harry G, Hashimoto K, Porter CJ, Andrade MA, Thebaud B, Michelakis ED: A mitochondria-K+ channel axis is suppressed in cancer and its normalization promotes apoptosis and inhibits cancer growth. Cancer Cell. 2007, 11: 37-51. 10.1016/j.ccr.2006.10.020.
19 Michelakis ED, Sutendra G, Dromparis P, Webster L, Haromy A, Niven E, Maguire C, Gammer TL, Mackey JR, Fulton D, Abdulkarim B, McMurtry MS, Petruk KC: Metabolic modulation of glioblastoma with dichloroacetate. Sci Transl Med. 2010, 2: 31-34. 10.1126/scitranslmed.3000677.
20 Wong JY, Huggins GS, Debidda M, Munshi NC, De VI: Dichloroacetate induces apoptosis in endometrial cancer cells. Gynecol Oncol. 2008, 109: 394-402. 10.1016/j.ygyno.2008.01.038.
21 Cao W, Yacoub S, Shiverick KT, Namiki K, Sakai Y, Porvasnik S, Urbanek C, Rosser CJ: Dichloroacetate (DCA) sensitizes both wild-type and over expressing Bcl-2 prostate cancer cells in vitro to radiation. Prostate. 2008, 68: 1223-1231. 10.1002/pros.20788.
22 Robinson J, Roberts CH, Dodi IA, Madrigal JA, Pawelec G, Wedel L, Marsh SG: The European searchable tumour line database. Cancer Immunol Immunother. 2009, 58: 1501-1506. 10.1007/s00262-008-0656-5.
23 Wu M, Neilson A, Swift AL, Moran R, Tamagnine J, Parslow D, Armistead S, Lemire K, Orrell J, Teich J, Chomicz S, Ferrick DA: Multiparameter metabolic analysis reveals a close link between attenuated mitochondrial bioenergetic function and enhanced glycolysis dependency in human tumor cells. Am J Physiol Cell Physiol. 2007, 292: C125-C136. 10.1152/ajpcell.00247.2006.
24 Brand MD, Nicholls DG: Assessing mitochondrial dysfunction in cells. Biochem J. 2011, 435: 297-312. 10.1042/BJ20110162.
25 Dahl C, Christensen C, Jonsson G, Lorentzen A, Skjodt ML, Borg A, Pawelec G, Guldberg P: Mutual Exclusivity Analysis of Genetic and Epigenetic Drivers in Melanoma Identifies a Link Between p14ARF and RARbeta Signaling. Mol Cancer Res. 2013, 11: 1166-1178. 10.1158/1541-7786.MCR-13-0006.
26 Sondergaard JN, Nazarian R, Wang Q, Guo D, Hsueh T, Mok S, Sazegar H, MacConaill LE, Barretina JG, Kehoe SM, Attar N, von EE, Zuckerman JE, Chmielowski B, Comin-Anduix B, Koya RC, Mischel PS, Lo RS, Ribas A: Differential sensitivity of melanoma cell lines with BRAFV600E mutation to the specific Raf inhibitor PLX4032.J Transl Med 2010, 8:39.,
27 Nazarian R, Shi H, Wang Q, Kong X, Koya RC, Lee H, Chen Z, Lee MK, Attar N, Sazegar H, Chodon T, Nelson SF, McArthur G, Sosman JA, Ribas A, Lo RS: Melanomas acquire resistance to B-RAF(V600E) inhibition by RTK or N-RAS upregulation. Nature. 2010, 468: 973-977. 10.1038/nature09626.
28 Shi H, Moriceau G, Kong X, Lee MK, Lee H, Koya RC, Ng C, Chodon T, Scolyer RA, Dahlman KB, Sosman JA, Kefford RF, Long GV, Nelson SF, Ribas A, Lo RS: Melanoma whole-exome sequencing identifies (V600E)B-RAF amplification-mediated acquired B-RAF inhibitor resistance.Nat Commun 2012, 3:724.,
29 de Barbi MM, Vincent G, Fayewicz SL, Bateman NW, Hood BL, Sun M, Suhan J, Duensing S, Yin Y, Sander C, Kirkwood JM, Becker D, Conrads TP, Van HB, Moschos SJ: Mitochondrial respiration–an important therapeutic target in melanoma.PLoS One 2012, 7:e40690.,
30 Gilbert RJ, Klein RA:Pyruvate kinase: a carnitine-regulated site of ATP production in Trypanosoma brucei brucei. Comp Biochem Physiol B. 1984, 78: 595-599.
31 Kobayashi K, Neely JR:Mechanism of pyruvate dehydrogenase activation by increased cardiac work. J Mol Cell Cardiol. 1983, 15: 369-382. 10.1016/0022-2828(83)90321-8.
32 Stockwin LH, Yu SX, Borgel S, Hancock C, Wolfe TL, Phillips LR, Hollingshead MG, Newton DL:Sodium dichloroacetate selectively targets cells with defects in the mitochondrial ETC. Int J Cancer. 2010, 127: 2510-2519. 10.1002/ijc.25499.
33 Madhok BM, Yeluri S, Perry SL, Hughes TA, Jayne DG:Dichloroacetate induces apoptosis and cell-cycle arrest in colorectal cancer cells. Br J Cancer. 2010, 102: 1746-1752. 10.1038/sj.bjc.6605701.
34 Tong J, Xie G, He J, Li J, Pan F, Liang H: Synergistic antitumor effect of dichloroacetate in combination with 5-fluorouracil in colorectal cancer.J Biomed Biotechnol 2011, 2011:740564.,
35 Niehr F, Von EE, Attar N, Guo D, Matsunaga D, Sazegar H, Ng C, Glaspy JA, Recio JA, Lo RS, Mischel PS, Comin-Anduix B, Ribas A: Combination therapy with vemurafenib (PLX4032/RG7204) and metformin in melanoma cell lines with distinct driver mutations.J Transl Med 2011, 9:76.,
36 Yuan P, Ito K, Perez-Lorenzo R, Del GC, Lee JH, Shen CH, Bosenberg MW, McMahon M, Cantley LC, Zheng B:Phenformin enhances the therapeutic benefit of BRAFV600E inhibition in melanoma. Proc Natl Acad Sci U S A. 2013, 110: 18226-18231. 10.1073/pnas.1317577110.
37 Buzzai M, Jones RG, Amaravadi RK, Lum JJ, DeBerardinis RJ, Zhao F, Viollet B, Thompson CB:Systemic treatment with the antidiabetic drug metformin selectively impairs p53-deficient tumor cell growth. Cancer Res. 2007, 67: 6745-6752. 10.1158/0008-5472.CAN-06-4447.
38 Choi YW, Lim IK:Sensitization of metformin-cytotoxicity by dichloroacetate via reprogramming glucose metabolism in cancer cells. Cancer Lett. 2014, 346: 300-8. 10.1016/j.canlet.2014.01.015.
39 Bergman U, Boman G, Wiholm BE:Epidemiology of adverse drug reactions to phenformin and metformin. Br Med J. 1978, 2: 464-466. 10.1136/bmj.2.6135.464.
40Dunbar EM, Coats BS, Shroads AL, Langaee T, Lew A, Forder JR, Shuster JJ, Wagner DA, Stacpoole PW:Phase 1 trial of dichloroacetate (DCA) in adults with recurrent malignant brain tumors. Invest New Drugs. 2013, 32: 452-64. 10.1007/s10637-013-0047-4.
41 Corazao-Rozas P, Guerreschi P, Jendoubi M, Andre F, Jonneaux A, Scalbert C, Garcon G, Malet-Martino M, Balayssac S, Rocchi S, Savina A, Formstecher P, Mortier L, Kluza J, Marchetti P:Mitochondrial oxidative stress is the Achille’s heel of melanoma cells resistant to Braf-mutant inhibitor. Oncotarget. 2013, 4: 1986-1998.