Parts of this article's content have been modified or rectified in:
Erratum: Glutathione levels are associated with methotrexate resistance in acute lymphoblastic leukemia cell lines
Rafael Renatino Canevarolo1†
Carolina Pereira de Souza Melo1Nathalia Moreno Cury1
Leonardo Luiz Artico1
Juliana Ronchi Corrêa1Yanca Tonhasca Lau1Samara Sousa Mariano1
Praneeth Reddy Sudalagunta2
Silvia Regina Brandalise1Ana Carolina de Mattos Zeri3
José Andrés Yunes1,4*- 1Centro de Pesquisa Boldrini, Centro Infantil Boldrini, Campinas, SP, Brazil
- 2Department of Cancer Physiology, H. Lee Moffitt Cancer Center & Research Institute, Tampa, FL, United States
- 3Brazilian Biosciences National Laboratory (LNBio), Brazilian Center for Research in Energy and Materials (CNPEM), Campinas, SP, Brazil
- 4Medical Genetics Department, Faculty of Medical Sciences, State University of Campinas, Campinas, SP, Brazil
Introduction: Methotrexate (MTX), a folic acid antagonist and nucleotide synthesis inhibitor, is a cornerstone drug used against acute lymphoblastic leukemia (ALL), but its mechanism of action and resistance continues to be unraveled even after decades of clinical use.
Methods: To better understand the mechanisms of this drug, we accessed the intracellular metabolic content of 13 ALL cell lines treated with MTX by 1H-NMR, and correlated metabolome data with cell proliferation and gene expression. Further, we validated these findings by inhibiting the cellular antioxidant system of the cells in vitro and in vivo in the presence of MTX.
Results: MTX altered the concentration of 31 out of 70 metabolites analyzed, suggesting inhibition of the glycine cleavage system, the pentose phosphate pathway, purine and pyrimidine synthesis, phospholipid metabolism, and bile acid uptake. We found that glutathione (GSH) levels were associated with MTX resistance in both treated and untreated cells, suggesting a new constitutive metabolic-based mechanism of resistance to the drug. Gene expression analyses showed that eight genes involved in GSH metabolism were correlated to GSH concentrations, 2 of which (gamma-glutamyltransferase 1 [GGT1] and thioredoxin reductase 3 [TXNRD3]) were also correlated to MTX resistance. Gene set enrichment analysis (GSEA) confirmed the association between GSH metabolism and MTX resistance. Pharmacological inhibition or stimulation of the main antioxidant systems of the cell, GSH and thioredoxin, confirmed their importance in MTX resistance. Arsenic trioxide (ATO), a thioredoxin inhibitor used against acute promyelocytic leukemia, potentiated MTX cytotoxicity in vitro in some of the ALL cell lines tested. Likewise, the ATO+MTX combination decreased tumor burden and extended the survival of NOD scid gamma (NSG) mice transplanted with patient-derived ALL xenograft, but only in one of four ALLs tested.
Conclusion: Altogether, our results show that the cellular antioxidant defense systems contribute to leukemia resistance to MTX, and targeting these pathways, especially the thioredoxin antioxidant system, may be a promising strategy for resensitizing ALL to MTX.
Introduction
Acute lymphoblastic leukemia (ALL) is a subtype of leukemia that affects B- or T-cell precursor cells. It is the most common type of cancer in children, accounting for one quarter of all childhood cancer cases (1). Although the cure rates of childhood ALL have reached 85% in developed countries (2, 3), drug resistance is considered the major cause of treatment failure; as so, it is an issue of continuous and intense research (4–7).
Methotrexate (MTX), a folic acid competitor (antifolate) first synthesized in 1947, has been widely used as a chemotherapeutic agent and immune system suppressor. MTX is a cornerstone drug in ALL therapy. Once internalized by the cell through the folate receptor (RFC1), MTX undergoes sequential glutamylation by folyl-polyglutamate synthase (FPGS), thus becoming MTX polyglutamates (MTX-PG2-6). The molecule with a 2-6 glutamic residues tail is the active version of the drug; although its chemical activity is similar to the original drug, MTX-PG2-6 is preferentially kept in the intracellular environment [whereas non-glutamylated MTX is quickly exported by ABC transporters (8)], and cellular sensitivity to MTX is directly related to the intracellular pool of MTX-PGs (9). Recently, it was shown that ALL subtypes are associated with differential accumulation of MTX-PG levels: ALL harboring TCF3-PBX1 or ETV6-RUNX1 fusions presented lower MTX-PG levels, whereas hyperdiploid and BCR-ABL-like ALL had higher MTX-PG levels than other subtypes (10).
MTX primarily acts by inhibiting dihydrofolate reductase (DHFR), which catalyzes the sequential conversion of folic acid to dihydrofolate (DHF) and tetrahydrofolate (THF). Once folic acid is a precursor in de novo synthesis of thymidine, MTX imbalances nucleotide pools, ultimately leading to incorrect DNA synthesis and halting cell division. A high percentage of cells in S-phase is associated with a better clinical response to MTX (11) and is accompanied by increased expression of genes involved in nucleotide biosynthesis, such as TYMS, CPTS, and an MTX target enzyme, DHFR (11). MTX also inhibits the bifunctional purine biosynthesis protein PURH (ATIC), blocking the conversion of the aminoimidazole carboxamide ribonucleotide (AICAR) into inosine monophosphate (IMP), thus inhibiting purine synthesis and the production of adenosine. AICAR accumulation leads to decreased one-carbon metabolism gene expression, cell proliferation, and increased global bioenergetic capacity via AMPK activation (12). Alternatively, MTX’s mechanism of resistance has been attributed to decreased MTX uptake due to deficiencies in membrane transport (13–15); increased expression of MTX efflux transporters (16–18); deficient polyglutamylation of MTX by FPGS (18–20); altered levels or structures of target enzymes (21, 22); and increased hydrolysis of MTXPGs by lysosomal gamma-glutamyl hydrolases (GGH) (23–25), although the influence of GGH on MTX efficacy remains controversial (26–28).
In parallel with MTX’s most well-known targets, several NAD+-dependent enzymes were inhibited by the antifolate in vitro: 2-oxoglutarate, isocitrate, malate, pyruvate, succinate, 6-phosphogluconate and glucose-6-phosphate dehydrogenases, glutamate-cysteine ligase, glutathione reductase, and glutathione peroxidase (29). MTX was shown to inhibit methionine S-adenosyltransferase (MAT) gene expression, protein levels amounts and activity in vitro and in vivo, resulting in decreased levels of S-adenosylmethionine (SAM) (30), the main methyl donor of metabolism.
Despite the number and importance of such enzymes to cellular metabolism, few studies have been conducted to better understand the effects of MTX on cellular metabolic homeostasis or find metabolic fragilities that could be exploited clinically. For instance, the depletion of multiple genes belonging to the histidine catabolism pathway increases resistance to MTX (31). This is caused by histidine catabolism draining the cellular pool of THF, which is already low in MTX-treated cells. The same study showed that in vivo dietary supplementation of histidine increased THF consumption by the histidine degradation pathway and enhanced the sensitivity of leukemia xenografts to MTX.
The antimetabolite nature of MTX and its inhibitory effect on many important metabolic enzymes prompted us to study the consequences of MTX treatment on cellular metabolism. We aimed to discover novel metabolic vulnerabilities that could be exploited therapeutically. Using a panel of 13 ALL cell lines (6 B-cell precursors and 7 T-cell ALL precursors) and untargeted hydrogen nuclear magnetic resonance (1H-NMR)-based metabolomics, we characterized MTX-triggered metabolic disturbance in detail. Though short term (48 hours) resistance to MTX seemed to be directly related to cellular proliferation, resistance in longer exposure times (96 hours) was associated with glutathione (GSH) levels. Gene expression analysis and pharmacological modulation confirmed the relevance of the cellular antioxidant system to MTX-mediated cell death susceptibility. Finally, we showed that pharmacological inhibition of the thioredoxin system by the chemotherapeutic drug arsenic trioxide (ATO) enhanced MTX cytotoxicity in vivo against one out of four primary human ALL cells, suggesting a novel therapeutic rationale to be further explored against ALL.
Material and Methods
Cell lines
Six B-cell precursor (Nalm6, Nalm16, Nalm30, REH, RS4;11, and 697) and 7 T-ALL (ALL-SIL, CCRF-CEM, HPB-ALL, Jurkat, Molt-4, P12-ICHIKAWA and TALL-1) cell lines were used in the study. The cell lines REH, 697, and RS4;11 were kindly provided by Dr. Sheila A. Shurtleff from St. Jude Children’s Research Hospital, Memphis, TN, USA; the cell lines Nalm16 and Nalm30 were provided by Dr. Akira Harashima from Hayashibara Biochemical Lab, Fujisaki, Japan; the cell lines ALL-SIL, HPB-ALL, P12-ICHIKAWA, and TALL-1 were kindly provided by Dr. João Barata from the Instituto de Medicina Molecular, Faculdade de Medicina da Universidade de Lisboa, Lisbon, Portugal; and the cell lines Jurkat, Nalm-6, CCRF-CEM (CEM) and Molt-4 were kindly provided by Dr. Alexandre E. Nowill from the State University of Campinas, Campinas, Brazil. Cell lines were cultured in RPMI-1640 culture medium (Cultilab, Campinas, Brazil); supplemented with 10% fetal bovine serum (Cultilab), 100 IU/ml of penicillin, and 100 μg/ml of streptomycin (Sigma-Aldrich, Saint Louis, MO, USA); and maintained at 37°C with a 5% CO2 atmosphere in all experiments.
Reagents
All compounds were purchased from Sigma-Aldrich, resuspended according to the manufacturers’ instructions, and stored at –20°C. Subsequent dilutions used in the experiments were prepared at the moment of use.
Cell viability assays
Eighty microliters of cell suspension (4 × 105 cells/ml) were seeded in a 96-well cell culture plate, followed by the addition of 20 μl of drug solution in culture medium. Each dose was tested in triplicate. Negative controls received only vehicle, whereas positive (death) controls were added with 2 μl of an 18% paraformaldehyde solution immediately before viability determination. After the treatment period, the culture medium of each well was replaced by 0.2 ml of calcein AM solution (Sigma-Aldrich) at 2 μM in phosphate buffer [PBS 1X] and the culture plates were put into the incubator for 30 minutes before fluorescence reading using a Synergy H1 microplate reader (Biotek, Winooski, VT, USA; excitation/emission: 492/518 nm). To account for heterogeneous cell distribution within each well, fluorescence was measured in 25 points (a 5 × 5 matrix layout), with values then integrated (summed) at the end. Alternatively, the MTT (thiazolyl blue tetrazolium bromide; Sigma-Aldrich) reduction test was also used. In those cases, after 4 hours of incubation, 20 μl of MTT (5 mg/ml) were added to each well for cell metabolization. After overnight incubation, 0.1 ml of a dodecyl sulfate solution (10%) containing HCl (10 mM) was added to the wells. Absorbance was measured at 570 nm. Cell viability was expressed in relation to controls. Given that almost identical IC50 values for MTX were obtained by MTT and calcein AM at 48 and 96 hours, either one or the other method was used to access cell viability, depending on experimental and laboratory logistics.
Determination of the doubling time
Two hundred microliters of a 2.5 × 105 cell/ml cell suspension were seeded in a 96-well culture plate in triplicates. Fifteen microliters of the cell suspension were collected daily for cell count under the microscope (exclusion by trypan blue). The doubling time of each cell line was determined from the proliferation curves obtained.
Cell treatment and metabolites extraction for 1H-NMR
Two hundred million cells in the exponential growth phase were suspended in culture medium (200 ml) and equally divided into 2 cell culture flasks. After 12 hours of acclimatization, either MTX (25 nM) or vehicle were added to the flasks. The doses were representative of the average MTX IC50 (96 hours) across all cell lines (18.4 nM). Following 24 hours of treatment, cells were washed with phosphate buffered saline (PBS 1X) and stored at –80°C. Both conditions (treatment and control) were analyzed in 3 biological triplicates for each cell line. The intracellular polar metabolic content was obtained following the methanol/chloroform extraction protocol, adapted from La Belle (32). Briefly, cells were disrupted (Sonics Vibra Cell, Sonics & Materials Inc., Newtown, CT, USA) in an ice-cold mixture of methanol (1.66 ml) and chloroform (0.83 ml) for 3 minutes. Next, ice-cold Milli Q water (1.25 ml) and chloroform (1.25 ml) were added to the sonicated mixture, followed by rapid vortexing. After centrifugation (3.4 × 103 g for 20 minutes at 4°C), 2.6 ml of the supernatant phase (containing water, methanol, and polar metabolites) were collected and freeze-dried.
Acquisition of 1H-NMR spectra and metabolic profiling
The lyophilized metabolic extracts were suspended in 0.6 ml of deuterated water (Cambridge Isotope Laboratories Inc., Tewksbury, MA, USA) containing phosphate buffer (100 mM) and TSP (3-3-(trimethylsilyl)-2,2,3,3-tetradeuteropropionic acid or TMSP-d4, Cambridge Isotope Laboratories) (0.5 mM) and inserted into a 5 mm NMR tube for immediate acquisition. Spectra were acquired on a Varian Inova NMR spectrometer (Agilent Technologies Inc., Santa Clara, CA, USA), equipped with a triple resonance probe and operating at a resonant frequency of 1H at 500 MHz and a constant temperature of 25°C. Two hundred and fifty-six scans were performed with delays of 1.5 seconds, a reading window of 16 ppm, and an acquisition time of 4 seconds. Phase and baseline corrections, as well as metabolic identification and quantification, were performed in Chenomx NMR Suite 7.1 (Chenomx Inc., Edmonton, Canada). In total, 78 1H-NMR spectra were analyzed. Metabolite set enrichment analyses were performed with the platform MetaboAnalyst (33) (http://www.metaboanalyst.ca).
Gene expression analysis
Cell lines (5 to 10 × 106 cells) in the exponential growth phase were maintained overnight in fresh culture medium and then had their RNA extracted using the Illustra RNAspin Mini Kit (GE Healthcare Life Sciences, Pittsburgh, PA, USA). Samples were processed with the One-Cycle Target Labeling and Control Reagents Kit (Affymetrix, Santa Clara, CA, USA) and hybridized on HG-U133 Plus 2.0 Arrays (Affymetrix). Expression values were obtained with the iterPLIER+16 algorithm and expressed in a log2 scale. For additional details, please see Silveira et al. (34). Gene set enrichment analyses (GSEA) (http://www.broadinstitute.org/gsea/) were performed with 1,000 gene set permutation. Phenotypic labels (eg, GSH levels or MTX IC50) were used as continuous labels. The GSEA platform used Pearson’s correlation to determine the degree of linear relationship between gene sets and phenotypes. Only probe sets/transcript clusters annotated with a Gene Symbol were used in the analyses. Gene sets were considered significantly enriched when P ≤ 0.05 and FDR ≤ 0.05.
High-throughput in vitro drug and combination screening assay
A high-throughput in vitro drug screening assay was used to quantify the chemosensitivity of leukemia cell lines as described previously (35, 36). Briefly, leukemia cells were seeded in 384-well, black, optically clear flat-bottom PhenoPlate (PerkinElmer, Waltham, MA) with collagen I (PureCol, Advanced Biomatrix, Carlsbad, CA) to a total volume of 9 μl, containing approximately 3,000 leukemia cells in each well. After collagen polymerization (1 hour at 5% CO2 and 37°C), each well was filled with 73 μl of culture media and left overnight in the incubator for cell acclimatization. In the next day, 9 μl of the drug dilutions (BSO, Auranofin, BSO+Auranofin) were added to the wells so that every drug/combination was tested at 5 (fixed concentration ratio, for combinations) concentrations (1:3 serial dilution) in 2 replicates. Negative controls containing vehicle (DMSO, maximum concentration of 0.5%) were included. Plates were placed in a motorized stage microscope (EVOS Auto FL 2, Thermo Fisher Scientific, Waltham, MA) equipped with an on-stage incubator and maintained at 5% CO2 and 37°C. Each well was imaged every 30 minutes for a total duration of 96 hours. A digital image analysis algorithm (37) was implemented to determine changes in viability of each well longitudinally across the 96-hour interval. This algorithm computes differences in sequential brightfield images and identifies live cells with continuous membrane deformations resulting from their interaction with the surrounding extracellular matrix. These interactions cease upon cell death. By applying this operation to all images acquired for each well, we were able to quantify nondestructively the effect of the drugs and their combinations as a function of concentration and exposure time – ie, by area under the curve (AUC) and lethal dose for median effect (LD50). Cell line-specific doses were chosen based on the MTT reduction assay and on previous standardization experiments. Synergy was determined by using the method described by Sudalagunta et al (38), in which percent live cells across time and the serially diluted drug(s) concentrations are used to compute additive response using the Bliss Independence Model. The additive response serves as a reference to determine the extent of synergy observed in each cell line by comparing it with the percent live cells measured when treated with the combination (at a fixed ratio of the 2 constituent single agents). The metric used to quantify in vitro drug synergy is given by the log2 ratio of the additive LD50 (computed from the two single agent responses as described above) and the combination LD50 for each cell line. This in vitro drug sensitivity characterization and combination effect framework was developed and pioneered by Dr Ariosto S. Silva’s group (Moffitt Cancer Center), which resulted in several significant contributions in the field of hematological malignancies over the past five years (39–42).
Calcein AM for cell viability analysis by flow cytometry
Cells were seeded in 96 well plates at 3.0e5 cells/well in 100 µl of RPMI-1640 with 10% FBS. Methotrexate was added at different concentrations in triplicates. After 48 hours of incubation, a fixed dose of Arsenic Trioxide was added and 48 hours later cell viability was measured by the MTT method. To confirm the MTT results, a replicate plate run in parallel was developed by adding 5 µl of 2.5 µM Calcein-AM (ThermoFisher Scientific) to selected wells and incubated 15 min at room temperature. The number of viable cells was measured by flow cytometry after addition of 5 µl spike in fluorescently labeled bead (CountBright™ Plus Absolute Counting Beads, ThermoFisher Scientific) for normalization. The number of Calcein-positive cells per 1.0e3 beads was computed. Some of the cell lines that were more sensitive to the drug combination were assayed in 48-hour experiments instead of 96 hours. In this case, ATO was added after the first 24 hours of treatment with MTX and the number of viable cells was measured 24 hour later. In this case, cells were cultured in 24-well plates at 1.2e5 cells/well in 500 µl of RPMI-1640 with 10% FBS. At the end of the experiment (48 hours), 100 µl were transferred to a cytometer tube, 0.5 µl of 25 µM Calcein-AM were added and, after 15 min of incubation at room temperature, 5 µl of fluorescently labeled beads plus 195 µl of PBS 1X were added and cells were analyzed by flow cytometry. Flow cytometry was performed on a LSRFortessa X-20 (Becton Dickinson, Franklin Lakes, NJ) and data were analyzed using FlowJo (v7.1.6, TreeStar).
In vivo efficacy of ATO+MTX against primary ALL
Cryopreserved primary ALL cells obtained from a mouse of the IL7R<CPT> knock-in model (43) or from 4 pediatric ALL patients (2 BCP- and 2 T-ALL) were used in this study. Cells were thawed, washed in PBS 1X buffer, and injected into nonirradiated mice via the tail vein. Sick animals were killed and had their spleen macerated, and ALL cells were obtained by Ficoll gradient centrifugation. Fresh ALL cells were then transplanted into more mice for the experiments. The murine ALL was transplanted into C57BL/6 mice. Treatment started 2 days after transplantation. Mice (n = 12) were randomly distributed into 4 treatment groups: vehicle (PBS 1X with 1% NaOH), arsenic trioxide (ATO; 2.5 mg/kg), methotrexate (MTX; 1, 5, or 10 mg/kg), and ATO+MTX. Drugs were intraperitoneally administered once a day, 5× per week. When in combination, MTX was given 6 hours before ATO (based on Szymanska et al (44), who used this interval for other drug combinations). Overall survival was measured from after 2 days the date of leukemia injection to the date of animal death or when humane endpoint criteria were reached. Differences in survival were determined by Log-rank (Mantel-Cox) test using the GraphPad Prism v7.0 (GraphPad, La Jolla, CA) software. The study was approved by the Animal Experimentation Ethics Committees of University of Campinas (CEUA/UNIMCAP, protocol 4557-1(A)/2018). Experiments with the human ALLs were carried out in nonirradiated NSG (NOD.Cg-PrkdcscidIl2rgtm1Wjl/SzJ) mice (Jackson Laboratory). Leukemia progression was monitored weekly in peripheral blood post-Ficoll mononuclear cells by flow cytometry using anti-human CD45-APC (clone HI30, ImmunoTools) and anti-mouse CD45-BV605 (clone 30-F11, BD Horizon). Flow cytometry was performed on a LSRFortessa X-20 and data were analyzed using FlowJo. When hCD45+ cells reached 0.1% of total CD45+ peripheral blood cells in half of the mice (n = 24), animals were randomly distributed among the 4 treatment groups. T-ALL PDXs were treated with vehicle, ATO (2.5 mg/kg), MTX (5 mg/kg for BCP-ALL; 10 mg/kg for T-ALL) or ATO+MTX, as described above. Treatment lasted for 4 weeks and was administered in consecutive weeks for T-ALL xenografts or every other week for BCP-ALL (the BCP-ALL cells used were more sensitive to MTX than T-ALL). Death by acute toxicity was censored. Overall survival was measured from the date of leukemia engraftment to either the animal death, when the proportion of hCD45+ cells in the peripheral blood > 25% (45), or when humane endpoint criteria were reached. The study was approved by the Animal Experimentation Ethics Committees of Centro Infantil Boldrini (CEUA/BOLDRINI, protocol 0016-2020) and Human Ethics Committees of Centro Infantil Boldrini (CAAE 34601120.7.0000.5376).
Statistical analyses
Paired nonparametric Wilcoxon signed rank tests and other analyses and graphs (dose-response curves, estimates of IC50 values, correlation plots, Kaplan-Meier curves, and log-rank tests for survival analyses) were made in Prism v7.0 (GraphPad, La Jolla, CA, USA). P ≤ 0.05 was considered significant.
Results
MTX impact on cellular metabolism
Initially, cells were treated with increasing doses of MTX for 48 or 96 hours to determine their IC50. IC50 at 48 hours ranged from 10.8 to 85.5 nM in T-ALL and from 14.7 to 101.1 nM in BCP-ALL. IC50 at 96 hours ranged from 5.5 to 38.4 nM in T-ALL and from 9.1 to 31.7 n in BCP-ALL (Supplementary Figure 1).
To determine the impact of MTX on cell metabolism, cell lines underwent unsupervised global metabolic profiling by 1H-NMR in either the presence or in the absence of the drug (a representative 1H-NMR spectrum is shown in Supplementary Figure 2). All cell lines were treated with a fixed dose of MTX (25 nM) for 24 hours, as this dose was representative of MTX’s average IC50 at 96 hours across all cell lines (ie, 18.4 nM). By treating the cells with the same dose of MTX, we intended to find cell line-specific metabolic differences under the same drug challenge. We reasoned that using cell-specific IC50 doses of MTX would have overthrown any metabolic differences among cell lines. This approach is in line with previous metabolomic studies that have adopted a similar rationale (46, 47). Additionally, this condition elicited the greatest metabolic variance, compared to high-dose MTX and shorter treatment time, in a pilot study on a MTX-resistant vs. a MTX-sensitive cell line (Supplementary Figure 3).
In total, 70 metabolites were identified and quantified. Metabolites’ identifiers, concentrations, and fold changes across all samples can be found in Supplementary Tables 1–3, respectively. We found that 19 metabolites were significantly increased by MTX, whereas 12 were decreased in response to the drug treatment (Figure 1A and Supplementary Table 4). Several amino acids had their concentrations increased following MTX treatment: glycine, alanine, phenylalanine, threonine, histidine, valine, arginine, leucine, isoleucine, tyrosine, and proline. 2-oxobutyrate, a byproduct of the conversion of cystathionine into cysteine and first intermediate in threonine degradation, was also increased with drug treatment. MTX treatment produced effects consistent with alterations in the folate cycle (formate); pyrimidine (uridine and dCTP) and purine (AMP, ADP, and NAD+) nucleotide synthesis; the pentose phosphate pathway (UDP-glucuronate, UDP-glucose, and myo-inositol); phospholipid synthesis (phosphocholine, sn-glycero-3-phosphocholine, phosphoethanolamine, and trimethylamine); the citric acid cycle (2-oxoglutarate); and the uptake of bile acids (cholate and glycocholate).
FIGURE 1
Figure 1 MTX-driven metabolic perturbations in ALL. (A) Metabolic disturbance caused by MTX treatment. Each dot is the mean of three biological replicates for each cell line. Metabolites’ concentrations fold change (in MTX-treated/untreated cells) above 1 indicate increase caused by MTX, whereas values below 1 indicate antifolate-induced decrease (P-values, Wilcoxon signed ranked test, after comparison of the distributions with the hypothetical mean value of 1). (B) Enrichment score for metabolites in (A) indicate the metabolic pathways and processes most impacted by MTX. (C) Correlation between metabolite concentration in untreated cells and MTX resistance (IC50 at 96 hours). Black and white beads represent BCP-ALL and T-ALL cell lines, respectively. P-values and r are from Pearson correlation.
A metabolite set enrichment analysis made from the 31 metabolites significantly modulated by MTX confirmed the major impact of the drug on glycine and serine metabolism, the urea cycle, alanine metabolism, and arginine and proline metabolism (Figure 1B). A scheme integrating all MTX-modulated metabolites into a single metabolic panel can be found in Supplementary Figure 4. Interestingly, a few cases of differential metabolic changes (ie, fold change of MTX treatment over basal levels) were associated with MTX resistance (at 96 hours); for example, the greatest fold changes for alanine were observed among the most sensitive cells to MTX, whereas 2-oxoglutarate (α-ketoglutarate) greatest fold changes were found among the most resistant cells to the antifolate (Figure 1C). Although a causal relationship cannot be established, this finding indicates that MTX resistance in vitro is associated with specific metabolic responses by the cells under drug action.
Metabolic levels associated with MTX resistance
In addition to the association between metabolic fold change and MTX resistance, we searched for correlations between metabolic concentration and resistance to MTX (expressed as IC50) at 48 and 96 hours in either the presence or absence of the drug (Table 1). Given that MTX is more effective against dividing cells because of its effect on the nucleotide synthesis, we also examined correlations between metabolites concentrations and cell proliferation rate, as this could be a confounding factor (proliferation curves for doubling time determination are in Supplemental Figure 5). In fact, MTX resistance at 48 hours and proliferation rate were tightly correlated (Pearson r = 0.829, P < 0.001) (Figure 2A); however, the association was lost for drug resistance at 96 hours (Pearson r = 0.111, P = 0.718) (Figure 2B), and MTX resistance at 48 hours was not associated with resistance at 96 hours (Pearson r = 0.268, P = 0.399) (Figure 2C), indicating that factors unrelated to cell proliferation were contributing to cellular resistance to MTX during longer exposure times.
TABLE 1
Table 1 Correlation between the concentration of selected metabolites and cell proliferation (expressed in terms of doubling time, in days) or MTX resistance (at 48 or 96 hours) in cells either treated or not with MTX.
FIGURE 2
Figure 2 Correlation between cellular proliferation, MTX resistance and GSH levels. (A) Correlation between proliferation rate (cell doubling time) and MTX resistance (IC50) at 48 hours or (B) 96 hours, and (C) across MTX resistance over time (48 × 96 hours). (D) Correlation between constitutive GSH and MTX resistance at 48 hours or (E) 96 hours, and (F) between GSH levels in MTX-treated cells and MTX IC50 at 96 hours. Black and white beads represent BCP-ALL and T-ALL cell lines, respectively. P-values and r are from Pearson correlation.
Several metabolites were negatively associated with MTX resistance at 48 hours: acetate, carnitine, cholate, formate, hypoxanthine, malate, and succinate (Table 1). Acetate is the product of glycolysis, malate and succinate are citric acid cycle elements, carnitine is important for lipid beta-oxidation, and hypoxanthine is a purine derivative, suggesting that, at shorter exposure times (48 hours), MTX sensitivity is associated with increased cellular metabolic activity. Alternatively, MTX resistance at 96 hours was positively associated with asparagine, guanosine, and glutathione (GSH) levels (Table 1).
GSH is the main antioxidant of the cell and is important in preventing damage caused by reactive oxygen species (ROS). Interestingly, although not associated with MTX resistance at 48 hours (Pearson r = 0.287, P = 0.366) (Figure 2D), GSH was the only metabolite associated with MTX resistance in both untreated (Pearson r = 0.764, P = 0.003) (Figure 2E) and treated cells (Pearson r = 0.780, P = 0.003) (Figure 2F) at 96 hours, indicating a constitutive metabolic trait predictive of cellular resistance to the MTX insult. Its association with MTX resistance is in line with the fact that MTX is a ROS-generating drug (48, 49). Correlation coefficients between all metabolites concentrations, cell proliferation rate, and MTX resistance (at 48 and 96 hours) in treated and untreated samples are in Supplementary Table 5.
GSH levels associated with the expression of GSH metabolism genes
We then searched for genetic traits that could explain the metabolite variability among the cell lines. For this purpose, the concentration of asparagine, guanosine, and GSH was correlated to the expression of genes belonging to their respective metabolic pathways. Though no associations for asparagine or guanosine-related genes were found, 8 genes that code for enzymes that participate in the GSH metabolism were positively correlated to GSH levels: CBS, GCLC, GSTM3, GSTA4, GGT1, GGT1/GGT2/GGT3P, GGT7, and TXNRD3 (Figure 3A). The expression of all genes involved in asparagine, guanosine, and GSH metabolisms that did not correlate to their respective metabolite’s levels are in Supplementary Figure 6. Interestingly, 2 of the GSH metabolism genes that correlated with GSH levels (GGT1 and TXNRD3) were also associated with MTX resistance (Figure 3B). Moreover, a Gene Set Enrichment Analysis (GSEA) confirmed that the expression of KEGG’s “Glutathione Metabolism” gene set (Supplementary Table 6) was enriched in MTX-resistant samples (Figure 3C). These findings offer a genetic basis for the differential concentration of GSH across the samples with distinct MTX sensitivity.
FIGURE 3
Figure 3 GSH-related transcriptome associates with metabolic levels and MTX resistance. (A) Association between glutathione (GSH) metabolism genes’ expression and GSH concentration; (B) Gene set enrichment analysis score plot showing KEGG’s Glutathione Metabolism gene set enrichment in high-GSH ALL cell lines. (C) Expression of two GSH-related genes also correlated to MTX resistance. Black and white beads represent BCP-ALL and T-ALL cell lines, respectively. P-values and r are from Pearson correlation.
Promoting GSH metabolism induces proliferation and MTX resistance
Given that both metabolic and gene expression data pointed to GSH metabolism as an important factor contributing to MTX resistance in ALL cells, we evaluated the effect of either enhancing or inhibiting GSH metabolism on cellular resistance to MTX. Initially, we treated the cell lines with increasing doses of MTX in the presence of N-acetylcysteine (NAC), a precursor of the de novo synthesis of GSH that is converted into cysteine by glutamate-cysteine ligase (GCL). We observed that NAC not only protected cells from MTX cytotoxicity but also promoted cell proliferation, including in RS4;11, 697, ALL-SIL and P12-ICHIKAWA (Figure 4A) — 4 cell lines with relatively low constitutive GSH levels — whereas cell lines with relatively high GSH (eg, HPB-ALL and Jurkat) tended to present less pronounced effects. The graphs for all cell lines are in Supplementary Figure 7.
FIGURE 4
Figure 4 In vitro combination of GSH pathway modulators with MTX or a thioredoxin reductase inhibitor. (A) Survival curves to increasing doses of MTX in either the presence or absence of NAC in four representative cell lines, accessed by the calcein AM method after 96 hours of treatment. A.U. = arbitrary units. (B) BCP-ALL RS4;11 survival in response to auranofin (AURA), buthionine sulphoximine (BSO), BSO + AURA and vehicle (CONTROL) accessed by bright field imaging (high-throughput drug combination screening assay) over 96 hours. (C) Idem (B), for TALL-1 cell line. In the BSO + AURA plots, curve colors indicate the corresponding doses of the single agents that were combined.
Inhibiting GSH metabolism does not sensitize cells to MTX
We also treated leukemia cells with increasing doses of 4 interferents of GSH metabolism: phenethylisothiocyanate (PEITC) and piperlongumine (PL), which deplete GSH and increase ROS (50–52); hydrogen peroxide (H2O2), a ROS itself; and buthionine sulphoximine (BSO), a GCL inhibitor and, thus, a pharmacological inhibitor of GSH synthesis. The IC50 obtained for all these compounds across the cell lines are listed in Supplementary Table 7. Surprisingly, these ROS-promoting compounds did not enhance MTX cytotoxicity significantly (Supplementary Figure 8), even though we found a correlation between PL or H2O2 resistance and MTX resistance (Supplementary Figure 9). In sum, these results show that, while promoting GSH metabolism consistently increases leukemia resistance to MTX — sometimes even accompanied by simultaneous increase in proliferation — its inhibition does result in sensitization to the antifolate.
Inhibiting the thioredoxin system potentiates MTX cytotoxicity in vitro
The intriguing asymmetrical contribution of GSH metabolism on MTX resistance suggested that a compensatory mechanism could be counterbalancing GSH metabolism inhibition and, thus, protecting cells from MTX insult.
As showed previously, thioredoxin reductases 3 (TXNRD3) expression was correlated to both intracellular GSH levels and MTX resistance (Figures 3A, B). TXNRD3 is a member of the cellular antioxidant arsenal whose main role is to chemically reduce thioredoxins, antioxidant proteins that facilitate the reduction of oxidized client proteins by cysteine thiol-disulfide exchange. Oxidized thioredoxin reductases, in turn, are reduced back to their functional state by NADPH. GSH and thioredoxin, therefore, are the 2 known antioxidant systems whose co-occurrence is crucial for the maintenance of the cellular redox homeostasis. We used a high-throughput single agent and drug combination screening assay based on bright-field images to determine cell survival under treatment with auranofin (Aura), a clinically approved thioredoxin reductase inhibitor for rheumatoid arthritis, and the GSH synthesis blocker BSO (Figures 4B, C and Supplementary Figure 10). A synergistic interaction between the drugs was observed in 5 out 6 cell lines tested using the method described by Sudalagunta et al (38), in which percent live cells across time and the serially diluted drug(s) concentrations are used to compute additive response using the Bliss Independence model (Supplementary Table 8). We confirmed the co-dependence of ALL cells on these antioxidant systems in 3 cell lines by the MTT viability assay (Supplementary Figure 11).
Next, we investigated the role of the antioxidant thioredoxin system on MTX sensitivity in vitro. For this purpose, we tested the effect of MTX in combination with 2 thioredoxin reductase inhibitors, arsenic trioxide (53) (ATO), and auranofin. We observed a mild potentiation of the MTX effect by ATO (Figure 5A) and Auranofin (Supplementary Figure 12) but only for few of the cell lines tested.
FIGURE 5
Figure 5 Combination of ATO and MTX in vitro and in PDX models of ALL. (A) Cell survival accessed by flow cytometry analysis of calcein-positive cells. The number of viable calcein+ cells per 1,000 fluorescent beads (spike in) was determined after 48 or 96 hours, as indicated. MTX was added at time zero while ATO was added at half-time. Results are the mean ± standard error (SE) of triplicate samples. P-values for Tukey’s Test following ANOVA. *P<0.05; ***P<0.001. Doses used: CCRF-CEM: ATO = 1.5 µM, MTX = 12.5 nM; P12-ICHIKAWA: ATO = 3 µM, MTX = 13 nM; Jurkat: ATO = 1.5 µM, MTX = 36 nM; RS4;11: ATO = 1.4 µM, MTX = 18 nM; 697: ATO = 0.75 µM, MTX = 6 nM. (B) Human CD45-positive (hCD45+) leukemia cells in peripheral blood of PDX-transplanted NSG mice, accessed by flow cytometry. Treatment started one week after engraftment when hCD45+ cells > 0.5% in half of the mice (Day 7) and was administered for 4 consecutive weeks against T-ALL or every other week against BCP-ALL. Doses used: MTX = 10 mg/kg daily (T-ALL) or 5 mg/kg daily (BCP-ALL); ATO = 2.5 mg/kg daily; ATO was administered 6 hours after MTX in the co-treatment cohort. L.E. = leukemia engraftment period. P-value for Log-rank test. P-value for unpaired t test between MTX and ATO+MTX. *P≤0.05; **P≤0.01; ****P≤0.0001; ns, not significant.
Altogether, our results indicate that the thioredoxin system exerts a compensatory effect over the GSH pathway that is not reciprocal and suggest that the thioredoxin system might be a potential target for MTX re-sensitization in leukemia cells.
In vivo antileukemic effects of ATO and MTX in murine and patient-derived leukemia
Although mild, the antileukemic effects of the combinations between thioredoxin reductase inhibitors and MTX in vitro motivated us to test them in vivo. Here we chose the ATO+MTX combination given ATO’s current use in the clinic as a chemotherapeutic agent against acute promyelocytic leukemia. Our first investigation was using an immunocompetent mouse model transplanted with a primary murine B-ALL obtained from our Il7r knock-in mouse model (43). We tested different doses of MTX (1, 5, and 10 mg/kg) and found a dose-response pattern, with 10 mg/kg showing best results. ATO+MTX significantly increased overall survival compared to vehicle or single agents, though the difference was modest (Supplementary Figure 13).
Next, we tested the efficacy of the drug combination in 4 patient-derived xenografts (PDXs; 2 BCP-ALL and 2 T-ALL) in vivo. The PDXs represented a diverse genetic background: T-ALL 1 carries an activating IL7R mutation that confers a particular aggressive phenotype (54), and BCP-ALL 1 expresses the TCF3-PBX1 fusion protein as a result of the t(1;19) chromosomal translocation. Mice transplanted with T-ALL PDXs were treated with the same dose (10 mg/kg) and schedule as mice transplanted with murine ALL cells, whereas BCP-ALL PDX-transplanted mice received MTX at a lower dose (5 mg/kg) and were treated every other week to minimize the contribution of MTX to disease control and thus maximize ATO’s effect. However, this approach resulted in longer treatment periods and increased gastrointestinal toxicity with some fatalities that were censored in the analysis (Supplementary Table 9). Of the 4 PDXs tested, only one T-ALL was resistant to MTX (remission lasted for 20 days compared to ≥ 50 days in the other 3 PDXs). In this case, the ATO+MTX combination significantly decreased tumor burden in the peripheral blood (Figure 5B) and extended overall survival (Supplementary Figure 14) compared to vehicle alone or either drug used as a monotherapy. Altogether, these results suggest that MTX-resistant ALL patients could benefit from therapies that include thioredoxin reductase inhibitors in combination with MTX.
Discussion
Seventy years after MTX was first synthesized, our comprehension of its mechanism of action and impact on cellular metabolism continues to expand. In this study, we explored the metabolic-wide changes caused by MTX in a panel of T- and B-cell precursors ALL cell lines, aiming to characterize the perturbations caused by the drug on the intracellular metabolome and provide new insights about metabolic vulnerabilities for future therapeutic exploitation.
We were able to identify and quantify 70 metabolites from diverse metabolic pathways; 19 were increased with drug treatment and 12 were decreased. Enrichment analyses showed that the metabolism of several amino acids was particularly impacted, highlighting the glycine and serine metabolism as the topmost altered pathway. Tedeschi et al (55) reported a 2-fold increase in glycine concentration in breast cancer cell lines treated with MTX, which is similar to what we observed in our study. We propose that glycine accumulation is the consequence of an MTX-induced halt in the glycine cleavage system (GCS). The GCS is main mechanism for glycine degradation and requires THF and NAD+ for glycine decarboxylation. MTX leads to THF shortage by inhibiting DHFR, and NAD+ was also decreased following MTX treatment (as seen in Figure 1A). Glycine can be derived from threonine, whose concentration was also increased in MTX-treated cells.
The pentose phosphate pathway (PPP) is an anabolic pathway diverted from glycolysis through glucose-6-phosphate dehydrogenase (G6PD), producing NADPH as well as ribose-5-phosphate, a precursor of purines biosynthesis. MTX was shown to inhibit G6PD in vitro (29). In our study, UDP-glucose, UDP-glucuronate, and myo-inositol, all intermediates of the PPP, were decreased after drug treatment, suggesting that G6PD inhibition in situ is also occurring and is an important feature of MTX action. Additionally, decreased levels of adenine, hypoxanthine, NAD+, AMP, ADP, and ATP also point towards ATIC inhibition — MTX’s main mechanism of action as an anti-inflammatory drug (56). Decreases in NAD+, AMP, ATP, and hypoxanthine levels were also observed in the erythroblastic cell line K562 following MTX treatment (57). Pyrimidine biosynthesis appears to be impaired by the antifolate as well. Thymidylate synthase (TS) was shown to be inhibited in vitro by MTX (58, 59); with TS inhibition, its substrate dUMP would be diverged from dTMP to dCMP synthesis, a precursor of uridine and dCTP — both of which were increased in our data set.
Decreases of several metabolites related to phospholipid metabolism (o-phosphocholine, sn-glycero-3-phosphocholine, phosphoethanolamine, and trimethylamine) indicated that MTX interferes with membrane synthesis. Lipid and membrane synthesis impairment has been reported in MTX-treated osteosarcoma cells (60, 61). In fact, cell cycle progression is a key determinant of phospholipid metabolism and membrane homeostasis: G1/S blockade by RNAi during cellular growth stimulation decreases phosphatidylcholine levels, the main component of cell membranes (62). Coupling lipid metabolism with cell cycle progression is thought to be an evolutionary strategy to coordinate proliferation with cell and organelle growth. The nucleotide shortage provoked by MTX, with the subsequent halt of cell division in the S phase of the cycle, seems to trigger a signal to restrain phospholipid metabolism, thus preventing decompensation of cellular components.
MTX also decreased 2 bile acids, cholate and glycocholate. Bile acids are produced in the liver from cholesterol and might have been interiorized by the cells from the fetal serum present in the culture medium. MRP3, MRP8, and SLC19A1 (folate receptor) are membrane transporters capable of interiorizing MTX, cholate, and glycocholate (63, 64). Polymorphisms and expression levels of these transporters have been associated with prognosis in childhood ALL (65–68) due to impairment on chemotherapy internalization by leukemia cells. In our gene expression data, SLC19A1 presented a negative correlation to MTX resistance, confirming that cellular sensitivity to the drug is related to enhanced drug uptake, whereas cholate and glycocholate concentrations were marginally associated with SLC19A1 expression (P < 0.1) (data not shown). Bile acids such as taurocholate and cholate have proved to be efficient competitors of MTX absorption in hepatocytes (69, 70). Further investigation is required to determine whether seric bile acids could predict MTX uptake by primary leukemia cells and offer therapeutic guidance in clinics.
We found a significant inverse correlation between formate levels in either MTX-treated or untreated cells and cell doubling time or MTX resistance at 48 hours. Formate is a one-carbon molecule that can be synthesized and metabolized in folate-dependent or -independent reactions (71). Metabolic formate elimination appears to occur via its oxidation to carbon dioxide. Formate, THF, and ATP are combined to form 10-formyl-THF, which is then metabolized into THF, NADPH, and CO2. We found that formate levels increased upon MTX treatment (Figure 1A), possibly due to its accumulation by the lack of sufficient THF caused by the drug treatment, as formate production was shown to be markedly reduced in folate deficiency (72). The fact that we found formate levels in both untreated and treated cells to be negatively correlated with cell proliferation and MTX IC50 at 48 hours suggests that high formate level is a constitutive trait of MTX-sensitive cells. Further investigation will detail the mechanistic relationship between formate levels and MTX resistance, as well as show its applicability as an intracellular metabolic MTX marker.
Contrary to the conspicuous associations between formate and cell proliferation and MTX resistance, S-adenosylhomocysteine (SAH) levels were negatively correlated to these parameters only in MTX-treated cells. SAH is produced from methionine S-adenosylmethionine (SAM; the main methyl donor of the metabolism) demethylation and degraded to adenosine and L-homocysteine by the enzyme L-homocysteine hydrolase (AHCY). L-homocysteine is then converted to cystathionine (a precursor of glutathione) or is salvaged to methionine in a reaction that requires 5-methyl-THF (Supplementary Figure 15). We hypothesize that the lower levels of THF and its correlates (including 5-methyl-THF) in MTX-sensitive cells after drug treatment would contribute to SAH accumulation. Our findings showed that methionine cycle intermediates in MTX-treated cells correlate to antifolate resistance, suggesting their potential use as a metabolic biomarker of MTX response that deserves further investigation.
Asparagine regulates mTORC1 activity (73) and can act as an exchange factor, influencing serine metabolism and nucleotide synthesis by regulating serine uptake (74). We speculate that constitutive higher levels of asparagine would predispose cells to resist MTX insult by conferring an extra pool of molecules to be used as exchange factors for serine uptake, leading to 5,10-methylene-THF production and, ultimately, to pyrimidine biosynthesis even in an MTX-mediated DHFR inhibition scenario.
Guanosine is a purine nucleoside and indirect precursor of GTP; it is plausible to assume that higher constitutive levels of this intermediate would confer a natural resistance to MTX, although, to the best of our knowledge, these findings (ie, a positive correlation between basal asparagine and guanosine levels and MTX resistance) have never been described in the literature.
MTX’s ability to induce the formation of ROS has already been described extensively in the literature (75–77). Our work showed that GSH levels were positively associated with MTX resistance in both untreated and treated cells, meaning that cell’s constitutive antioxidant defense is predictive of MTX resistance and remains stable even after drug insult. We found the expression of 8 genes coding for enzymes involved in the GSH metabolism to be positively associated with the metabolite’s level. Among these are the genes GSTM3 and GSTA4, which code for GSH S-transferases capable of MTX glutathionylation (78). Intriguingly, GSTM3 expression in ALL lymphoblasts was positively associated with favorable prognosis (79). Some GSH hydrolases (or gamma-glutamyltransferases [GGTs]) were also correlated with GSH levels. GGTs participate in the gamma-glutamyl cycle, a route of synthesis/degradation of GSH and detoxification of drugs and xenobiotics (80), which explains the positive association we observed between GGT1 expression with both GSH content and MTX resistance. GGT’s expression decrease during ALL treatment returns to normal levels during the chemotherapy-induced remission stage (81), though circulating enzyme levels did not show prognostic value (82).
Parallel to the GSH detox metabolism, thioredoxins are small oxidoreductase enzymes with a dithiol disulfide site that are essential for mammalian life and involved in cellular redox balance and signaling (83). When thioredoxins are oxidized, thioredoxin reductases (TrxR’s) convert them back to their reduced state at the expense of NADPH, which is oxidized to NADP+. Interestingly, among the 3 thioredoxin reductase isozymes found in mammals, only thioredoxin reductase 3 (TXNRD3) contains an additional N-terminal glutaredoxin domain, which makes it a participant of both thioredoxin and glutathione systems. In this study, we found TXNRD3 expression to be positively associated with both GSH content and MTX resistance. TrxR’s enzymes are overexpressed or constitutively expressed in various cancers, contributing to tumor growth and chemoresistance (84), so its inhibition has drawn increasing attention. There are several studies demonstrating that pharmacological inhibitors of TrxR’s can be effective as single agents or combined with chemo/radiotherapy in a variety of cancer types (85, 86). In a study of 14 cancer cell lines from diverse tissues treated with a panel of 19 chemotherapeutic agents, MTX resistance was directly related to the enzymatic activity of glutathione reductase and thioredoxin reductase but not to the intracellular concentration of GSH, which, in turn, was associated with ThioTEPA (an alkylating agent) and doxorubicin resistance (87). Recently, by using CRISPR-Cas9, Oshima et al (88) showed that TXNRD3 (and TXNRD2) knockout in REH cells resulted in increased sensitivity to MTX, confirming the role of the thioredoxin system on MTX detoxification.
In this study, the GSH synthesis promoter N-acetylcysteine was shown to increase MTX resistance and proliferation of ALL cell lines. NAC’s protective effect against diverse drugs in varied cancer types has been extensively described in the literature. However, we noticed an inverse correlation between GSH levels and NAC effect, a phenomenon that we attribute to the feedback loop that coordinates GSH synthesis, through which GSH inhibits GCL activity (89) and thus decreases metabolite synthesis when its cellular pool is already full.
If, by one hand, NAC supplementation boosts MTX resistance in ALL cell lines, co-treatment with GSH scavengers had a marginal effect on sensitizing cells to the antifolate. We hypothesized that this was due to the thioredoxin metabolism, the other antioxidant system of the cell which can compensate an impaired or blocked GSH metabolism (90, 91). We demonstrated the interplay between both antioxidant systems by treating ALL cells with auranofin, a thioredoxin reductase, and BSO, a GSH synthesis inhibitors. The remarkable combination effect observed has been previously reported in head and neck cancer (92, 93), breast cancer (94), and lymphoma and multiple myeloma (95) cells. In some cases, the thioredoxin system may be even more important than GSH to protect cells from oxidative stress. For instance, in prostate and breast cancer, the synergistic combination of 2-deoxy-glucose (2DG) with DHEA (a phosphate-pentose pathway inhibitor) was potentiated by auranofin but not BSO, whereas the cytotoxic effect was completely abrogated by the addition of NAC (96). Our metabolomic results suggest that MTX also inhibits the pentose-phosphate pathway — a theoretically similar effect to that caused by DHEA. Therefore, like the duo MTX+DHEA, the MTX+BSO+auranofin triple combination may deserve further investigation in primary ALL.
Studies in fibroblasts showing that folic acid supplementation abrogates the cytotoxic effects of arsenite have led to the suggestion that MTX could synergize with ATO (97). By causing a relative methyl insufficiency, MTX could potentiate the effects of ATO in the presence of excess folate. We investigated the effects of the thioredoxin system on MTX resistance in vitro by treating ALL cell lines with thioredoxin reductase inhibitors in the presence of the antifolate. Both auranofin and ATO — a rheumatologic and a chemotherapy agent, respectively — potentiated MTX cytotoxicity slightly in a few cell lines, indicating that the thioredoxin antioxidant system could be protecting the cells from the oxidative stress caused by MTX.
ATO has been used against acute promyelocytic leukemia (APL) for the last 2 decades (98), usually in combination with all-trans retinoic acid (ATRA) (99, 100) during both induction and consolidation phases. Prolonged maintenance with continuous low-dose chemotherapy (which includes MTX) and intermittent ATRA was beneficial for high-risk APL (101–103). Aligned with this fact and our in vitro data, we showed that ATO potentiated MTX cytotoxicity in an in vivo model of an aggressive primary murine ALL and in 1 of 4 human ALL PDXs tested, for which tumor burden in the peripheral blood was significantly decreased by the drug combination as well. Interestingly, this ATO-responsive PDX presented a particular MTX-resistant phenotype, whereas in the other 3 cases, single-agent MTX was sufficient to control leukemia progression, thus hampering any possible additive effect by ATO. Both high-dose MTX and ATO are known to cause gastrointestinal toxicity (104, 105), and although no mouse that received either MTX or ATO single-agent experienced severe gastrointestinal toxicity, this was so for 5 out of 24 animals treated with ATO+MTX, indicating that future studies are required to better evaluate the safety of the ATO+MTX combination for the treatment of ALL patients.
Although relapse/refractory ALL did not respond to single-agent ATO in a previous clinical trial (106), which is in line with our PDX results, few studies have investigated its benefits in ALL in combination with other drugs. For instance, Szymanska et al (44) tested the addition of ATO (in the same dose and regimen as ours) with vincristine, dexamethasone, and asparaginase (VXL) in 4 ALL PDXs. Although ATO, in combination with VXL, statistically improved the progression delay of 2 xenografts tested, the authors concluded that the benefit appears unlikely to be of biological significance. We believe that our study brings a new perspective by proposing the use of thioredoxin reductase inhibitors, specially ATO, in combination with MTX. Further investigation should aim to test the Auranofin+MTX combination, the benefits of a third drug inclusion (such as BSO), and the molecular features of the ALL that would most benefit from this new drug combination.
MTX’s mechanisms of action and resistance are still being rewritten even after several decades of research and therapeutic use. Our study showed that the antioxidant systems of the cell are an important component of leukemia resistance to MTX, and targeting these pathways, especially the thioredoxin antioxidant system, may contribute to re-sensitize ALL to MTX. In a broad sense, this work illustrates how metabolomics can be employed as an alternative approach to unravel unexplored mechanisms of action and/or resistance of drugs, leveraging insights about new drug combinations with translational potential for life-threatening diseases such as cancer.
Data availability statement
The data presented in the study are deposited in the Gene Expression Omnibus (GEO) repository, accession number GSE218348.
Ethics statement
The studies involving human participants were reviewed and approved by Comitê de Ética em Pesquisa, Centro Infantil Boldrini, CAAE 34601120.7.0000.5376. Written informed consent to participate in this study was provided by the participants’ legal guardian/next of kin. The animal study was reviewed and approved by Comissão de Ética no Uso de Animais (CEUA), Instituto de Biologia, Universidade Estadual de Campinas (UNICAMP), protocol 4557-1(A)/2018; and Comissão de Ética no Uso de Animais (CEUA), Centro Infantil Boldrini, protocol 0016-2020.
Author contributions
RRC, CPSM, SRB, ACMZ and JAY conceived the study. RRC, CPSM, NMR, LLA, JRC, YTL, SSM, and PRS performed the experiments and analyzed the data. ACMZ, SRB and JAY supervised the study. RC and JY wrote the manuscript. All authors contributed to the article and approved the submitted version.
Funding
This work was supported by São Paulo Research Foundation (FAPESP) grants (08/10034–1 to JAY and 07/00952-0 to ACMZ). RRC (2009/04167-1 and 2012/11952-0), NMC (2014/08247-8), LLA (2017/03239-5), and JRC (2017/02400-7) received FAPESP fellowships. JAY receives a productivity fellowship from the National Counsel of Technological and Scientific Development (CNPq, processes 301596.2017-4 and 308399/2021-8).
Acknowledgments
We acknowledge the Brazilian Biosciences National Laboratory (LNBio) for the use of its facilities and Dr. André B. Silveira for kindly sharing the microarray data of the cell lines.
Conflict of interest
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.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fonc.2022.1032336/full#supplementary-material
References
1. Smith OP, Hann IM. Clinical Features and Therapy of Lymphoblastic Leukemia. In: Arceci RJ, Hann IM, Smith OP, editors. Pediatric Hematology, 3rd ed. (India: Blackwell Publishing) (2006). p. 450–81.
2. Pui C-H, Mullighan CG, Evans WE, Relling MV. Pediatric acute lymphoblastic leukemia: where are we going and how do we get there? Blood (2012) 120(6):1165–74. doi: 10.1182/blood-2012-05-378943
3. Hunger SP, Lu X, Devidas M, Camitta BM, Gaynon PS, Winick NJ, et al. Improved survival for children and adolescents with acute lymphoblastic leukemia between 1990 and 2005: a report from the children's oncology group. J Clin Oncol Off J Am Soc Clin Oncol (2012) 30(14):1663–9. doi: 10.1200/JCO.2011.37.8018
4. Inukai T. Mechanisms of drug resistance in acute lymphoblastic leukemia. [Rinsho ketsueki] Japanese J Clin hematology (2016) 57(10):2022–8. doi: 10.11406/rinketsu.57.2022
5. Zahreddine H, Borden KLB. Mechanisms and insights into drug resistance in cancer. Front Pharmacol (2013) 4:28. doi: 10.3389/fphar.2013.00028
6. Ross DD. Novel mechanisms of drug resistance in leukemia. Leukemia. Mar (2000) 14(3):467–73. doi: 10.1038/sj.leu.2401694
7. Housman G, Byler S, Heerboth S, Lapinska K, Longacre M, Snyder N, et al. Drug resistance in cancer: an overview. Cancers (2014) 6(3):1769–92. doi: 10.3390/cancers6031769
8. Chen ZS, Lee K, Walther S, Raftogianis RB, Kuwano M, Zeng H, et al. Analysis of methotrexate and folate transport by multidrug resistance protein 4 (ABCC4):MRP4 is a component of the methotrexate efflux system. Cancer Res (2002) 62(11):3144–50.
9. Masson E, Relling MV, Synold TW, Liu Q, Schuetz JD, Sandlund JT, et al. Accumulation of methotrexate polyglutamates in lymphoblasts is a determinant of antileukemic effects in vivo: A rationale for high-dose methotrexate. J Clin Invest (1996) 97(1):73–80. doi: 10.1172/JCI118409
10. Lopez-Lopez E, Autry RJ, Smith C, Yang W, Paugh SW, Panetta JC, et al. Pharmacogenomics of intracellular methotrexate polyglutamates in patients' leukemia cells in vivo. J Clin Invest (2020) 130(12):6600–15. doi: 10.1172/JCI140797
11. Sorich MJ, Pottier N, Pei D, Yang W, Kager L, Stocco G, et al. In vivo response to methotrexate forecasts outcome of acute lymphoblastic leukemia and has a distinct gene expression profile. PloS Med (2008) 5(4):0646–56. doi: 10.1371/journal.pmed.0050083
12. Papadopoli DJ, Ma EH, Roy D, Russo M, Bridon G, Avizonis D, et al. Methotrexate elicits pro-respiratory and anti-growth effects by promoting AMPK signaling. Sci Rep May 12 (2020) 10(1):7838. doi: 10.1038/s41598-020-64460-z
13. Leyva-Vázquez MA, Organista-Nava J, Gómez-Gómez Y, Contreras-Quiroz A, Flores-Alfaro E, Illades-Aguiar B. Polymorphism G80A in the reduced folate carrier gene and its relationship to survival and risk of relapse in acute lymphoblastic leukemia. J Invest Med Off Publ Am Fed Clin Res (2012) 60(7):1064–7. doi: 10.2310/JIM.0b013e31826803c1
14. Wojtuszkiewicz A, Peters GJ, van Woerden NL, Dubbelman B, Escherich G, Schmiegelow K, et al. Methotrexate resistance in relation to treatment outcome in childhood acute lymphoblastic leukemia. J Hematol Oncol 2015/12// (2015) 8(1):61–1. doi: 10.1186/s13045-015-0158-9
15. Schuetz JD, Westin EH, Matherly LH, Pincus R, Swerdlow PS, Goldman ID. Membrane protein changes in an L1210 leukemia cell line with a translocation defect in the methotrexate-tetrahydrofolate cofactor transport carrier. J Biol Chem (1989) 264(27):16261–7. doi: 10.1016/S0021-9258(18)71616-5
16. Kremer JM. Toward a better understanding of methotrexate. Arthritis rheumatism (2004) 50(5):1370–82. doi: 10.1002/art.20278
17. Ranganathan P, McLeod HL. Methotrexate pharmacogenetics - The first step toward individualized therapy in rheumatoid arthritis. Arthritis Rheum-Us (2006) 54(5):1366–77. doi: 10.1002/art.21762
18. Kager L, Cheok M, Yang W, Zaza G, Cheng Q, Panetta JC, et al. Folate pathway gene expression differs in subtypes of acute lymphoblastic leukemia and influences methotrexate pharmacodynamics. J Clin Invest (2005) 115(1):110–7. doi: 10.1172/JCI22477
19. Fabre I, Fabre G, Goldman ID. Polyglutamylation, an important element in methotrexate cytotoxicity and selectivity in tumor versus murine granulocytic progenitor cells in vitro. Cancer Res (1984) 44(8):3190–5.
20. Li WW, Lin JT, Tong WP, Trippett TM, Brennan MF, Bertino JR. Mechanisms of natural resistance to antifolates in human soft tissue sarcomas. Cancer Res (1992) 52(6):1434–8.
21. Matherly LH, Taub JW, Ravindranath Y, Proefke SA, Wong SC, Gimotty P, et al. Elevated dihydrofolate reductase and impaired methotrexate transport as elements in methotrexate resistance in childhood acute lymphoblastic leukemia. Blood (1995) 85(2):500–9. doi: 10.1182/blood.V85.2.500.500
22. Srimatkandada S, Schweitzer BI, Moroson BA, Dube S, Bertino JR. Amplification of a polymorphic dihydrofolate reductase gene expressing an enzyme with decreased binding to methotrexate in a human colon carcinoma cell line, HCT-8R4, resistant to this drug. J Biol Chem (1989) 264(6):3524–8. doi: 10.1016/S0021-9258(18)94097-4
23. Rhee MS, Wang Y, Nair MG, Galivan J. Acquisition of resistance to antifolates caused by enhanced gamma-glutamyl hydrolase activity. Cancer Res (1993) 53(10 Suppl):2227–30.
24. Wang Y, Dias JA, Nimec Z, Rotundo R, O'Connor BM, Freisheim J, et al. The properties and function of gamma-glutamyl hydrolase and poly-gamma-glutamate. Adv Enzyme Regul (1993) 33:207–18. doi: 10.1016/0065-2571(93)90019-a
25. Organista-Nava J, Gómez-Gómez Y, Del Moral-Hernandez O, Illades-Aguiar B, Gómez-Santamaria J, Rivera-Ramírez AB, et al. Deregulation of folate pathway gene expression correlates with poor prognosis in acute leukemia. Oncol letters (2019) 18(3):3115–27. doi: 10.3892/ol.2019.10650
26. Rots MG, Pieters R, Peters GJ, van Zantwijk CH, Mauritz R, Noordhuis P, et al. Circumvention of methotrexate resistance in childhood leukemia subtypes by rationally designed antifolates. Blood (1999) 94(9):3121–8. doi: 10.1182/blood.V94.9.3121
27. Panetta JC, Yanishevski Y, Pui C-H, Sandlund JT, Rubnitz J, Rivera GK, et al. A mathematical model of in vivo methotrexate accumulation in acute lymphoblastic leukemia. Cancer chemotherapy Pharmacol (2002) 50(5):419–28. doi: 10.1007/s00280-002-0511-x
28. Waltham MC, Li WW, Gritsman H, Tong WP, Bertino JR. gamma-Glutamyl hydrolase from human sarcoma HT-1080 cells: characterization and inhibition by glutamine antagonists. Mol Pharmacol (1997) 51(5):825–32. doi: 10.1124/mol.51.5.825
29. Babiak RMV, Campello AP, Carnieri EGS, Oliveira MBM. Methotrexate: pentose cycle and oxidative stress. Cell Biochem Funct (1998) 16(4):283–93. doi: 10.1002/(SICI)1099-0844(1998120)16:4<283::AID-CBF801>3.0.CO;2-E
30. Wang YC, Chiang EP. Low-dose methotrexate inhibits methionine S-adenosyltransferase in vitro and in vivo. Mol Med (2012) 18:423–32. doi: 10.2119/molmed.2011.00048
31. Kanarek N, Keys HR, Cantor JR, Lewis CA, Chan SH, Kunchok T, et al. Histidine catabolism is a major determinant of methotrexate sensitivity (vol 559 2021). Nature (2022) 602(7895):pg 632:E17-E18. doi: 10.1038/s41586-021-03487-2
32. Le Belle JE, Harris NG, Williams SR, Bhakoo KK. A comparison of cell and tissue extraction techniques using high-resolution 1H-NMR spectroscopy (vol 15, pg 37, 2002). Nmr Biomed Apr (2016) 29(4):527–7. doi: 10.1002/nbm.3503
33. Xia J, Wishart DS. Using metaboanalyst 3.0 for comprehensive metabolomics data analysis. Curr Protoc Bioinf (2016) 2016:14.10.11–14.10.91. doi: 10.1002/cpbi.11
34. Silveira AB, Laranjeira ABA, Rodrigues GOL, Leal R, Cardoso BA, Barata JT, et al. PI3K inhibition synergizes with glucocorticoids but antagonizes with methotrexate in T-cell acute lymphoblastic leukemia. Oncotarget (2015) 6(15):13105–18. doi: 10.18632/oncotarget.3524
35. Khin ZP, Ribeiro MLC, Jacobson T, Hazlehurst L, Perez L, Baz R, et al. A Preclinical Assay for Chemosensitivity in Multiple Myeloma. Cancer Res (2014) 74(1):56–67. doi: 10.1158/0008-5472.CAN-13-2397
36. Silva A, Silva MC, Sudalagunta P, Distler A, Jacobson T, Collins A, et al. An Ex Vivo Platform for the Prediction of Clinical Response in Multiple Myeloma. Cancer Res (2017) 77(12):3336–51. doi: 10.1158/0008-5472.CAN-17-0502
37. Silva A, Jacobson T, Meads M, Distler A, Shain K. An Organotypic High Throughput System for Characterization of Drug Sensitivity of Primary Multiple Myeloma Cells. J Vis Exp (2015) 101):e53070. doi: 10.3791/53070
38. Sudalagunta P, Silva MC, Canevarolo RR, Alugubelli RR, DeAvila G, Tungesvik A, et al. A pharmacodynamic model of clinical synergy in multiple myeloma. EBioMedicine (2020) 54:102716. doi: 10.1016/j.ebiom.2020.102716
39. Zhao X, Ren Y, Lawlor M, Shah BD, Park PMC, Lwin T, et al. BCL2 Amplicon Loss and Transcriptional Remodeling Drives ABT-199 Resistance in B Cell Lymphoma Models. Cancer Cell (2019) 35(5):752–766.e759. doi: 10.1016/j.ccell.2019.04.005
40. Mostofa A, Distler A, Meads MB, Sahakian E, Powers JJ, Achille A, et al. Plasma cell dependence on histone/protein deacetylase 11 reveals a therapeutic target in multiple myeloma. JCI Insight (2021) 6(24). doi: 10.1172/jci.insight.151713
41. Ren Y, Bi C, Zhao X, Lwin T, Wang C, Yuan J, et al. PLK1 stabilizes a MYC-dependent kinase network in aggressive B cell lymphomas. J Clin Invest (2018) 128(12):5517–30. doi: 10.1172/JCI122533
42. Zhou L, Zhang Y, Meads MB, Dai Y, Ning Y, Hu X, et al. IAP and HDAC inhibitors interact synergistically in myeloma cells through noncanonical NF-κB– and caspase-8–dependent mechanisms. Blood Advances (2021) 5(19):3776–88. doi: 10.1182/bloodadvances.2020003597
43. Almeida ARM, Neto JL, Cachucho A, Euzebio M, Meng X, Kim R, et al. Interleukin-7 receptor alpha mutational activation can initiate precursor B-cell acute lymphoblastic leukemia. Nat Commun (2021) 12(1):7268. doi: 10.1038/s41467-021-27197-5
44. Szymanska B, Wilczynska-Kalak U, Kang MH, Liem NL, Carol H, Boehm I, et al. Pharmacokinetic modeling of an induction regimen for in vivo combined testing of novel drugs against pediatric acute lymphoblastic leukemia xenografts. PloS One (2012) 7(3)::e33894. doi: 10.1371/journal.pone.0033894
45. Lock RB, Liem N, Farnsworth ML, Milross CG, Xue C, Tajbakhsh M, et al. The nonobese diabetic/severe combined immunodeficient (NOD/SCID) mouse model of childhood acute lymphoblastic leukemia reveals intrinsic differences in biologic characteristics at diagnosis and relapse. Blood. Jun 1 (2002) 99(11):4100–8. doi: 10.1182/blood.v99.11.4100
46. Tiziani S, Lodi A, Khanim FL, Viant MR, Bunce CM, Gunther UL. Metabolomic profiling of drug responses in acute myeloid leukaemia cell lines. PloS One (2009) 4(1):e4251. doi: 10.1371/journal.pone.0004251
47. Sharaf BM, Giddey AD, Alniss H, Al-Hroub HM, El-Awady R, Mousa M, et al. Untargeted Metabolomics of Breast Cancer Cells MCF-7 and SkBr3 Treated With Tamoxifen/Trastuzumab. Cancer Genomics Proteomics (2022) 19(1):79–93. doi: 10.21873/cgp.20305
48. Paul M, Hemshekhar M, Thushara RM, Sundaram MS, NaveenKumar SK, Naveen S, et al. Methotrexate Promotes Platelet Apoptosis via JNK-Mediated Mitochondrial Damage: Alleviation by N-Acetylcysteine and N-Acetylcysteine Amide. PloS One (2015) 10(6):e0127558. doi: 10.1371/journal.pone.0127558
49. Kubota M, Nakata R, Adachi S, Watanabe K, Heike T, Takeshita Y, et al. Plasma homocysteine, methionine and S-adenosylhomocysteine levels following high-dose methotrexate treatment in pediatric patients with acute lymphoblastic leukemia or Burkitt lymphoma: association with hepatotoxicity. Leuk Lymphoma (2014) 55(7):1591–5. doi: 10.3109/10428194.2013.850684
50. Hoffman JD, Ward WM, Loo G. Effect of antioxidants on the genotoxicity of phenethyl isothiocyanate. Mutagenesis (2015) 30(3):421–30. doi: 10.1093/mutage/gev003
51. Hong Y-H, Uddin MH, Jo U, Kim B, Song J, Suh DH, et al. ROS Accumulation by PEITC Selectively Kills Ovarian Cancer Cells via UPR-Mediated Apoptosis. Front Oncol (2015) 5:167. doi: 10.3389/fonc.2015.00167
52. Raj L, Ide T, Gurkar AU, Foley M, Schenone M, Li X, et al. Selective killing of cancer cells by a small molecule targeting the stress response to ROS. Nature (2011) 475(7355):231–4. doi: 10.1038/nature10167
53. Lu J, Chew EH, Holmgren A. Targeting thioredoxin reductase is a basis for cancer therapy by arsenic trioxide. P Natl Acad Sci USA (2007) 104(30):12288–93. doi: 10.1073/pnas.0701549104
54. Zenatti PP, Ribeiro D, Li W, Zuurbier L, Silva MC, Paganin M, et al. Oncogenic IL7R gain-of-function mutations in childhood T-cell acute lymphoblastic leukemia. Nat Genet (2011) 43(10):932–9. doi: 10.1038/ng.924
55. Tedeschi PMPM, Johnson-Farley N, Lin H, Shelton LM, Ooga T, Mackay GMG, et al. Quantification of folate metabolism using transient metabolic flux analysis. Cancer Metab (2015) 3(1):6–6. doi: 10.1186/s40170-015-0132-6
56. Budzik GP, Colletti LM, Faltynek CR. Effects of methotrexate on nucleotide pools in normal human T cells and the CEM T cell line. Life Sci (2000) 66(23):2297–307. doi: 10.1016/S0024-3205(00)00559-2
57. Funk RS, Singh RK, Becker ML. Metabolomic Profiling to Identify Molecular Biomarkers of Cellular Response to Methotrexate In Vitro. Clin Transl Sci (2020) 13(1):137–46. doi: 10.1111/cts.12694
58. Hornung N, Stengaard-Pedersen K, Ehrnrooth E, Ellingsen T, Poulsen JH. The effects of low-dose methotrexate on thymidylate synthetase activity in human peripheral blood mononuclear cells. Clin Exp Rheumatol (2000) 18(6):691–8.
59. Chu E, Drake JC, Boarman D, Baram J, Allegra CJ. Mechanism of thymidylate synthase inhibition by methotrexate in human neoplastic cell lines and normal human myeloid progenitor cells. J Biol Chem (1990) 265(15):8470–8. doi: 10.1016/S0021-9258(19)38912-4
60. Lamego I, Duarte IF, Marques MPM, Gil AM. Metabolic Markers of MG-63 Osteosarcoma Cell Line Response to Doxorubicin and Methotrexate Treatment: Comparison to Cisplatin. J Proteome Res (2014) 13(12):6033–45. doi: 10.1021/pr500907d
61. Lamego I, Marques MP, Duarte IF, Martins AS, Oliveira H, Gil AM. Impact of the Pd2Spermine Chelate on Osteosarcoma Metabolism: An NMR Metabolomics Study. J Proteome Res (2017) 16(4):1773–83. doi: 10.1021/acs.jproteome.7b00035
62. Sanchez-Alvarez M, Zhang Q, Finger F, Wakelam MJO, Bakal C. Cell cycle progression is an essential regulatory component of phospholipid metabolism and membrane homeostasis. Open Biol (2015) 5(9):150093–3. doi: 10.1098/rsob.150093
63. Kruh GD, Zeng H, Rea PA, Liu G, Chen ZS, Lee K, et al. MRP subfamily transporters and resistance to anticancer agents. J Bioenerg Biomembr (2001) 33(6):493–501. doi: 10.1023/a:1012827221844
64. Chen Z-S, Guo Y, Belinsky MG, Kotova E, Kruh GD. Transport of bile acids, sulfated steroids, estradiol 17-beta-D-glucuronide, and leukotriene C4 by human multidrug resistance protein 8 (ABCC11). Mol Pharmacol (2005) 67(2):545–57. doi: 10.1124/mol.104.007138
65. Ansari M, Sauty G, Labuda M, Gagné V, Rousseau J, Moghrabi A, et al. Polymorphism in multidrug resistance-associated protein gene 3 is associated with outcomes in childhood acute lymphoblastic leukemia. pharmacogenomics J (2012) 12(5):386–94. doi: 10.1038/tpj.2011.17
66. Steinbach D, Wittig S, Cario G, Viehmann S, Mueller A, Gruhn B, et al. The multidrug resistance-associated protein 3 (MRP3) is associated with a poor outcome in childhood ALL and may account for the worse prognosis in male patients and T-cell immunophenotype. Blood (2003) 102(13):4493–8. doi: 10.1182/blood-2002-11-3461
67. Laverdière C, Chiasson S, Costea I, Moghrabi A, Krajinovic M, Laverdiere C, et al. Polymorphism G80A in the reduced folate carrier gene and its relationship to methotrexate plasma levels and outcome of childhood acute lymphoblastic leukemia. Blood (2002) 100(10):3832–4. doi: 10.1182/blood.V100.10.3832
68. Dulucq S, Laverdière C, Sinnett D, Krajinovic M. Pharmacogenetic considerations for acute lymphoblastic leukemia therapies. Expert Opin Drug Metab Toxicology (2014) 10(5):699–719. doi: 10.1517/17425255.2014.893294
69. Said HM, Strum WB, Hollander D. Inhibitory effect of unconjugated bile acids on the enterohepatic circulation of methotrexate. J Pharmacol Exp Ther Dec (1984) 231(3):660–4.
70. Honscha W, Petzinger E. Characterization of the bile acid sensitive methotrexate carrier of rat liver cells. Naunyn Schmiedebergs Arch Pharmacol May (1999) 359(5):411–9. doi: 10.1007/pl00005369
71. Lamarre SG, Morrow G, Macmillan L, Brosnan ME, Brosnan JT. Formate: an essential metabolite, a biomarker, or more? Clin Chem Lab Med (2013) 51(3):571–8. doi: 10.1515/cclm-2012-0552
72. Morrow GP, MacMillan L, Lamarre SG, Young SK, MacFarlane AJ, Brosnan ME, et al. In vivo kinetics of formate metabolism in folate-deficient and folate-replete rats. J Biol Chem (2015) 290(4):2244–50. doi: 10.1074/jbc.M114.600718
73. Meng D, Yang Q, Wang H, Melick CH, Navlani R, Frank AR, et al. Glutamine and asparagine activate mTORC1 independently of Rag GTPases. J Biol Chem (2020) 295(10):2890–9. doi: 10.1074/jbc.AC119.011578
74. Krall AS, Xu S, Graeber TG, Braas D, Christofk HR. Asparagine promotes cancer cell proliferation through use as an amino acid exchange factor. Nat Commun (2016) 7:11457. doi: 10.1038/ncomms11457
75. Yiang GT, Chou PL, Hung YT, Chen JN, Chang WJ, Yu YL, et al. Vitamin C enhances anticancer activity in methotrexatetreated Hep3B hepatocellular carcinoma cells. Oncol Rep (2014) 32(3):1057–63. doi: 10.3892/or.2014.3289
76. Baquero JM, Benitez-Buelga C, Rajagopal V, Zhenjun Z, Torres-Ruiz R, Muller S, et al. Small molecule inhibitor of OGG1 blocks oxidative DNA damage repair at telomeres and potentiates methotrexate anticancer effects. Sci Rep (2021) 11(1):3490. doi: 10.1038/s41598-021-82917-7
77. Phillips DC, Woollard KJ, Griffiths HR. The anti-inflammatory actions of methotrexate are critically dependent upon the production of reactive oxygen species. Br J Pharmacol (2003) 138(3):501–11. doi: 10.1038/sj.bjp.0705054
78. Barros S, Mencia N, Rodríguez L, Oleaga C, Santos C, Noé V, et al. The Redox State of Cytochrome C Modulates Resistance to Methotrexate in Human MCF7 Breast Cancer Cells. PloS One (2013) 8(5):e63276. doi: 10.1371/journal.pone.0063276
79. Kearns PR, Chrzanowska-Lightowlers ZM, Pieters R, Veerman A, Hall AG. Mu class glutathione S-transferase mRNA isoform expression in acute lymphoblastic leukaemia. Br J Haematol (2003) 120(1):80–8. doi: 10.1046/j.1365-2141.2003.04039.x
80. Courtay C, Oster T, Michelet F, Visvikis A, Diederich M, Wellman M, et al. γ-Glutamyltransferase: Nucleotide sequence of the human pancreatic cDNA: Evidence for a ubiquitous γ-glutamyltransferase polypeptide in human tissues. Biochem Pharmacol (1992) 43(12):2527–33. doi: 10.1016/0006-2952(92)90140-E
81. Russo SA, Harris MB, Greengard O. Diminished lymphocyte and granulocyte gamma-glutamyltranspeptidase activity in acute lymphocytic leukemia and response to chemotherapy. Am J Hematol (1987) 26(1):67–75. doi: 10.1002/ajh.2830260108
82. Roberts BE, Child JA, Cooper EH, Turner R, Stone J. Evaluation of the usefulness of serum gamma-glutamyl transpeptidase levels in the management of haematological neoplasia. Acta Haematol (1978) 59(2):65–72. doi: 10.1159/000207747
83. Nordberg J, Arnér ESJ. Reactive oxygen species, antioxidants, and the mammalian thioredoxin system. Free Radical Biol Med (2001) 31(11):1287–312. doi: 10.1016/S0891-5849(01)00724-9
84. Liu JJ, Liu Q, Wei HL, Yi J, Zhao HS, Gao LP. Inhibition of thioredoxin reductase by auranofin induces apoptosis in adriamycin-resistant human K562 chronic myeloid leukemia cells. Pharmazie (2011) 66(6):440–4. doi: 10.1691/ph.2011.0835
85. Wang H, Bouzakoura S, de Mey S, Jiang H, Law K, Dufait I, et al. Auranofin radiosensitizes tumor cells through targeting thioredoxin reductase and resulting overproduction of reactive oxygen species. Oncotarget (2017) 8(22):35728–42. doi: 10.18632/oncotarget.16113
86. Zhang J, Li X, Han X, Liu R, Fang J. Targeting the Thioredoxin System for Cancer Therapy. Trends Pharmacol Sci Sep (2017) 38(9):794–808. doi: 10.1016/j.tips.2017.06.001
87. Bracht K, Liebeke M, Ritter CA, Gr??nert R, Bednarski PJ. Correlations between the activities of 19 standard anticancer agents, antioxidative enzyme activities and the expression of ATP-binding cassette transporters: comparison with the National Cancer Institute data. Anti-Cancer Drugs (2007) 18(4):389–404. doi: 10.1097/CAD.0b013e3280140001
88. Oshima K, Zhao J, Perez-Duran P, Brown JA, Patino-Galindo JA, Chu T, et al. Mutational and functional genetics mapping of chemotherapy resistance mechanisms in relapsed acute lymphoblastic leukemia. Nat Cancer (2020) 1(11):1113–27. doi: 10.1038/s43018-020-00124-1
89. Griffith OW. Biologic and pharmacologic regulation of mammalian glutathione synthesis. Free Radic Biol Med (1999) 27(9-10):922–35. doi: 10.1016/s0891-5849(99)00176-8
90. Benhar M, Shytaj IL, Stamler JS, Savarino A. Dual targeting of the thioredoxin and glutathione systems in cancer and HIV. J Clin Invest (2016) 126(5):1630–9. doi: 10.1172/JCI85339
91. Scarbrough PM, Mapuskar KA, Mattson DM, Gius D, Watson WH, Spitz DR. Simultaneous inhibition of glutathione- and thioredoxin-dependent metabolism is necessary to potentiate 17AAG-induced cancer cell killing via oxidative stress. Free Radical Biol Med (2012) 52(2):436–43. doi: 10.1016/j.freeradbiomed.2011.10.493
92. Sobhakumari A, Love-Homan L, Fletcher EVM, Martin SM, Parsons AD, Spitz DR, et al. Susceptibility of Human Head and Neck Cancer Cells to Combined Inhibition of Glutathione and Thioredoxin Metabolism. PloS One (2012) 7(10):e48175. doi: 10.1371/journal.pone.0048175
93. Roh JL, Jang H, Kim EH, Shin D. Targeting of the Glutathione, Thioredoxin, and Nrf2 Antioxidant Systems in Head and Neck Cancer. Antioxid Redox Sign (2017) 27(2):106–14. doi: 10.1089/ars.2016.6841
94. Harris Isaac S, Treloar Aislinn E, Inoue S, Sasaki M, Gorrini C, Lee Kim C, et al. Glutathione and Thioredoxin Antioxidant Pathways Synergize to Drive Cancer Initiation and Progression. Cancer Cell (2015) 27(2):211–22. doi: 10.1016/j.ccell.2014.11.019
95. Kiebala M, Skalska J, Casulo C, Brookes PS, Peterson DR, Hilchey SP, et al. Dual targeting of the thioredoxin and glutathione antioxidant systems in malignant B cells: A novel synergistic therapeutic approach. Exp Hematology (2015) 43(2):89–99. doi: 10.1016/j.exphem.2014.10.004
96. Li L, Fath MA, Scarbrough PM, Watson WH, Spitz DR. Combined inhibition of glycolysis, the pentose cycle, and thioredoxin metabolism selectively increases cytotoxicity and oxidative stress in human breast and prostate cancer. Redox Biol (2015) 4:127–35. doi: 10.1016/j.redox.2014.12.001
97. Ruan Y, Peterson MH, Wauson EM, Waes JG, Finnell RH, Vorce RL. Folic acid protects SWV/Fnn embryo fibroblasts against arsenic toxicity. Toxicol Lett (2000) 117(3):129–37. doi: 10.1016/s0378-4274(00)00254-x
98. Lengfelder E, Hofmann WK, Nowak D. Impact of arsenic trioxide in the treatment of acute promyelocytic leukemia. Leukemia (2012) 26(3):433–42. doi: 10.1038/leu.2011.245
99. Platzbecker U, Avvisati G, Cicconi L, Thiede C, Paoloni F, Vignetti M, et al. Improved Outcomes With Retinoic Acid and Arsenic Trioxide Compared With Retinoic Acid and Chemotherapy in Non-High-Risk Acute Promyelocytic Leukemia: Final Results of the Randomized Italian-German APL0406 Trial. J Clin Oncol (2017) 35(6):605–12. doi: 10.1200/JCO.2016.67.1982
100. Lo-Coco F, Avvisati G, Vignetti M, Thiede C, Orlando SM, Iacobelli S, et al. Retinoic acid and arsenic trioxide for acute promyelocytic leukemia. N Engl J Med (2013) 369(2):111–21. doi: 10.1056/NEJMoa1300874
101. Ades L, Guerci A, Raffoux E, Sanz M, Chevallier P, Lapusan S, et al. Very long-term outcome of acute promyelocytic leukemia after treatment with all-trans retinoic acid and chemotherapy: the European APL Group experience. Blood (2010) 115(9):1690–6. doi: 10.1182/blood-2009-07-233387
102. Muchtar E, Vidal L, Ram R, Gafter-Gvili A, Shpilberg O, Raanani P. The role of maintenance therapy in acute promyelocytic leukemia in the first complete remission. Cochrane Database Syst Rev (2013) 3):CD009594. doi: 10.1002/14651858.CD009594.pub2
103. Huang H, Qin Y, Xu R, You X, Teng R, Yang L, et al. Combination therapy with arsenic trioxide, all-trans retinoic acid, and chemotherapy in acute promyelocytic leukemia patients with various relapse risks. Leuk Res (2012) 36(7):841–5. doi: 10.1016/j.leukres.2012.03.027
104. Howard SC, McCormick J, Pui CH, Buddington RK, Harvey RD. Preventing and Managing Toxicities of High-Dose Methotrexate. Oncologist (2016) 21(12):1471–82. doi: 10.1634/theoncologist.2015-0164
105. Zhang P. On arsenic trioxide in the clinical treatment of acute promyelocytic leukemia. Leuk Res Rep (2017) 7:29–32. doi: 10.1016/j.lrr.2017.03.001
Keywords: acute lymphoblastic leukemia, methotrexate, glutathione, metabolomics, drug resistance, arsenic trioxide, thioredoxin reductase
Citation: Canevarolo RR, Melo CPS, Cury NM, Artico LL, Corrêa JR, Lau YT, Mariano SS, Sudalagunta PR, Brandalise SR, Zeri ACM and Yunes JA (2022) Glutathione levels are associated with methotrexate resistance in acute lymphoblastic leukemia cell lines. Front. Oncol. 12:1032336. doi: 10.3389/fonc.2022.1032336
Received: 30 August 2022; Accepted: 02 November 2022;
Published: 01 December 2022.
Edited by:
Giuseppe Maurizio Campo, University of Messina, ItalyCopyright © 2022 Canevarolo, Melo, Cury, Artico, Corrêa, Lau, Mariano, Sudalagunta, Brandalise, Zeri and Yunes. 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(s) 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: José Andrés Yunes, YW5kcmVzQGJvbGRyaW5pLm9yZy5icg==
†Present address: Rafael Renatino Canevarolo, Department of Cancer Physiology, H. Lee Moffitt Cancer Center & Research Institute, Tampa, FL, United States