Allosteric activation of the metabolic enzyme GPD1 inhibits bladder cancer growth via the lysoPC-PAFR-TRPV2 axis
Journal of Hematology & Oncology volume 15, Article number: 93 (2022)
Bladder cancer is the most common malignant tumor of the urinary system. Surgical resection and chemotherapy are the two mainstream treatments for bladder cancer. However, the outcomes are not satisfactory for patients with advanced bladder cancer. There is a need to further explore more effective targeted therapeutic strategies.
Proteomics were performed to compare protein expression differences between human bladder cancer tissues and adjacent normal tissues. The function of GPD1 on bladder cancer cells were confirmed through in vivo and in vitro assays. Transcriptomics and metabolomics were performed to reveal the underlying mechanisms of GPD1. Virtual screening was used to identify allosteric activator of GPD1.
Here, we used proteomics to find that GPD1 expression was at low levels in bladder cancer tissues. Further investigation showed that GPD1 overexpression significantly promoted apoptosis in bladder cancer cells. Based on transcriptomics and metabolomics, GPD1 promotes Ca2+ influx and apoptosis of tumor cells via the lysoPC-PAFR-TRPV2 axis. Finally, we performed a virtual screening to obtain the GPD1 allosteric activator wedelolactone and demonstrated its ability to inhibit bladder tumor growth in vitro and in vivo.
This study suggests that GPD1 may act as a novel tumor suppressor in bladder cancer. Pharmacological activation of GPD1 is a potential therapeutic approach for bladder cancer.
Bladder cancer, as one of the most common malignant tumors of the urinary system, ranks 4th in incidence among male malignant tumors . More than 400,000 new cases occur each year worldwide, with over 165,000 deaths associated with bladder cancer [2, 3]. The latest data for 2020 show that 573,000 new bladder cancer cases and 213,000 deaths occurred worldwide . In China, the incidence and mortality rate of bladder cancer are increasing gradually . Based on the depth of bladder wall infiltration, bladder cancer is clinically classified into muscle-invasive bladder cancer (MIBC) and nonmuscle-invasive bladder cancer (NMIBC). Of these, NMIBC is usually multifocal, accounting for approximately 80%, while MIBC accounts for approximately 20% and is characterized by a high incidence of distant metastases even after total cystectomy and systemic chemotherapy, with a 5-year survival rate of less than 5% . Clinical management of bladder cancer has barely progressed in the last 30 years. Surgical resection and chemotherapy are the two mainstream treatments for bladder cancer. Gemcitabine in combination with cisplatin is the first-line chemotherapy regimen for bladder cancer . However, patients with advanced bladder cancer treated with this chemotherapy regimen also had only a 40–60% objective remission rate and only a 5% increase in overall survival [8, 9]. Therefore, understanding the molecular mechanisms underlying bladder cancer progression is essential for exploring more effective targeted therapies.
Glycerol 3-phosphate dehydrogenase 1 (GPD1) catalyzes the conversion of dihydroxyacetone phosphate (DHAP) from glucose and NADH to glycerol-3-phosphate (G3P) and NAD+ . The product glycerol-3-phosphate and glycerol is the backbone for lipid biosynthesis . GPD1, together with an isoenzyme in mitochondria named GPD2, facilitates the glycerol-3-phosphate shuttle pathway, which plays an important role in transporting reducing equivalents from the cytoplasm to the mitochondria . GPD1 is widely distributed in tissues, with the highest levels in subcutaneous fat, duodenum and mesenteric fat [13, 14]. The abnormal expression of GPD1 may negatively affect human growth. Mutant GPD1 overexpression causes transient hypertriglyceridemia and fatty liver in infancy . It has been reported that GPD1-deficient mice exhibit enhanced locomotor activity by increasing lipid oxidation , while they have reduced obesity and body weight . GPD1 has been identified as a tumor suppressor in breast cancer [17, 18] and exerts antitumor effects in lung and prostate cancers [19, 20]. It has been demonstrated that GPD1 exhibits extremely low expression levels or absence in some tumor tissues [21, 22]. However, GPD1 acts as a brain tumor stem cell marker to promote the progression of glioblastoma . The underlying mechanisms of GPD1 in cancer are still not well studied, especially in human bladder cancer.
Here, we performed proteomics to compare protein expression differences between human bladder cancer tissues and adjacent normal tissues. Many differentially expressed proteins were identified, among which GPD1 expression was significantly decreased in bladder cancer and further validated that GPD1 overexpression showed a suppressive effect on tumor growth in vivo and in vitro, which suggests that GPD1 may be a novel therapeutic target for bladder cancer. Based on virtual screening, we report the identification of wedelolactone (WE), the first reported selective and cell-active activator of GPD1 in our knowledge. WE is a well-characterized allosteric GPD1 activator and a promising small molecule to inhibit bladder cancer growth.
Cell lines and cell culture
5637 and T24 cells were cultured in RPMI 1640 medium, J82 and UMUC3 cells were cultured in DMEM, and SV-HUC-1 cells were cultured in F12K medium. All media were supplemented with 10% fetal bovine serum (FBS), 2 mM l-glutamine, and 100 units/ml penicillin with 100 units/ml streptomycin. All cells were cultured in a humidified atmosphere containing 5% CO2.
Label-free proteomic profiling analysis
All tissues were collected from patients undergoing surgery at the Affiliated Drum Tower Hospital of Medical School of Nanjing University (Nanjing, China). The study was approved by Ethical Committee of Nanjing Drum Tower Hospital, Medical School of Nanjing University with written consent provided by all patients (project number:2018-015-01). Combine the diagnostic imaging features to determine the location of the tumor site. Select intact, erythematous tumors for tumor samples. Take the tissue at least 2 cm away from the tumor site as the normal samples, and the pathologist confirmed that the sample was cancer-free by fast freezing pathology. Tumor samples and para-tumor samples were paired. Additional file 1: Table S1 summarizes the characteristics of the patients. The whole process of label-free proteomics was performed by Jingjie PTM Biolab Co., Ltd. (Hangzhou, China).
Five-micrometer-thick sections were obtained from paraffin-embedded tissue. Sections were deparaffinized, rehydrated, and endogenous peroxidase activity was eliminated. After washing in PBS buffer (pH 7.4), sections were blocked with 10% goat serum and then incubated with the primary antibody for GPD1(1:200 dilution, 13451-1-AP, Proteintech, China) overnight at 4 °C, followed by incubation with secondary antibodies for 45 min at room temperature according to the manufacturer's instructions. Staining was visualized, and each section was reacted with diaminobenzidine and then counterstained with hematoxylin. Negative control was performed by staining with only secondary antibodies.
Apoptotic cells in tumor tissues were analyzed by terminal deoxynucleotidyl transferase-mediated UTP nick end labeling (TUNEL) staining. Apoptotic cell nuclei were stained with green fluorescence according to the manufacturer's (TUNEL FITC Apoptosis Detection Kit, Vazyme, Jiangsu, China) protocol. Cell nuclei were stained with DAPI.
Cells were collected, and proteins were extracted by whole cell lysis containing protease and phosphatase inhibitors. Cell debris was removed by centrifugation at 4 °C, and the protein concentration was determined by Pierce BCA assay. Protein content was electrophoresed on 10% SDS–PAGE gels followed by immunoblotting on polyvinylidene difluoride membranes (Biosciences, USA). Antibodies against GPD1 (13451-1-AP, Proteintech) and β-actin (700068, ZenBio) were used.
Total RNA from cell or tissue samples was prepared using TRIzol reagent (Invitrogen). Total RNA (1 μg) was reverse transcribed using HiScript Q-RT SuperMix (TaKaRa Biotech). β-actin was used as a normalization gene. Q-PCR assays were performed using a Q-PCR kit (Vazyme Biotech) with QuantStudio 6-Flex (Applied Biosystems, Foster City, CA, USA). The results were expressed as fold change using the comparative threshold method.
The cDNA library was constructed by Biomarker Technologies Corporation (Beijing, China). The constructed cDNA libraries were sequenced using an Illumina HiSeq™ sequencing platform. The constructed cDNA libraries of were sequenced on a flow cell using an Illumina HiSeq™ sequencing platform. DESeq and Q value were employed and used to evaluate differential gene expression. The false discovery rate (FDR) control method was used to identify the threshold of the P value in multiple tests in order to compute the significance of the differences. Here, only gene with an absolute value of log2 ratio ≥ 2 and FDR significance score < 0.01 was used for subsequent analysis. The raw reads were deposited into NCBI Sequence Read Archive (SRA) database with accession number of PRJNA849215.
Untargeted metabolomics based on ultra-high-performance liquid chromatography quadrupole time of flight mass spectrometry (UHPLC-QTOF-MS)
The LC/MS system was used for metabolomic analysis. Positive ion mode: 0.1% formic acid aqueous solution and 0.1% formic acid acetonitrile. Negative ion mode: 0.1% formic acid aqueous solution and 0.1% formic acid acetonitrile. The raw data were collected using MassLynx V4.2 and then processed by Progenesis QI software based on the Progenesis QI software online METLIN database and Biomark’s self-built library for identification, and at the same time, theoretical fragment identification and mass deviation were all within 100 ppm.
Flow cytometry staining
5637 cells or T24 cells were incubated with anti-TRPV2 antibody (1:100 dilution, ACC-139, Alomone Labs) at 4 °C for 30 min. Subsequently, cells were washed and incubated with Alexa Fluor™ 594-labeled goat anti-rabbit IgG (H + L) (1:100 dilution, A-11037, Invitrogen, CA, USA) for 1 h. After staining, cells were washed and suspended in PBS for flow cytometric analysis. For detection of intracellular Ca2+ concentration in 5637 cells or T24 cells by flow cytometry, Fluo-3 AM (S1056, Beyotime, Jiangsu, China) was used according to the manufacturer's instructions.
Cell counting kit-8 (CCK-8) assay
Cell viability was analyzed using the CCK-8 assay. A total of 5637 or T24 cells (5000 cells per well) were inoculated in 96-well plates. At the indicated time points, the cells were incubated with 10 μL CCK-8 solution for 2 h at 37 °C, followed by measurement of absorbance at 452 nm. All experiments were performed in triplicate.
Measurement of G3P and NAD+
G3P levels were measured using Glycerol-3-phosphate (G3P) Assay Kit (ab174094, Abcam, Cambridge, MA, USA) according to the manufacturer's instructions. NAD+ levels were measured using NAD+/NADH Assay Kit (Beyotime, Haimen, Jiangsu, China) according to the manufacturer's instructions.
Cell apoptosis detection
Apoptotic cells were measured using the membrane-linked protein V/propidium iodide (PI) double staining assay. According to the manufacturer's instructions, 5637 or T24 cells were resuspended in 100 µL of 1× binding buffer and stained with both 5 µL of APC-coupled Annexin V and 5 µL of PI. Cells were then incubated in the dark for 15 min and subjected to flow cytometry analysis.
Transwell migration assay
The assay was performed using the method described in the previous study .
Tumor sphere formation assay
For the sphere formation assay of 5637 and T24 cells, 5000 cells/well were cultured in serum-free DMEM/F-12 (1:1 ratio) medium supplemented with B-27 (Gibco, USA), 10 ng/mL bFGF, 20 ng/mL EGF, and N2 supplement (Thermo Fisher Scientific, USA) in ultralow attachment 6-well plates to form tumor spheres. Images were collected with light microscopy.
Colony formation assay
T24 or 5637 cells were seeded into 12-well plates (300 cells per well) and cultured for 10 days. The cells were then fixed in 10% formaldehyde for 15 min and subsequently stained with 1.0% crystal violet for 20 min. The number of colonies formed was counted by ImageJ.
Cell cycle detection
Cells were collected and fixed in 70% cold ethanol at 4 °C overnight. The ethanol-fixed cells were washed again with PBS, incubated with 50 μg/mL PI containing 20 μg/mL RNase A for 30 min at room temperature, and subsequently analyzed using flow cytometry.
TRPV2 siRNA interference assay
TRPV2 siRNA sequences were synthesized as follows by GenePharma (Shanghai, China). Synthetic scrambled siRNAs were used as negative controls. Transfections were performed using Lipofectamine 3000 (Invitrogen). Total RNA and proteins were collected 48 h post-transfection.
Tumor models and studies
Six-week-old BALB/c Nude mice were purchased from GemPharmatech Co., Ltd. (Nanjing, China). All animals were assessed to be healthy and free of disease before tumor implantation. Animal study was approved by the Institutional Animal Care Committee of Jiangsu Province and the Ethics Committee of Nanjing Drum Tower Hospital, Medical School of Nanjing University (project number: IACUC-2007015). A total of 5 × 106 T24 cells or 5637 cells were injected subcutaneously into nude mice. Mice were killed after 4 weeks. Tumors were collected, and tumor weights were measured. For WE treatment experiments, vehicle or 50 mg/kg WE was administered intraperitoneally every 3 days starting on day 7 after injection.
Prediction of potential allosteric sites of GPD1 and virtual screening
The potential allosteric sites were predicted based on the structure of GPD1 (PDB ID:6E8Y) using the Allosite server (http://mdl.shsmu.edu.cn/AST/). Grid-based ligand docking from energetics (GLIDE) software (Schrödinger Maestro 11.4) was used to virtually screen three commercial compound libraries containing over 65,000 compounds to target the predicted sites. The 10 highest scoring compounds were purchased from MedChemExpress for experimental use.
Microscale thermophoresis (MST)
To determine the binding affinity of GPD1 with WE, microscale thermophoresis (MST) was performed using Monolith NT.115 (NanoTemper Technologies GmbH, Munich, Germany) and standard capillaries. GPD1 recombinant protein was purchased from Proteintech (Cat. Ag4471), and wedelolactone was purchased from MedChemExpress (HY-N0551). The GPD1 recombinant protein was labeled with the fluorescent dye Red-NHS using NanoTemper’s labeling kit MO-L011 Monolith™ according to the manufacturer's instructions. A fluorophore/aptamer ratio of 0.5 was determined by measuring absorbance at 260 and 546 nm using a Nanophotometer™ (Implen GmbH, Munich, Germany). Incubate 200 nM labeled GPD1 protein with 1 µM to 62.5 µM WE. Analysis was performed in binding buffer including 0.05% Tween 20 at 25 °C. The fluorescence signal was normalized, and the Hill equation was fitted using MO Affinity Analysis v2.1.3 software (NanoTemper).
GPD1 enzyme activity assay
For GPD1 enzyme activity assays, the 100-μL reaction system contained TEA buffer (20 mM), 0.1 mg/mL BSA, 200 μM NADH, 8 mM DHAP, 1 μM GPD1, and compound/DMSO at I = 0.12 (NaCl). The initial rate of DHAP reduction was determined by measuring the change in NADH content via 340 nM absorbance. The absorbance was measured using Tecan infinite 200Pro plate reader (Tecan Ltd., Switzerland).
Gene expression of GPD1 in 31 tumors and their paired normal tissues in the TCGA database were obtained from GEPIA2 (website: gepia2.cancerpku.cn/#analysis).
Data are expressed as the mean ± SEM. When only two value sets were compared, statistical analysis was performed by Student’s t test. When the data involved three or more groups, two-way ANOVA was used. Survival analysis was calculated with Kaplan–Meier method, and differences were compared by the log-rank test. P < 0.05, P < 0.01 or P < 0.001 were considered statistically significant and indicated by *, ** or ***, respectively.
Proteomics reveals low expression of GPD1 in human bladder cancer tissues
To obtain the protein expression profiles of bladder cancer tissues, we collected 17 pairs of bladder cancer tissues and matched adjacent normal tissues to perform differential expression protein analysis using proteomics. A total of 98,211.0 peptides were identified by profiling, of which 92,295.0 were specific peptides. We identified a total of 8286.0 proteins, of which 6700.0 were quantifiable. With 1.5-fold as the differential expression change threshold and t test P value < 0.05 as the significance threshold by statistical test, we then identified 1584 proteins with upregulated expression and 518 proteins with downregulated expression in the cancer tissue among the quantified proteins (Fig. 1A and Additional file 2: Table S2). Among these proteins, GPD1 was the most significantly downregulated (Fig. 1B). We detected the expression of GPD1 in 4 pair normal and tumor tissues by western blotting, and the results showed that GPD1 expression was also significantly reduced in tumor tissue (Fig. 1C, D). In the TCGA database, GPD1 mRNA was shown to be at low levels in multiple tumor types, including breast cancer, lung cancer, and prostate cancer. However, no significant differences in GPD1 mRNA levels were observed between normal and tumor tissues in bladder cancer cohort (Fig. 1E). To determine the tissue expression pattern of GPD1 in human bladder tumors, we generated tissue microarrays containing 76 human bladder cancer samples (T1 stage = 24, T2 stage = 28, T3 stage = 14, T4 stage = 10). Staining of these samples using a validated antibody for GPD1 revealed reduced GPD1 expression in tumor tissue compared to normal tissue, and GPD1 expression gradually decreased with bladder cancer progression (Fig. 1F, G). Additionally, GPD1 expression was lower in the MIBC than in the NMIBC (Fig. 1H). In the MIBC cohort, low mRNA levels of GPD1 were also significantly correlated with decreased survival time for bladder cancer patients (Fig. 1I). We further examined the expression level of GPD1 in human bladder cancer cell lines (5637, T24, J82, and UMUC3) and a human urothelial epithelial cell line (SV-HUC-1). Consistent with our analysis of human bladder cancer tissues, the protein levels but not mRNA levels of GPD1 expression were significantly lower in bladder cancer cells than in SV-HUC-1 normal cells (Fig. 1J–L). Together, these results suggest that GPD1 is expressed at low levels in bladder cancer tissues and may have antitumor functions, which prompted us to further explore the role of GPD1 in the development of bladder cancer.
Overexpression of GPD1 leads to bladder cancer cell apoptosis and inhibits tumor growth
To test our hypothesis that GPD1 may have a cancer cell inhibitory effect, we established bladder cancer cell lines (5637, T24) stably expressing GPD1 via lentivectors transduction. Both constructed cell lines were demonstrated by western blotting (Fig. 2A, B). We observed that GPD1 overexpression could significantly reduce proliferation in both 5637 and T24 cells in vitro by CCK-8 cell assays (Fig. 2C, D) and did not affect the cell cycle (Additional file 5: Fig. S1A, B). The colony formation assay indicated that GPD1 overexpression in 5637 cells and T24 cells significantly decreased the colony number (Fig. 2E, F), additionally supporting the growth inhibitory effect of GPD1 on bladder tumor cells. We further found that GPD1 overexpression significantly promoted apoptosis in 5637 and T24 cells (Fig. 2G, H). Transwell migration and tumor sphere formation assays indicated that the migration (Additional file 5: Fig. S1C, D) and sphere-forming capacities (Additional file 5: Fig. S1E, F) of 5637 and T24 cells were also suppressed by overexpression of GPD1.
To investigate whether overexpression of GPD1 inhibits tumor growth in vivo, 5637 cells and T24 cells stably transfected with GPD1 or control vector were injected subcutaneously into BALB/c Nude mice. Compared with the control group, the GPD1 overexpression group had significantly reduced tumor weight (Fig. 2I, J). We further examined whether GPD1 overexpression was associated with enhanced necrosis or apoptosis of tumor cells in vivo. H&E staining and TUNEL assays showed more tumor necrosis and TUNEL-positive cells in tumor tissue from the GPD1 overexpression group (Fig. 2K, L). These results demonstrate that overexpression of GPD1 efficiently inhibits bladder cancer growth.
GPD1-induced inhibition of bladder cancer depends on its enzymatic activity
GPD1 is able to catalyze DHAP and NADH into G3P and NAD+. In GPD1-overexpressing bladder cancer cell lines (T24 and 5637), we observed an increase in the levels of G3P and NAD+ (Fig. 3A–D). Next, we examined whether increasing the concentration of G3P affects bladder cancer cell growth by CCK-8 assay. 5637 and T24 cells were treated with different dosages of G3P. As shown in Fig. 3E and F, the proliferation of 5637 and T24 cells was inhibited by G3P treatment at increasing concentrations. We found that G3P significantly promoted apoptosis in 5637 cells and T24 cells (Fig. 3G, H). Interestingly, the combination of G3P and NAD+ further enhanced the pro-apoptotic and anti-proliferative effects (Fig. 3I, J). To further determine that the suppression of tumor growth by GPD1 is dependent on its catalytic activity, which can promote the accumulation of G3P and NAD+, we next constructed 5637 and T24 cells overexpressing the GPD1 mutant (K120A variant) (Fig. 3K, L). Mutation of Lys120 in GPD1 led to the loss of catalytic activity , which was further evidenced by the fact that K120A GPD1 overexpression did not promote upregulation of G3P and NAD+ levels (Fig. 3M–P). As expected, the proliferation and apoptosis of 5637 and T24 cells were not affected by K120A GPD1 overexpression (Fig. 3Q–T). An in vivo study showed that loss of GPD1 catalytic activity did not inhibit bladder cancer growth (Fig. 3U, V). These findings revealed that GPD1 overexpression led to overgeneration of total cellular G3P and NAD+, which inhibited bladder cancer growth.
GPD1 upregulates TRPV2 expression to promote Ca2+ influx, leading to apoptosis of bladder tumor cells
To reveal the mechanism underlying GPD1-induced apoptosis of bladder cancer cells, we performed RNA-seq analysis of control 5637 cells, GPD1-OE 5637 cells and K120A GPD1-OE 5637 cells to identify the genes that are affected by GPD1. Over 2144 genes were differentially expressed (fold change > 1.5) between control and GPD1-OE cells, while over 1646 genes were differentially expressed (fold change > 1.5) between K120A GPD1-OE and GPD1-OE cells (Fig. 4A, B and Additional file 3: Table S3). A total of 1270 of the differentially expressed genes (DEGs) were shared by CON versus GPD1-OE and K120A GPD1-OE versus GPD1-OE. Based on KEGG functional annotation and gene set enrichment analysis (GSEA), the 1270 differentially expressed genes were significantly enriched in the NOD-like receptor signaling pathway (Fig. 4C and Additional file 4: Table S4). In this context, the expression of transient receptor potential vanilloid 2 (TRPV2), one of the key genes in the NOD-like receptor signaling pathway, was significantly increased in bladder cancer cells after GPD1 overexpression (Fig. 4D, E). Flow cytometry showed that GPD1 overexpression (Fig. 4F, G) or G3P/NAD+ (Fig. 4H, I) upregulated TRPV2 expression on the membranes of 5637 and T24 cells. TRPV2-mediated calcium influx can inhibit cell proliferation and induce apoptosis of bladder cancer cells . We further evaluated the intracellular calcium concentration with a Fluo-3AM probe. The results showed that GPD1 overexpression or G3P/NAD+ significantly promoted an increase in intracellular calcium concentration (Fig. 4J–M). To further determine whether GPD1 promotes tumor cell apoptosis via TRPV2-mediated Ca2+ influx, we silenced the expression of TRPV2 with siRNA (Fig. 4N). The knockdown of TRPV2 significantly inhibited the elevation of Ca2+ concentration (Fig. 4O, P) and reduced apoptosis induced by G3P/NAD+ in bladder cancer cells (Fig. 4Q, R). Similar results were observed by using the TRPV2 inhibitor tranilast (Additional file 5: Fig. S2). These data suggest that GPD1 promotes apoptosis in tumor cells through TRPV2-mediated Ca2+ influx.
GPD1 promotes TRPV2 upregulation via lysoPC-PAFR axis
To understand how GPD1, as a metabolic enzyme, upregulates the expression of TRPV2, we performed untargeted metabolomics of control 5637 cells, GPD1-OE 5637 cells and K120A GPD1-OE 5637 cells. A total of 261 differential metabolites relevant to GPD1 activity were identified in POS mode, while 38 differential metabolites relevant to GPD1 activity were identified in NEG mode (Fig. 5A). KEGG was used to explore relevant metabolic pathways, among which glycerophospholipid metabolism was the top metabolic pathway regulated in GPD1-OE 5637 cells in both POS and NEG modes. In the glycerophospholipid metabolism pathway, we found that GPD1 overexpression increased the level of lysophosphatidylcholine (lysoPC), which can be converted from G3P (Fig. 5A). We further observed that lysoPC significantly promoted the apoptosis of 5637 and T24 cells (Fig. 5B, C). In addition, similar to GPD1 overexpression, lysoPC upregulated the expression of TRPV2 on the membrane and intracellular Ca2+ levels in 5637 and T24 cells (Fig. 5D–G).
To identify receptors for the response to lysoPC, we performed qPCR for all known lysoPC receptors and found the expression of PAFR shared by 5637 and T24 cells (Fig. 5H, I). Immunohistochemical staining indicated that PAFR protein expression in the bladder tumor tissue was elevated compared with the adjacent normal tissues (Additional file 5: Fig. S3A, B). Additionally, GPD1 overexpression had no effect on PAFR expression in 5637 and T24 cells (Additional file 5: Fig. S3C, D). We used ginkgolide B, an inhibitor of PAFR, to confirm that PAFR triggers the response of bladder cancer cells to lysoPC. The results showed that ginkgolide B significantly reduced lysoPC-induced apoptosis (Fig. 5J, K) and TRPV2 expression upregulation (Fig. 5L, M). To further determine that the lysoPC-PAFR axis is necessary for GPD1-induced apoptosis, we evaluated the effect of a PAFR antagonist on apoptosis and TRPV2 expression in 5637 and T24 cells in the presence of G3P/NAD+. The results showed that ginkgolide B abrogated G3P/NAD+-induced apoptosis (Fig. 5N) and TRPV2 expression upregulation (Fig. 5O). In summary, these data demonstrate that GPD1 regulates TRPV2 expression through the lysoPC-PAFR axis.
Discovery of a cellularly active GPD1 allosteric activator
Given the antitumor function of GPD1 in bladder cancer, there is considerable interest in discovering activators of the enzyme. First, we hypothesized that functionally critical allosteric sites for GPD1 activation might exist. For identifying these sites, we predicted allosteric sites in GPD1 by using Allosite (http://mdl.shsmu.edu.cn/AST) (Fig. 6A). The predictions indicated that a pocket on the surface of GPD1 may function as an allosteric site, and key residues of the pocket are G10/G12/W39/L50/N60/V92/K100. Next, we virtually docked more than 65,000 compounds into the predicted site (Fig. 6B). Based on the top-ranked GPD1-compound binding models, we selected and purchased 10 compounds (Fig. 6C and Additional file 5: Fig. S4). To evaluate the activity of the compounds on the catalytic activity of GPD1, enzyme activity assays were performed spectrophotometrically at 340 nm. Among the 10 compounds, wedelolactone (WE) from Bioactive Compound Library Plus in MedChemExpress LLC had substantial activity in activating the catalytic activity of GPD1 (Fig. 6D, E). We then determined the binding affinity (Kd) between WE and the GPD1 protein by using microscale thermophoresis (MST). The results showed that the binding affinity for WE had a dissociation constant (Kd) of 505 nM (Fig. 6F). To evaluate the activation effect of WE on the catalytic activity of GPD1 in bladder cell lines, we monitored G3P and NAD+ levels in the 5637 and T24 cell lines after WE treatment. We found that WE increased both G3P and NAD+ at a concentration of 10 µM (Fig. 6G, H). Therefore, WE activates endogenous GPD1 catalytic activity in bladder cancer cells.
Wedelolactone inhibits bladder cancer growth in vitro and in vivo
To further investigate the function of WE in bladder cancer, we assessed the growth and death of bladder cancer cells after WE treatment. Our results showed that WE reduced the proliferation ability of 5637 and T24 cells by CCK-8 cell assay, while other compounds found in virtual screening did not affect proliferation (Fig. 7A, B). The colony formation assay indicated that WE significantly decreased the colony number of 5637 cells and T24 cells (Fig. 7C, D). Additionally, WE induced apoptosis in 5637 cells and T24 cells (Fig. 7E, F). Transwell migration and tumor sphere formation assays indicated that the migration and sphere-forming capacities of bladder cancer cell lines were also suppressed by WE (Additional file 5: Fig. S5). Similar to the effect of GPD1 overexpression, WE upregulated TRPV2 expression (Fig. 7G, H) and intracellular Ca2+ concentration (Fig. 7I, J) in 5637 cells and T24 cells. WE is also known as an inhibitor of IKK that is critical for activation of NF-κB. We selected IKK-16, another IKK inhibitor, to confirm that WE promotes apoptosis through GPD1 activation rather than NF-κB inhibition. The results showed that IKK-16 did not promote apoptosis in 5637 and T24 cells and had no effect on G3P and NAD+ levels (Additional file 5: Fig. S6). Furthermore, to determine whether WE might suppress bladder cells in vivo, we implanted 5637 cells or T24 cells into BALB/c Nude mice. Compared with the control group, WE treatment significantly reduced tumor weight (Fig. 7K–N). H&E staining and TUNEL assays showed more tumor necrosis (Fig. 7O, P) and TUNEL-positive cells (Fig. 7Q, R) in tumor tissue from the WE treatment group. These results demonstrate that WE efficiently inhibits the growth of bladder cancer in vivo.
In this study, our data show that GPD1 expression is at a low level in bladder cancer. We demonstrate for the first time that GPD1 overexpression can inhibit bladder cancer growth. We found that GPD1 promotes cellular Ca2+ influx, leading to apoptosis through the lysoPC-PAFR-TRPV2 signaling axis in bladder cancer cells. Finally, wedelolactone, an allosteric GPD1 activator, was identified for the first time based on virtual screening and was shown to effectively inhibit bladder cancer growth in vitro and in vivo (Additional file 5: Fig. S7).
Over the past 30 years, little progress has been made in the clinical management of bladder cancer. With large-scale gene sequencing efforts, a deeper understanding of bladder cancer biology has been achieved, which has led to the emergence of more favorable targeted therapies in the clinical context. In this study, we compared protein profiles in bladder cancer tissues and adjacent normal tissues by proteomics, from which we found that GPD1 protein levels were significantly downregulated in bladder cancer tissues. Little is known about the role of GPD1 in cancer. GPD1 expression has been found to be low in prostate, lung, and breast cancers [17, 18], and GPD1 was shown to have antitumor effects, which is consistent with our findings and predicts that GPD1 may be a potential therapeutic target in bladder cancer. However, specific expression of GPD1 was observed in brain tumor stem cells and was involved in tumor initiation . Our RNA-seq and metabolomics results suggest that GPD1 can regulate many cellular signaling pathways and metabolic pathways, but these pathways may differ between cancer types. Therefore, the applicability of GPD1 as a target in different tumor types needs to be further investigated.
GPD1 mRNA levels in tumor tissues were not significantly different from those in normal tissues in bladder cancer cohort from TCGA datasets (Fig. 1E). RT-qPCR analysis showed that GPD1 mRNA levels in bladder cancer cells were not lower than in SV-HUC-1 normal cells (Fig. 1L). These data suggest that downregulation of GPD1 protein levels in bladder cancer cells is independent of mRNA transcriptional level alteration. Post-transcriptional regulation is a complex process that encompasses multiple mechanisms. It has been reported that the RNA-binding protein Pub1p reduces GPD1 protein levels and enzymatic activity, but does not alter mRNA levels . We also observed upregulation of several RNA-binding proteins expression in bladder cancer tissues in our proteomics data (Additional file 2: Table S2). This requires further studies to investigate whether these proteins are involved in the downregulation of GPD1 in bladder cancer.
There are few works on the mechanism underlying the role of GPD1 in tumors. Upregulation of G3P levels, an intracellular GPD1 catalytic product, has been found to inhibit tumor cell proliferation in vitro . We have confirmed that the antitumor effect of GPD1 is dependent on the production of G3P by using amino acid mutations to lose GPD1 catalytic activity. Moreover, NAD+, another GPD1 catalytic product, was able to enhance the pro-apoptotic effect of G3P. However, the role of G3P in cancer has not been fully investigated. Here, we found that GPD1 activates the NOD-like receptor signaling pathway based on RNA-seq data. Furthermore, we found that an increase in intracellular G3P concentration significantly upregulated TRPV2 expression. Transient receptor potential vanilloid 2 (TRPV2), which belongs to the thermo TRPV subfamily (TRPV1-TRPV4) of transient receptor potential (TRP) channels [28, 29], is a calcium-permeable cation channel. Cannabidiol, an activator of TRPV2, triggers Ca2+ influx and apoptosis in the T24 cell line , which is consistent with our observations.
We further performed metabolomics and found that overexpression of GPD1 promotes increased levels of lysoPC. A prospective metabolomics study showed that high levels of lysoPC were associated with a reduced incidence of common cancers . In a study of 59 patients with various tumors, reduced lysoPC levels were found to be associated with an increase in inflammatory process parameters (CRP, albumin reduction) and severe weight loss . In our study, lysoPC promoted Ca2+ influx through TRPV2 and apoptosis in bladder cancer cells, which was similar to the effect of GPD1 overexpression. A similar study showed that lysoPC induces Ca2+ influx and apoptosis through activation of TRPC1/TRPC3 channels in human coronary artery smooth muscle cells . The signaling cell cascade activated by lysoPC is triggered by receptors. The receptors that have been reported are G2A, GPR17, GPR4, and PAFR. We found that lysoPC promotes TRPV2 expression and Ca2+ influx via PAFR in 5637 and T24 cells. Overall, this study provides new information that GPD1 induces Ca2+ influx and apoptosis through the lysoPC-PAFR-TRPV2 axis.
Based on the potential antitumor function of GPD1, the development of drugs that modulate GPD1 activity may provide a new therapeutic strategy for bladder cancer. Allosteric regulation is a direct, rapid, and widely effective way to regulate protein function, and its dysfunction or abnormal regulation leads to the development of many major human diseases. After more than 50 years of development, the number of allosteric proteins and allosteric regulatory molecules has been increasing, and research on allosteric mechanisms has deepened and gradually translated into drug development and other applications . With the natural advantages of higher target selectivity, fewer side effects, and lower toxicity, the pharmaceutical industry and academia are investing increasing efforts in the development of allosteric drugs [34, 35]. We used the computational chemical biology tool developed in the AlloSteric Database (ASD) to analyze the GPD1 structure and identify potential allosteric regulatory sites. Based on this site, we combined computer-assisted virtual screening and GPD1 catalytic activity assays in vitro to screen and obtain the small-molecule agonist wedelolactone, which can modulate GPD1 catalytic activity, and preliminarily demonstrated its antitumor functions in vitro and in vivo. Wedelolactone, a medicinal plant-derived natural compound, suppresses LPS-induced caspase-11 expression by directly inhibiting the IKK complex . Some studies have also reported that wedelolactone induces caspase-dependent apoptosis in prostate cancer cells via downregulation of PKCε without inhibiting Akt . These studies suggest that wedelolactone may be a multitargeted small-molecule compound beyond targeting GPD1. To improve the target selectivity of wedelolactone, medicinal chemistry optimization is necessary. The combination of treatments may be relatively more effective in eliminating cancer cells with minimal side effects. Metformin, a GPD2 inhibitor, leads to an increase in the concentration of G3P in cancer cells. It has been reported that GPD1 overexpression enhances the anticancer effect of metformin through synergistically increasing G3P levels . This information gives us an insight that the allosteric activators of GPD1 we found in combination with metformin may be a better therapeutic strategy for bladder cancer.
Taken together, our findings highlight GPD1 as a novel tumor suppressor in bladder cancer. Our data showed that lysoPC converted from the GPD1 catalytic product G3P induced TRPV2 expression via PAFR to promote Ca2+ influx, which led to apoptosis of tumor cells. Finally, we screened and obtained the allosteric agonist wedelolactone for GPD1, demonstrating that pharmacological activation of GPD1 is a potential therapeutic approach for bladder cancer.
Availability of data and materials
The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request.
Glycerol 3-phosphate dehydrogenase 1
Muscle-invasive bladder cancer
Nonmuscle-invasive bladder cancer
Platelet-activating factor receptor
Transient receptor potential vanilloid 2
Terminal deoxynucleotidyl transferase-mediated UTP nick end labeling
Ultra-high-performance liquid chromatography quadrupole time of flight mass spectrometry
The cancer genome atlas
Cell counting kit-8
Siegel RL, Miller KD, Fuchs HE, Jemal A. Cancer statistics, 2021. CA Cancer J Clin. 2021;71(1):7–33. https://doi.org/10.3322/caac.21654.
Cattrini C, Boccardo F. Atezolizumab and bladder cancer: facing a complex disease. Lancet. 2018;391(10118):305–6. https://doi.org/10.1016/S0140-6736(18)30095-3.
Hurst CD, Knowles MA. Bladder cancer: multi-omic profiling refines the molecular view. Nat Rev Clin Oncol. 2018;15(4):203–4. https://doi.org/10.1038/nrclinonc.2017.195.
Sung H, Ferlay J, Siegel RL, Laversanne M, Soerjomataram I, Jemal A, et al. Global cancer statistics 2020: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J Clin. 2021;71(3):209–49. https://doi.org/10.3322/caac.21660.
Chen W, Zheng R, Baade PD, Zhang S, Zeng H, Bray F, et al. Cancer statistics in China, 2015. CA Cancer J Clin. 2016;66(2):115–32. https://doi.org/10.3322/caac.21338.
Gakis G. Management of muscle-invasive bladder cancer in the 2020s: challenges and perspectives. Eur Urol Focus. 2020;6(4):632–8. https://doi.org/10.1016/j.euf.2020.01.007.
Wong YNS, Joshi K, Pule M, Peggs KS, Swanton C, Quezada SA, et al. Evolving adoptive cellular therapies in urological malignancies. Lancet Oncol. 2017;18(6):e341–53. https://doi.org/10.1016/S1470-2045(17)30327-3.
Funt SA, Rosenberg JE. Systemic, perioperative management of muscle-invasive bladder cancer and future horizons. Nat Rev Clin Oncol. 2017;14(4):221–34. https://doi.org/10.1038/nrclinonc.2016.188.
Meeks JJ, Bellmunt J, Bochner BH, Clarke NW, Daneshmand S, Galsky MD, et al. A systematic review of neoadjuvant and adjuvant chemotherapy for muscle-invasive bladder cancer. Eur Urol. 2012;62(3):523–33. https://doi.org/10.1016/j.eururo.2012.05.048.
Menaya J, Gonzalez-Manchon C, Parrilla R, Ayuso MS. Molecular cloning, sequencing and expression of a cDNA encoding a human liver NAD-dependent alpha-glycerol-3-phosphate dehydrogenase. Biochim Biophys Acta. 1995;1262(1):91–4. https://doi.org/10.1016/0167-4781(95)00069-s.
Reshef L, Olswang Y, Cassuto H, Blum B, Croniger CM, Kalhan SC, et al. Glyceroneogenesis and the triglyceride/fatty acid cycle. J Biol Chem. 2003;278(33):30413–6. https://doi.org/10.1074/jbc.R300017200.
Joshi M, Eagan J, Desai NK, Newton SA, Towne MC, Marinakis NS, et al. A compound heterozygous mutation in GPD1 causes hepatomegaly, steatohepatitis, and hypertriglyceridemia. Eur J Hum Genet. 2014;22(10):1229–32. https://doi.org/10.1038/ejhg.2014.8.
Basel-Vanagaite L, Zevit N, Har Zahav A, Guo L, Parathath S, Pasmanik-Chor M, et al. Transient infantile hypertriglyceridemia, fatty liver, and hepatic fibrosis caused by mutated GPD1, encoding glycerol-3-phosphate dehydrogenase 1. Am J Hum Genet. 2012;90(1):49–60. https://doi.org/10.1016/j.ajhg.2011.11.028.
Ratner PL, Fisher M, Burkart D, Cook JR, Kozak LP. The role of mRNA levels and cellular localization in controlling sn-glycerol-3-phosphate dehydrogenase expression in tissues of the mouse. J Biol Chem. 1981;256(7):3576–9.
Sato T, Morita A, Mori N, Miura S. Glycerol 3-phosphate dehydrogenase 1 deficiency enhances exercise capacity due to increased lipid oxidation during strenuous exercise. Biochem Biophys Res Commun. 2015;457(4):653–8. https://doi.org/10.1016/j.bbrc.2015.01.043.
Brown LJ, Koza RA, Marshall L, Kozak LP, MacDonald MJ. Lethal hypoglycemic ketosis and glyceroluria in mice lacking both the mitochondrial and the cytosolic glycerol phosphate dehydrogenases. J Biol Chem. 2002;277(36):32899–904. https://doi.org/10.1074/jbc.M202409200.
Descotes F, Jezequel P, Spyratos F, Campion L, Grenot C, Lerebours F, et al. Identification of potential prognostic biomarkers for node-negative breast tumours by proteomic analysis: a multicentric 2004 national PHRC study. Int J Oncol. 2012;41(1):92–104. https://doi.org/10.3892/ijo.2012.1456.
Zhou C, Yu J, Wang M, Yang J, Xiong H, Huang H, et al. Identification of glycerol-3-phosphate dehydrogenase 1 as a tumour suppressor in human breast cancer. Oncotarget. 2017;8(60):101309–24. https://doi.org/10.18632/oncotarget.21087.
Xie J, Ye J, Cai Z, Luo Y, Zhu X, Deng Y, et al. GPD1 enhances the anticancer effects of metformin by synergistically increasing total cellular glycerol-3-phosphate. Cancer Res. 2020;80(11):2150–62. https://doi.org/10.1158/0008-5472.CAN-19-2852.
Liu R, Feng Y, Deng Y, Zou Z, Ye J, Cai Z, et al. A HIF1alpha-GPD1 feedforward loop inhibits the progression of renal clear cell carcinoma via mitochondrial function and lipid metabolism. J Exp Clin Cancer Res. 2021;40(1):188. https://doi.org/10.1186/s13046-021-01996-6.
Boxer GE, Shonk CE. Low levels of soluble DPN-linked alpha-glycerophosphate dehydrogenase in tumors. Cancer Res. 1960;20:85–91.
Shonk CE, Arison RN, Koven BJ, Majima H, Boxer GE. Enzyme patterns in human tissues. 3. Glycolytic enzymes in normal and malignant tissues of the colon and rectum. Cancer Res. 1965;25:206–13.
Rusu P, Shao C, Neuerburg A, Acikgoz AA, Wu Y, Zou P, et al. GPD1 specifically marks dormant glioma stem cells with a distinct metabolic profile. Cell Stem Cell. 2019;25(2):241–57. https://doi.org/10.1016/j.stem.2019.06.004.
Lu Q, Liu T, Feng H, Yang R, Zhao X, Chen W, et al. Circular RNA circSLC8A1 acts as a sponge of miR-130b/miR-494 in suppressing bladder cancer progression via regulating PTEN. Mol Cancer. 2019;18(1):111. https://doi.org/10.1186/s12943-019-1040-0.
He R, Cristobal JR, Gong NJ, Richard JP. Hydride transfer catalyzed by glycerol phosphate dehydrogenase: recruitment of an acidic amino acid side chain to rescue a damaged enzyme. Biochemistry. 2020;59(51):4856–63. https://doi.org/10.1021/acs.biochem.0c00801.
Yamada T, Ueda T, Shibata Y, Ikegami Y, Saito M, Ishida Y, et al. TRPV2 activation induces apoptotic cell death in human T24 bladder cancer cells: a potential therapeutic target for bladder cancer. Urology. 2010;76(2):509. https://doi.org/10.1016/j.urology.2010.03.029.
Orozco H, Sepulveda A, Picazo C, Matallana E, Aranda A. RNA binding protein Pub1p regulates glycerol production and stress tolerance by controlling Gpd1p activity during winemaking. Appl Microbiol Biotechnol. 2016;100(11):5017–27. https://doi.org/10.1007/s00253-016-7340-z.
Park U, Vastani N, Guan Y, Raja SN, Koltzenburg M, Caterina MJ. TRP vanilloid 2 knock-out mice are susceptible to perinatal lethality but display normal thermal and mechanical nociception. J Neurosci. 2011;31(32):11425–36. https://doi.org/10.1523/JNEUROSCI.1384-09.2011.
Peralvarez-Marin A, Donate-Macian P, Gaudet R. What do we know about the transient receptor potential vanilloid 2 (TRPV2) ion channel? FEBS J. 2013;280(21):5471–87. https://doi.org/10.1111/febs.12302.
Kuhn T, Floegel A, Sookthai D, Johnson T, Rolle-Kampczyk U, Otto W, et al. Higher plasma levels of lysophosphatidylcholine 18:0 are related to a lower risk of common cancers in a prospective metabolomics study. BMC Med. 2016;14:13. https://doi.org/10.1186/s12916-016-0552-3.
Taylor LA, Arends J, Hodina AK, Unger C, Massing U. Plasma lyso-phosphatidylcholine concentration is decreased in cancer patients with weight loss and activated inflammatory status. Lipids Health Dis. 2007;6:17. https://doi.org/10.1186/1476-511X-6-17.
Wang Y, Wang Y, Li GR. TRPC1/TRPC3 channels mediate lysophosphatidylcholine-induced apoptosis in cultured human coronary artery smooth muscles cells. Oncotarget. 2016;7(32):50937–51. https://doi.org/10.18632/oncotarget.10853.
Nussinov R, Tsai CJ. Allostery in disease and in drug discovery. Cell. 2013;153(2):293–305. https://doi.org/10.1016/j.cell.2013.03.034.
Digby GJ, Conn PJ, Lindsley CW. Orthosteric- and allosteric-induced ligand-directed trafficking at GPCRs. Curr Opin Drug Discov Devel. 2010;13(5):587–94.
Christopoulos A, May LT, Avlani VA, Sexton PM. G-protein-coupled receptor allosterism: the promise and the problem(s). Biochem Soc Trans. 2004;32(Pt 5):873–7. https://doi.org/10.1042/BST0320873.
Kobori M, Yang Z, Gong D, Heissmeyer V, Zhu H, Jung YK, et al. Wedelolactone suppresses LPS-induced caspase-11 expression by directly inhibiting the IKK complex. Cell Death Differ. 2004;11(1):123–30. https://doi.org/10.1038/sj.cdd.4401325.
Sarveswaran S, Gautam SC, Ghosh J. Wedelolactone, a medicinal plant-derived coumestan, induces caspase-dependent apoptosis in prostate cancer cells via downregulation of PKCepsilon without inhibiting Akt. Int J Oncol. 2012;41(6):2191–9. https://doi.org/10.3892/ijo.2012.1664.
This study was funded by the National Natural Science Foundation of China (82072822, 81802535, 81802532, 81772710, 81972388), China postdoctoral fund (223427), and Nanjing Medical Science and technique Development Foundation (YKK 18064). This study was also supported by the Project of Invigorating Health Care through Science, Technology and Education Jiangsu Provincial Key Medical Discipline (ZDXKB2016014).
Ethics approval and consent to participate
The study was approved by the Ethical Committee of Nanjing Drum Tower Hospital, Medical School of Nanjing University, and all included patients provided written informed consent. All methods were carried out in accordance with relevant guidelines and regulations. All animal experiments were carried out in accordance with the Guide for the Care and Use of Laboratory Animals (The Ministry of Science and Technology of China) and the related ethical regulations of Nanjing Drum Tower Hospital, Medical school of Nanjing University.
Consent for publication
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Summary of the characteristics of 17 patients.
The proteomic data per sample.
The list of DEG identified in RNA-seq.
The pathways base on gene set enrichment analysis.
The effect of GPD1 overexpression on tumor cell phenotype. Figure S2. Flow cytometry analysis of apoptosis in 5637 cells and T24 cells treated with G3P/NAD+ in presence of tranilast or not. Figure S3. GPD1 did not affect the expression of PAFR in 5637 and T24 cells. Figure S4. The structure and docking score of compounds on the basis of the top-ranked GPD1–compound binding models. Figure S5. The effect of Wedelolactone on tumor cell phenotype. Figure S6. The effect of IKK-16 on apoptosis and GPD1 activation in bladder cancer cells. Figure S7. Schematic summary.
About this article
Cite this article
Zhang, W., He, X., Yin, H. et al. Allosteric activation of the metabolic enzyme GPD1 inhibits bladder cancer growth via the lysoPC-PAFR-TRPV2 axis. J Hematol Oncol 15, 93 (2022). https://doi.org/10.1186/s13045-022-01312-5
- Bladder cancer
- Metabolic enzyme
- Tumor suppressor
- Allosteric activator