Targeting of Evolutionarily Acquired Cancer Cell Phenotype by Exploiting pHi-Metabolic Vulnerabilities

Evolutionary dynamics can be used to control cancers when a cure is not clinically considered to be achievable. Understanding Darwinian intratumoral interactions of microenvironmental selection forces can be used to steer tumor progression towards a less invasive trajectory. Here, we approach intratumoral heterogeneity and evolution as a dynamic interaction among subpopulations through the application of small, but selective biological forces such as intracellular pH (pHi) and/or extracellular pH (pHe) vulnerabilities. Increased glycolysis is a prominent phenotype of cancer cells under hypoxia or normoxia (Warburg effect). Glycolysis leads to an important aspect of cancer metabolism: reduced pHe and higher pHi. We recently showed that decreasing pHi and targeting pHi sensitive enzymes can reverse the Warburg effect (WE) phenotype and inhibit tumor progression. Herein, we used diclofenac (DIC) repurposed to control MCT activity, and Koningic acid (KA) that is a GAPDH partial inhibitor, and observed that we can control the subpopulation of cancer cells with WE phenotype within a tumor in favor of a less aggressive phenotype without a WE to control progression and metastasis. In a 3D spheroid co-cultures, we showed that our strategy can control the growth of more aggressive MDA-MB-231 cells, while sparing the less aggressive MCF7 cells. In an animal model, we show that our approach can reduce tumor growth and metastasis. We thus propose that evolutionary dynamics can be used to control tumor cells’ clonal or sub-clonal populations in favor of slower growth and less damage to patients. We propose that this can result in cancer control for tumors where cure is not an option.


Introduction
The initiation and development of cancers is associated with major metabolic alterations in response to dynamically changing microenvironmental conditions such as hypoxia and acidosis [1]. These forces will select for the fittest phenotype in the context of the current microenvironment. Cancer cells that are more plastic and adaptable can acclimate to an increasing array of emergent microenvironmental selective forces. Combined, these factors define tumors as highly dynamic ecosystems in which many different cancer cell subpopulations compete for space and resources [2]. The microenvironment around the cells can alter the local fitness of cancer cell subpopulations in a tumor; leading to possible dramatically different evolutionary trajectories selecting for cells with different phenotypic and genotypic properties [3]. We have previously reported how acid-producing and acid-resistant phenotypes can engineer the tumor ecosystem in order to increase

Diclofenac and Koningic Acid Reduce the Glycolytic Activity of Cancer Cells
H + -monocarboxylate transporters (MCTs) play a major role in transporting organic anions such as lactate across cell membranes. Since the transport cycle is coupled to H + ions, the activity of MCTs also affects intracellular pH (pHi) within cells [22]. It has been shown that MCT inhibition decreases cancer cells' aggressive phenotypes [9,23,24] by reducing the lactate efflux, leading to acidification of pHi and buildup of intracellular [lactate] [6]. It has been shown that diclofenac can inhibit MCTs [7] and that koningic acid (KA) can inhibit GAPDH [10,25,26], which is very sensitive to pHi [6] ( Figure 1A). As proof of principle, we measured the effect of diclofenac and KA on glucose consumption and lactate production of MCF7 and MDA-MB-231 breast cancer cell lines. In both cell lines, diclofenac and KA reduced the lactate production ( Figure 1B) and glucose consumption ( Figure 1C) significantly within 72 h of treatment. Notably the effect was larger in the more highly glycolytic MDA-mb-231 cells. We then measured the activity of both drugs in the 3D culture of MCF7 and MDA-MB-231 spheroids over time. We observed that both drugs showed effects as early as 24 h ( Figure S1A,B). These results indicated that both drugs are effective in inhibiting glycolysis as proposed and can work both in 2D and 3D experiments. We also measured the glutamine and glutamate concentration in the media following treatments of both cell lines with both drugs under normoxia and hypoxia and found no differences among the groups ( Figure S1C).

Diclofenac Inhibits MCT Activity
To test whether the effect of diclofenac treatment on glycolysis was due to inhibition of lactate export, we measured MCT activity before and after treatment with diclofenac using pHi as a reporter of transmembrane flux. To rule out if the effect on pHi change is because of NSAID activity of diclofenac, we also measured the MCT activity on the same cells treated with aspirin. Using a superfusion microscopy system and cells loaded with the pH reporter dye SNARF-1, we measured the time course of pHi during superfusate maneuvers that trigger net transmembrane flux of lactate ( Figure 1D). When presented with extracellular lactate, cells acidified, and when extracellular lactate was withdrawn, cells alkalinized. The rate of these pHi changes are measures of MCT activity. In the presence of diclofenac, however, efflux of lactate was completely inhibited in both cell lines. Controls included another NSAID, aspirin, which was without effect. As these are acute effects, cell death was not expected to play a role (vide infra), and we further ruled it out, due to the inability of dead cells to retain the SNARF-1 dye. To confirm that diclofenac was acting on a specific transport process, i.e., MCT, rather than causing a generalized tightening of membrane permeability, experiments were repeated with acetate in place of lactate. Compared to lactate, acetate can permeate the cell membrane in an MCT-independent manner by partitioning across the lipid matrix as acetic acid, and its efflux was unaffected by diclofenac ( Figure 1E). This suggests that the action of diclofenac is targeted to MCT, the main route for lactate traffic in and out of cells. The results of this assay are not indicative of what is expected to be seen under physiological conditions. The large changes in pH observed in these experiments are not sustained for long periods of time, and are induced by superphysiological lactate levels. The purpose of these experiments was solely to determine the effects of Diclofenac on a cells ability to export lactate; not to demonstrate specific pHi values under drug treatment. cells were treated with diclofenac or Aspirin to measure the effect of drug on their lactate transport abilities. Line (a) represents the switch from a 30 mM lactate environment to a 0 mM Lactate environment. Line (b) represents the switch from a 0 mM lactate environment to a 30 mM lactate environment. Line (c) represents the switch from a 30 mM lactate environment − drug treatment to a 30 mM lactate environment + drug treatment. Line (d) represents the switch from a 30 mM lactate environment + drug treatment to a 0 mM lactate environment + drug treatment. (E) MCF7 and MDA-MB-231 cells were treated with diclofenac or Aspirin to measure the effect of drug on their pHi altering effects. Line (a) represents the switch from a 30 mM acetate environment to a 0 mM acetate environment. Line (b) represents the switch from a 0 mM acetate environment to a 30 mM acetate environment. Line (c) represents the switch from a 30 mM acetate environment − drug treatment to a 30 mM acetate environment + drug treatment. Line (d) represents the switch from a 30 mM acetate environment + drug treatment to a 0 mM acetate environment + drug treatment. p-values are represented as follows: *** p < 0.001, **** p < 0.0001.

Diclofenac and Koningic Acid Decrease the Viability of Cancer Cells
To test if the effects on glycolysis translate into the effective killing of cancer cells, we assessed the survival of cells using CCK8 viability assays and CelltiterGlo (see Methods) in both hypoxia and normoxia. We experimented with both normoxic and hypoxic conditions as they both commonly occur in the microenvironment of solid tumors such as breast cancer, and cell viability will be expected to be much more sensitivity to glycolysis inhibition under hypoxic conditions. In solid tumors, cells that are further than 160 micrometers from the vasculature they experience a hypoxic condition that makes the cells glycolysis dependent [27,28]. Many cancer cells overexpress or over-activate MCTs to deplete the lactate produced under the hypoxic conditions, although this response is not universal as shown in supplemental Figure S2. Under normoxia ( Figure 2A) and hypoxia ( Figure 2B) cell survival was measured after treatment with diclofenac, KA, or both in MCF7 and MDA-MB-231 cells. As shown in Figure 2A, diclofenac alone reduced the viability of MDA-mb-231, and MCF10/DCIS cells, which are highly glycolytic and tumorigenic, as well as MCF-7 cells which are tumorigenic and only moderately glycolytic. Notably there was no effect on viability of MCF10AT cells, which are neither tumorigenic nor glycolytic. More pronounced effects with a similar distribution were observed with KA and the combination of DIC + KA. Surprisingly, the effect of DIC on viability under hypoxia ( Figure 2B) was less significant cf. normoxia (Figure 2A), whereas the effects of KA and the DIC + KA combination were no different. To test the activity of the drug on viability and H + export flux in the same experiment we performed a viability assay while monitoring the acidification of conditioned media. Inhibition of MCTs or GAPDH should reduce rate of media acidification. After treatment and before viability assays we measured the pHe of the media using the Five Easy TM /FiveGo TM pH meter (Mettler Toledo, Columbus, OH, USA see Methods). The acidification was reduced by both drugs and the pH of media was neutral, validating the effect of both diclofenac and KA on lowering glycolytic activity of cells ( Figure S3).
To determine the effect of diclofenac and KA on the pHi of cancer cells, we measured the pHi of a panel of breast cancer cell lines (MDA-MB-231, MCF-DCIS, MCF7, 4T1, H605, and T47D), loaded with SNARF-1-AME ( Figure 2C). For 16 h, cells were treated with diclofenac (100 µM), KA (1 µM), or aspirin (300 µM), and then imaged in buffered RPMI media. All cell lines exhibited a significantly lower pHi with the treatment of diclofenac, and 4 of the 6 diclofenac treated cell lines had significantly lower pH than the aspirin group.
Considering the different microenvironment in solid tumors such as hypoxia and acidosis and their combination we wanted to test if inhibiting MCTs and GAPDH have an effect on the viability of cancer cells. These conditions are also very important to be targeted because they are specific for tumors and not experienced by normal cells. Therefore, we performed the survival assays on both MCF7 and MDA-MB-231 cells under four possible combinations of change in oxygen level and pH as following: (i) Normal pH-Normoxia, (ii) Low pH-Normoxia, (iii) Normal pH-Hypoxia, and (iv) Low pH-Hypoxia. Targeting MCTs and GAPDH under different microenvironments reduces the cancer cell's viability ( Figure  S4). The effect was greatest with the combination of both drugs and under the most unique condition of solid tumors -hypoxia, and acidic pH ( Figure S4).

Diclofenac and Koningic Acid Reduce the Warburg Phenotype of Cancer Cells
The Warburg effect (WE) phenotype is defined by a high glycolytic rate even in the presence of oxygen (aerobic glycolysis) and is associated with the progression and aggressiveness of cancers [29][30][31][32]. We recently showed that inhibiting proton pumps such as MCTs and pH-sensitive glycolytic enzymes such as GAPDH or GPI, reduces the WE phenotype of cancer cells [6]. Here we used the same strategy to reduce the WE phenotype by repurposing diclofenac from a NSAID to an MCT inhibitor and a natural product produced by fungi, koningic acid to target GAPDH. Above, we already showed that these compounds can reduce the glycolytic phenotype of cancer cells, and here we are using Seahorse measurements to confirm the inhibitor effect of these products on aerobic glycoysis. Seahorse can measure the glycolytic condition through simultaneous of proton efflux (ECAR) and mitochondrial respiration (OCR), which can be used to assess WE phenotype [6,33]. In the first experiment we interrogated at the real-time effect of diclofenac of KA on ECAR and OCR of MCF7 cancer cells ( Figure 3A,B). KA and diclofenac monotonically increased the OCR compared to control and the combination of the two had the greatest effect ( Figure 3A). KA and diclofenac both caused an initial drop of ECAR that lasted for 150 min that slowly returned to control levels. The recovery of ECAR to KA was more rapid compared to either conditions that contained diclofinac, which may be due to differences between GAPDH and MCT regulation ( Figure 3B). This could imply that even if the cells over-activate or overexpress the GAPDH, because the pHi is too low the phenotype can't be survived. We also measured the WE phenotype (ECAR/OCR) of these cells in three different time points: (i) right before drug injection, (ii) 5 h after injections, and (iii) 12 h after injections (black arrows in Figure 3A). These results showed that a combination of KA and diclofenac has the highest effect on decreasing the Warburg phenotype of cancer cells in mentioned time points ( Figure 3C). To confirm this finding, we also performed a Seahorse experiment based on glycolytic rate assay that measured the maximum glycolytic capacity of cells by, first, shutting down mitochondrial respiration and then glycolysis (Methods). This assay also showed the highest inhibition of WE phenotype by the combination of KA and diclofenac in both MCF7 and MDA-MB-231 breast cancer cells ( Figure 3D).  Figure 4C shows that, in coculture conditions, the green signal that belongs to MCF7-GFP was increased in response to both drug and the combination, compared to the control group. These results show that it may be possible to control population dynamics in tumors based on their metabolic vulnerability.

Evolutionary Designed Therapy Controls Tumor Growth and Metastasis In Vivo
To translate our in vitro 2D and 3D results to in vivo we then performed animal experiments in female NSG mice. The experimental design was similar to the 3D spheroid experiment co-culture -MCF7 and MDA-MB-231 cells were mixed in 1:1 ratio and inoculated into the mice mammary fat poad. Mice were randomized into four groups of six mice per group and treated with diclofenac, KA, both, or none (DMSO). Primary tumor growth was unaffected by either monotherapy cf. controls, but there was a significant decrease caused by the combination of diclofenac and KA ( Figure 5A).
At the end of the experiment we extracted the tumors and all the vital organs from each mouse and scored for metastasis. Notably, the diclofenac + KA combination group was metastasis-free, while all other groups had measurable metastases across multiple sites ( Figure 5B). These results indicate that evolutionary therapeutic approaches based on metabolic targeting can be used to control tumor growth and even change aggressive phenotypes such as metastasis.

Reducing Warburg Phenotype Can Change Population Dynamics
Cancer evolution is a result of genetic diversification and clonal selection within the adaptive landscapes of tissue ecosystems. Therapeutic approaches can destroy some cancer clones, but it may also create a new space and different selective pressure resulting in expansion of resistant and most probably more aggressive clones aka "competitive release" in ecology and evolution. Evolutionary based therapeutic design can be used to control cancers when a cure is not possible. For that reason, understanding Darwinian intratumoral dynamics and their interactions with microenvironmental selection forces is critical to steer a tumor into a less invasive phenotype.
In 3D spheroid experiments we investigated the competitive release strategy to develop a new treatment design. Two distinct cell types were used to provide insight into the dynamic interactions among tumor cell subpopulations, a crucial element of intratumor heterogeneity. Mixed-culture competition experiments in spheroids can successfully track competitive outcomes amongst different cell types [34]. In 3D spheroid experiments differences in growth rates, carrying capacities and competition coefficients can be determined accurately [34]. For this experiment, we mono-and co-cultured MCF7-GFP and MDA-MB-231-RFP cells (1:1 ratio) and treated them with DMSO (control), diclofenac, KA, or combination of KA and diclofenac ( Figure 4A-C). As a control experiment we grew spheroids in different microenvironmental conditions and with different ratios of MCF7 to MDA-MB-231 cells. In any microenvironmental condition and with any ratio of the two cancer cells, MDA-MB-231 cells always eventually dominated the MCF7 cells and took over the whole population ( Figure S5). After treating the spheroids with diclofenac, KA, or the combination, we showed that these metabolically designed drugs can steer the population dynamics towards the less aggressive phenotype cells (MCF7-GFP). Figure 4A,B are controls showing effect of the agents on monocultures. Figure 4C shows that, in coculture conditions, the green signal that belongs to MCF7-GFP was increased in response to both drug and the combination, compared to the control group. These results show that it may be possible to control population dynamics in tumors based on their metabolic vulnerability.

Evolutionary Designed Therapy Controls Tumor Growth and Metastasis In Vivo
To translate our in vitro 2D and 3D results to in vivo we then performed animal experiments in female NSG mice. The experimental design was similar to the 3D spheroid experiment co-culture -MCF7 and MDA-MB-231 cells were mixed in 1:1 ratio and inoculated into the mice mammary fat poad. Mice were randomized into four groups of six mice per group and treated with diclofenac, KA, both, or none (DMSO). Primary tumor growth was unaffected by either monotherapy cf. controls, but there was a significant decrease caused by the combination of diclofenac and KA ( Figure 5A).   At the end of the experiment we extracted the tumors and all the vital organs from each mouse and scored for metastasis. Notably, the diclofenac + KA combination group was metastasis-free, while all other groups had measurable metastases across multiple sites ( Figure 5B). These results indicate that evolutionary therapeutic approaches based on metabolic targeting can be used to control tumor growth and even change aggressive phenotypes such as metastasis.

Discussion
Tumor evolution follows the Darwinian principles; meaning nature selects for phenotype not genotype [5]. When cancer initiates and develops, many clones and sub-clones emerge that are genetically heterogeneous with many of them presenting the same phenotype such as the Warburg effect (WE) phenotype [35]. The emergence of one phenotype from many clones with different genetic backgrounds implies the fitness of the selected phenotype to the selection forces. The WE phenotype is associated with progression and invasion of many cancers and is defined by a high glycolytic rate in the absence or presence of oxygen (aerobic glycolysis). Most cancer cells reprogram their metabolism in favor of aerobic glycolysis despite the presence of plentiful oxygen in their microenvironment, implying higher fitness of WE phenotype.
During the evolution of tumor cells and through adaptation to constantly variable microenvironment tumor cells adapt and acquire phenotypes helping them survive and grow. These acquired characteristics can be used as a new vulnerability against them. As mentioned previously, we recently showed that certain targets for inhibition of glycolysis have increased effects of inhibition at a lower intracellular pH [6]. To test this principle, we coupled Diclofenac, a compound we demonstrated to be a potent inhibitor of MCT's, with Koningic Acid, a known inhibitor of GAPDH. The MCT inhibition by Diclofenac provided the decrease in phi needed to enhance the effectiveness of KA. In this paper we showed that targeting the WE phenotype in more aggressive cancer cells can take their competitiveness from them so the less malignant cells can compete with them in a tumor niche. We used a simple model of co-culture of two cell lines with extremely opposing characteristics: MDA-MB-231 cells that are triple-negative fast proliferating, tumorigenic and metastatic in mouse xenograft models that are also highly glycolytic with WE phenotype and MCF7 cells that are ER/PR positive, slow-growing and tumorigenic in mouse xenografts with the help of estrogen pellets but not metastatic and very low glycolytic and not WE phenotype. The co-culture of these two cell lines was used as a simple population model of tumor heterogeneity to monitor the population dynamics of the tumor. In the future, these findings will be confirmed and expanded in a more complex, transgenic tumor model, to observe the effect of real intratumor heterogeneity in cell subpopulation. Using our metabolic targeted strategy, we showed that it is possible to control the population dynamic and heterogeneity. This population dynamic adjustment was able to successfully control tumor growth and metastasis in the mouse xenograft model. Therefore, we propose a novel strategy of cancer control for the tumor where a cure is not an option. We propose evolutionary principles of tumor growth that can be used to control tumor cells' clonal or sub-clonal population in favor of slower growth and less damage to patients.

Hypoxia Cell Culture
For hypoxia conditions we used a Biospherix X-Vivo Hypoxia Chamber. The conditions in each hypoxia chamber kept at 37 • C, 5% CO 2 , with 0.1% O 2 and 94.9% N 2 or 1.0% O 2 and 95% N 2 . The chamber was quality controlled and calibrated according to the manufacturer's specifications (Biospherix, Parish, NY, USA).

Mettler Toledo Five Easy TM /FiveGo TM pH Meter
The pH meter was first calibrated by using standard pre-made calibration buffer solutions for pH 4, 7, and 10. Then using the pH probe the pH of the control media used to grow the spheroids was measured. This pH was used as a baseline to compare the pH of the pHe media of the spheroids.

CCK8
Cell viability was measured after treatment with drugs (diclofenac and koningic acid or both) and DMSO as control using Cell Counting Kit-8 (CCK-8) under different microenvironmental conditions such as normoxia, hypoxia, normal pH, or low pH. CCK8 is a sensitive colorimetric-based viability assay based on Dojindo's highly water-soluble tetrazolium salt, with WST-8 as its active agent. CCK8 was used to measure viability as it is not pH sensitive and can be added to the cells directly in their niche, without fixation or change of media. For measuring viability, cells were seeded in a 96-well plate (with triplicate of the same samples), and viability was measured at the indicated intervals. The experiments were repeated at least two times.

CellTiter-Glo ®
This 3D Cell Viability Assay is designed to determine the number of viable cells in 3D cell culture based on quantitation of the ATP through based on Luminescent Cell Viability Assay chemistry. The lytic capability of the reagents is much higher than the 2D assay and the assay is compatible with 96 well-plate formats of our spheroid experiments that make it ideal for high-throughput screening. In our assay 72 h after treatment of spheroid with inhibitors, the lysis buffer directly added to each well containing one sphere of MDA-MB-231 or MCF7 cells. The plates were incubated in RT for 10 min on a rocking shaker. Then the assay mixture containing luciferin and luciferase was added and the luminescent was measured after 5 min incubation in RT.

Glycolytic Rate Measurements (Seahorse)
Glycolytic rate of MCF7 and MDA-MB-231 cancer cells was measured using Seahorse XF96 extracellular flux analyzer and a glycolysis rate kit (Seahorse Biosciences, Billerica, MA, USA). All the seahorse experiment has been performed in the absence of CO 2 /HCO 3− . Oxygen consumption rate (OCR) and extracellular acidification rate (ECAR) of cancer cells were determined by seeding them on XF96 microplates in their growth medium until they reached over 90% confluence. In these studies, seeding started with 20,000 cells (80% of the well area). Measurements were determined 24 h later when the cells reached 90% confluence. One hour before the seahorse measurements culture media were removed and cells were washed 3 times with PBS followed by the addition of base medium (nonbuffered DMEM supplemented with 25 mM glucose). Finally, data were normalized for total protein content of each well using the Bradford protein assay (ThermoFisher, Waltham, MA, USA). Seahorse measurements were performed with 4-6 technical replicates and these experiments were repeated at least 2 times. The WE phenotype ("Warburgness") can be expressed as the ratio of glycolysis (ECAR) to oxidative phosphorylation (expressed as the oxygen consumption rate, OCR) from the GST.

Metabolic Profiling
Cells were seeded in a regular 96-well plate for 2D culture and in U bottom 96 well plates for 3D spheroid culture in their growth medium containing 10% FBS under standard culture condition. Once cells reach 90% confluence for 2D or caring capacity for 3D (usually when the growth of spheres is steady and they don't get bigger), the growth media were removed and placed in a new 96 well plate with the same order as the original plate for metabolic profiling using YSI machine (YSI 2900 multi-analyte system (YSI, Yellow Springs, OH, USA)). We used at least 150 µL of media in each well to make sure there was enough media for each assay in the machine. Sensors in YSI were changed weekly and calibrated before each experiment. Final data of the lactate, glucose, glutamine, and glutamate in conditioned media were normalized to the protein amount per well or the confluency of each well imaged right before metabolite measurements.

MCT Activity
Cells were plated on four chamber Lab-Tek chambered cover glass slides 1 day prior to experimenting, and cultured in RPMI media containing 10% FBS and 1% P/S. On day of experimentation, media was replaced with RPMI containing 10% FBS, 1% P/S, 25 mM HEPES, and 25 mM PIPES, to stabilize pH while outside of the 5% CO 2 incubator. Once the assay was ready to begin for a given group, 5-(and-6)-Carboxy SNARF TM -1, Acetoxymethyl Ester, Acetate was added 8 min prior to experimentation at a concentration of 10 µg/mL. Prior to experimentation, the pH calibration curve of the 5-(and-6)-Carboxy SNARF TM -1, Acetoxymethyl Ester, Acetate dye was conducted using the nigericin method as described previously [36]. Once the 8 min incubation period was up, imaging began and superfusion with the 30 mM lactate solution (solutions listed in extended data file) and the SNARF-1 ratio was allowed to stabilize. The superfusion setup used in our experiments consists of a peristaltic to flow the media, a custom manufactured fluidic switch system which allows for the rapid changing of solutions, and a vacuum pump to aspirate media from the opposite side of the chamber. Once stabilized, cells were instantaneously switched to the 0 mM lactate solution and maintained in such until the velocity of the SNARF-1 ratio began to level off. Cells were then switched back to the 30 mM lactate condition and allowed to stabilize again. Once stabilized, the cells were switched to the same 30 mM lactate solution, this time either containing (X concentration of Diclofenac or Aspirin). Cells were maintained in the drug treatment solutions for 20 min in order for the drugs to take effect, and then switched to the 0 mM lactate solution containing the same concentration of drug they were treated with.

pHi Measurement
Cells were plated on four chamber Lab-Tek chambered cover glass slides 1 day prior to experimentation and cultured in RPMI media containing 10% FBS and 1% P/S. On day of experimentation, media was replaced with RPMI containing 10% FBS, 1% P/S, 25 mM HEPES, and 25 mM PIPES, to stabilize pH while outside of the 5% CO 2 incubator. Once the assay was ready to begin for a given group, 5-(and-6)-Carboxy SNARF TM -1, Acetoxymethyl Ester, Acetate was added 10 min prior to experimentation at a concentration of 10 µg/mL. After the 10 min, cells were washed 2 times with drug containing media and fresh drug containing media was resupplied to the cells. Change in pHi caused by each drug was measured by converting the SNARF-1 ratio to pHi and comparing pre-and post-treatment pHi's.

Spheroid Mono-and Co-Culture
Perfecta3 ® 96-Well Hanging Drop Plates or non-adhesive U shape bottom 96 well plates were used to grow the primary spheres containing 10,000 total cells. 10% Matrigel was used for the primary spheroid construction as follows. After cells were counted they were diluted as 200 cells per µL of media. Then the cells mixture was cooled down in ice for 10 min and 10% Matrigel (melted in ice in the 4-degree Celsius walk-in fridge overnight) was added and directly seeded in the U plates. The plates were centrifuged for 1 min 500 RPM. For each experimental condition, the ratio between MCF-7-GFP and MDA-MB-231-RFP was 0, 50, or 100 percent as follows: 0/100, 50/50, 100/0. Each ratio had 4-8 replicates for each experiment. An Incucyte microscope kept at 37 • C and 5% CO 2 was used to image the spheroid growth every 6 h over approximately 30 days. Images and fluorescent intensity were analyzed in Incucyte built-in software and image J. Relative fluorescent units (RFUs) were normalized by averaging red/green fluorescence when the corresponding cell type was absent then subtracting this value from other fluorescent values of that cell type. Drug treatments were directly added to the growth media of the spheroids and renewed every 3 days.

Animal Experiments
MCF7 and MDA-MB-231 breast cancer cell lines were grown in T-225 flask and harvested at 70-80% confluency. Both cells were authenticated by short tandem repeat analysis and tested for mycoplasma. 10 million cells in 200 µL cold PBS and 10% Matrigel was inoculated into cleared mammary fat pads of SCID mice (eight-to ten-week old females). Tumor measurements were done by calipers every other day and ultrasound once a week. Mice had free access to water and food for the whole duration of the experiment.
For mixed cultures, a 1:1 mixture of 5 million MDA-MB-231 and 5 million MCF7 cells were injected. One week before cell injection, an estrogen pellet (0.72 mg slow-release, Innovative Research of America) was implanted to allow for the growth of ER-positive MCF7 tumors. The concentration of diclofenac given to each mouse was 40 mg/kg. As for the KA the concentration was 1 mg/kg. These concentrations remained the same when the mice were treated with both drugs. The control mice were injected with DMSO in PBS. The mice were treated by intraperitoneal injections once on Monday, Wednesday, and Friday weekly. At the end of the experiment tumors were extracted, then the size and weight of them were measured. Vital organs were also extracted and looked for metastasis.
Animal studies were performed by the guidelines of the IACUC of the H. Lee Moffitt Cancer Center (that was approved by the University of South Florida IACUC committee: IACUC 5331R).

Conclusions
We propose a novel strategy of cancer control for the tumor when a cure is not an option. We propose evolutionary principles of tumor growth can be used to control tumor cells' clonal or sub-clonal population in favor of slower growth and less damage to patients. We showed that slight metabolic adaptations acquired during the evolution of cancer cells such as increased pHi can be used as vulnerability against them.
Supplementary Materials: The following are available online at https://www.mdpi.com/2072 -6694/13/1/64/s1, Figure S1: Glutamine and glutamate concentration in conditioned media of MCF7 and MDA-MB-231 cells treated with diclofenac, Koningic acid, and both compared to nontreated control in hypoxia (0.1% oxygen) and normoxia, Figure S2: ICC Analysis of MCT1 and MCT4 expression in breast cancer cell lines. Figure S3: Extracellualr pH measurements of cells exposed to various treatment conditions. Treatment with KA, DIC, or combination of both results in an increase in extracellular pH. Figure S4: Survival assay of breast cancer cell lines in different microenvironmental conditions. Targeting MCTs and GAPDH under different microenvironments reduces the cancer cells' viability. The effect is the most with the combination of both drug and under the most unique condition of solid tumors, hypoxia, and acidic pH. Data is normalized to non-treated for each cell line in each condition and presented as mean with SD as error bars. Each experiment has 4 replicates, Figure S5: Mono-and co-culture of fluorescently tagged breast cancer cell line. In different microenvironmental conditions and with any initial ration of MCF7: MDA-MB-231, the more aggressive cells, MDA-MB-231 always win. The spheroids are grown in normoxia but considering the size of spheres that reaches to 1-2 mm we know that oxygen concentration can vary from normoxia at the surface of spheres to hypoxia in the center.

Data Availability Statement:
All data presented in this study and not included within the supplementary materials are available by request submitted to the corresponding author.

Conflicts of Interest:
The authors declare no competing interest.