- Open Access
Cysteine catabolism and the serine biosynthesis pathway support pyruvate production during pyruvate kinase knockdown in pancreatic cancer cells
Cancer & Metabolism volume 7, Article number: 13 (2019)
Pancreatic ductal adenocarcinoma (PDAC) is an aggressive cancer with limited treatment options. Pyruvate kinase, especially the M2 isoform (PKM2), is highly expressed in PDAC cells, but its role in pancreatic cancer remains controversial. To investigate the role of pyruvate kinase in pancreatic cancer, we knocked down PKM2 individually as well as both PKM1 and PKM2 concurrently (PKM1/2) in cell lines derived from a KrasG12D/-; p53-/- pancreatic mouse model.
We used liquid chromatography tandem mass spectrometry (LC-MS/MS) to determine metabolic profiles of wildtype and PKM1/2 knockdown PDAC cells. We further used stable isotope-labeled metabolic precursors and LC-MS/MS to determine metabolic pathways upregulated in PKM1/2 knockdown cells. We then targeted metabolic pathways upregulated in PKM1/2 knockdown cells using CRISPR/Cas9 gene editing technology.
PDAC cells are able to proliferate and continue to produce pyruvate despite PKM1/2 knockdown. The serine biosynthesis pathway partially contributed to pyruvate production during PKM1/2 knockdown: knockout of phosphoglycerate dehydrogenase in this pathway decreased pyruvate production from glucose. In addition, cysteine catabolism generated ~ 20% of intracellular pyruvate in PDAC cells. Other potential sources of pyruvate include the sialic acid pathway and catabolism of glutamine, serine, tryptophan, and threonine. However, these sources did not provide significant levels of pyruvate in PKM1/2 knockdown cells.
PKM1/2 knockdown does not impact the proliferation of pancreatic cancer cells. The serine biosynthesis pathway supports conversion of glucose to pyruvate during pyruvate kinase knockdown. However, direct conversion of serine to pyruvate was not observed during PKM1/2 knockdown. Investigating several alternative sources of pyruvate identified cysteine catabolism for pyruvate production during PKM1/2 knockdown. Surprisingly, we find that a large percentage of intracellular pyruvate comes from cysteine. Our results highlight the ability of PDAC cells to adaptively rewire their metabolic pathways during knockdown of a key metabolic enzyme.
Pancreatic ductal adenocarcinoma (PDAC) is a devastating disease with < 10% of patients surviving beyond 5 years after diagnosis . It is currently the fourth leading cause of cancer-related death in western countries and is expected to be the second leading cause by 2030 . Currently, treatment options for pancreatic cancer patients include surgical resection, radiation therapy, and/or systemic chemotherapy [3, 4]. However, current therapies often fail because most patients are diagnosed at advanced stages. Given the poor outlook for PDAC patients, it is critical to improve our understanding of pancreatic cancer cells to design improved treatment strategies.
Reprogramming cancer metabolism is recognized as a hallmark of cancer . Many tumor cells exhibit the “Warburg effect,” fermenting glucose to lactate even in the presence of abundant of oxygen [6, 7]. Indeed, pancreatic cancer cells rewire metabolism to utilize a wide range of nutrients including glucose, extracellular proteins, and various amino acids to support survival and proliferation [8,9,10]. Targeting metabolic pathways specifically upregulated in PDAC cells may be a promising direction for therapy .
Altered metabolism in cancer cells is supported in part by the expression of a specific isoform of pyruvate kinase (PK), a glycolytic enzyme which catalyzes the conversion of phosphoenolpyruvate (PEP) and ADP into pyruvate and ATP . Pyruvate kinase has four isoforms (L, R, M1, and M2) encoded by two genes . The PKLR gene encodes PKL and PKR, and the PKM gene encodes PKM1 and PKM2 through alternative splicing of exons 9 and 10, respectively . PKL is mainly expressed in the liver, kidney, and small intestine and PKR in erythrocytes [15, 16]. PKM1 is mainly expressed in differentiated tissues such as muscle, heart, and brain, whereas PKM2 is expressed in various adult tissues and many proliferating cells, including embryonic and tumor cells [13, 15,16,17]. Generally, expression of PKM1 and PKM2 is mutually exclusive in a given cell type, and loss of PKM2 leads to compensatory expression of PKM1 . PKM2 is highly expressed in a variety of human cancer cells, including pancreatic cancer cells . PKM2 has been reported to promote proliferation, migration, invasion, and tumorigenesis in pancreatic cancer [20,21,22,23]. However, there are conflicting reports regarding the expression of PKM2 and overall patient survival: some studies show that PKM2 expression is associated with worse overall survival of pancreatic cancer patients [20, 24, 25], while others show improved overall survival , and still others show no effect on overall survival [27, 28]. A recent study demonstrated that expression of PKM1, but not PKM2, promotes small-cell lung cancer cell growth . Other studies reported that PKM2 is dispensable for leukemia, liver cancer, colon cancer, lymphoma, lung cancer, and pancreatic cancer [30,31,32,33,34].
Given the controversial roles of pyruvate kinase in cancer, we investigated the function of PKM1/2 for pancreatic cancer in PDAC cells derived from a KrasG12D/-; p53-/- pancreatic mouse model. Our results demonstrate that knockdown of PKM2 results in expression of PKM1 and does not affect pancreatic cancer cell proliferation. Additionally, pancreatic cancer cells are able to proliferate even with concurrent knockdown of both PKM1 and PKM2 isoforms. Further, they are able to produce pyruvate from glucose with knockdown of both PKM1 and PKM2 isoforms. We explored the contributions of alternative pathways to pyruvate production, such as the serine biosynthesis pathway and the sialic acid pathway during PKM1/2 knockdown. Knockout of phosphoglycerate dehydrogenase (PHGDH), the rate-limiting enzyme in the serine biosynthesis pathway, decreased pyruvate production from glucose in PKM1/2 knockdown cells. However, knockout of N-acetylneuraminate pyruvate lyase (NPL), which can convert phosphoenolpyruvate to pyruvate in the sialic acid pathway, did not decrease pyruvate production from glucose. Using multiple isotopically labeled precursors, we discovered that glucose contributes to only ~ 40% of intracellular pyruvate, and amino acid cysteine contributes ~ 20% of intracellular pyruvate in these cells. The relatively low contribution of glucose to pyruvate in combination with alternative sources for pyruvate generation may explain the minimal impact of PKM1/2 knockdown on pancreatic cancer cell proliferation.
A13M2-1 and A13M13 cell lines were derived from KrasG12D/-; p53-/- pancreatic mouse tumor and contain doxycycline-inducible hairpins that target PKM2 or both PKM1 and PKM2, respectively [35, 36]. Cells were cultured in DMEM (Fisher Scientific, MT10017CV) without sodium pyruvate, supplemented with 10% fetal bovine serum (FBS), 1% penicillin and streptomycin (P/S), and 1% glutamine and cultured in a humidified incubator with 5% CO2 at 37 °C. For cells infected with inducible vectors, doxycycline aqueous solution (1 mg/ml as stock solution) was added at the time of plating at a final concentration of 1 μg/ml. Cells were passaged in doxycycline-containing media for 7 days to maximize knockdown of PKM1/2 before experiments were conducted.
Cell proliferation assay
The cells were seeded at 25,000 per well in 6-well plates in normal cell growth medium. Cell counts were obtained using a Cellometer Auto T4 Cell Counter (Nexcelom).
Gene knockout by CRISPR/Cas9 gene editing
We used CRISPR/Cas9-mediated genome editing to achieve gene knockout with lentivirus-mediated gene expression . Dual-guide RNAs targeting PHGDH or NPL gene were designed by CRISPR DESIGN (http://crispr.mit.edu/) and set just before the protospacer adjacent motif (PAM), a DNA sequence immediately following the Cas9-targeted DNA sequence. All the specific target sequences were amplified and cloned into lentiviral vectors and verified by DNA sequencing. CRISPR gene editing plasmid vectors with gRNA and Cas9 co-expression were acquired from VectorBuilder. The VSVG plasmid was a gift from Bob Weinberg (Addgene plasmid # 8454; http://n2t.net/addgene:8454; RRID:Addgene 8454). The psPAX2 plasmid was a gift from Didier Trono (Addgene plasmid # 12260; http://n2t.net/addgene:12260; RRID:Addgene 12260). To produce lentivirus, HEK293T cells seeded in 10-cm plates were transfected with 10.0 μg lentivirus plasmids, 0.5 μg VSVG, and 5.0 μg psPAX2 plasmids. The following morning, fresh DMEM with 15% FBS and 1% P/S was added, and cells were grown for another 48 h to generate virus. For transduction with lentivirus, the A13M13 cells (1 × 105 cells) were seeded in 10-cm plates and the supernatant of transfected HEK293T was collected and passed through 0.45 micron PVDF syringe filter. Five milliliters of the viral supernatant and 5 ml of fresh media were added to recipient A13M13 cell plates with polybrene (Fisher Scientific, TR1003G) at final concentration of 4 μg/ml. The cells were cultured for 24 h followed by adding fresh DMEM medium supplemented with 10% FBS and treated for 10 days with 10 μg/ml blasticidin (Fisher Scientific, A1113903) for selection. The blasticidin selected cells were then resuspended to a concentration of 5 cells/ml and seeded 1 cell/well on 96-well plates. Surviving clones were expanded and analyzed for successful gene knockout. Genomic DNA was extracted using DNeasy Blood and Tissue Kit (Qiagen) to check for successful gene editing. Dual gRNAs for PHGDH knockout plasmid vector are 5′-CGGGCTCAGCCCTCCGACCC-3′ and 5′-GGTGCTCCCTACCAAGCCGT-3′. Dual gRNAs for NPL knockout plasmid vector are 5′- GAGCGTCTCTGAACGTCGCC-3′ and 5′-CGTGGGAGCACTAAACGTGA-3′. To address PHGDH/NPL dual knockout, we designed a plasmid vector containing gRNAs for both PHGDH and NPL. The gRNAs are 5′-CGTGGGAGCACTAAACGTGA-3′ (targeting PHGDH) and 5′-CGGGCTCAGCCCTCCGACCC-3′ (targeting NPL). The sequence for scramble control plasmid with blasticidin selection maker is 5′-GCACTACCAGAGCTAACTCA-3′.
Western blot analysis
Cell lysis and Western blot analysis were carried out according to standard protocols. The following dilutions of primary commercial antibodies were used as probes: 1:1000 dilution of anti-PKM1 (Cell Signaling Technology, 7067S), 1:1000 dilution of anti-PKM2 (Cell Signaling 4053S), 1:1000 dilution of β-actin (13E5) (Cell Signaling Technology, 4970S), 1:10000 dilution of anti-vinculin (E1E9V) (Cell Signaling Technology, 13901S), 1:1000 dilution of anti-GAPDH (Cell Signaling Technology, 5174S), and 1:1000 dilution of anti-PHGDH (Cell Signaling Technology, 13428S). Primary antibodies were diluted in 5% non-fat milk and incubated overnight at 4 °C. Secondary antibodies (Cell Signaling Technology, 7074S) were diluted in 5% non-fat milk at a dilution of 1:1000 and incubated at room temperature for 1 h.
Metabolomic profiling and stable isotope labeling
For metabolite quantification, cells were seeded in triplicates (n = 3) in 6-well plates with DMEM supplemented with 10% FBS and 1% PS. For stable isotope labeling, media was switched to labeling media containing appropriate tracer, 25 mM [U-13C6]-glucose, 5 mM [U-13C5]-glutamine, 2 mM [U-13C3]-serine, 2 mM [U-13C3]-cysteine, 3.5 mM [U-13C11]-tryptophan, 7.2 mM [U-13C4]-threonine, 25 mM [1,2-13C2]-glucose, or 5 mM [5-13C1]-glutamine (all from Cambridge Isotopes Laboratories, Inc). Samples collected at T = 0 (unlabeled), 5 min, 30 min, 60 min, 120 min, and 24 h after starting the experiment. Metabolite extraction was performed as described previously . Protein left from the extraction was dissolved in 0.2 M potassium hydroxide aqueous solution overnight, then quantified using Pierce BCA Protein Assay Kit (Fisher Scientific, PI23225). Dried metabolite extracts were resuspended in HPLC-grade water containing 1 μM 1,4-piperazinediethanesulfonic acid (PIPES; Sigma-Aldrich, P6757) as an internal standard. To normalize sample concentrations, samples were resuspended at volumes corresponding to their protein quantification values. For amino acid analysis, 20 μl of resuspended sample was added to 80 μl methanol and derivatized with 10 μl triethylamine and 2 μl benzylchloroformate. For pyruvate, lactate, and citrate analysis, 20 μl of resuspended sample was added to a mixture of 20 μl of 250 mM 3-nitrophenylhydrazine (3-NPH) in 50% methanol, 20 μl of 150 mM 1-ethyl-3-(3-dimethylaminopropyl) carbodiimide HCl (EDC) in methanol, and 20 μl of 7.5% pyridine in methanol and allowed to react at 30 °C for 30 min . After this reaction, 16 μl of 2 mg/ml butylated hydroxytoluene (BHT) in methanol was quickly added to these solutions, which were then diluted with 104 μl of water. Samples with and without derivatization were transferred to HPLC vials for analysis.
LC-MS/MS analysis was performed with ion-pairing reverse phase chromatography using an Ascentis Express column (C18, 5 cm × 2.1 mm, 2.7 μm, Sigma-Aldrich) for separation and a Waters Xevo TQ-S triple quadrupole mass spectrometer operated in negative mode as mass analyzer. LC parameters were described previously . Peak processing was performed in MAVEN . Isotope labeling data was corrected for the natural abundance of different isotopes using IsoCor . Heat maps were generated using Cluster .
Statistical analysis was performed using a two-tailed Student’s t test.
PKM2 and PKM1/2 knockdown do not decrease pancreatic cancer cell proliferation
To study the role of pyruvate kinase in pancreatic cancer proliferation, we characterized PDAC cell lines (A13M2-1 and A13M13) derived from a KrasG12D/-; p53-/- mouse pancreatic tumor . A13M2-1 contains a doxycycline-inducible hairpin that knocks down ~ 85% of PKM2 expression after 7 days of doxycycline treatment. Knockdown of PKM2 does not arrest proliferation of PDAC cells (Fig. 1a). Consistent with previous reports , knockdown of PKM2 in PDACs induced PKM1 expression (Fig. 1a). Since PKM1 expression may compensate for PKM2 knockdown, we further knocked down both PKM1 and PKM2 expression in the A13M13 cell line, which contains a doxycycline-inducible short hairpin RNA (shRNA) that knocks down ~ 85% of both PKM1 and PKM2 after 7-day doxycycline treatment (Fig. 1b and Additional file 1: Figure S1A). Unexpectedly, knockdown of both M1 and M2 isoforms still does not affect PDAC cell proliferation (Fig. 1b). To ensure that PKL or PKR expression was not induced after PKM1/2 knockdown, we probed the expression of PKL/R in PDAC cells. Western blot results indicated that there is no PKL/R expressed in A13M13 cells (Additional file 1: Figure S1B). We also measured intracellular metabolites levels following knockdown of both M1 and M2 isoforms in A13M13 PDAC cells (Fig. 1c). Knockdown of PKM1/2 greatly elevated upstream intermediates of glycolysis, such as 2/3-phosphoglycerate (2/3-PG) and phosphoenolpyruvate (Fig. 1c). This is consistent with decreased glycolytic flux through PK, indicating that PK is the rate-limiting step in glycolysis during PKM knockdown. Pyruvate levels were unchanged, but downstream intermediate lactate was significantly decreased with PKM1/2 knockdown. We further performed 13C-glucose labeling experiments to determine how glucose is processed through different metabolic pathways. Surprisingly, PDACs with PKM1/2 knockdown are able to generate proportionally similar amounts of 13C-labeled pyruvate (Fig. 1d). Approximately 60% of pyruvate and 90% of lactate and alanine were labeled in PDAC cells with PKM1/2 knockdown, indicating that pyruvate can still be made from glucose despite targeting a key glycolytic enzyme.
Serine biosynthesis pathway is upregulated in PDACs following PKM1/2 knockdown
While PKM1/2 knockdown does not decrease PDAC cell proliferation (Fig. 1b), it does impact cellular metabolism. PDACs displayed different intracellular metabolic profiles upon PKM1/2 knockdown (Fig. 2a). As expected, PKM1/2 knockdown causes an accumulation of glycolytic intermediates upstream of pyruvate, including fructose 1,6-bisphosphate (FBP), 2/3-PG, and PEP. In addition, serine biosynthesis pathway intermediates phosphoserine and serine were vastly elevated. In contrast, many other amino acids were decreased in PKM1/2 knockdown cells (Fig. 2a). Further investigation with 13C-glucose labeling studies showed that glucose flux to serine is upregulated in PDACs with PKM1/2 knockdown (Fig. 2b). Since serine can be converted to pyruvate via the action of serine dehydratase , we postulated that the serine biosynthesis pathway may be one way to circumvent PKM1/2 knockdown in pancreatic cancer cells (Fig. 2c).
PHGDH knockout depletes intracellular serine and decreases pyruvate production
It is possible that cells bypass PKM1/2 knockdown and convert glucose to pyruvate through the serine biosynthesis pathway (Fig. 2c). To explore this possibility, we targeted the serine biosynthesis pathway by deleting phosphoglycerate dehydrogenase (PHGDH), the rate-limiting enzyme in this pathway. PHGDH converts NAD+ and 3-phosphoglycerate to NADH and 3-phosphohydroxypyruvate [44,45,46]. We deleted PHGDH using CRISPR/Cas9 gene editing technology. We successfully generated PHGDH knockout (KO) PDAC populations and single cell clones, as confirmed by Western blotting (Additional file 1: Figure S2F) and DNA sequencing (Additional file 1: Figure S3). We find that PHGDH KO does not affect PDAC population cell proliferation following PKM1/2 knockdown (Fig. 3a and Additional file 1: Figure S2A–E). We also measured the levels of serine in PHGDH KO clone cells. As expected, knocking out PHGDH significantly decreased serine levels (Fig. 3b). We also find that PHGDH KO indeed led to reduced pyruvate levels (Fig. 3b). Interestingly, PHGDH KO resulted in decreased levels of upstream glycolytic intermediates 2/3-PG and PEP (Fig. 3b), suggesting that targeting the serine biosynthesis pathway may downregulate glycolysis. We further performed 13C-glucose labeling and confirmed that knockout of PHGDH abolished serine and glycine production from 13C-glucose as indicated by loss of M3 and M2 labeling in serine and glycine, respectively (Fig. 3c). Additionally, the fraction of M3-labeled pyruvate from 13C-glucose was slightly decreased, indicating that the serine biosynthesis pathway indeed contributes to pyruvate production during PKM1/2 knockdown.
PHGDH/NPL dual knockout does not impact cell growth rate and pyruvate generation
Since 13C-glucose still contributes to a significant fraction of labeled pyruvate in PHGDH KO cells with PKM1/2 knockdown (Fig. 3c), we investigated whether the sialic acid pathway contributes to producing pyruvate from upstream glycolytic intermediates (Fig. 2c). N-acetylneuraminate synthase (NANS) produces N-acetylneuraminic acid (sialic acid) from N-acetylmannosamine and PEP, while N-acetylneuraminate pyruvate lyase (NPL) catalyzes the degradation of sialic acid into N-acetylmannosamine and pyruvate [47, 48]; hence, this pathway could potentially bypass pyruvate kinase to convert glucose-derived phosphoenolpyruvate into pyruvate. To test this possibility, we generated NPL knockout PDAC cells by using CRISPR/Cas9. NPL knockout PDAC clones can proliferate as well as wild-type PDAC cells following PKM1/2 knockdown (Additional file 1: Figure S4A). Since it is possible that cells use the serine biosynthesis pathway for pyruvate generation during NPL knockout, we further generated PHGDH/NPL dual knockout (dual KO) PDAC cells by using CRISPR/Cas9 (Additional file 1: Figure S5), as confirmed by sequencing (Additional file 1: Figure S6). We find that dual KO cells do not suffer loss of proliferative ability with PKM1/2 knockdown (Fig. 4a). We further investigated the metabolic impacts in PHGDH/NPL dual KO cells using 13C-glucose labeling. Glucose can still be converted into pyruvate in wild-type and dual KO cells following PKM1/2 knockdown (Fig. 4b). Although N-acetylneuraminate is labeled by 13C-glucose after 24 h (Additional file 1: Figure S1C), it was not labeled in either wild-type or PHGDH/NPL dual KO cells after 60 min, indicating that 13C-glucose did not go through the sialic acid pathway to generate labeled pyruvate in this experiment.
Other sources of pyruvate
A significant fraction of pyruvate remains unlabeled even after incubation with 13C-glucose for 24 h in both wild-type and PHGDH KO cell lines with or without PKM1/2 knockdown (Fig. 3c). This suggests that pyruvate can be generated from sources other than glucose in these cell lines. Potential pathways that can generate pyruvate during PKM1/2 knockdown are shown in Fig. 5a. Studies of pyruvate metabolism demonstrate that anaplerotic pyruvate entry by pyruvate carboxylase (PC) in the mitochondria will generate oxaloacetate, which can produce phosphoenolpyruvate through mitochondrial phosphoenolpyruvate carboxykinase (PEPCK-M). PEP is then transported out of the mitochondrial matrix by an anion transporter into the cytosol. Cytosolic malate also can be converted to pyruvate through malic enzyme. To determine whether PC, PEPCK-M, or malic enzyme actively contribute to pyruvate production in PDAC cells, we characterized 13C-glutamine labeling in dual KO cells. While a large proportion of citrate (~ 60%) and malate (~ 80%) were labeled from 13C-glutamine, the downstream metabolites of pyruvate—lactate and alanine—remained virtually unlabeled (Fig. 5b). These results suggest that glutamine-derived TCA cycle intermediates are not converted into pyruvate in these cells. We further extended the labeling time and showed that pyruvate, alanine, and lactate are not labeled by 13C-glutamine even after 24 h (Additional file 1: Figure S7A). Pyruvate can also be produced from several other amino acids. Serine can be converted to pyruvate by serine dehydratase. Threonine can be cleaved to yield glycine , which is converted to serine by serine hydroxymethyltransferase (SHMT) and then converted to pyruvate by serine dehydratase. Tryptophan catabolism via the kynurenine pathway yields alanine , which can be converted to pyruvate via transamination. Cysteine may also be catabolized to yield pyruvate and inorganic bisulfite . We further investigated pyruvate generation from amino acids by using 13C-serine, 13C-tryptophan, 13C-threonine, and 13C-cysteine tracers. Surprisingly, we find that 13C-serine has negligible contribution to pyruvate and lactate labeling (Fig. 5c and Additional file 1: Figure S7D). On the other hand, cysteine contributes to an unexpectedly large proportion (~ 20%) of pyruvate in dual KO cells (Fig. 5d). Cysteine has not been previously reported as a major source of pyruvate production in PDAC cells. WT PDAC cells also produce labeled pyruvate (~ 10%) from 13C-cysteine (Additional file 1: Figure S8), indicating that these cells generally produce pyruvate from cysteine. The lower production of pyruvate from cysteine in WT cells (~ 10%) compared to PHGDH/NPL dual KO cells (~ 20%) is likely due to the dual KO cells having impaired pyruvate production through other pathways (e.g., serine biosynthesis pathway) and requiring a larger proportion of pyruvate production from cysteine. Other amino acids did not significantly contribute to the production of pyruvate (Fig. 5d).
To better understand the sources of pyruvate, we performed a multi-tracer experiment in PHGDH/NPL dual KO cells. Cells were cultured with 1,2-13C-glucose, 5-13C-glutamine, and U-13C-cysteine for 60 min (Fig. 6). Approximately 40% of pyruvate is generated from glucose, as 20% of labeled pyruvate (M2) and 20% of unlabeled pyruvate (M0) comes from 1,2-13C-glucose. Consistent with the lack of alanine and lactate labeling from 13C-glutamine (Fig. 5b), pyruvate, alanine, and lactate were not labeled by 13C-glutamine in this multi-tracer experiment. Interestingly, 20% of labeled pyruvate (M3) comes from 13C-cysteine.
In this study, we demonstrate that pancreatic cancer cells can rewire their metabolism and continue to proliferate during PKM1/2 knockdown. PKM2 in cancer has been studied since the early twenty-first century, but its roles remain controversial. High expression of PKM2 is correlated with poor prognosis in pancreatic cancer patients [20, 24, 25], and suppression of PKM2 expression resulted in decreased cell survival [20, 52]. On the contrary, our results, like other published studies investigating PKM2 function in pancreatic cancer [53, 54], showed that knockdown of PKM2 expression has no effect on PDAC cell proliferation. Similar results were observed in other cancer models, including breast cancer , hepatocellular carcinoma , colon cancer , and leukemia . Differences between our results and published reports may be due to incomplete knockout of PKM1/2 in our cells; despite ~ 85% depletion in PKM1/2 expression, a small fraction of PKM1/2 may be enough for maintaining pyruvate production. Complete knockout of PKM1/2 did not yield any surviving clones in multiple experiments, indicating a small amount of PKM1/2 is likely essential for in vitro survival of these cells. The impact of pyruvate kinase on cancer cell proliferation is likely dependent on the context of each cancer cell, including genetic makeup, mutations, and the microenvironment.
We find that downregulation of PKM1/2 expression does not impact pyruvate production in PDAC cells. However, we do observe rewiring of metabolism following PKM1/2 knockdown: there is accumulation of upstream glycolytic intermediates, consistent with previous reports [30, 52, 56]. We further find that PKM1/2 knockdown causes an increase in serine biosynthesis intermediates. This may be due to an accumulation of upstream glycolytic intermediates causing increased flux through the serine biosynthesis pathway. Knockout of PHGDH, the rate-limiting enzyme in the serine biosynthesis pathway, decreases pyruvate production as indicated by both pyruvate levels (Fig. 3b) and pyruvate labeling from 13C-glucose (Fig. 3c). Thus, our data suggests that the serine biosynthesis pathway contributes to the production of pyruvate from glucose during PKM1/2 knockdown. However, since addition of 13C-serine to the media does not result in pyruvate labeling (Fig. 5c), serine imported from the media does not appear to be directly converted into pyruvate. It is possible that flux through the serine biosynthesis pathway may have a regulatory effect on pyruvate production, since serine is an allosteric activator of PKM2 . Despite reports that suppression of PHGDH can impair cancer cell proliferation [45, 46, 58], we find that PHGDH deletion does not attenuate cell proliferation in PDAC cells. This may be related to the fact that 13C-glucose-derived serine comprises only a small fraction of overall serine; most serine remained unlabeled after 24 h (Fig. 3c) and presumably originated from the culture media. This data agrees with another study that shows the majority of intracellular serine is obtained extracellularly in T cells . Finally, we note that PHGDH KO results in accumulation of nucleotide intermediates (PRPP, UMP, UTP, UDP, dUMP, dTMP, dTTP, and ATP), TCA cycle intermediates (α-ketoglutarate, fumarate, and citrate/isocitrate), and many amino acids (Additional file 1: Figure S9), indicating that PHGDH deletion causes glucose flux to be shunted to nucleotides, the TCA cycle, and amino acid biosynthesis.
During PKM1/2 knockdown, cell growth is supported in part by alternative pathways that generate pyruvate, including breakdown of exogenous amino acids. We explored the breakdown of cysteine, serine, glutamine, tryptophan, and threonine into pyruvate. Cysteine is a key sulfur-containing semi-essential amino acid which plays important functions in redox homeostasis, protein function, and metabolism . Interestingly, our results show that about 20% of pyruvate came from 13C-cysteine (Fig. 6). This is an unexpectedly high fraction, since cysteine has not been previously reported as a major source of intracellular pyruvate in PDAC cells. WT PDAC cells also generate pyruvate from 13C-cysteine (Additional file 1: Figure S8), indicating that this is a general pathway of cysteine production in these cells. The higher fraction of cysteine-derived pyruvate in PHGDH/NPL dual KO cells (~20%; Fig. 6) compared to WT cells (~ 10%; Additional file 1: Figure S8) may be due to the dual KO cells having impaired pyruvate production through other pathways (e.g., serine biosynthesis pathway) and hence requiring a larger proportion of pyruvate generation from cysteine. We also observed a trend toward decreased cysteine pools in the WT cells upon PKM knockdown (Fig. 2a) which may be indicative of cysteine utilization to produce pyruvate, supporting the notion that cysteine is a major source of pyruvate in PDAC cells. The labeled pyruvate from 13C-cysteine is not further converted to lactate or alanine, as indicated by the absence of M3 lactate labeling in either metabolite (Fig. 6). Since lactate is produced in the cytosol, and alanine is predominantly produced in the mitochondria, this data may suggest that cysteine is converted to pyruvate within another cellular compartment that is distinct from either the cytosol or mitochondria. The dependence on cysteine is a potential metabolic adaption for pyruvate generation in PDAC cells. Glutamine is another potential source of pyruvate that we investigated. As the most abundant amino acid in plasma, upregulated glutamine metabolism has been found in several cancer types . Once transported into mitochondria, glutamine serves as a carbon source for the production of fatty acids, nucleotides, and TCA cycle intermediates, and up to 60% of glutamine has been reported to be used for generating lactate and alanine. In our study, glutamine does not contribute to pyruvate generation in PDACs (Fig. 5b and Additional file 1: Figure S7A), despite their dependence on glutamine for proliferation (Additional file 1: Figure S10). Tryptophan and threonine can also be converted to pyruvate, but this was not observed in our cells (Fig. 5d, Additional file 1: Figure S7B and S7C). Thus, our results show that ~ 40% of pyruvate is generated from glucose and ~ 20% of pyruvate is generated from cysteine. For our labeling experiments, we used DMEM formulated without pyruvate and dialyzed FBS that does not contain low molecular weight compounds like amino acids. We further checked the media by LCMS and confirmed that both dialyzed FBS and DMEM contain negligible amounts of pyruvate (Additional file 1: Figure S11). It is possible that the remaining unlabeled portion of pyruvate is produced from amino acids released by degradation of serum proteins in the media, as macropinocytosis and extracellular protein scavenging has been observed in pancreatic cancer cells [61,62,63]. Since the cells produce pyruvate from cysteine, cystine in the media—which can be imported via xCT and subsequently converted to cysteine intracellularly—may also contribute to the unlabeled portion of pyruvate.
Pancreatic cancer cells rewire their metabolism during pyruvate kinase knockdown. The serine biosynthesis pathway enables some conversion of glucose to pyruvate during pyruvate kinase knockdown; however, since direct conversion of serine to pyruvate was not observed, it appears that the serine pathway contributes to pyruvate generation through a potentially regulatory mechanism. We also find that a surprisingly large percentage of intracellular pyruvate comes from cysteine. Our study reveals the resilience of PDAC cells to pyruvate kinase knockdown, underscoring the metabolic flexibility of these cells to overcome environmental perturbation and maintain proliferation.
Availability of data and materials
All data generated or analyzed during this study are included in this article.
1-Ethyl-3-(3-dimethylaminopropyl) carbodiimide HCl
Liquid chromatography tandem mass spectrometry
N-Acetylneuraminate pyruvate lyase
Pancreatic ductal adenocarcinoma
Mitochondrial phosphoenolpyruvate carboxykinase
Piperazine-N,N-bis (2-ethanesulfonic acid)
Pyruvate kinase M2 isoform
Penicillin and streptomycin
5-year relative survival, 2008-2014, cancer statistic center, American Cancer Society. https://cancerstatisticscenter.cancer.org/#/
Siegel RL, Miller KD, Jemal A. Cancer statistics, 2017. CA Cancer J Clin. 2017;67(1):7–30.
Conroy T, Bachet JB, Ayav A, Huguet F, Lambert A, Caramella C, Marechal R, Van Laethem JL, Ducreux M. Current standards and new innovative approaches for treatment of pancreatic cancer. Eur J Cancer. 2016;57:10–22.
Kleeff J, Korc M, Apte M, La Vecchia C, Johnson CD, Biankin AV, Neale RE, Tempero M, Tuveson DA, Hruban RH, Neoptolemos JP. Pancreatic cancer. Nat Rev Dis Primers. 2016;2:16022.
Hanahan D, Weinberg RA. Hallmarks of cancer: the next generation. Cell. 2011;144(5):646–74.
Otto W. On the origin of cancer cells. Science. 1956;123(3191):309–14.
Lunt SY, Vander Heiden MG. Aerobic glycolysis: meeting the metabolic requirements of cell proliferation. Annu Rev Cell Dev Biol. 2011;27:441–64.
Falasca M, Kim M, Casari I. Pancreatic cancer: current research and future directions. Biochim Biophys Acta. 2016;1865(2):123–32.
Perera RM, Bardeesy N. Pancreatic cancer metabolism: breaking it down to build it back up. Cancer Discov. 2015;5(12):124712–61.
Cohen RNC, Tijeras-Raballand A, Faivre S, de Gramont A, Raymond E. Targeting cancer cell metabolism in pancreatic adenocarcinoma. Oncotarget. 2015;6(19):16832–47.
Hosein AN, Beg MS. Pancreatic cancer metabolism: molecular mechanisms and clinical applications. Curr Oncol Rep. 2018;20(7):56.
Jurica MSMA, Heath PJ, Shi W, Nowak T, Stoddard BL. The allosteric regulation of pyruvate kinase by fructose-1,6-bisphosphate. Structure. 1998;6(2):195–210.
Dayton TL, Jacks T, Vander Heiden MG. PKM2, cancer metabolism, and the road ahead. EMBO Rep. 2016;17(12):1721–30.
Takenaka MYK, Lu T, Kang R, Tanaka T, Noguchi T. Alternative splicing of the pyruvate kinase M gene in a minigene system. Eur J Biochem. 1996;235(1-2):366–71.
Osterman JFP. Pyruvate kinase isozymes- a comparative study in tissues of various mammalian species. Comp Biochem Physiol B. 1973;44(4):1077–85.
Marie JLM, Simon MP, Kahn A. Genetic and epigenetic control of the pyruvate kinase isozymes in mammals. Isozymes Curr Top Biol Med Res. 1983;7:221–40.
Sato T, Morita M, Nomura M, Tanuma N. Revisiting glucose metabolism in cancer: lessons from a PKM knock-in model. Mol Cell Oncol. 2018;5(4):e1472054.
Lunt SY, Muralidhar V, Hosios AM, Israelsen WJ, Gui DY, Newhouse L, Ogrodzinski M, Hecht V, Xu K, Acevedo PN, Hollern DP, Bellinger G, Dayton TL, Christen S, Elia I, Dinh AT, Stephanopoulos G, Manalis SR, Yaffe MB, Andrechek ER, Fendt SM, Vander Heiden MG. Pyruvate kinase isoform expression alters nucleotide synthesis to impact cell proliferation. Mol Cell. 2015;57(1):95–107.
Bandara IA, Baltatzis M, Sanyal S, Siriwardena AK. Evaluation of tumor M2-pyruvate kinase (Tumor M2-PK) as a biomarker for pancreatic cancer. World J Surg Oncol. 2018;16(1):56.
Li C, Zhao Z, Zhou Z, Liu R. PKM2 promotes cell survival and invasion under metabolic stress by enhancing warburg effect in pancreatic ductal adenocarcinoma. Dig Dis Sci. 2016;61(3):767–73.
Feng J, Ma T, Ge Z, Lin J, Ding W, Chen H, Zhu W, Zhou S, Tan Y. PKM2 gene regulates the behavior of pancreatic cancer cells via mitogen-activated protein kinase pathways. Mol Med Rep. 2015;11(3):2111–7.
Azoitei N, Becher A, Steinestel K, Rouhi A, Diepold K, Genze F, Simmet T, Seufferlein T. PKM2 promotes tumor angiogenesis by regulating HIF-1alpha through NF-kappaB activation. Mol Cancer. 2016;15:3.
Cheng TY, Yang YC, Wang HP, Tien YW, Shun CT, Huang HY, Hsiao M, Hua KT. Pyruvate kinase M2 promotes pancreatic ductal adenocarcinoma invasion and metastasis through phosphorylation and stabilization of PAK2 protein. Oncogene. 2018;37(13):1730–42.
Ogawa H, Nagano H, Konno M, Eguchi H, Koseki J, Kawamoto K, Nishida N, Colvin H, Tomokuni A, Tomimaru Y, Hama N, Wada H, Marubashi S, Kobayashi S, Mori M, Doki Y, Ishii H. The combination of the expression of hexokinase 2 and pyruvate kinase M2 is a prognostic marker in patients with pancreatic cancer. Mol Clin Oncol. 2015;3(3):563–71.
Mohammad GH, Olde Damink SW, Malago M, Dhar DK, Pereira SP. Pyruvate kinase M2 and lactate dehydrogenase A are overexpressed in pancreatic cancer and correlate with poor outcome. PLoS One. 2016;11(3):e0151635.
Lockney NA, Zhang M, Lu Y, Sopha SC, Washington MK, Merchant N, Zhao Z, Shyr Y, Chakravarthy AB, Xia F. Pyruvate kinase muscle isoenzyme 2 (PKM2) expression is associated with overall survival in pancreatic ductal adenocarcinoma. J Gastrointest Cancer. 2015;46(4):390–8.
Zhu HLH, Zhu X, Hu X, Zheng L, Zhu X. Pyruvate kinase M2 (PKM2) expression correlate with prognosis in solid cancers: a meta-analysis. Oncotarget. 2017;8(1):1628–40.
Wu J, Hu L, Chen M, Cao W, Chen H, He T. Pyruvate kinase M2 overexpression and poor prognosis in solid tumors of digestive system: evidence from 16 cohort studies. Onco Targets Ther. 2016;9:4277–88.
Morita M, Sato T, Nomura M, Sakamoto Y, Inoue Y, Tanaka R, Ito S, Kurosawa K, Yamaguchi K, Sugiura Y, Takizaki H, Yamashita Y, Katakura R, Sato I, Kawai M, Okada Y, Watanabe H, Kondoh G, Matsumoto S, Kishimoto A, Obata M, Matsumoto M, Fukuhara T, Motohashi H, Suematsu M, Komatsu M, Nakayama KI, Watanabe T, Soga T, Shima H, Maemondo M, Tanuma N. PKM1 confers metabolic advantages and promotes cell-autonomous tumor cell growth. Cancer Cell. 2018;33(3):355–67 e7.
Cortes-Cros M, Hemmerlin C, Ferretti S, Zhang J, Gounarides JS, Yin H, Muller A, Haberkorn A, Chene P, Sellers WR, Hofmann F. M2 isoform of pyruvate kinase is dispensable for tumor maintenance and growth. Proc Natl Acad Sci USA. 2013;110(2):489–94.
Wang YH, Israelsen WJ, Lee D, Yu VWC, Jeanson NT, Clish CB, Cantley LC, Vander Heiden MG, Scadden DT. Cell-state-specific metabolic dependency in hematopoiesis and leukemogenesis. Cell. 2014;158(6):1309–23.
Dayton TL, Gocheva V, Miller KM, Israelsen WJ, Bhutkar A, Clish CB, Davidson SM, Luengo A, Bronson RT, Jacks T, Vander Heiden MG. Germline loss of PKM2 promotes metabolic distress and hepatocellular carcinoma. Genes Dev. 2016;30(9):1020–33.
Lau AN, Israelsen WJ, Roper J, Sinnamon MJ, Georgeon L, Dayton TL, Hillis AL, Yilmaz OH, Di Vizio D, Hung KE, Vander Heiden MG. PKM2 is not required for colon cancer initiated by APC loss. Cancer Metab. 2017;5:10.
Hillis AL, Lau AN, Devoe CX, Dayton TL, Danai LV, Di Vizio D, Vander Heiden MG. PKM2 is not required for pancreatic ductal adenocarcinoma. Cancer Metab. 2018;6:17.
Gidekel Friedlander SY, Chu GC, Snyder EL, Girnius N, Dibelius G, Crowley D, Vasile E, DePinho RA, Jacks T. Context-dependent transformation of adult pancreatic cells by oncogenic K-Ras. Cancer Cell. 2009;16(5):379–89.
Muzumdar MD, Chen PY, Dorans KJ, Chung KM, Bhutkar A, Hong E, Noll EM, Sprick MR, Trumpp A, Jacks T. Survival of pancreatic cancer cells lacking KRAS function. Nat Commun. 2017;8(1):1090.
Ran FA, Hsu PD, Wright J, Agarwala V, Scott DA, Zhang F. Genome engineering using the CRISPR-Cas9 system. Nat Protoc. 2013;8(11):2281–308.
Ogrodzinski MP, Teoh ST, Yu L, Broadwater D, Ensink E, Lunt SY. Measuring the nutrient metabolism of adherent cells in culture. In: Fendt SM, Lunt S, editors. Methods Mol Biol, vol. 1862. New York: Metabolic Signaling. Methods Mol Biol: Humana Press; 2019. p. 37–52.
Han J, Gagnon S, Eckle T, Borchers CH. Metabolomic analysis of key central carbon metabolism carboxylic acids as their 3-nitrophenylhydrazones by UPLC/ESI-MS. Electrophoresis. 2013;34(19):2891–900.
Eugene Melamud LV, Joshua D. Rabinowitz. Metabolomic analysis and visualization engine for LC−MS data. Anal Chem. 2010;82(23):9819–26.
Millard P, Letisse F, Sokol S, Portais JC. IsoCor: correcting MS data in isotope labeling experiments. Bioinformatics. 2012;28(9):1294–6.
Yamada TKJ, Takata Y, Ogawa H, Pitot HC, Takusagawa F. Crystal structure of serine dehydratase from rat liver. Biochemistry. 2003;42(44):12854–65.
Rowsell EV, Carnie JA, Wahbi SD, Al-Tai AH, Rowsell KV. l-serine dehydratase and l-serine-pyruvate aminotransferase activities in different animal species. Comp Biochem Physiol B. 1979;63(4):543–55.
Zogg CK. Phosphoglycerate dehydrogenase: potential therapeutic target and putative metabolic oncogene. J Oncol. 2014;2014:524101.
Yoshino H, Nohata N, Miyamoto K, Yonemori M, Sakaguchi T, Sugita S, Itesako T, Kofuji S, Nakagawa M, Dahiya R, Enokida H. PHGDH as a key enzyme for serine biosynthesis in HIF2alpha-targeting therapy for renal cell carcinoma. Cancer Res. 2017;77(22):6321–9.
Zhang B, Zheng A, Hydbring P, Ambroise G, Ouchida AT, Goiny M, Vakifahmetoglu-Norberg H, Norberg E. PHGDH defines a metabolic subtype in lung adenocarcinomas with poor prognosis. Cell Rep. 2017;19(11):2289–303.
Comb DG, Roseman S. The sialic acids. I. The structure and enzymatic synthesis of N-acetylneuraminic acid. J Biol Chem. 1960;235:2529–37.
Ji W, Sun W, Feng J, Song T, Zhang D, Ouyang P, Gu Z, Xie J. Characterization of a novel N-acetylneuraminic acid lyase favoring N-acetylneuraminic acid synthesis. Sci Rep. 2015;5:9341.
Schirch LV, Gross T. Serine transhydroxymethylase identification as the threonine and allothreonine aldolases. J Biol Chem. 1968;243(21):5651–5.
Lima S, Kristoforov R, Momany C, Phillips RS. Crystal structure of homo sapiens kynureninase. Biochemistry. 2007;46(10):2735–44.
Stipanuk MH, Dominy JE, Lee J-I, Coloso RM. Mammalian cysteine metabolism: new insights into regulation of cysteine metabolism. J Nutr. 2006;136(6):1652S–9S.
Yokoyama M, Tanuma N, Shibuya R, Shiroki T, Abue M, Yamamoto K, Miura K, Yamaguchi K, Sato I, Tamai K, Satoh K. Pyruvate kinase type M2 contributes to the development of pancreatic ductal adenocarcinoma by regulating the production of metabolites and reactive oxygen species. Int J Oncol. 2018;52(3):881–91.
Hillis AL, Lau AN, Devoe CX, Dayton TL, Danai LV, Di Vizio D, Vander Heiden MG. PKM2 is not required for pancreatic ductal adenocarcinoma. Cancer Metabol. 2018;6:1.
Li X, Deng S, Liu M, Jin Y, Zhu S, Deng S, Chen J, He C, Qin Q, Wang C, Zhao G. The responsively decreased PKM2 facilitates the survival of pancreatic cancer cells in hypoglucose. Cell Death Dis. 2018;9(2):133.
Israelsen WJ, Dayton TL, Davidson SM, Fiske BP, Hosios AM, Bellinger G, Li J, Yu Y, Sasaki M, Horner JW, Burga LN, Xie J, Jurczak MJ, DePinho RA, Clish CB, Jacks T, Kibbey RG, Wulf GM, Di Vizio D, Mills GB, Cantley LC, Vander Heiden MG. PKM2 isoform-specific deletion reveals a differential requirement for pyruvate kinase in tumor cells. Cell. 2013;155(2):397–409.
Christofk HR, Vander Heiden MG, Harris MH, Ramanathan A, Gerszten RE, Wei R, Fleming MD, Schreiber SL, Cantley LC. The M2 splice isoform of pyruvate kinase is important for cancer metabolism and tumour growth. Nature. 2008;452(7184):230–3.
Chaneton B, Hillmann P, Zheng L, Martin ACL, Maddocks ODK, Chokkathukalam A, Coyle JE, Jankevics A, Holding FP, Vousden KH, Frezza C, O’Reilly M, Gottlieb E. Serine is a natural ligand and allosteric activator of pyruvate kinase M2. Nature. 2012;491(7424):458–62.
Ma EH, Bantug G, Griss T, Condotta S, Johnson RM, Samborska B, Mainolfi N, Suri V, Guak H, Balmer ML, Verway MJ, Raissi TC, Tsui H, Boukhaled G, Henriques da Costa S, Frezza C, Krawczyk CM, Friedman A, Manfredi M, Richer MJ, Hess C, Jones RG. Serine is an essential metabolite for effector T cell expansion. Cell Metab. 2017;25(2):345–57.
Locasale JW, Grassian AR, Melman T, Lyssiotis CA, Mattaini KR, Bass AJ, Heffron G, Metallo CM, Muranen T, Sharfi H, Sasaki AT, Anastasiou D, Mullarky E, Vokes NI, Sasaki M, Beroukhim R, Stephanopoulos G, Ligon AH, Meyerson M, Richardson AL, Chin L, Wagner G, Asara JM, Brugge JS, Cantley LC, Vander Heiden MG. Phosphoglycerate dehydrogenase diverts glycolytic flux and contributes to oncogenesis. Nat Genet. 2011;43(9):869–74.
Kohl JB, Mellis AT, Schwarz G. Homeostatic impact of sulfite and hydrogen sulfide on cysteine catabolism. Br J Pharmacol. 2019;176(4):554–70.
Yao W, Rose JL, Wang W, Seth S, Jiang H, Taguchi A, Liu J, Yan L, Kapoor A, Hou P, Chen Z, Wang Q, Nezi L, Xu Z, Yao J, Hu B, Pettazzoni PF, Ho IL, Feng N, Ramamoorthy V, Jiang S, Deng P, Ma GJ, Den P, Tan Z, Zhang SX, Wang H, Wang YA, Deem AK, Fleming JB, Carugo A, Heffernan TP, Maitra A, Viale A, Ying H, Hanash S, DePinho RA, Draetta GF. Syndecan 1 is a critical mediator of macropinocytosis in pancreatic cancer. Nature. 2019;568:410–4.
Kamphorst JJ, Nofal M, Commisso C, Hackett SR, Lu W, Grabocka E, Vander Heiden MG, Miller G, Drebin JA, Bar-Sagi D, et al. Human pancreatic cancer tumors are nutrient poor and tumor cells actively scavenge extracellular protein. Cancer Res. 2015;75(3):544–53.
Commisso C, Davidson SM, Soydaner-Azeloglu RG, Parker SJ, Kamphorst JJ, Hackett S, Grabocka E, Nofal M, Drebin JA, Thompson CB, et al. Macropinocytosis of protein is an amino acid supply route in Ras-transformed cells. Nature. 2013;497(7451):633–7.
We thank Talya Dayton and Tyler Jacks for sharing the A13M2-1 and M13 cell lines. We thank the Michigan State University RTSF Mass Spectrometry and Metabolomics Core for helpful technical assistance. We also thank Deanna Broadwater for critical reading of this manuscript.
This work is supported by the AACR-Incyte Corporation NextGen Grant for Transformative Cancer Research, Grant Number 16-20-46- LUNT, and Mayo Clinic Metabolomics Resource Core through grant number U24DK100469 from the National Institute of Diabetes and Digestive and Kidney Diseases and originates from the National Institutes of Health Director's Common Fund.
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
About this article
Cite this article
Yu, L., Teoh, S.T., Ensink, E. et al. Cysteine catabolism and the serine biosynthesis pathway support pyruvate production during pyruvate kinase knockdown in pancreatic cancer cells. Cancer Metab 7, 13 (2019). https://doi.org/10.1186/s40170-019-0205-z
- Pyruvate kinase
- Pancreatic cancer
- Liquid chromatography mass spectrometry