Activating mutations in NOTCH1 are common in T cell acute lymphoblastic leukemia (T-ALL). Here we identify glutaminolysis as a critical pathway for leukemia cell growth downstream of NOTCH1 and a key determinant of the response to anti-NOTCH1 therapies in vivo. Mechanistically, inhibition of NOTCH1 signaling in T-ALL induces a metabolic shutdown, with prominent inhibition of glutaminolysis and triggers autophagy as a salvage pathway supporting leukemia cell metabolism. Consequently, inhibition of glutaminolysis and inhibition of autophagy strongly and synergistically enhance the antileukemic effects of anti-NOTCH1 therapy in mice harboring T-ALL. Moreover, we demonstrate that Pten loss upregulates glycolysis and consequently rescues leukemic cell metabolism, thereby abrogating the antileukemic effects of NOTCH1 inhibition. Overall, these results identify glutaminolysis as a major node in cancer metabolism controlled by NOTCH1 and as therapeutic target for the treatment of T-ALL.
Activating mutations in NOTCH1 are common in T cell acute lymphoblastic leukemia (T-ALL). Here we identify glutaminolysis as a critical pathway for leukemia cell growth downstream of NOTCH1 and a key determinant of the response to anti-NOTCH1 therapies in vivo. Mechanistically, inhibition of NOTCH1 signaling in T-ALL induces a metabolic shutdown, with prominent inhibition of glutaminolysis and triggers autophagy as a salvage pathway supporting leukemia cell metabolism. Consequently, inhibition of glutaminolysis and inhibition of autophagy strongly and synergistically enhance the antileukemic effects of anti-NOTCH1 therapy in mice harboring T-ALL. Moreover, we demonstrate that Pten loss upregulates glycolysis and consequently rescues leukemic cell metabolism, thereby abrogating the antileukemic effects of NOTCH1 inhibition. Overall, these results identify glutaminolysis as a major node in cancer metabolism controlled by NOTCH1 and as therapeutic target for the treatment of T-ALL.
NOTCH signaling is a conserved signal transduction pathway with a prominent role in cell differentiation and tissue patterning during development[1]. In the hematopoietic system, NOTCH1 has been implicated in stem cell homeostasis and most prominently as a major driver of T-cell lineage specification in lymphoid progenitors and a master regulator of thymocyte development[2-4]. In addition, aberrant NOTCH1 signaling plays a major role in the pathogenesis of over 60% of T-ALLs harboring activating mutations in the NOTCH1 gene[5]. Most notably, oncogenic NOTCH1 has been proposed as a therapeutic target in NOTCH1-mutant leukemias and small molecule γ-secretase inhibitors (GSIs), which effectively block NOTCH activation via inhibition of a critical intramembrane proteolytic cleavage required for NOTCH signaling[6], are now in clinical development for the treatment of relapsed and refractory T-ALL. However, the clinical development of anti-NOTCH1 therapies in T-ALL has been hampered by the limited and delayed therapeutic response to these drugs, making the identification of highly effective and synergistic drug combinations capable of delivering strong antileukemic responses a top priority in the field. In addition, and most troubling, most T-ALL cell lines harboring activating mutations in NOTCH1 fail to respond to GSI therapy, a phenotype strictly associated with mutational loss of the Phophatase and tensin homolog (PTEN) tumor suppressor gene[7], making essential to establish the specific role and mechanisms of PTEN inactivation as driver of resistance to anti-NOTCH1 therapies.
RESULTS
Pten loss confers resistance to NOTCH inhibition in T-ALL
To analyze the effects of Pten inactivation in the response of primary NOTCH1-induced leukemia cells to GSI therapy in vivo we generated a mouse model of NOTCH1 induced T-ALL with conditional and inducible loss of Pten. Towards this goal we infected bone marrow hematopoietic progenitors from tamoxifen-inducible conditional Pten knockout mice (Rosa26) with retroviruses expressing a constitutively active oncogenic mutant form of the NOTCH1 receptor (NOTCH1L1601P Δ-PEST)[8] and transplanted them into isogenic recipients, which consequently developed NOTCH1-induced T-ALLs. We then injected these NOTCH1-induced Pten-conditional tumor cells into secondary recipients and treated them with vehicle only or tamoxifen in order to generate Pten-positive (non-deleted) and Pten-deleted isogenic tumors, respectively. Treatment of Pten-positive leukemia-bearing mice with DBZ, a highly active GSI[9], induced a marked response to therapy by in vivo bioimaging (Fig. 1a) and a significant improvement in survival compared with vehicle-only treated controls (P < 0.005) (Fig. 1b and Supplementary Fig. 1). In contrast, all mice harboring isogenic Pten-deleted tumors showed increased resistance to anti-NOTCH inhibition therapy and progressed under treatment (Fig. 1a,b and Supplementary Fig. 1). Analysis of cell proliferation in T-ALL lymphoblasts isolated from Pten-positive and Pten-deleted leukemia-bearing mice treated with DBZ showed decreased proliferation upon NOTCH inhibition in Pten-positive tumors, a phenotype rescued by genetic loss of Pten (Fig. 1c). Importantly, analysis of NOTCH1 signaling showed complete clearance of activated NOTCH1 protein (ICN1) both in Pten-positive and Pten-deleted tumors treated with DBZ, supporting that Pten loss does not impair the uptake or intrinsic activity of this GSI (Fig. 1d). Moreover, Myc, a critical downstream effector of the oncogenic effects of NOTCH1 was effectively downregulated in Pten-positive and Pten-deleted leukemias treated with DBZ (Fig. 1d and Supplementary Fig. 1), ruling out increased Myc expression secondary to Pten loss as a potential mechanism of escape from the antileukemic effects of NOTCH1 inhibition. Next, and to assess the effects of isogenic PTEN loss in human cells, we infected a human primary xenograft (PDTALL#19) with lentiviruses expressing a shRNA targeting PTEN (shPTEN) or a shRNA control (shLUC), and confirmed the knockdown of PTEN levels in cells expressing shPTEN (Supplementary Fig. 2). Expression of the shLUC did not alter the response to GSI (Supplementary Fig. 2). In contrast, and most notably, PTEN knockdown restored leukemia cell growth in the context of GSI treatment (Supplementary Fig. 2). Overall, these results show that Pten loss and consequent constitutive activation of the PI3K-AKT pathway can confer resistance to anti-NOTCH1GSI therapy in vivo.
Figure 1
Pten loss induces resistance to GSI treatment in vivo
(a) Changes in tumor load by bioimaging in mice allografted with NOTCH1-induced Pten-positive or Pten-deleted isogenic leukemias treated with vehicle or DBZ (n = 2 per group). (b) Kaplan-Meier survival curve of mice harboring Pten-positive and Pten-deleted isogenic leukemias treated with 4 cycles of vehicle or DBZ (5 mg Kg−1) on a 4 days ON (red blocks) and 3 days OFF schedule (Pten+/+/DMSO, n = 11; Pten+/+/DBZ, n = 9; Pten−/−/DMSO and Pten−/−/DBZ, n = 10) (c) Cell cycle analyses of leukemias acutely treated with vehicle or DBZ. (d) Western blot analyses of intracellular activated NOTCH1, Myc, Akt and Pten expression in leukemias acutely treated with vehicle or DBZ. (e) Relative mRNA expression of the NOTCH1 target Hes1 in leukemias acutely treated with vehicle or DBZ. (f) Volcano plot representations of gene expression changes induced by GSI treatment in Pten-positive and Pten-deleted leukemias. Downregulated or upregulated genes by DBZ in Pten-positive cells are marked in green or red, respectively. Dashed boxes indicate differentially expressed genes. Solid line boxes contain non-significantly differentially expressed genes. (g) Heat map representation of the top differentially expressed genes upon DBZ treatment in Pten-positive and Pten-deleted leukemias. (h) Heat map representation of the top differentially expressed genes upon DBZ treatment in Pten-positive lymphoblasts that are rescued by Pten loss. P values (c,e) were calculated using two-tailed Student’s t-test. Bar graphs indicate mean ± s.d. (n = 3 for Pten−/− groups and n = 6 for Pten+/+ groups).
To investigate the underlying mechanisms mediating resistance to NOTCH inhibition in Pten-null T-ALL tumor cells we performed gene expression profiling of isogenic Pten-positive and Pten-deleted leukemia lymphoblasts after acute treatment with DBZ in vivo. This analysis revealed that, while direct NOTCH1 target genes (such as Hes1, Dtx1, Ptcra, Heyl and Notch3) are effectively downregulated in both Pten-positive and Pten-deleted tumors (Fig 1e–g and Supplementary Fig. 1), genetic ablation of Pten elicits a global reversal of much of the transcriptional effects of NOTCH inhibition (Fig. 1f,h and Supplementary Fig. 1). Functional annotation of genes downregulated by NOTCH inhibition whose expression is restored upon Pten loss revealed a marked enrichment in pathways associated with cell anabolism, such as ribosomal RNA processing and amino acid and nucleobase biosynthesis (Fig. 1f and Supplementary Table 1). Conversely, genes selectively upregulated by GSI treatment in Pten-positive tumors only are predominantly implicated in apoptosis (e.g. regulation of cell death and caspase regulator activity), autophagy (e.g. lytic vacuole, regulation of autophagy and lysosome) and catabolism (e.g. lipid catabolism, ubiquitin conjugation and phagocytosis)(Fig. 1f and Supplementary Table 1).
Impaired metabolism by NOTCH inhibition in T-ALL
Following up on these results, we decided to explore the metabolic effects of NOTCH inhibition and Pten loss by performing a broad-based metabolomic analysis by LC-MS/MS of isogenic Pten-positive and Pten-deleted NOTCH1-induced leukemia cells treated with the DBZGSI in vivo. These analyses showed that inhibition of NOTCH signaling by DBZ in NOTCH1-induced Pten-positive tumors resulted in the accumulation of glucose and proximal (glucose-6-phosphate, fructose-6-phosphate and fructose-1,6-biphosphate) and distal glycolytic intermediates (3-phosphoglycerate and 2,3-diphosphoglycerate) (Fig. 2a), coupled with higher levels of ribose-5-phosphate and ribulose-5-phosphate/xylulose-5-phosphate in the pentose phosphate pathway. Moreover, inhibition of NOTCH signaling was associated with elevated free amino acid levels (Supplementary Fig. 3), potentially linked with increased autophagy; and increased levels of glutamine, but not glutamate (Fig. 2a). Notably, loss of Pten resulted in increased lactate levels (Fig. 2a) and reversed the accumulation of glycolytic intermediates induced by NOTCH1 inhibition in Pten-positive T-ALL cells (Fig. 2a).
Figure 2
Metabolic profiling analysis and metabolic rescue of NOTCH1 inhibition in T-ALL. (a) GC/MS and LC/MS/MS metabolic profiles (mass spectrometry scaled intensity) of DBZ-treated mice harboring Pten-positive or Pten-deleted isogenic leukemias. Changes in glycolysis, pentose phospate pathway, glutaminolysis and TCA cycle intermediates are shown (scaled intensity arbitrary units). Box plots represent the upper quartile to lower quartile distribution. + Sign indicates mean value and horizontal line the median value. Whiskers indicate the maximum and minimum values of the distribution. Open circles indicate extreme data points. (n = 3 for Pten−/− groups and n = 6 for Pten+/+ groups). (b,c) Representative flow cytometry histograms (b) and quantitative analyses (c) showing cell diameter changes induced by DBZ treatment in DND41 cells treated with DMSO only or DBZ, in the presence or absence of methyl pyruvate. (d) Cell cycle analysis of DND41 cells treated with vehicle only or DBZ, in the presence or absence of methyl pyruvate. (e,f) Representative flow cytometry histograms (e) and quantitative analyses (f) showing cell diameter changes induced by DBZ treatment in DND41 cells treated with DMSO only or DBZ, in the presence or absence of dimethyl ketoglutarate. (g) Cell cycle analysis of DND41 cells treated with DMSO only or DBZ, in the presence or absence of dimethyl ketoglutarate. MP= methyl pyruvate; DMK= dimethyl ketoglurate. P values were calculated using two-tailed Student’s t-test. Bar graphs indicate mean ± s.d of biological triplicates.
To directly assess the role of impaired carbon metabolism in mediating the antileukemic effects of NOTCH1 inhibition with GSIs, we evaluated the capacity of methyl pyruvate, a membrane soluble metabolite that bypasses glycolysis and can be incorporated directly into the tricarboxylic acid cycle (TCA cycle)[10], to rescue the effects of NOTCH inhibition in DND41, a NOTCH1-mutated and PTEN-positive T-ALL cell line[7]. Consistent with previous reports[11], inhibition of NOTCH1 signaling with DBZ induced decreased cell size and proliferation with cell cycle arrest in G1 (Fig. 2b–d). Methyl pyruvate treatment induced a slight increase in cell size, and abrogated the antileukemic effects of NOTCH inhibition on cell growth (10.5% reduction in size by DBZ in vehicle control cells vs. 2.6% decrease in cell diameters in DBZ treated cells grown in media supplemented with methyl pyruvate, P < 0.001) and proliferation (Fig. 2b–d). Similarly, bypass of glutaminolysis with membrane-soluble dimethyl α-ketoglutarate[12], effectively antagonized the inhibitory effects of NOTCH1 inhibition in cell size (7.7% reduction in size by DBZ in vehicle control cells vs. 2.6% decrease in cell diameters in DBZ treated cells grown in media supplemented with dimethyl α-ketoglutarate, P < 0.001) and proliferation (Fig. 2e–g), further supporting a major role for inhibition of carbon metabolism as a key effector of the antileukemic effects of NOTCH1 inhibition in T-ALL. We obtained similar results in a second NOTCH1-mutated PTEN-positive cell line (HPB-ALL) (Supplementary Fig. 4).
Pten-positive T-ALLs rely on autophagy upon NOTCH inhibition
The transcriptional and metabolic changes associated with NOTCH1 inhibition in Pten-positive T-ALL tumor cells are indicative of decreased anabolism and a concomitant increase in catabolic processes. In this context, electron microscopy analysis showed increased autophagy in vivo upon DBZ treatment in Pten-positive NOTCH-induced leukemias, which was effectively rescued in Pten-deleted tumors (Fig. 3a). In addition, western blot analyses showed increased levels of the LC3a II autophagic marker in Pten-positive leukemias treated with DBZ, whereas Pten loss efficiently rescued this phenotype (Fig. 3b and Supplementary Fig. 5). Based on these observations we hypothesized that induction of autophagy can contribute to sustain cell survival during NOTCH1 inhibition by recycling essential metabolites required for leukemia cell metabolism. To test this model, we generated NOTCH1-induced T-ALLs from tamoxifen-inducible conditional Atg7 knockout (Rosa267)[13] hematopoietic progenitors, which were subsequently treated with vehicle only or tamoxifen in order to generate Atg7-positive (non-deleted) and Atg7-deleted isogenic tumors, respectively (Fig. 3c). Treatment of NOTCH1-induced Atg7-positive leukemia-bearing mice with DBZ induced a marked in vivo antileukemic response and significantly improved survival (P < 0.005) (Fig. 3d), which was markedly increased after abrogation of autophagy via tamoxifen-induced deletion of Atg7 (P < 0.005) (Fig. 3d).
Figure 3
Autophagy supports leukemia cell growth in response to NOTCH1 inhibition. (a) Representative images of electron microscopy micrographs of Pten-positive and Pten-deleted lymphoblasts acutely treated with vehicle only or DBZ in vivo. Quantification of the percentage of cells with autophagosomes is shown on the right (n = 3 per group). (b) Quantification of LC3a II / LC3a I protein ratio after acute treatment of leukemic mice with vehicle only or DBZ. P values were calculated using two-tailed Student’s t-test. Bar graphs indicate mean ± s.d. (n = 3 per group) (c) Western blot analyses showing Atg7 expression in Rosa267 leukemias treated with tamoxifen. (d) Kaplan-Meier survival curve of mice harboring Atg7-positive and Atg7-deleted isogenic leukemias treated with 4 cycles of vehicle or DBZ (5 mg Kg−1) on a 4 days ON (red blocks) and 3 days OFF schedule (n = 5 per group). Scale bars, 500 nm (a).
Overall, these results support that inhibition of NOTCH1 signaling drives a metabolic crisis in Pten-positive T-ALL tumor cells, rendering leukemia lymphoblasts dependent on autophagy for growth and survival.
NOTCH-induced glutaminolysis is a key carbon source in T-ALL
Next, and to further explore the anabolic effects of NOTCH inhibition and Pten deletion in glycolysis and glutaminolysis we performed metabolic tracing studies with 13C isotope-labeled glucose (U-13C6) and glutamine (U-13C5). These studies revealed that NOTCH1-induced leukemias made prominent use of glutamine as a source of carbon. Thus, glutamine-derived 13C was readily converted to glutamate and effectively incorporated into all TCA cycle intermediates (~80% labeling), supporting a major role for glutaminolysis as a carbon source in the metabolism of NOTCH1-induced leukemia (Fig. 4a,b). Glucose labeled the lactate pool (~80% labeling) and contributed to 40–60% of the TCA cycle intermediates citrate, fumarate and malate (Fig. 4a,b). Notably, Pten deletion resulted in increased lactate production suggestive of increased glycolysis (Fig. 4c). Moreover, kinetic profiling analysis of glucose and glutamine-derived 13C incorporation into the TCA cycle showed that NOTCH1 inhibition with DBZ decreased the fractional contribution of glucose to lactate, but not to TCA cycle intermediates (Fig. 4d). In contrast, DBZ treatment decreased glutamine labeling of glutamate and TCA cycle intermediates (Fig. 4e). Importantly, the block in glycolysis observed in NOTCH1-induced Pten-positive leukemia lymphoblasts treated with DBZ was decreased by isogenic deletion of Pten (Fig. 4d). In addition, Pten loss reduced the DBZ-induced block in glutamine to glutamate conversion but had no clear effects in the incorporation of glutamine derived carbon to the TCA cycle (Fig. 4e).
Figure 4
Glucose and glutamine metabolic flux analysis of T-ALL cells upon NOTCH1 inhibition and Pten loss. (a) Percentage of 13C incorporation into lactate and glutamate after incubation of primary NOTCH1-induced T-ALL cells with 13C-glucose and 13C-glutamine, respectively. (b) Percentage of 13C incorporation into TCA cycle intermediates after incubation of primary NOTCH1-induced T-ALL cells with 13C-glucose and 13C-glutamine, respectively. (c) Extracellular levels of lactate produced by Pten-positive and Pten-deleted leukemia cells after 6 hours of culture in vitro using dialyzed medium. (d) Kinetic analysis of glucose incorporation in primary lymphoblasts incubated with 13C-glucose. (e) Kinetic analysis of glutamine incorporation in primary lymphoblasts incubated with 13C-glutamine. M+2, +3, +4 or +5 labeled compounds indicate molecules of those compounds that contain 2, 3, 4 or 5 13C atoms, respectively. Bar graphs and kinetic curves indicate mean ± s.d. P values were calculated using Student’s t-test between control and DBZ treated samples across all time points (***P < 0.005; **P < 0.01; *P < 0.05; n = 3 per group).
Next, we investigated the expression levels of Pten, p-AKT and glutaminase (Gls) in Pten wild-type NOTCH1-induced leukemias at the time of disease progression after the four cycles of treatment with DBZ. Notably, these leukemic cells show decreased levels of Pten, with a concomitant increase in p-AKT levels, as well as increased Gls (Fig. 5a). These results suggest that T-ALL lymphoblasts treated with GSI might be overcoming the therapeutic effects of NOTCH1 inhibition by downregulating Pten and upregulating Gls expression in vivo.
Figure 5
Rescue of GSI antileukemic effects in Pten-positive leukemias. (a) Western blot analyses showing intracellular activated NOTCH1, p(T308) Akt, Akt total, Pten and Gls expression in Pten-positive leukemias progressing after 4 cycles of treatment with vehicle or DBZ (5 mg Kg−1) on a 4 days ON (red blocks) and 3 days OFF schedule. (b-g) Peripheral blood leukemia infiltration in mice harboring NOTCH1-induced Pten-positive T-ALL cells expressing (b) MyrAKT plus mCherry, (c) Glutaminase plus mCherry, (d) PKM2 plus mCherry, (e) PHGDH plus mCherry (f) IL7R plus mCherry or (g) preTCRalpha (Ptcra) plus mCherry, upon treatment with DBZ. Changes in leukemia cell counts of non-infected (mCherry negative) cells are shown as internal control. Graphs indicate median ± s.d. P values were calculated using two-tailed Student’s t-test (n = 5 mice per group).
Following on these results we further explored the mechanisms downstream of NOTCH1 and PTEN in the control of leukemia cell growth. Expression of critical signaling factors controlled by NOTCH1 and downregulated by GSI treatment in both Pten-positive and Pten-deleted T-ALL cells such as IL7R and Ptcra had no effect on the response of Pten-positive cells to either vehicle or GSI therapy (Fig. 5b–g and Supplementary Fig. 6). In contrast, mice transplanted with Pten-positive T-ALL tumor cells infected with constitutively active AKT (myristoylated-AKT) reversed the anabolic (downregulated) and catabolic (upregulated) transcriptional changes associated with NOTCH1 inhibition (Supplementary Fig. 7), as well as the glycolytic block, the block in glutaminolysis and the increase in free amino acids observed by LC-MS/MS (Supplementary Fig. 8), and induced overt resistance to NOTCH1 inhibition and progression under treatment (Fig. 5b). In addition, analysis of metabolic regulators showed that expression of GLS significantly restored leukemia cell growth in the context of inhibition of NOTCH1 signaling (Fig. 5c). Similarly, expression of PKM2, a cancer associated isoform of the glycolytic enzyme pyruvate kinase resulted in enhanced leukemia cell growth in the context of inhibition of NOTCH1 signaling (Fig. 5d). In contrast, expression of 3-Phosphoglycerate dehydrogenase (PHGDH), the first and rate-limiting step in the phosphorylated pathway of serine biosynthesis also associated with tumor cell growth[14], had limited or no effect in the response of tumor cells to NOTCH1 inhibition with DBZ (Fig. 5e). Cells overexpressing GLS increase their glutamine utilization (Supplementary Fig. 9) but still show downregulation of glutaminolysis upon DBZ treatment in glutamine-derived 13C kinetic labeling analyses (Supplementary Fig. 9), suggesting that increased GLS expression may confer resistance to DBZ treatment in this experiment by sheer mass effect. Overall, these results demonstrate that metabolic manipulation can modify the in vivo response to GSI treatment in NOTCH1-induced leukemias.
Therapeutic targeting of NOTCH1 and glutaminase in T-ALL
The identification of glutaminolysis as a critical effector of the antileukemic response to NOTCH1 inhibition in T-ALL supports that pharmacologic inhibition of glutaminase may impair leukemic cell growth and sensitize T-ALL cells to NOTCH1 inhibition therapies. To test this hypothesis, we analyzed the effects of glutaminase inhibition on leukemia cell growth and survival. Treatment of HPB-ALL and DND41 T-ALL cells with BPTES, a potent and specific glutaminase inhibitor[15], impaired leukemia cell growth (Fig. 6a and Supplementary Fig. 10), and showed strong and significantly synergistic antileukemic effects in combination with NOTCH1 inhibition with DBZ (Combination Index at ED50 = 0.012 in HPB-ALL and 0.666 in DND-41 cells) (Fig. 6b and Supplementary Fig. 10). Notably, this phenotype was primarily driven by increased cytotoxicity in cells treated with BPTES and DBZ in combination (Fig. 6c,d and Supplementary Fig. 10). In contrast, treatment with 2-deoxyglucose, a non-metabolizable glucose analog that inhibits glycolysis[16], showed additive antileukemic effects in combination with DBZ, further supporting a dominant role for glutaminolysis over glycolysis downstream of NOTCH1 signaling in T-ALL (Supplementary Fig. 11). Importantly, Crispr-Cas9-mediated PTEN inactivation in the HPB-ALL cell line conferred resistance to treatment with DBZ alone and in combination with BPTES, further supporting a role for PTEN loss as driver of resistance to GSI therapy in human cells (Supplementary Figure 12). To further explore the clinical relevance of glutaminase inhibition therapies in combination with GSIs we next tested the efficacy of the NOTCH1 inhibition with DBZ, glutaminase inhibition with BPTES and the combination of DBZ plus BPTES compared to vehicle only (DMSO) treatment in mice xenografted with two independent NOTCH1-mutant (PDTALL#10: NOTCH1L1601P, R1609H; PDTALL#19: NOTCH1 P1606LV, V2412fs) and PTEN wild type human-derived primary T-ALL xenografts. In this experiment, mice treated with BPTES showed progressive tumor growth in vivo similar to that observed in vehicle-treated controls, while NOTCH inhibition with DBZ induced significant antitumor responses (Fig. 6e,f), which were markedly and significantly increased by glutaminase inhibition in the DBZ plus BPTES treatment group (Fig. 6e,f).
Figure 6
Synergistic antileukemic effects of glutaminase inhibition and GSI treatment in T-ALL. (a) Differential growth of HPB-ALL cells treated with DMSO, DBZ, BPTES, or DBZ plus BPTES combination in vitro. (b) Isobologram analysis of DBZ plus BPTES combination in HPB-ALL cells. X shape mark in red shows value for combination index at ED50. (c,d) Representative flow cytometry plots of Annexin V/7AAD staining (c) and quantitation of apoptosis (d) of HPB-ALL cells treated with the different drug combinations. (e,f) Representative images (left) and quantitation (right) of tumor burden changes assessed by bioimaging in mice xenografted with human primary T-ALL cells PDTALL#19 (e) or PDTALL#10 (f) treated with the different drug combinations (n = 5 per group except PDTALL#19 DBZ+BPTES, n = 4). (g,h) Representative images (left) and quantitation (right) of tumor burden changes assessed by bioimaging in mice allografted with NOTCH1-induced Pten-positive (g) or Pten-deleted (h) mouse leukemia lymphoblasts treated with the different drug combinations (n = 5 per group except: Pten−/−/Vehicle, n = 4; Pten−/−/BPTES, n = 3). (i) Western blot analyses showing Gls expression in Rosa26 leukemias treated with vehicle or tamoxifen. (j) Kaplan-Meier survival curve of mice harboring Gls-positive and Gls-deleted isogenic leukemias treated with 4 cycles of vehicle or DBZ (5 mg Kg−1) on a 4 days ON (red blocks) and 3 days OFF schedule (***P < 0.005; n = 10 per group). P values (a,d) were calculated using two-tailed Student’s t-test. Bar graphs indicate mean ± s.d. of biological triplicates. Scale bars, 1 cm (e–h).
Next, and to better assess the interaction between NOTCH signaling and Pten loss in the response to anti-NOTCH1 and glutaminase inhibition therapies, we analyzed the effects of vehicle only (DMSO), BPTES, DBZ and BPTES plus DBZ in mice transplanted with our NOTCH1-induced Pten-positive or NOTCH1-induced Pten-deleted mouse isogenic tumor cells. In this experiment, inhibition of glutaminolysis significantly increased the antileukemic effects of NOTCH1 inhibition in NOTCH1-induced Pten-positive T-ALL bearing animals treated with DBZ plus BPTES (Fig. 6g and Supplementary Fig. 13), which translated into increased survival (Supplementary Fig. 14). In contrast, and most notably, Pten-deleted T-ALL cells showed not only impaired response to GSI therapy, but also to treatment with DBZ plus BPTES in combination (Fig. 6h). These results are consistent with our metabolic studies showing that Pten loss induces a hyperglycolytic phenotype, which would render T-ALL cells not only resistant to NOTCH inhibition, but also less sensitive to inhibition of glutaminolysis as result of increased glucose derived carbon input to the TCA cycle.Finally, to genetically test the interaction between NOTCH signaling and glutaminase in vivo and rule out a potential off-target effect of BPTES, we generated NOTCH1-induced T-ALLs from tamoxifen-inducible conditional Gls knockout (Rosa26) hematopoietic progenitors, which were subsequently transplanted into secondary recipients and treated with vehicle only or tamoxifen in order to generate Gls-positive (non-deleted) and Gls-deleted isogenic tumors, respectively (Fig. 6i). Treatment of NOTCH1-induced Gls-positive leukemia-bearing mice with DBZ induced a marked in vivo antileukemic response and significantly improved survival (P < 0.005) (Fig. 6j). In addition, tamoxifen-induced deletion of glutaminase resulted in overt anti-leukemic effects, which were markedly and significantly enhanced (P < 0.005; 40% complete remission) upon GSI treatment (Fig. 6j). This finding fully validates our pharmacological results and highlights glutaminolysis as a key player and therapeutic target in NOTCH1-induced T-ALL.
DISCUSSION
Co-occurring activating mutations in NOTCH1 and mutations and deletions in PTEN are characteristic of GSI-resistant T-ALL cell lines[7] and can also be found in primary patient samples[7,17]. However, PTEN mutant leukemias may be still at least partially responsive to NOTCH1 inhibition[18] and epigenetic reprogramming of NOTCH1-driven T-ALL cell lines has been recently implicated in tumor escape from anti-NOTCH1 therapies after prolonged in vitro treatment with GSIs[19], supporting that additional genetic and epigenetic factors may coordinately contribute to drive resistance to GSI based anti-NOTCH1 therapies. To formally address the specific role of Pten inactivation as driver of resistance to GSI therapy we engineered a model of NOTCH1 induced leukemia with conditional inducible loss of Pten, which allowed the direct comparison of Pten-positive and Pten-deleted isogenic tumors in immune-competent isogenic secondary recipients in vivo. In this model, NOTCH1-induced leukemias showed effective but variable antileukemic responses to GSI therapy even in the context of tumors generated with the same activated NOTCH1 allele. These results substantiate the value of our analyses across heterogeneous tumors with different levels of primary NOTCH1 oncogene addiction. Most notably, although Pten deletion did not render the cells completely insensitive to NOTCH inhibition with GSI, in each case mutational loss of Pten effectively impaired the response to NOTCH1 inhibition therapy and all mice harboring Pten-deleted leukemias showed overt progression under therapy, a hallmark criteria of therapy resistance in the clinic, with 100% mortality by the end of the fourth cycle of treatment. Similarly, PTEN knockdown in a primary humanT-ALL xenograft also impaired the therapeutic efficacy of GSI treatment. These results conclusively demonstrate that mutational loss of Pten can induce resistance to GSI therapy in primary tumor cells in vivo.Mechanistically, gene expression profiling analyses revealed that inhibition of the γ-secretase complex effectively suppressed NOTCH1 activation and that key genes controlled by NOTCH1 were effectively downregulated in both Pten-wild type and Pten-deleted tumors. In addition pathway analysis of gene expression programs controlled by NOTCH1 in Pten-wild type and Pten-deleted T-ALL tumor cells revealed a dominant effect NOTCH1 inhibition in the suppression of anabolic routes with concomitant upregulation of catabolic genes. The broad transcriptional effects of NOTCH1 inhibition across different metabolic pathways are reminiscent of those induced by KrasG12D inactivation in pancreatic ductal adenocarcinoma tumor cells[20] and support that the antileukemic effects of NOTCH1 inhibition in T-ALL may be mediated at least in part by inhibition of cell metabolism. Notably, this model is consistent with the decreased cell growth and cell size phenotypes typically induced by NOTCH1 inhibition in PTEN wild type GSI-sensitive human T-ALL cell lines[5,7]. Unexpectedly, outside NOTCH1 direct target genes, Pten deletion had profound effects in the transcriptional response of NOTCH1-induced leukemias to anti-NOTCH1 therapy and induced effective reversal of GSI-induced anabolic gene downregulation and abrogation of the catabolic signature associated with NOTCH1 inhibition in Pten-wild type tumors. These results support a significant transcriptional input of the PI3K-AKT signaling pathway in the control of cell metabolism that is antagonistic to that of suppression of oncogenic NOTCH1 in T-ALL. Consistently, we observed increased autophagy in Pten-wild type tumors, but not in Pten-deleted leukemias upon NOTCH1 inhibition. Activation of autophagy upon GSI treatment and increased GSI antileukemic effects in autophagy-deficient Atg7 null tumors support that NOTCH1 inhibition forces tumor cells to resource to this catabolic mechanism to obtain critical metabolic intermediates for sustained cell growth in Pten-wild type T-ALLs[21].Integration of these results with global metabolic profiling and metabolic isotope tracing analyses revealed a negative effect of NOTCH1 inhibition in glycolysis and glutaminolysis in Pten-wild type leukemias. Notably, genetic loss of Pten or activated AKT expression resulted in increased intracellular lactate levels supporting increased glycolytic efficiency in Pten-deleted tumors, a finding consistent with the well-established role of the PI3K-AKT pathway in inducing aerobic glycolysis[10,22-25]. Moreover, and perhaps most surprisingly, metabolic isotope labeling analyses unveiled a most prominent role for glutaminolysis in feeding carbon into the TCA cycle. Consistently, glutaminase inhibition or genetic deletion of glutaminase were highly synergistic with GSI treatment and induced marked therapeutic responses in vivo. Notably, the therapeutic response to combined NOTCH and glutaminase inhibition is also abrogated by loss of Pten. This is not surprising since the switch to a hypergycolytic phenotype induced by Pten loss would make T-ALL cells less dependent on glutaminolysis because of increased glucose derived carbon input to the TCA cycle.Overall, our results highlight the fundamental importance of NOTCH1 signaling in the control of leukemia cell metabolism, extend our understanding of the prominent role of PTEN and the PI3K pathway in the control of oncogenic cell growth and provide the basis for the design of new therapeutic strategies targeting cell metabolism for the treatment of T-ALL.
ONLINE METHODS
Mice and animal procedures
We maintained all animals in specific pathogen-free facilities at the Irving Cancer Research Center at Columbia University Medical Campus. The Columbia University Institutional Animal Care and Use Committee (IACUC) approved all animal procedures. Rosa26mice expressing a tamoxifen-inducible form of the Cre recombinase from the ubiquitous Rosa26 locus[26], Pten conditional knockout (Pten)[27] and Atg7 conditional knockout (Atg7)[13] mice have been previously described. To generate NOTCH1-induced T-ALL tumors in mice, we performed retroviral transduction of bone marrow cells enriched in Lineage negative cells (using Miltenyi kit #130-090-858, following manufacturer’s guidelines) with activated forms of the NOTCH1 oncogene (NOTCH1L1601P ΔPEST in Pten cells or ΔE-NOTCH1 in Atg7 and Gls cells) and transplanted them via intravenous injection into lethally irradiated recipients as previously described[8].For all subsequent studies, we bought secondary recipient C57BL6 females (6–8 week old) from Taconic Farms. Animals were randomly assigned to the different treatment groups and no blinding was done. For survival studies, we transplanted leukemia cells from primary recipients into secondary recipients sub-lethally irradiated (4 Gy); two days after the transplant, we treated the micetamoxifen (5mg per mouse by intra-peritoneal injection) to induce deletion of either the Pten, the Atg7 or the Gls locus, or with vehicle only (Corn oil) for the control groups; 4 days later, we subdivided the mice in groups treated with vehicle only (2.3% DMSO in 0.00.5% Methylcellulose, 0.1% Tween-80) or with DBZ (5 mg kg−1) on a 4 days-on and 3 days-off basis.For acute treatment analyses, we transplanted leukemia cells from primary recipients into secondary recipients, which were treated with vehicle or tamoxifen (to induce Pten deletion), as described before. We monitored the mice until they showed more than 70% of GFP positive cells in peripheral blood and then, mice were treated twice, 16 hours apart, with vehicle only or DBZ (5 mg kg−1). 4 hours after the second treatment, mice were sacrificed and splenic samples were collected and subsequently analyzed (13C-labeling) or snap-frozen for further analyses (microarrays, RT-PCR, Western Blots, metabolomics).For the cherry rescue experiment in vivo, we infected leukemia cells from the mousetumors previously generated with retroviral particles expressing the red cherry fluorescent protein (MSCV empty-mCherry) or a fusion protein of the red cherry fluorescent protein together with either IL7R (MSCVIL7R-IRES-mCherry), preTCR-alpha (MSCV Ptcra-IRES-mCherry), myristoylated-AKT (MSCV Myr-AKT-IRES-mCherry), Glutaminase (MSCVGLS-IRES-mCherry), PKM2 (MSCVPKM2-IRES-mCherry) or PHGDH (MSCVPHGDH-IRES-mCherry) and re-injected them in sub-lethally irradiated C57/BL6 mice (4 Gy).For drug synergism studies in vivo, we infected leukemia cells from the mousetumors previously generated with retroviral particles expressing a fusion protein between the red cherry fluorescent protein and luciferase (MigR1 mCherryLuc) and re-injected them in sub-lethally irradiated C57/BL6 mice (4 Gy). For the human primary leukemia xenograft experiment, we infected PDTALL#19 and PDTALL#10 cells with lentiviral particles expressing the red cherry fluorescent protein and luciferase (FUW-mCherry-Luc) and re-injected them in NRG mice (6 – 8 week old, both males and females). We analyzed the efficacy of glutaminase inhibition with BPTES in combination with DBZ in secondary NRG recipient mice transplanted with these mCherry-Luciferase expressing cells. In these experiments, we treated the mice with daily intraperitoneal doses of vehicle (DMSO), BPTES (25 mg kg−1), DBZ (5 mg kg−1) or BPTES (25 mg kg−1) plus DBZ (5 mg kg−1) for 6 days. We evaluated disease progression and therapy response by in vivo bioimaging with the In Vivo Imaging System (IVIS, Xenogen).
Microarray gene expression profiling
We isolated RNA from leukemic cells from spleen samples collected from mice acutely treated with Vehicle or DBZ using RNeasy plus mini kit (Qiagen) according to manufacturer’s protocol. RNA (1 μg) from Pten-positive and Pten-deleted leukemic cells treated with vehicle or DBZ was amplified and labeled using 3′ IVT Express Kit (Affymetrix) and hybridized on GeneChip Mouse Genome 430 2.0 Arrays (Affymetrix) according to the manufacturer’s protocol. RNA (1 μg) from Pten-positive mCherry and Pten-positive myr-AKTleukemic cells treated with vehicle or DBZ was amplified and labeled with the Illumina TotalprepRNA kit and hybridized on Illumina MouseWG-6 v2.0 BeadChip arrays according to the manufacturer’s protocol. For the Affymetrix arrays, we performed inter-array normalization with the GC-RMA algorithm using open-source Bioconductor software[28]. For the Illumina arrays, expression values were log2 transformed and quantile normalized. We evaluated group differences using t-test and fold change.For the Pathway Enrichment analyses, gene sets of interest, including differentially expressed genes by GSI that were rescued by Pten deletion, were tested for functional annotations enrichment using the web-based DAVID bioinformatics tools available at http://david.abcc.ncifcrf.gov.We compared the enrichment of genes differentially expressed by GSI in Pten-positive and Pten-deleted tumors with the genes differentially expressed by GSI in Pten-positive mCherry expressing and in Pten-positive mCherry-myrAKT expressing leukemias by GSEA[29] using the t-test metric and 10,000 permutations of the gene list.
Metabolomic analyses
We analyzed leukemic spleen samples from mice acutely treated with Vehicle or DBZ by GC/MS and LC/MS/MS. Briefly, we collected splenic samples and snap froze them 4 hours after the second round of in vivo treatment with Vehicle or DBZ. We then extracted the samples and prepared them for analysis using standard solvent extraction methods. The extracted samples were split into equal parts for analysis on the GC/MS and LC/MS/MS platforms in Metabolon facilities.Dataset comprises a total 267 named biochemicals. Following log transformation and imputation with minimum observed values for each compound, we used ANOVA contrasts for pairwise comparisons to identify biochemicals that differed significantly between experimental groups. Analysis by two-way ANOVA identified biochemicals exhibiting significant interaction and main effects for experimental parameters of Pten expression and DBZ treatment. We calculated an estimate of the false discovery rate (q-value) to take into account the multiple comparisons that normally occur in metabolomic-based studies. For example, when analyzing 200 compounds, we would expect to see about 10 compounds meeting the p ≤ 0.05 cut-off by random chance.
13C labeling studies
We cultured leukemia cells collected from the spleens of acutely treated mice (as described before) in RPMI media supplemented with 10% dialyzed fetal bovine serum and 10 mM 13C-Glucose or 4 mM 13C-Glutamine (Cambridge Isotope Laboratories). Cells were cultured for 6 hours for the steady-state glucose and glutamine experiments and for 30 minutes, 1 hour, 2 hours and 6 hours for the glucose and glutamine kinetics experiment. After labeling, we rinsed the cells in ice-cold normal saline and lysed them with three freeze-thaw cycles in cold 50% methanol. The lysates were centrifuged to remove precipitated protein, a standard (50 nmols of sodium 2-oxobutyrate) was added, and the samples were evaporated and derivatized by trimethylsilylation (Tri-Sil HTP reagent, Thermo Scientific). We then injected three microliters of the derivatized material onto an Agilent 6970 gas chromatograph equipped with a fused silica capillary GC column (30 m length, 0.25 mm diameter) and networked to an Agilent 5973 mass selective detector. Retention times of all metabolites of interest were validated using pure standards. The measured distribution of mass isotopomers was corrected for natural abundance of 13C.Extracellular lactate was measured in 100ul of media from cells that had been cultured for 6 hrs in media made with dialyzed serum. Removal of lipids and proteins from the media was accomplished with a 1:1:1 chloroform:methanol:water extraction. Prior to evaporation, a uniformly labeled 13C lactate standard (Cambridge) was added to the aqueous phase. After being dried, the sample was derivatized with 100ul of Tri-Sil (Thermo Scientific) for 10min at 75C then 10min at room temperature. 5ul of the derivatized sample was injected onto an Agilent 6890 GC networked to an Agilent 5975 Mass Selective Detector. The comparison of the labeled lactate standard to the unlabeled lactate pool derived from cellular metabolism allowed for the quantitation of lactate abundance in the sample. We measured ammonia with a spectrophotometric assay (Megazyme).
Quantitative real time PCR
We generated cDNA with the ThermoScript RT-PCR system (Invitrogen) and analyzed it by quantitative real-time PCR (FastStart Universal SYBR Green Master Mix (Roche) using a 7300 Real-Time PCR System (Applied Biosystems). Relative expression levels were based on Gapdh as a reference control.
Western Blotting
We performed Western blots using standard procedures. Antibodies against Cleaved Notch1 (#4147), GAPDH (#5174), LC3a (#4599), p(S473)AKT (#9271), p(T308)AKT (#9275), AKT total (#9272) and PTEN (#9188) were from Cell Signaling technologies (1:1000 dilution). Glutaminase antibody was from proteintech (1:1000, #20170-1-AP) or a purified rabbit polyclonal antibody against GLS[30] (1:1000 dilution). ATG7 antibody was from Sigma-Aldrich (1:1000, #A2856). Myc antibody (1:200, N-262) was obtained from Santa Cruz. GFP antibody (1:1000, #11814460001) was obtained from Roche.
Inhibitors and drugs
DBZ ((S)-2-(2-(3,5-Difluorophenyl)acetamido)-N-((S)-5-methyl-6-oxo-6,7-dihydro-5H-dibenzo[b,d]azepin-7-yl)propanamide) was obtained from Syncom (#SIC-020042). BPTES (Bis-2-(5-phenylacetamido-1,3,4-thiadiazol-2-yl)ethyl sulfide; #SML0601), 2-deoxyglucose (#D6134), tamoxifen (#T5648), dimethyl ketoglutarate (#349631) and methyl pyruvate (#371173) were obtained from Sigma-Aldrich.
Cell lines
(HEK) 293T (purchased from ATCC), HPB-ALL and DND41 (obtained from DSMZ) cell lines were cultured in standard conditions in RPMI media supplemented with 10% FBS and 1% Penicillin/Streptomycin. Cell lines were regularly authenticated and tested for mycoplasma. To CRISPR-out PTEN in HPB-ALL cells, we designed a gRNA against PTEN using the E-CRISP software (http://www.e-crisp.org/E-CRISP/designcrispr.html). This gRNA was subsequently cloned into pL-CRISPR.EFS.GFP vector (Addgene plasmid #57818)[31]. We then infected HPB-ALL cells with this construct and the GFP positive population was sorted and further analyzed.
Human primary xenografts
T-ALL samples were provided by Columbia Presbyterian Hospital and University of Padova with informed consent and analyzed under the supervision of the Columbia University Medical Center Institutional Review Board committee.
Cell viability, cell size and flow cytometric analysis
We analyzed cell line viability upon treatment with BPTES (10 μM) and DBZ (250 nM) alone and in combination using the Cell Proliferation Kit I (Roche). We analyzed apoptosis by flow cytometry with APC AnnexinV Apoptosis Kit I (BD Biosciences). We used Propidium Iodide (Sigma) DNA staining to analyze cell cycle distribution.For the metabolic rescue experiment, cells were maintained on RPMI buffered with 40 mM HEPES, and treated with methyl pyruvate (10 mM) or dimethyl ketoglutarate (8 mM) and/or DBZ (250 nM). We added new media every two days.For the measurement of PTEN intracellular levels in PDTALL 19 human primary xenografted samples, we removed red cells from peripheral blood samples by incubation with red blood cell lysis buffer (155 mM NH4Cl, 12 mM KHCO3 and 0.1 mM EDTA) for 5 min at room temperature. We then stained the samples with anti-human-CD45 (eBioscience, Cat# 17-0459-42, clone HI30, dilution 1:200). Next, we performed fixation/permeabilization using Foxp3/Transcription Factor Staining Buffer Set (eBioscience Cat# 00-5523), following manufacturer’s instructions. Finally, we performed intracellular staining with anti-PTEN-PE (BD Phosflow, Cat# 560002, clone A2B1, dilution 1:5) and we used PE-labeled mouseIgG1 isotype as a control (BD Pharmingen, Cat# 551436, dilution 1:5).
Statistical analyses
We performed statistical analysis by Student’s t-test. We considered results with P < 0.05 as statistically significant. We analyzed drug synergism using the median-effect method of Choy and Talay[32] and used the CalcuSyn software (Biosoft, Great Shelford, Cambridge, UK) to calculate the combination index (CI) and perform isobologram analysis of drug interactions. Survival in mouse experiments was represented with Kaplan-Meier curves and significance was estimated with the log-rank test (Prism GraphPad).
Authors: Andrew P Weng; Adolfo A Ferrando; Woojoong Lee; John P Morris; Lewis B Silverman; Cheryll Sanchez-Irizarry; Stephen C Blacklow; A Thomas Look; Jon C Aster Journal: Science Date: 2004-10-08 Impact factor: 47.728
Authors: Teresa Palomero; Wei Keat Lim; Duncan T Odom; Maria Luisa Sulis; Pedro J Real; Adam Margolin; Kelly C Barnes; Jennifer O'Neil; Donna Neuberg; Andrew P Weng; Jon C Aster; Francois Sigaux; Jean Soulier; A Thomas Look; Richard A Young; Andrea Califano; Adolfo A Ferrando Journal: Proc Natl Acad Sci U S A Date: 2006-11-17 Impact factor: 11.205
Authors: Monica Buzzai; Daniel E Bauer; Russell G Jones; Ralph J Deberardinis; Georgia Hatzivassiliou; Rebecca L Elstrom; Craig B Thompson Journal: Oncogene Date: 2005-06-16 Impact factor: 9.867
Authors: Aravind Subramanian; Pablo Tamayo; Vamsi K Mootha; Sayan Mukherjee; Benjamin L Ebert; Michael A Gillette; Amanda Paulovich; Scott L Pomeroy; Todd R Golub; Eric S Lander; Jill P Mesirov Journal: Proc Natl Acad Sci U S A Date: 2005-09-30 Impact factor: 11.205
Authors: José A Campos-Sandoval; Amada R López de la Oliva; Carolina Lobo; Juan A Segura; José M Matés; Francisco J Alonso; Javier Márquez Journal: Int J Biochem Cell Biol Date: 2006-12-21 Impact factor: 5.085
Authors: Georgia Hatzivassiliou; Fangping Zhao; Daniel E Bauer; Charalambos Andreadis; Anthony N Shaw; Dashyant Dhanak; Sunil R Hingorani; David A Tuveson; Craig B Thompson Journal: Cancer Cell Date: 2005-10 Impact factor: 31.743
Authors: Teresa Palomero; Maria Luisa Sulis; Maria Cortina; Pedro J Real; Kelly Barnes; Maria Ciofani; Esther Caparros; Jean Buteau; Kristy Brown; Sherrie L Perkins; Govind Bhagat; Archana M Agarwal; Giuseppe Basso; Mireia Castillo; Satoru Nagase; Carlos Cordon-Cardo; Ramon Parsons; Juan Carlos Zúñiga-Pflücker; Maria Dominguez; Adolfo A Ferrando Journal: Nat Med Date: 2007-09-16 Impact factor: 53.440
Authors: Mary M Robinson; Steven J McBryant; Takashi Tsukamoto; Camilo Rojas; Dana V Ferraris; Sean K Hamilton; Jeffrey C Hansen; Norman P Curthoys Journal: Biochem J Date: 2007-09-15 Impact factor: 3.857
Authors: Evgeni Efimenko; Utpal P Davé; Irina V Lebedeva; Yao Shen; Maria J Sanchez-Quintero; Daniel Diolaiti; Andrew Kung; Brian J Lannutti; Jianchung Chen; Ronald Realubit; Zoya Niatsetskaya; Vadim Ten; Charles Karan; Xi Chen; Andrea Califano; Thomas G Diacovo Journal: Mol Cancer Ther Date: 2017-07-17 Impact factor: 6.261
Authors: Nina Fenouille; Christopher F Bassil; Issam Ben-Sahra; Lina Benajiba; Gabriela Alexe; Azucena Ramos; Yana Pikman; Amy S Conway; Michael R Burgess; Qing Li; Frédéric Luciano; Patrick Auberger; Ilene Galinsky; Daniel J DeAngelo; Richard M Stone; Yi Zhang; Archibald S Perkins; Kevin Shannon; Michael T Hemann; Alexandre Puissant; Kimberly Stegmaier Journal: Nat Med Date: 2017-02-13 Impact factor: 53.440
Authors: Brett M Noel; Steven B Ouellette; Laura Marholz; Deborah Dickey; Connor Navis; Tzu-Yi Yang; Vinh Nguyen; Sarah J Parker; David Bernlohr; Zohar Sachs; Laurie L Parker Journal: J Proteome Res Date: 2019-02-21 Impact factor: 4.466