Protein Palmitoylation Regulates Cell Survival by Modulating XBP1 Activity in Glioblastoma Multiforme

Glioblastoma multiforme (GBM) almost invariably acquires an invasive phenotype, resulting in limited therapeutic options. Protein palmitoylation markedly affects tumorigenesis and malignant progression in GBM. The role of protein palmitoylation in GBM, however, has not been systematically reported. This study aimed to investigate the effect of protein palmitoylation on GBM cell survival and the cell cycle. In this study, most palmitoyltransferases were upregulated in GBM and its cell lines, and protein palmitoylation participated in signaling pathways controlling cell survival and the GBM cell cycle. Inhibition of protein palmitoylation with substrate-analog inhibitors, that is, 2-bromopalmitate, cerulenin, and tunicamycin, induced G2 cell cycle arrest and cell death in GBM cells through enhanced endoplasmic reticulum (ER) stress. These effects are primarily attributed to the palmitoylation inhibitors activating pro-apoptotic pathways and ER stress signals. Further analysis revealed was the accumulation of SUMOylated XBP1 (X-box binding protein 1) and its transcriptional repression, along with a reduction in XBP1 palmitoylation. Taken together, the present results indicate that protein palmitoylation plays an important role in the survival of GBM cells, further providing a potential therapeutic strategy for GBM.


INTRODUCTION
Glioblastoma multiforme (GBM) is considered the most common type of intracranial tumor. 1,2 Despite improvements in conventional surgical, radiotherapy, and chemotherapy, the median survival of GBM patients is only approximately 14 months. 3,4 Thus, new therapeutic strategies or drugs targeting GBM are urgently needed.
Disturbances in endoplasmic reticulum (ER) homeostasis activate the ER stress response, which is mediated by three parallel signaling pathways initiated by protein kinase R-like ER kinase (PERK), inositolrequiring enzyme 1a (IRE1a), and activating transcription factor-6 (ATF6). 5,6 X-box binding protein 1 (XBP1) is a key transcription factor in the IRE1a signaling branch. 7,8 Upon ER stress, the XBP1 transcript is spliced by the endoribonuclease, IRE1a, to generate the active form, spliced XBP1 (XBP1s), which regulates genes implicated in protein folding, trafficking, and secretion, thus enhancing cell survival and ER homeostasis under ER stress. In contrast, unspliced XBP1 (XBP1u) functions as the dominant negative form antagonizing the function of XBP1s. 9,10 Furthermore, IRE1a serves as a kinase to activate pro-apoptotic c-Jun N-terminal kinase (JNK) signaling to induce cell death. 11,12 The seemingly paradoxical dual functions of IREa in inducing both pro-survival and pro-apoptotic signaling upon ER stress is challenging to understand. Protein functions are influenced by their expression level, localization, interaction with other proteins, and posttranslational modifications. [13][14][15] Numerous proteins can be modified through palmitoylation at cysteine residues by a family of enzymes containing a unique DHHC domain. 16,17 At least 23 members of this palmitoyltransferase (PAT) family have been reported in the mammalian genome. Protein localization, trafficking, and stability are regulated by this modification. 18,19 Palmitoylation markedly influences tumorigenesis and tumor progression through different substrates, particularly in glioma development and malignant progression. 20,21 For example, ZDHHC5 is upregulated in p53 mutant glioma cells and enhances their invasiveness and tumorigenicity. 22 Moreover, ZDHHC23 and ZDHHC18 contribute to the transition of glioma stem cells and cell survival in a stressful tumor microenvironment. 23 Thus, palmitoylation inhibitors would potentially be effective for treating GBM.
2-Bromopalmitate (2BP), which irreversibly inhibits PATs, is the most commonly used palmitoylation inhibitor. 24 Two additional important lipid-based DHHC inhibitors are cerulenin (Cer) and tunicamycin (Tun). 25 In the present study, we initially analyzed the expression pattern of DHHC members in GBM and investigated the effect of protein palmitoylation on GBM cell survival and the cell cycle.

Protein Palmitoylation Occurs in GBM and Cell Lines
To determine the role of protein palmitoylation in GBM pathogenesis and malignant progression, we initially evaluated the expression of the 23 potential PATs in GBM based on the Human Protein Atlas ( Figure 1A). Most of the DHHC members were upregulated in GBM, especially for ZDHHC21, ZDHHC12, ZDHHC17, and ZDHHC18. We then quantified total protein palmitoylation via acyl-biotin exchange analyses. Compared to normal brain tissues, GBM displayed several protein bands with covalent palmitoyl modifications ( Figure 1B). These results suggest that GBM involves changes in protein palmitoylation, and aberrant protein palmitoylation may be a potent biomarker for GBM. Furthermore, we analyzed the expression of all DHHCs in six human GBM cell lines via RT-PCR (see Table S1 for primer sequences). Of the 23 PATs, 20 were detected in these six cell lines ( Figure 1C). ZDHHC7, ZDHHC11, and ZDHHC22 were undetectable, suggesting that they were downregulated in GBM cells. In subsequent experiments, we used SF126 cells to investigate the roles of protein palmitoylation during GBM pathogenesis because the 20 PATs were most upregulated in this cell line. Consistent with the GBM tissues, some palmitoylated proteins were detected in SF126 cells ( Figure 1D). that is, 2BP, Cer, and Tun. Consequently, 50 mM 2BP, 25 mM Cer, and 2.5 mM Tun adequately repressed protein palmitoylation in glioma cells by 85%, 75%, and 65%, respectively ( Figure 1D). EZH2 can be palmitoylated; 22 however, in this study, EZH2 palmitoylation levels were reduced after 2BP, Cer, or Tun treatment to different degrees ( Figure S1A).
The effect of palmitoylation on cell survival was determined by using 2BP at 10-100 mM, Cer at 5-75 mM, and Tun at 1.0-7.5 mM in SF126 cells. Consequently, treatment with 50 mM 2BP, 25 mM Cer, and 2.5 mM Tun adequately decreased the cell survival rate by 50% (Figure 2A) and adequately decreased cell viability ( Figure 2B), colony formation ( Figure 2C), and cell proliferation ( Figure 2D) among GBM cells. Thus, we used palmitoylation inhibitors at these concentrations in subsequent experiments.
To investigate the effects of protein palmitoylation on cell cycle progression, GBM cells were exposed to 2BP, Cer, and Tun for 48 h, and the cell cycle phase distribution was analyzed using fluorescence-acti-vated cell sorting ( Figures 3A-3C). Quantification of untreated SF126 cells revealed that 76.5% (76.5% for negative control; 75.1% for dimethyl sulfoxide [DMSO]) of the cells were in the G 0 /G 1 phase, 7.27% (7.27% for negative control; 8.17% for DMSO) were in the S phase, and 14.8% (14.8% for negative control; 17.2% for DMSO) were in the G 2 /M phase of the cell cycle 48 h after plating. Treatment of SF126 cells with 50 mM 2BP for 48 h increased the percentage of cells in the G 2 /M phase to 62.8% and reduced cells in the G 0 /G 1 and S phases to 11.8% and 4.96%, respectively. The sub-G 1 cell population exceeded 20.6%. The percentage of SF126 cells treated with 25 mM Cer for 48 h in the G 2 /M phase increased to 52.0% and the sub-G 1 population of cells exceeded 16.1%. Tun (2.5 mM) increased the percentage of cells in the G 2 /M phase to 25.0% and the sub-G 1 population of cells exceeded 29.0%. These data indicate that SF126 GBM cells underwent cell cycle arrest in the G 2 /M phase upon treatment with 2BP, Cer, or Tun.
Furthermore, we performed annexin V-propidium iodide (PI) double staining with fluorescence-activated cell sorting to verify that SF126 GBM cells treated with palmitoylation inhibitors underwent apoptosis ( Figure 3D). Control GBM cells exhibited a very low apoptotic rate (1.36% for negative control; 1.67% for DMSO). The apoptotic population increased significantly to 17.6% with 50 mM 2BP, 12.3% with 25 mM Cer, and 24.7% with 2.5 mM Tun ( Figure 3E). These data indicate that apoptosis plays a role in palmitoylation inhibitor-induced cell death.

Transcriptomes Are Altered after Palmitoylation Inhibition
We performed RNA sequencing (RNA-seq) to investigate alterations in the transcriptome upon 2BP, Cer, and Tun treatment. On treating the SF126 GBM cell line with 50 mM 2BP, 433 genes were upregulated and 441 were downregulated in comparison with the control DMSO-treated cells. These upregulated genes were associated with biological processes, including responses to organic substances and carboxylic acid metabolic processes (Figure 4A). Upon treatment with 25 mM Cer, 488 genes were downregulated and 271 were upregulated in comparison with the DMSO control group. These upregulated genes were involved in responses to organic substances and ER stress ( Figure 4B). Upon treatment with Tun (2.5 mM versus DMSO), 590 genes were upregulated and 833 genes were downregulated. These upregulated genes were involved in responses to ER stress and unfolded proteins ( Figure 4C). In total, 148 transcripts were differentially expressed in the 2BP versus DMSO, Cer versus DMSO, and Tun versus DMSO datasets. We performed a Gene Ontology analysis with respect to the biological process to determine the biological relevance of the 148 differentially expressed gene transcripts obtained from the aforementioned three RNA-seq datasets. Several biological processes were enriched in palmitoylation inhibitor-treated samples as opposed to untreated samples, primarily including responses to stress, ER stress, and chemical stimuli ( Figure 4D). Taken together, these results indicate that inhibition of palmitoylation in GBM cells activates responses to ER stress, increases the levels of unfolded proteins, and causes other changes associated with ER stress.

Inhibition of Palmitoylation Triggers the Pro-apoptotic ER Stress Response
Most PATs were localized at the ER, Golgi, and in cytoplasmic vesicles ( Figure 5A). PDI (protein disulfide isomerase) is an ER enzyme that catalyzes thiol-disulfide exchange, thus facilitating disulfide bond formation and rearrangement reactions. 26 Immunostaining revealed that PDI accumulated in cells treated with palmitoylation inhibitors, suggesting the occurrence of ER stress ( Figure 5B). Consistent with these results, electron microscopy revealed progressive morphological changes in the ER after inhibition of palmitoylation, including swollen and abnormal areas ( Figure 5C), directly indicating ER stress.
Because ER stress is one of the major cell stress responses associated with apoptosis, we investigated whether palmitoylation inhibitors activated the ER stress response, as well as its proapoptotic effects in GBM cells. As shown in Figure 6, treatment of SF126 cells with 2BP, Cer, or Tun increased IRE1a phosphorylation and, concurrently, IRE1a levels. Since IRE1a activation results in JNK activation, 11 herein JNK phosphorylation showed a similar time course pattern as IRE1a. Consistent with IRE1a phosphorylation, XBP1s protein levels significantly increased. On treatment with the anti-XBP1s antibody, XBP1u protein levels decreased.
The ER stress response of SF126 GBM cells exposed to 2BP, Cer, or Tun was also evidenced by the increase in CCAAT-enhancer binding protein homologous protein (CHOP) levels. Apoptosis in these cells was reflected by increased cleavage of caspase-3 ( Figure 6). These results show that inhibition of palmitoylation induces the ER stress response, which correlates with palmitoylation inhibitor-induced apoptosis.

Inhibition of Palmitoylation Represses Pro-survival XBP1 Signaling
Having observed that 2BP, Cer, and Tun activate the ER stress response, we further assessed the effect of inhibition of palmitoylation on target mRNA (XBP1, VEGFA, GADD34, and CHOP) levels in these signaling pathways in GBM cells. Among these target genes, both CHOP and GADD34 are regulated by non-IRE1a (PERK and ATF6) branches of the ER stress response. 27 Human VEGFA contains putative XBP1s that binds at sites in its promoter, and these sites are conserved across species, including mice, rats, and humans. 28 XBP1 is a downstream target of the XBP1s gene in positive feedback loops. 28 In the present study, 2BP, Cer, and Tun suppressed XBP1 and VEGFA mRNAs in SF126 cells and potently induced CHOP and GADD34 mRNAs ( Figure 7A).
Furthermore, we used a 5Â unfolded protein response element (UPRE) luciferase reporter construct to assess the transcriptional activity of XBP1s. Consequently, 2BP, Cer, and Tun treatment inhibited the relative luciferase activity of 5Â UPRE that was induced by cotransfected XBP1s ( Figure 7B). This finding indicates that palmitoylation inhibitors potentially inhibit the transcriptional activity of XBP1s selectively while enhancing the mRNA expression of downstream factors, as reflected via upregulation of GADD34 and CHOP mRNAs (target genes of the non-IRE1a cascade ER stress response, such as the PERK/eIF2a signaling pathway) in GBM cells.
XBP1s levels were increased, whereas its transcriptional activity was repressed, after treatment with palmitoylation inhibitors. SUMOylation potentially suppresses the transcriptional activity of XBP1s during ER stress. 29 To test this possibility, we investigated whether SUMOylation of XBPs is increased in palmitoylation inhibitortreated cells. As shown in Figure 7C, XBP1s SUMOylation was increased after 2BP, Cer, or Tun treatment in comparison with untreated cells, suggesting that reduction of protein palmitoylation specifically results in the accumulation of SUMOylated XBP1s. Indeed, an increase in SUMOylated XBP1s levels is potentially associated with the inhibition of palmitoylated XBP1s because palmitoylation levels of XBP1s were discernibly decreased upon treatment with 2BP, Cer, or Tun. As predicted using CSS-Palm 4.0, XBP1s has three potential palmitoylation sites at its C terminus: Cys325, Cys331, and Cys339. These potential palmitoylation sites are proximal to SUMOylation sites at XBP1s, that is, Lys281 and Lys302. Palmitoylation of XBP1s may hinder XBP1s SUMOylation. As assumed, mutations at the potential palmitoylation sites in XBP1s decreased XBP1s palmitoylation levels and increased the XBP1s SUMOylation levels (Figure 7D). Considering the interaction between a PAT and its substrate, we tested some DHHC members upregulated in GBM that may bind to XBP1s in SF126 cells. Co-immunoprecipitation analysis revealed the physical interaction between XBP1s and three DHHC members, that is, ZDHHC1, ZDHHC6, and ZDHHC17 ( Figure 7E). Concurrently, knockdown of ZDHHC1, ZDHHC6, or ZDHHC17 resulted in a reduction in XBP1s palmitoylation but an increase in XBP1s SUMOylation ( Figure 7F). These results suggest that XBP1s palmitoylation specifically inhibits XBP1s SUMOylation, and palmitoylation inhibitors can increase XBP1s SUMOylation and repress its transcriptional activity.

Treatment Prospects of Palmitoylation Inhibitors for GBM In Vivo
Recently, radiotherapy and temozolomide (TMZ) chemotherapy treatments have become the gold standard treatment methods for GBM after surgery. We initially analyzed the expression of the 23 potential PATs in SF126 cells after treatment with X-ray irradiation (IR) and TMZ via RT-PCR ( Figure S2). More interestingly, ZDHHC5 and ZDHHC8 were upregulated only in the IR-treated group ( Figure S2A), suggesting that these PATs potentially contribute to this phenomenon. ZDHHC4/9/11/14/17/21/22/24 were upregulated only in the TMZ-treated group ( Figure S2B), suggesting that these PATs potentially corroborate the effects of TMZ treatment. ZDHHC7 was upregulated in both IR-and TMZ-treated cells ( Figure S2C). ZDHHC2 was upregulated after IR treatment but downregulated after TMZ treatment ( Figure S2D). However, ZDHHC12/13/16/18/23 were slightly downregulated in the IR-or TMZ-treated cells rather than the untreated cells ( Figure S2E), while ZDHHC1/3/6/15/19/20 levels remained unchanged ( Figure S2F). These results indicate that PATs markedly participate in IR and TMZ treatments among GBM patients.
Moreover, palmitoylation inhibitors (2BP, Cer, or Tun) displayed cooperative effects with IR and TMZ and significantly repressed GBM cell proliferation ( Figure 8A), migration ( Figure 8B), and invasion ( Figure 8C). A xenograft model was used to evaluate the in vivo anti-tumor activity of palmitoylation inhibitors. Approximately 1 Â 10 6 SF126 cells were injected at the dorsa of nude mice. After tumors grew to approximately 100 mm 3 , the mice (six animals per group) were subcutaneously administered 3 mg/kg 2BP, 5 mg/kg Cer, or 2 mg/kg Tun daily for 10 days, along with IR and TMZ treatment. Consequently, all three agents significantly inhibited tumor growth in comparison with untreated mice or the TMZ/IR treatment group ( Figure 8D). Immunohistochemical assessment revealed that Ki-67 expression levels in vivo were lower in inhibitor-treated versus untreated mice or TMZ/IR-treated mice ( Figure 8E). As expected, palmitoylation levels of XBP1s were decreased to different degrees (Figure S1B). The in vivo analysis highlighted treatment prospects via IR and TMX treatment for GBM patients.

DISCUSSION
The present results reveal a novel mechanism through which palmitoylation inhibitors (2BP, Cer, and Tun) induce cell death in GBM cells. These inhibitors activate ER stress signaling and its proapoptotic effects, while promoting SUMOylated XBP1 accumulation and inhibiting XBP1s-mediated transcription and pro-survival signaling during ER stress. These results suggest that palmitoylation inhibitors are potential pharmacotherapeutic candidates for GBM by subjecting GBM cells to double jeopardy, thereby skewing the ER stress response survival/apoptosis balance toward apoptosis.
Almost all DHHC family members are associated with human diseases, spanning neurological disorders to some cancers, through different substrates such as oncogenic Ras, 30 Hedgehog, 31 Wnt, 32 and epidermal growth factor receptor 33 proteins. Hence, palmitoylation inhibitors are potentially beneficial for treating some cancers. Well-established lipid-based palmitoylation inhibitors include 2BP, Tun, and Cer. Within cells, 2BP is metabolically incorporated as 2BP-coenzyme A (CoA) and inhibits PAT activity of all DHHCs that have been tested. 34,35 The mode of action for Tun is still unclear; Figure 6. Inhibition of Palmitoylation Activates Inositol-Requiring Enzyme 1a (IRE1a) and Causes Death in SF126 Cells (A-C) SF126 cells were treated with 2BP (50 mM; A), Cer (25 mM; B), or Tun (2.5 mM; C) for the indicated periods. Whole-cell lysates were harvested and analyzed for phosphorylated and total IRE1a, phosphorylated and total c-Jun N-terminal kinase (JNK), spliced (XBP1s) and unspliced X-box binding protein 1 (XBP1)u, CCAAT-enhancer binding protein homologous protein (CHOP), pro-caspase-3, cleaved caspase-3, and the loading control, b-tubulin. ***p < 0.001, unpaired t test.
www.moleculartherapy.org Molecular Therapy: Oncolytics Vol. 17 June 2020 however, it may compete with the palmitoyl-CoA substrate for binding to DHHC-PATs. 36 Cer is proposed to react irreversibly with the catalytic cysteine of DHHC, or cysteine residues of the substrate proteins, through its epoxy carboxamide group. 37 Protein palmitoylation significantly influences tumorigenesis and progression in GBM. Thus, we used these lipid-based palmitoylation inhibitors to assess the role of palmitoylation in GBM pathogenesis. Our results show that inhibition of protein palmitoylation with 2BP, Cer, or Tun induced GBM cell G 2 cell cycle arrest, increased cell death, and inhibited glioma growth. Regarding the regulation of GBM cell survival by protein palmitoylation, we found that inhibition of palmitoylation induces ER stress through the IRE1a/JNK/XBP1 signaling axis. IRE1a/JNK/ XBP1 signals play both pro-apoptotic and pro-survival roles under stress conditions mediated by JNK and XBP1s, respectively. 38,39 XBP1s activates proteins regulating protein trafficking, folding, and quality control through the ER-associated protein degradation pathway. 40 GBM is a solid tumor in cells that can survive hypoxia, nutrient deprivation, and low pH. Furthermore, Bip and PERK are upregulated in GBM cells, 41,42 suggesting an intrinsic dependence on the ER stress pathway for survival. The ER stress sensor IRE1a contributes to GBM progression through XBP1 mRNA that splices and regulates IRE1-dependent mRNA decay. 42 Interference with the ER stress pathway may affect cell survival. For example, pharmacologic inhibition of IRE1a endoribonuclease activity inhibits XBP1 mRNA splicing and decreases cellular sensitivity to ER stress-induced death. Bortezomib increases intrinsic ER stress in GBM cells by inducing the accumulation of numerous unfolded proteins and/or and blocking the degradation of misfolded proteins. 43 Bortezomib, A 5Â unfolded pathway response element-luciferase reporter and XBP1s expression construct were used to determine the transcriptional activity of XBP1s. The firefly luciferase value was divided by the Renilla luciferase value to normalize each sample. Data are expressed as means ± SD (n = 3). (C) Increased accumulation of SUMOylated XBP1 in 2BP (50 mM), Cer (25 mM), or Tun (2.5 mM) treatment versus that in control group. XBP1 was immunoprecipitated with anti-XBP1 antibody (IP) from these cell lysates. Bound proteins were blotted with anti-XBP1 or anti-SUMO1 antibody (IB). (D) Cys325, Cys331, and Cys339 were determined as XBP1 palmitoylation sites. Palmitoylation sites Cys325, Cys331, and Cys339 of XBP1 were predicted using CSS-Palm 4.0 and mutated to Ala, respectively, the palmitoylation level of XBP1 was detected via the ABE method, and the SUMOylated XBP1 was also analyzed. (E) SF126 was transfected with hemagglutinin (HA)-tagged DHHC members upregulated in GBM, and subjected to immunoprecipitation (IP) of HA. (F) Palmitoylation and SUMOylation levels of XBP1 in SF126 cells transfected with siRNAs for different DHHCs. *p < 0.05, **p < 0.01, ***p < 0.001, unpaired t test. ns, not significant. however, represses the function of XBP1s via elevations in XBP1u levels, the unspliced and transcriptionally inactive form of XBP1 that antagonizes the functions of XBP1s. 44 Furthermore, our results indicate a distinct mechanism through which palmitoylation inhibitors affect XBP1s signaling in GBM cells. The present results show that although 2BP, Cer, and Tun activated IRE1a activity and enhanced mRNA splicing of XBP1, these agents exerted a specific inhibitory effect on XBP1s transcriptional activity, paradoxically without elevating XBP1u, leading to a defective autoregulation loop for XBP1, wherein XBP1s regulated the expression of total XBP1s. This concept is supported by the reduction in total XBP1 mRNA levels despite an increase in XBP1s protein levels in GBM cells in response to palmitoylation inhibitor treatment. Furthermore, these Data are expressed as means ± SD (n = 6). Con (PBS), control, phosphate-buffered saline. *p < 0.05, **p < 0.01, ***p < 0.001, unpaired t test. ns, not significant.
inhibitors repressed mRNA levels of another XBP1s target gene, VEGFA. Taken together, these results suggest that inhibition of protein palmitoylation induces ER stress and inhibits XBP1s-inducible mRNAs, including XBP1 and VEGFA, thus promoting apoptosis in GBM cells upon ER stress.
The present results show that accumulation of SUMOylated XBP1 and the repression of its transcriptional function were accompanied by a reduction in XBP1 palmitoylation. SUMOylation suppresses transcription through multiple mechanisms. 45 Co-repressors are recruited through transcription factors via SUMO-dependent interactions, serving as a potential mechanism contributing to transcriptional regulation via SUMO. SUMOylation of certain co-repressors affects their transcriptional repressor activity, thus serving as another potential mechanism underlying SUMOdependent transcriptional repression. 45,46 However, the effects of palmitoylation on XBP1 SUMOylation and its transcriptional activity warrant further investigation.
In summary, the present results show that most DHHC members are upregulated in GBM and inhibition of palmitoylation results in G 2 cell cycle arrest and cell death. These effects are primarily attributed to the potential of palmitoylation inhibitors to activate ER stress signals and consequent pro-apoptotic effects. Further analysis revealed the accumulation of XBP1 SUMOylation and the suppression of its transcriptional function after treatment with palmitoylation inhibitors. Taken together, these results show that palmitoylation inhibitors potentially help treat malignant GBM.

Cell Lines and Chemical Reagents
Human GBM cell lines U87, U251, SH-SY-5Y, A172, SF126, and SW1088 were purchased from the American Type Culture Collection (Manassas, VA, USA) and characterized via DNA fingerprinting and isozyme detection. These cell lines were cultured in Dulbecco's modified Eagle's medium containing 10% fetal bovine serum (FBS) and 100 ng/mL of both penicillin and streptomycin. Cells were maintained at 37 C in a humidified atmosphere of 5% CO 2 . All cultures were free of mycoplasma.

Cell Viability Assay
Cells were seeded in 96-well plates at 4,000 cells per well and incubated overnight. Following treatment for 24 h, WST-8 reagent from Cell Counting Kit-8 (CCK-8; Beyotime Biotechnology, Haimen, China) was added to each well (final volume ratio was 10%). Absorbance was measured at 450 nm. All experiments were performed in triplicate.

EdU Assay
Cells were seeded in serum-free media overnight to allow for cell cycle synchronization and incubated with 5-ethynyl-2 0 -deoxyuridine (EdU; Cell-Light EdU Apollo 488 in vitro kit, Guangzhou RiboBio) for 2 h before fixation in 2% paraformaldehyde (PFA) and subsequent EdU detection per the manufacturer's protocol.

Cell Cycle Assay
The cell cycle status was assayed via PI staining using a Cytomics FC500 flow cytometer. The cell cycle profiles were determined using CXP analysis software (Beckman Coulter, Brea, CA, USA).

Apoptosis Assay
Cells were harvested via centrifugation and stained with annexin V and PI (V13242, Invitrogen, Carlsbad, CA, USA) in accordance with the manufacturer's instructions. Apoptotic cells were quantified via flow cytometry.

Transwell Migration and Invasion Assays
Transwell chambers with 8-mm pores (Corning Life Sciences) were utilized to assess cell migration. The Transwell membrane was previ-ously coated with Matrigel matrix (30 mL) for the invasion assay of tumor cells (1:3 mixed with PBS; BD Biosciences). Cells (1 Â 10 5 ) were loaded to the top chamber of the Transwell dish (8-mm pore size; Corning Life Sciences). FBS (10%) was loaded in the bottom chamber as a chemoattractant. In the bottom chamber, the cells were fixed and then stained with 0.005% (w/v) crystal violet. The number of invaded or migrated cells was measured by counting those five random areas in every membrane.
RNA Isolation, Real-Time PCR, and Sequencing RNA was isolated from each cell sample using TRIzol (Invitrogen) in accordance with the manufacturer's instructions. Total RNA quality and quantity were verified using a NanoDrop 2000 spectrometer (Thermo Scientific, Rockford, IL, USA).
Transcript levels were determined via quantitative real-time PCR with iQ SYBR Green supermix (Bio-Rad, Hercules, CA, USA) and performed with a qTOWER 2.0 thermal cycler (Analytik Jena, Jena, Germany). Amplification conditions were as follows: 95 C for 3 min and then 95 C for 15 s, 55 C for 15 s, and 72 C for 30 s for 40 cycles. Primers used for this study were synthesized by Invitrogen and are enlisted in Table S1.
The Bioconductor RNA-seq workflow was followed to detect differential gene expression using DESeq 2 software and other Bioconductor packages in R.

Western Blot Analyses
Cells were rinsed twice with ice-cold phosphate-buffered saline and harvested by scraping into 100 mL of radio immunoprecipitation assay buffer, followed by rocking at 4 C for 30 min. The cell lysates were then centrifuged at 14,000 Â g for 20 min at 4 C. Protein concentrations in the supernatants were quantified using a bicinchoninic acid protein assay kit. Using SDS-PAGE, proteins (50 mg) were resolved and then electrotransferred onto a polyvinylidene fluoride membrane. The membrane was blocked with 5% non-fat milk for 1 h and probed with the respective primary antibodies overnight at 4 C. The membrane was then washed three times with Tris-buffered saline/0.1% Tween 20 and probed with appropriate secondary antibodies for 1 h. An enhanced chemiluminescence detection system was used in accordance with the manufacturer's instructions.

Palmitoylation Assay
Cell lysates were prepared in the presence of 10 mM N-ethylmaleimide (23030, Pierce, Carlsbad, CA, USA), followed by denaturation with chloroform/methanol. Samples were further incubated overnight with 10 mM N-ethylmaleimide. After extensive washing with chloroform/methanol, samples were treated with hydroxylamine (HAM). After each chloroform/methanol precipitation, protein pellets were dissolved in aqueous buffers via sonication. After washing, HPDP-biotin (A35390, Pierce) was used for sulfhydryl biotinylation. Biotin-exchanged samples were affinity purified using the NeutrAvidin resin (53151, Pierce). Purified proteins were eluted by boiling in nonreducing SDS-PAGE loading buffer and detected with horseradish peroxidase-conjugated NeutrAvidin, or with dithiothreitol for antibody-based detection.

Electron Microscopy
For electron microscopy, cells treated with 2BP, Cer, or Tun were harvested via centrifugation and fixed with 2% PFA and 3% glutaraldehyde for 1 h at room temperature. The fixed cells were harvested via centrifugation, washed three times with PBS, and incubated in PBS solution overnight at 4 C. Thereafter, the cells were post-fixed in 1% OsO 4 (pH 7.4) for 2 h at room temperature, dehydrated using an ascending ethanol series, and infiltrated with 100% acetone/resin (1:5) overnight at room temperature. Thereafter, cells were embedded in fresh resin and polymerized at 55 C for 1 h. Cut sections were supported on grids (150 nm), and the slides were stained with uranyl acetate for 15 min and lead citrate for 10 min at room temperature. Images were captured using JEM-1200EX transmission electron microscope (Japan Electron Optics Laboratory, Tokyo, Japan).

Tumor Xenograft Model
Experiments were carried out using 6-to 8-week-old C57BL/6 mice weighing 22-28 g. When the tumor volume approached 100 mm 3 , the animals were randomly segregated into four groups: control, TMZ/IR, TMZ/IR+2BP, TMZ/IR+Cer, and TMZ/IR+Tun. Each group comprised six mice on day 0. These mice were treated (IR) with 20 Gy (4.5-4.6 Gy/min) for tumor induction and administered by gastric infusion with TMZ (50 mg/kg/day). 2BP (3 mg/kg), Cer (5 mg/kg), and Tun (2 mg/kg), dissolved in DMSO, were subcutaneously administered daily (distal from the sites of tumor inoculation) from days 0 to 10. The control group was treated with DMSO (or PBS) only. All animal experiments were approved by the Institutional Animal Care and Use Committee at Hefei Institutes of Physical Science, Chinese Academy of Sciences.

Statistical Analyses
A Student's t test was performed for all quantitative data between different groups. Statistical significance is indicated as *p < 0.05, **p < 0.01, and ***p < 0.001. All quantitative data are presented as means ± SD from at least three independent experiments.

AUTHOR CONTRIBUTIONS
Conception and design: X.C., H.W., and Z.F. Acquisition of data: H.L., X.F., C.Z., Z.Z., and L.H. Analysis and interpretation of data: X.C., K.Y., H.W., and H.M. Writing, review, and/or revision of the manuscript: X.C., X.F., and Z.F. Administrative, technical, or material support: L.H. and H.M. Study supervision: X.C. and K.Y. All authors read and approved the final manuscript.

CONFLICTS OF INTEREST
The authors declare no competing interests.