Thymoquinone-induced conformational changes of PAK1 interrupt prosurvival MEK-ERK signaling in colorectal cancer
Molecular Cancer volume 13, Article number: 201 (2014)
Thymoquinone (TQ) was shown to reduce tumor growth in several cancer models both in vitro and in vivo. So far only a few targets of TQ, including protein kinases have been identified. Considering that kinases are promising candidates for targeted anticancer therapy, we studied the complex kinase network regulated by TQ.
Novel kinase targets influenced by TQ were revealed by in silico analysis of peptide array data obtained from TQ-treated HCT116wt cells. Western blotting and kinase activity assays were used to determine changes in kinase expression patterns in colorectal cancer cells (HCT116wt, DLD-1, HT29). To study the viability/apoptotic effects of combining the PAK1 inhibitor IPA-3 and TQ, crystal violet assay and AnnexinV/PI staining were employed. Interactions between PAK1 and ERK1/2 were investigated by co-immunoprecipitation and modeled by docking studies. Transfection with different PAK1 mutants unraveled the role of TQ-induced changes in PAK1 phosphorylation and TQ´s effects on PAK1 scaffold function.
Of the 104 proteins identified, 50 were upregulated ≥2 fold by TQ and included molecules in the AKT-MEK-ERK1/2 pathway. Oncogenic PAK1 emerged as an interesting TQ target. Time-dependent changes in two PAK1 phosphorylation sites generated a specific kinase profile with early increase in pPAKThr212 followed by late increase in pPAKThr423. TQ induced an increase of pERK1/2 and triggered the early formation of an ERK1/2-PAK1 complex. Modeling confirmed that TQ binds in the vicinity of Thr212 accompanied by conformational changes in ERK2-PAK1 binding. Transfecting the cells with the non-phosphorylatable mutant T212A revealed an increase of pPAKThr423 and enhanced apoptosis. Likewise, an increase in apoptosis was observed in cells transfected with both the kinase-dead K299R mutant and PAK1 siRNA. Using structural modeling we suggest that TQ interferes also with the kinase domain consequently disturbing its interaction with pPAKThr423, finally inhibiting MEK-ERK1/2 signaling and disrupting its prosurvival function. pERK1/2 loss was also validated in vivo.
Our study shows for the first time that the small molecule TQ directly binds to PAK1 changing its conformation and scaffold function. Because TQ affects the central RAF/MEK/ERK1/2 pathway, the combination of TQ with targeted therapies is worth considering for future anticancer treatments.
Colorectal cancer (CRC) affects yearly more than 1 million people worldwide , therefore there is constant need to achieve more effective cures. Over the recent years naturally occurring compounds have received increasing attention because of their anticancer effects . Thymoquinone (TQ), the active compound extracted from Nigella sativa, is a very promising anticancer drug, whether used separately or in combination with conventional medicines [2–7]. Interestingly, TQ was found to have only limited toxicity to normal intestinal cells in vitro and not to affect the survival or body weight of animals when used at doses up to 25 mg/kg in colorectal cancer animal models [9, 10]. Previous work documented TQ´s growth inhibitory and apoptosis triggering effects on colon and other solid tumors such as uterine sarcoma, breast, and pancreatic cancer in a dose- and time-dependent manner [9, 11–14]. Moreover, TQ was shown to reduce tumor growth and to induce apoptosis in various murine cancer models [5–7, 10, 15, 16]. So far, the anticancer mechanism of TQ is not fully understood; however, several modes of action have been described depending on the stimulus and the cellular context .
Protein kinases, the protagonists of phosphorylation, commonly work in complex networks and have become novel candidates for targeted therapy. There are a few studies describing the involvement of TQ in the regulation of kinases i.e. AKT1  , JAK2 , JNK , IKKb, ERK2 , CHEK1  andPlk1 . Most of these kinases are frequently deregulated in colon cancer. In the last 10 years several kinase inhibitors have been clinically tested or are currently undergoing clinical trials in targeted therapy . Such therapy blocks the growth and spread of cancer by specifically interfering with molecules affecting tumor cell proliferation and progression. In colorectal cancer we have shown that kinases such as c-Jun N-terminal kinase (JNK) and the prosurvival extracellular-signal-regulated kinase (ERK) were activated under TQ treatment before the onset of apoptosis . The TQ resistance observed in HCT116 p53−/− cells was associated with an upregulation of CHEK1, a serine/threonine kinase cell cycle checkpoint mediator, indicating an inefficient defense mechanism against TQ-induced DNA damage when p53 is absent .
So far despite a few reports, the interaction of TQ with the kinase network of colon cancer cells is only poorly understood. Kinase signaling has been extensively investigated in the last decade and new techniques of quantitative proteomics have been developed simplifying the study of complex systems. Thus we performed a quantitative phospho-proteomic analysis to identify new targets in TQ- treated HCT116wt colorectal cancer cells. For the first time we show a direct binding of the plant derived small molecule TQ to the PAK1/ERK1/2 kinase complex leading to conformational changes in protein structures and consequently to apoptosis.
Results and discussion
Identification of new TQ targets by peptide array analysis
Kinome analysis allows the characterization and quantification of the phosphorylation profile of a given kinase target under various experimental conditions. Using a peptide array, the treatment of colorectal HCT116wt cells with 40 μM TQ for 24 h led to the identification of 104 proteins with a significant phosphorylation upregulation, among which were 50 proteins and kinases upregulated by ≥2 fold (out of 1152 kinase substrate peptides). Many of these proteins were previously described as TQ targets (such as p21Cip1, p53, ERK or NFkB) verifying the feasibility of the array data (Additional file 1: Table S1). The analysis of the phosphorylation motifs of all 104 candidates (Figure 1A) showed that TQ had a propensity for serine phosphorylation as well as an 80% probability to induce phosphorylation of a neutral polar amino acid (S group). In silico analysis revealed that the most influenced pathways (Figure 1B) and networks (Figure 1C,1D) were those involved in cancer, cell cycle, cell death, and survival mechanisms. Moreover 24 of the top 50 candidate proteins were grouped into the cancer-related networks “cytoskeleton”, “PI3K/AKT” and “Wnt signaling” (Table 1). The Venn-diagram for the most relevant biological functions identified 11 TQ-modulated proteins that are common in apoptosis, proliferation, and inflammation pathways. Besides the epigenetic marker CEBPB, the analysis included key molecules in the AKT-MEK-ERK pathway, a central signaling network for current targeted therapies (Figure 1E). With a fold-change of 2.12, p21 protein (Cdc42/Rac)-activated kinase 1 (PAK1) was an appealing candidate considering its role in cell growth, invasion, cell migration, cell survival, mitosis and cytoskeletal remodeling of cancer cells [21, 22]. Recently PAK1 was considered a novel subject for targeted cancer therapeutic approaches . Furthermore for the first time, PAK1 has been identified as a TQ target. An extensive pathway mapping of genes differentially phosphorylated by TQ at 24 h and subsequent pathway annotation clustering resulted in 5 significant pathway based clusters. It is interesting that PAK1 was mostly associated with AKT1 and RAF1 (Table 2) in all the pathways in the Annotation cluster3. The top 4 pathways obtained through clustering approach (T-cell receptor signaling, angiogenesis, MAPK and chemokine signaling) also included PAK1 making it the target of our more detailed study of TQ action.
TQ effect on PAK1 phosphorylation
In HCT116wt cells, the TQ-induced activation of pPAK1Thr212 observed in the array analysis was verified by western blotting (Figure 2A). Although PAK1 has several different phosphorylation sites, only PAK1Thr212 was spotted on the PEPSCAN peptide array. In addition, we investigated two other well-studied phosphorylation sites: pPAK1Thr423 and pPAK1Ser144, by western blotting. Time-dependent changes in the expression patterns of pPAK1Thr212 and pPAK1Thr423 generated a specific kinase profile at early and late time points after TQ exposure. While pPAK1Thr212 showed a fast induction at 1 and 3 h, pPAK1Thr423, found within the activation loop of PAK1, increased at 24 h (Figure 2A). In contrast, pPAK1Ser144 and total PAK1 levels did not change over time (Figure 2A). Investigating another colorectal cancer cell line, DLD-1 showed the same profile of pPAK1Thr212 and pPAK1Thr423 regulation over time (Figure 2B). As expected, the TQ resistant HT29 cells  showed a different kinase profile where both pPAK1Thr212 and pPAK1Thr423 levels increased at the same time from 3 to 6 h and decreased at 24 h (Figure 2C).
TQ is a possible PAK1 inhibitor
To further understand the role of PAK1 in TQ-induced apoptosis, we used 1,1′-disulfanediyldinaphthalen-2-ol (IPA-3), the allosteric inhibitor of PAK1. IPA-3 prevents pPAK1Thr423 phosphorylation by targeting the auto-regulatory site and disrupting the functional interaction between PAK1 and Cdc42 in addition to inhibiting the activity of other kinases [23, 24]. Here we show that TQ did not change the activity of AKT2, GSKα, GSKβ and p38, suggesting that this molecule might be more specific than IPA-3 in interacting with PAK1 (Additional file 2: Figure S1). In addition, we performed a crystal violet viability assay and obtained an IC50 of 141 μM ±30.7 in normal intestinal epithelial HCEC cells in comparison to an IC50 of 50 μM ± 6 in HCT116 tumor cells (Additional file 3: Figure S2), thus confirming the limited toxicity of TQ to normal cells.
Since IPA-3 does not inhibit already activated PAK1, HCT116wt cells were pre-incubated with IPA-3 for 1 h and then stimulated with 40 μM TQ (Figure 3). Crystal violet staining showed that 30% of the cells were dead after 10 μM IPA-3 (Figure 3A). Interestingly, the combination of TQ and IPA-3 (10 μM) caused significantly more cell death and decreased cell viability by 70% (Figure 3A). AnnexinV/PI staining revealed that most dead cells were in late apoptosis (Figure 3B). TQ induced the “apoptotic” cleavage of PARP (89 kDa fragment) at 24 h, whereas in response to IPA-3 the cleavage was stronger and occurred earlier at 6 h (Figure 3C). Furthermore, transfection with kinase-dead dominant negative K299R mutant led to a larger increase in PARP cleavage showing more apoptosis when the kinase activity is lost (Figure 3D). Similarly, PAK1 siRNA transfection of cells followed by TQ treatment induced an increase in the pro-apoptotic response (Figure 3E).
In a next step, we determined changes in the expression of the different phosphorylation sites of PAK1 in IPA-3 and/or TQ treated cells (Figure 3F). As expected, pPAK1Thr423 protein level decreased when IPA-3 was combined with TQ at 3 and 6 h but this decrease was not sustained at 24 h. To further understand this finding, we modeled the interaction between IPA-3, PAK1 and TQ. Docking IPA-3 to PAK1 (Additional file 4: Figure S3A) revealed that IPA-3 binds to the CRIB motif (75–90) in the autoregulatory region of PAK1 dimer in a mode that enhances the interaction between both monomers (binding site first monomer: His83 and Gly98, second monomer: Met99). When TQ is docked into the IPA-3-bound conformation, IPA-3 still binds with both monomers. TQ also interacts with residues of both PAK1 monomers in the auto-inhibited dimer conformation (Additional file 4: Figure S3B and Additional file 5: Table S2). In this binding mode, TQ interaction with the kinase inhibitory segment of PAK1 (Asn38 residue) possibly causes the masking of the PAK1 activation loop and prevents the action of IPA-3 on Thr423 (Additional file 4: Figure S3).
The phosphorylation of pPAK1Thr212 in TQ and IPA-3 co-treated cells showed a remarkable increase in comparison to TQ alone (Figure 3F). This was surprising considering that IPA-3 only interrupts the interaction between Cdc42 and the autoregulatory site at Thr423 residue of PAK1 [24, 25]. Knowing that pPAK1Thr212 is a major target of ERK2 , we further studied the role of ERK1/2/PAK1 interaction in response to TQ.
TQ triggers conformational changes in PAK1 and induces ERK1/2-PAK1 binding
In a next step we investigated the effect of TQ on ERK1/2 in HCT116wt cells. Kinase activity assay revealed that ERK1/2 activity significantly increased after 1 h of TQ treatment (Figure 4A). This was confirmed by western blotting, where pERK1/2 level was elevated at the same time point followed by a later decrease to normal levels (Figure 4B). In addition to the previously published increase in TUNEL and cleaved caspase-3 , immunohistochemical staining of the prosurvival pERK1/2 was significantly (p < 0.05) lost in the nuclei of the TQ treated tumor xenografts (Figure 4C and Additional file 6: Figure S4).
Furthermore, time course analysis of pPAKThr212 and pERK1/2 levels at early time points revealed that pPAKThr212 increase can be observed as early as 30 min whereas pERK1/2 increase started after 45 min of TQ treatment (Additional file 7: Figure S5). This suggests an involvement of other kinases in the phosphorylation of PAK1 at the Thr212 site in response to TQ. Interestingly, co-immunoprecipitation studies showed that TQ induced the formation of a PAK1 and ERK1/2 complex (Figure 4D). To better understand the nature of ERK1/2/PAK1 interaction in response to TQ, we modeled the PAK1-TQ interaction by docking PAK1 and TQ and identified a possible ligand binding site at the vicinity of Thr212 (Figure 5A) with a docking score of −2.316 (Additional file 5: Table S2) [27–30]. The energy of ERK2 bound to PAK1 conformation changed in the presence of TQ from -21883 kcal/mol to −22076 kcal/mol suggesting that TQ strengthens ERK2 binding to PAK1 . It may also be noted from Figure 5B and Figure 5C that ERK2 binds in a different mode to PAK1 in the presence of TQ. ERK2 binding to a different conformation of PAK1 may prevent ERK2 from phosphorylating the Thr212 of PAK1. This could explain why pPAK1Thr212 levels decrease over time in response to TQ (Figure 2A). Furthermore, TQ-induced pERK1/2 levels were not reduced by IPA-3 treatment, instead ERK1/2 activity was further enhanced from 1 h to 6 h when TQ was combined with IPA-3 (Figure 4B), suggesting the involvement of other upstream kinases in phosphorylating ERK1/2. If the TQ-triggered closer binding between ERK1/2 and PAK1 is disrupted under IPA-3 this might explain the slight increase in pPAK1Thr212 levels when TQ is combined with IPA-3 (Figure 3F).
From these findings we propose a functional relation between both phosphorylation sites of PAK1. Indeed, T212A mutant induced a significant increase in pPAK1Thr423 levels whereas the hyperphosphorylated T212E mutant showed a lesser extent of increase (Figure 4E). As presented in Additional file 5: Table S2, TQ binds in the vicinity of Thr212 and no binding occurs next to the 423 site. Therefore we speculate that the chain reaction orchestrated by TQ starts at Thr212 site, which correlates with the western blot pattern (Figure 2A). Furthermore, T423E resulted in a decrease in ERK1/2 phosphorylation (Figure 4F) suggesting an impaired interference of Thr423 residue with the kinase domain of PAK1 under TQ treatment. We analysed the structural changes induced by TQ on PAK1 catalytic site and activation loop. Figure 4G and Figure 4H depict the structural conformation of the catalytic domain of PAK1 in the absence (Figure 4G) or presence of TQ (Figure 4H). The kinase domain shows a root-mean-square deviation of 0.38 Å accompanied by rearrangements in the main activation loop. Hydrogen bond analysis (Additional file 8: Table S3) showed that a higher number of hydrogen bonds are formed in the catalytic loop region in the presence of TQ involving residues Arg388 and Arg421 which are known to interact with Thr423 for the catalytic activity. Finally this results in a disturbed interaction between the Thr423 site and the catalytic kinase domain which inhibits the PAK1 kinase activity and its prosurvival signaling (Figure 6). There are three other lines of evidence for this hypothesis: first, early TQ-induced ERK1/2 activation is inhibited when PAK1Thr423 is maximally phosphorylated at 24 h. Second, when IPA-3 is combined with TQ there is a decrease in Thr423 phosphorylation at early time points accompanied by a significant upregulation of pERK1/2 levels that confirms an early activation of MEK-ERK signaling. The lack of inhibition of pPAK1Thr423 at 24 h is closely associated with a decrease in prosurvival ERK1/2 activation and enhanced apoptosis induction at 24 h. This can be explained by another structural modeling showing that TQ has the ability to effectively bind to the autoregulatory domain of PAK1 preventing IPA-3 ability to interrupt the interaction between Cdc42 and pPAK1Thr423. Third, the pPAK1Thr212 upregulation is not as dramatic as the pERK1/2 activation after combined TQ and IPA-3 treatment reinforcing the inhibitory loop for PAK1 activation. Although pERK1/2 should disappear completely by the higher level of pPAK1Thr423 at 24 h, it is worth mentioning that other secondary kinase loops could be at play. Consequently the phosphorylation status of Thr423 and PAK1 kinase activity is affected by TQ proving the experimental observations reported here.
PAK1 has been described as a scaffold for ERK1/2/MEK to recruit MEK to RAF at the membrane  thus facilitating signaling through the ERK1/2 pathway [26, 33]. We have shown previously a prosurvival function of ERK1/2 under TQ, since inhibition of ERK1/2 by the MEK inhibitor PD98059 resulted in enhanced apoptosis . In the present study and for the first time we document that a small natural molecule, TQ, directly binds to PAK1/ERK kinase complex, induces considerable conformational changes of PAK1 and interrupts its scaffold function. Upon TQ the interference of Thr423 phosphorylation with the kinase domain is disrupted inhibiting the PAK1 kinase activity and its prosurvival signaling . TQ-induced loss of ERK1/2 dependent phosphorylation at Thr212 of PAK1 provides another negative feedback of prosurvival PAK1/MEK/ERK1/2 signaling. The loss of nuclear pERK1/2 in TQ-treated mouse xenografts confirms the in vivo relevance of our finding. Overall, our approach allowed an identification of several novel TQ targets which could be used in the future to evaluate the therapeutic benefit of combination therapies.
Cell lines and reagents
Human colon carcinoma HCT116wt, DLD-1 and HT29 cells were grown in RPMI 1640 medium (PAA Laboratories GmbH, Pasching Austria) supplemented with 10% fetal bovine serum (FBS) and 1% penicillin-streptomycin (PANTM Biotech GmbH, Germany) and kept at 37°C in a humidified incubator (95% air, 5%CO2). Cells were treated at 50% confluency with 40 or 60 μM TQ (Sigma-Aldrich) dissolved in DMSO (final concentration less than 0.1%) and collected at different time points. PAK1 inhibition was performed by preincubating the cells with 10 μM IPA-3 (Sigma-Aldrich) for 1 hour.
Peptide Array PepChip Kinomics v2
HCT116wt cells (1.5x105/ml) were treated with 40 μM TQ and collected after 24 hours. Proteins were extracted and shipped to manufacturer (Pepscan, Netherland) according to manufacturer protocol. Briefly, cell pellets were lysed (20 mM Tris–HCl (pH 7.5), 150 mM NaCl, 1 mM EDTA, 1 mM EGTA, 1% Triton X-100, 2.5 mM sodium pyrophosphate, 1 mM beta-glycercophosphate, 1 mM sodium orthovanadate (Na3VO3), 1 mM NaF, Roche Complete protease inhibitor cocktail). After sonication, the lysates were spin down and snap frozen before being shipped to the manufacturer. The concentration of the samples was determined using Bradford Protein Assay (Bio-Rad, Hercules, CA). Samples were then hybridized on PepChip Kinomics v2 chip and the kinome profiling was performed by the Pepscan Company. The peptide array was done with samples obtained from two independent experiments and every sample was spotted two times on the array.
Array Analysis software Ingenuity pathway analysis (IPA)
Peptide array data was analyzed by Ingenuity pathway analysis software (IPA, version 9.0, Ingenuity © Systems, http://www.ingenuity.com, Mountain View, CA, USA). The data set contained upregulated proteins (fold-change ≥ 2 after 24 hours of 40 μM TQ treatment) identified by the pepscan peptide array and their protein ID number. Statistical significance of the different pathways obtained by IPA was calculated using a right-tailed Fisher’s Exact test . Pathways with a p–value ≤ 0.05 were selected. WebLogo v3 application (http://weblogo.threeplusone.com/) was used to analyze the phosphorylation consensus motifs of all the kinases and proteins identified by the array.
Pathway Annotation Clustering
The probable pathways by which the 50 candidate proteins, with fold-change ≥ 2 after 24 hours of 40 μM TQ treatment, could interfere were predicted using DAVID bioinformatics tool (http://david.abcc.ncifcrf.gov) [35, 36]. For pathway mapping, Uniprot IDs of differentially phosphorylated proteins were submitted to the DAVID web resource with colorectal cancer specific proteins from Human Protein Atlas (http://www.proteinatlas.org/) as background. Pathway predictions were made by including the subset pathway databases: KEGG, Panther, Biocarta, Reactome and BBID followed by pathway annotation clustering.
HCT116wt, DLD-1 and HT29 cells were grown to 50% confluency and treated with 40 μMTQ and harvested over different time points. The cells were then lysed with RIPA lysis buffer (150 mM NaCl, 0.5% DOC, 0.1% SDS, 1% NP40, 50 mM Tris pH8) containing a cocktail of protease inhibitors. The protein concentration of all samples was determined by by DC BioRad protein assay kit (BioRad Laboratories, Hercules, CA) using bovine serum albumin as standard. Total proteins (50 μg /sample) were separated on an 8–12% SDS polyacrylamide gel and transferred to nitrocellulose membrane by blotting. After blocking with 5% non-fat dry milk in TBST buffer, membranes were incubated with primary antibody at 4°C overnight, washed three times, with TBST buffer, and incubated again with the corresponding HRP-conjugated secondary antibody at room temperature for 1 h. The membranes were then washed with TBST buffer and protein bands were detected by enhanced chemiluminescence. PAK1, ERK1/2, pERK1/2 and PARP were purchased from Cell signaling; pPAKThr212 and pPAKSer144 from Abcam and pPAKThr423 from Abgent. Densitometric analysis was performed for all western blots using ImageJ 1.45 s software.
Cells treated with 60 μM TQ for 10, 45, 60 minutes and 24 hours were collected and lysed using RIPA buffer supplemented with proteases and phosphatases inhibitors. Immunoprecipitiation was performed using the Dynabeads Protein G magnetic separation KIT as per manufacturer protocol (Invitrogen). After incubating the beads with 600 μg of proteins, they were incubated with anti-PAK1 and anti-ERK1/2 antibodies. The precipitated lysates were then loaded into SDS-PAGE gels and immunoblotted against total PAK1 and ERK1/2.
Plasmids and siRNA transfections
Plasmids (PAK1wt, K299R, T212E, T212A, T423E) were a gift from Prof. Jonathan Chernoff (Fox Chase Cancer Center, 333 Cottman Avenue, Philadelphia). Cells seeded in 6-well plates were transfected at 90% confluency with 1 μg of the different plasmids for 6 hours, using Invitrogen Lipofectamine 2000 according to manufacturer instructions. Afterwards cells were treated with 100 μM TQ (due to the high cellular density) and collected after 24 hours. Cell lysates were then used for Western blotting analysis.
PAK1 and scrambled siRNA reagents were purchased from Thermo Fisher Scientific (Dharmacon RNAi Technologies). Lipofectamine RNAiMAX reagent was used to transfect HCT116 cells with 10 μM siRNA, according to manufacturer’s instructions. The next day the transfected cells were split and 50% confluent plates were treated with 60 μM TQ. Cell lysates were then used for Western blotting analysis.
Structural analysis by docking
To understand the interaction between PAK1 and TQ, we modeled the structure of PAK1 and performed molecular docking with TQ using the Schrodinger suite (Maestro, version 9.3, Schrodinger, Inc, New York, NY, 2012). A model of the possible active state conformation of PAK1 was obtained by integrating the available crystal structures of the autoinhibitory domain 78-147(PDB ID: 1F3M)  and the kinase domain 250-542(PDB ID: 3Q52)  of PAK1 and a threaded model of the region corresponding to 148–249 using ITASSER [29, 30]. The possible ligand binding sites in PAK1 was analyzed using the SiteMap module (SiteMap, version 2.6, Schrödinger, Inc., New York, NY, 2012) of the Schrodinger suite. In order to analyze the binding mode of TQ near to Thr212 site, grids were generated focusing on Thr212 site and the interacting residues were identified. To study the interaction details of PAK1 and ERK2, the crystal structure of ERK2 available with PDB ID:2ERK  was docked to the TQ bound and unbound conformations of PAK1 using ClusPro [37, 38]. The energies of the complexes were determined through energy minimization using AMBER11  . Crystal structure of the autoinhibited dimer of PAK1 (PDB ID: 1F3M) was used to understand the binding mode of IPA-3 on PAK1 . The structural coordinates of the residues (416–422), missing in the PDB, belong to the kinase activation segment of PAK1 and were modeled using Modeller9v12  . IPA-3 was then docked to the autoregulatory region of PAK1 (autoinhibited conformation). To analyze the combined effect of TQ and IPA-3 on PAK1, TQ was docked to the IPA-3 bound conformation of PAK1. All the dockings were performed using the extra precision mode of the Glide program of Schrodinger (Glide, version 5.8, Schrödinger, Inc., New York, NY, 2012). The structures of TQ and IPA-3 were obtained from the PubChem database. Protein and ligand structures for docking were prepared using Protein Preparation Wizard and LigPrep utilities of Schrodinger (LigPrep, version 2.5, Schrödinger, Inc., New York, NY, 2012). All renderings were done using Chimera 1.8 .
Kinase Activity assay
Multi-kinase ELISA array was performed as per the manufacturer protocol (Symansis). Briefly, Cells treated with 40 μM TQ were collected after 1, 3, 6, 24 hours and lysed with 1X denaturing cell lysis buffer (Symansis CLB001). 30 μg of proteins were loaded on pre-coated strips with antibodies corresponding to the investigated kinases. After several washing and hybridization steps, absorbance of each well was measured at 450 nm using VICTORTM X3 multilabel reader.
Cell Viability assay
Cellular viability was measured by crystal violet staining. HCT116wt (3.75 × 104/ml) cells were seeded in 96 well plate and treated with 40 μM TQ and/or 1-50 μM IPA-3 for 24 hours. The treated cells were washed once with PBS, fixed for 15 min in a crystal violet solution (0.5% crystal violet in 20% methanol) at room temperature, then washed twice with water and air-dried. The stained cells were solubilized with methanol for 15 min with mild agitation. Absorbance was measured at 595 nm using VICTORTM X3 multilabel reader. IC50 values were calculated using EXCEL 2010. Statistical analysis was performed using SPSS version 19. One tailed Student T-Test was performed by comparing 2 samples assuming that they have equal variances.
Annexin V/Propidium Iodide staining
Apoptosis was measured using Annexin V/PI co-staining. HCT116wt cells (1.5 × 105/ml) were pretreated with IPA-3 (10 μM) for 1 hour then treated with TQ (40 μM) for 24 hours. After collection, cells were centrifuged at 200 g for 5 min, 4°C and washed with 1X PBS. The pellet was resuspended in 100 μl Annexin-V-Fluos labeling solution (10 μl annexin reagent and 10 μl PI solution in 150 μl incubation buffer (according to manufacturer, Roche). The samples were incubated for 7 min in the dark, at room temperature then 100 μl incubation buffer was added. The cellular fluorescence was then measured using a Fluorescence Activated Cell Sorter (FACS) flow cytometer (BDFACS CantoTM II). Each sample was collected as 20,000 ungated events and the different cell populations were determined using FlowJo software (FlowJo7.6.5).
Immunohistochemistry (IHC) was performed on available paraffin fixed tissue blocks from a xenograft experiment  to detect the expression of pERK1/2. Rehydration of tissue sections was performed in descending concentrations solutions of ethanol (96% to 70%). Antigen was retrieved by heating in a pressure cooker (1 mmol/L Tris-EDTA buffer, 120°C, 5 min). The slices were incubated in blocking solution (Dako, Glostrup, Denmark) to prevent nonspecific binding sites. Next pERK1/2 (1:2000) primary antibody was added to the slices and incubated for 30 min at room temperature. The sections were then washed with washing buffer (Dako) and incubated with EnVision + System horseradish peroxidase-linked secondary antibody (goat anti-rabbit, Dako) at room temperature for 30 min. Positive immunoreactivity was detected using diaminobenzidine + (Dako). Positive and negative IHC controls were included in this study. Percentage and intensity of positively stained epithelial cells (cytoplasmic versus nuclei) was quantified and scored by an expert pathologist (T.T.R.). Statistical analysis was performed using SPSS. Two tailed student t-test was performed by comparing 2 samples assuming that they have equal variances.
Protein kinase B
- Plk1 PBD:
Serine/threonine kinase Polo-like kinase 1
p21 protein (Cdc42/Rac)-activated kinase 1
Cunningham D, Atkin W, Lenz HJ, Lynch HT, Minsky B, Nordlinger B, Starling N: Colorectal cancer. Lancet. 2010, 375: 1030-1047. 10.1016/S0140-6736(10)60353-4
Schneider-Stock R, Fakhoury IH, Zaki AM, El-Baba CO, Gali-Muhtasib HU: Thymoquinone: fifty years of success in the battle against cancer models. Drug Discov Today. 2014, 19: 18-30. 10.1016/j.drudis.2013.08.021
Nessa MU, Beale P, Chan C, Yu JQ, Huq F: Synergism from combinations of cisplatin and oxaliplatin with quercetin and thymoquinone in human ovarian tumour models. Anticancer Res. 2011, 31: 3789-3797.
Woo CC, Kumar AP, Sethi G, Tan KH: Thymoquinone: potential cure for inflammatory disorders and cancer. Biochem Pharmacol. 2012, 83: 443-451. 10.1016/j.bcp.2011.09.029
Jafri SH, Glass J, Shi R, Zhang S, Prince M, Kleiner-Hancock H: Thymoquinone and cisplatin as a therapeutic combination in lung cancer: In vitro and in vivo. J Exp Clin Cancer Res. 2010, 29: 87- 10.1186/1756-9966-29-87
Banerjee S, Kaseb AO, Wang Z, Kong D, Mohammad M, Padhye S, Sarkar FH, Mohammad RM: Antitumor activity of gemcitabine and oxaliplatin is augmented by thymoquinone in pancreatic cancer. Cancer Res. 2009, 69: 5575-5583. 10.1158/0008-5472.CAN-08-4235
Lei X, Lv X, Liu M, Yang Z, Ji M, Guo X, Dong W: Thymoquinone inhibits growth and augments 5-fluorouracil-induced apoptosis in gastric cancer cells both in vitro and in vivo. Biochem Biophys Res Commun. 2012, 417: 864-868. 10.1016/j.bbrc.2011.12.063
El-Najjar N, Chatila M, Moukadem H, Vuorela H, Ocker M, Gandesiri M, Schneider-Stock R, Gali-Muhtasib H: Reactive oxygen species mediate thymoquinone-induced apoptosis and activate ERK and JNK signaling. Apoptosis. 2010, 15: 183-195. 10.1007/s10495-009-0421-z
Gali-Muhtasib H, Kuester D, Mawrin C, Bajbouj K, Diestel A, Ocker M, Habold C, Foltzer-Jourdainne C, Schoenfeld P, Peters B: Thymoquinone triggers inactivation of the stress response pathway sensor CHEK1 and contributes to apoptosis in colorectal cancer cells. Cancer Res. 2008, 68: 5609-5618. 10.1158/0008-5472.CAN-08-0884
Gali-Muhtasib H, Ocker M, Kuester D, Krueger S, El-Hajj Z, Diestel A, Evert M, El-Najjar N, Peters B, Jurjus A: Thymoquinone reduces mouse colon tumor cell invasion and inhibits tumor growth in murine colon cancer models. J Cell Mol Med. 2008, 12: 330-342.
Gali-Muhtasib H, Diab-Assaf M, Boltze C, Al-Hmaira J, Hartig R, Roessner A, Schneider-Stock R: Thymoquinone extracted from black seed triggers apoptotic cell death in human colorectal cancer cells via a p53-dependent mechanism. Int J Oncol. 2004, 25: 857-866.
Gali-Muhtasib HU: Abou Kheir WG, Kheir LA, Darwiche N, Crooks PA: Molecular pathway for thymoquinone-induced cell-cycle arrest and apoptosis in neoplastic keratinocytes. Anticancer Drugs. 2004, 15: 389-399. 10.1097/00001813-200404000-00012
Yi T, Cho SG, Yi Z, Pang X, Rodriguez M, Wang Y, Sethi G, Aggarwal BB, Liu M: Thymoquinone inhibits tumor angiogenesis and tumor growth through suppressing AKT and extracellular signal-regulated kinase signaling pathways. Mol Cancer Ther. 2008, 7: 1789-1796. 10.1158/1535-7163.MCT-08-0124
Woo CC, Hsu A, Kumar AP, Sethi G, Tan KH: Thymoquinone inhibits tumor growth and induces apoptosis in a breast cancer xenograft mouse model: the role of p38 MAPK and ROS. PLoS One. 2013, 8: e75356- 10.1371/journal.pone.0075356
Attoub S, Sperandio O, Raza H, Arafat K, Al-Salam S: Al Sultan MA, Al Safi M, Takahashi T, Adem A: Thymoquinone as an anticancer agent: evidence from inhibition of cancer cells viability and invasion in vitro and tumor growth in vivo. Fundam Clin Pharmacol. 2013, 27: 557-569. 10.1111/j.1472-8206.2012.01056.x
Siveen KS, Mustafa N, Li F, Kannaiyan R, Ahn KS, Kumar AP, Chng WJ, Sethi G: Thymoquinone overcomes chemoresistance and enhances the anticancer effects of bortezomib through abrogation of NF-kappaB regulated gene products in multiple myeloma xenograft mouse model. Oncotarget. 2014, 5: 634-648.
Li F, Rajendran P, Sethi G: Thymoquinone inhibits proliferation, induces apoptosis and chemosensitizes human multiple myeloma cells through suppression of signal transducer and activator of transcription 3 activation pathway. Br J Pharmacol. 2010, 161: 541-554. 10.1111/j.1476-5381.2010.00874.x
Sethi G, Ahn KS, Aggarwal BB: Targeting nuclear factor-kappa B activation pathway by thymoquinone: role in suppression of antiapoptotic gene products and enhancement of apoptosis. Mol Cancer Res. 2008, 6: 1059-1070. 10.1158/1541-7786.MCR-07-2088
Yin Z, Song Y, Rehse PH: Thymoquinone blocks pSer/pThr recognition by Plk1 Polo-box domain as a phosphate mimic. ACS Chem Biol. 2013, 8: 303-308. 10.1021/cb3004379
Fabbro D, Cowan-Jacob SW, Mobitz H, Martiny-Baron G: Targeting cancer with small-molecular-weight kinase inhibitors. Methods Mol Biol. 2012, 795: 1-34. 10.1007/978-1-61779-337-0_1
Zhao ZS, Lim JP, Ng YW, Lim L, Manser E: The GIT-associated kinase PAK targets to the centrosome and regulates Aurora-A. Mol Cell. 2005, 20: 237-249. 10.1016/j.molcel.2005.08.035
Eswaran J, Li DQ, Shah A, Kumar R: Molecular pathways: targeting p21-activated kinase 1 signaling in cancer–opportunities, challenges, and limitations. Clin Cancer Res. 2012, 18: 3743-3749. 10.1158/1078-0432.CCR-11-1952
Kichina JV, Goc A, Al-Husein B, Somanath PR, Kandel ES: PAK1 as a therapeutic target. Expert Opin Ther Targets. 2010, 14: 703-725. 10.1517/14728222.2010.492779
Wong LL, Lam IP, Wong TY, Lai WL, Liu HF, Yeung LL, Ching YP: IPA-3 inhibits the growth of liver cancer cells by suppressing PAK1 and NF-kappaB activation. PLoS One. 2013, 8: e68843- 10.1371/journal.pone.0068843
Zenke FT, King CC, Bohl BP, Bokoch GM: Identification of a central phosphorylation site in p21-activated kinase regulating autoinhibition and kinase activity. J Biol Chem. 1999, 274: 32565-32573. 10.1074/jbc.274.46.32565
Sundberg-Smith LJ, Doherty JT, Mack CP, Taylor JM: Adhesion stimulates direct PAK1/ERK2 association and leads to ERK-dependent PAK1 Thr212 phosphorylation. J Biol Chem. 2005, 280: 2055-2064.
Lei M, Lu W, Meng W, Parrini MC, Eck MJ, Mayer BJ, Harrison SC: Structure of PAK1 in an autoinhibited conformation reveals a multistage activation switch. Cell. 2000, 102: 387-397. 10.1016/S0092-8674(00)00043-X
Wang J, Wu JW, Wang ZX: Structural insights into the autoactivation mechanism of p21-activated protein kinase. Structure. 2011, 19: 1752-1761. 10.1016/j.str.2011.10.013
Zhang Y: I-TASSER server for protein 3D structure prediction. BMC Bioinformatics. 2008, 9: 40- 10.1186/1471-2105-9-40
Roy A, Kucukural A, Zhang Y: I-TASSER: a unified platform for automated protein structure and function prediction. Nat Protoc. 2010, 5: 725-738. 10.1038/nprot.2010.5
Canagarajah BJ, Khokhlatchev A, Cobb MH, Goldsmith EJ: Activation mechanism of the MAP kinase ERK2 by dual phosphorylation. Cell. 1997, 90: 859-869. 10.1016/S0092-8674(00)80351-7
Wang Z, Fu M, Wang L, Liu J, Li Y, Brakebusch C, Mei Q: p21-activated kinase 1 (PAK1) can promote ERK activation in a kinase-independent manner. J Biol Chem. 2013, 288: 20093-20099. 10.1074/jbc.M112.426023
Howe AK, Aplin AE, Juliano RL: Anchorage-dependent ERK signaling–mechanisms and consequences. Curr Opin Genet Dev. 2002, 12: 30-35. 10.1016/S0959-437X(01)00260-X
Zhan X, Desiderio DM: Signaling pathway networks mined from human pituitary adenoma proteomics data. BMC Med Genomics. 2010, 3: 13- 10.1186/1755-8794-3-13
Da Huang W, Sherman BT, Lempicki RA: Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc. 2009, 4: 44-57.
Da Huang W, Sherman BT, Lempicki RA: Bioinformatics enrichment tools: paths toward the comprehensive functional analysis of large gene lists. Nucleic Acids Res. 2009, 37: 1-13. 10.1093/nar/gkn923
Comeau SR, Gatchell DW, Vajda S, Camacho CJ: ClusPro: a fully automated algorithm for protein-protein docking. Nucleic Acids Res. 2004, 32: W96-99. 10.1093/nar/gkh354
Comeau SR, Gatchell DW, Vajda S, Camacho CJ: ClusPro: an automated docking and discrimination method for the prediction of protein complexes. Bioinformatics. 2004, 20: 45-50. 10.1093/bioinformatics/btg371
Gohlke H, Kiel C, Case DA: Insights into protein-protein binding by binding free energy calculation and free energy decomposition for the Ras-Raf and Ras-RalGDS complexes. J Mol Biol. 2003, 330: 891-913. 10.1016/S0022-2836(03)00610-7
Sanchez R, Sali A: Comparative protein structure modeling. Introduction and practical examples with modeller. Methods Mol Biol. 2000, 143: 97-129.
Pettersen EF, Goddard TD, Huang CC, Couch GS, Greenblatt DM, Meng EC, Ferrin TE: UCSF Chimera–a visualization system for exploratory research and analysis. J Comput Chem. 2004, 25: 1605-1612. 10.1002/jcc.20084
We would like to thank Krista Winkelmann for immunohistochemical staining assistance and Mariana Bublik and Anna Zarling for technical assistance.
This study was supported by the Deutsche Forschungsgemeinschaft (SCHN477/12-1), to RSS and HGM.
The authors declare that they have no competing interests.
CEB performed all in vitro assays; she was involved in conception and design, development and methodology, acquisition of data, analysis and interpretation of data, writing and revision of the manuscript. VM performed the docking and clustering analysis and revised the manuscript. FBF performed the Ingenuity pathway analysis. SM performed docking and pathway clustering. TTR performed the IHC scoring, critically discussed and revised the manuscript. HGM was involved in data analysis and interpretation, writing and revision of the manuscript. RSS is the corresponding author as well as the study supervisor; she was involved in conception and design of this work, development and methodology, analysis and interpretation of data, writing and revision of the manuscript. All authors read and approved the final manuscript.
Electronic supplementary material
Additional file 2: Figure S1: TQ did not induce changes in the main kinases affected by IPA-3. Cells treated with 40 μM TQ were collected after 1, 3, 6 and 24 hours. Untreated cells were used as a control. 30 μg of proteins per well were used to assess the activity of AKT2, GSKα, GSKβ and p38. Data is presented as fold activity normalized to control activity. Each value is the mean ± SD of two independent experiments each done in duplicates. (TIFF 99 KB)
Additional file 3: Figure S2: TQ show low toxicity on normal intestinal cells. HCT116 (A) and HCEC (B) cells were treated with different TQ concentrations (0-100 μM) for 24 hours to define the IC50 value. Cell viability was assessed by crystal violet staining. Data are presented as percentage of control. Each value is the mean ± SD of three independent experiments done in quadruplicates. (TIFF 118 KB)
Additional file 4: Figure S3: IPA-3 and TQ bind to the autoinhibited dimer conformation of PAK1. A. Model of IPA-3 bound autoinhibited dimer conformation of PAK1 (PDB id:1kx5). Autoregulatory region (blue and green) and kinase domain (pink and blue) of PAK1 dimer are shown. IPA-3 is shown in red using a stick representation and the residues interacting with IPA-3 labeled in black. B. Model of IPA-3 and TQ bound autoinhibited dimer conformation of PAK1. TQ is shown in purple as a stick representation. The residues interacting with TQ are labeled in black. The rendering follows color legend used in A.(TIFF 5 MB)
Additional file 6: Figure S4: TQ induces down regulation of prosurvival pERK1/2 in vivo. Mouse xenograft experiment showing A. Immunohistochemical detection of pERK1/2 on tumors of control (left) and TQ treated (right) animals. B. Box plot analysis of the percentage of tumor cells expressing cytoplasmic pERK1/2 after TQ treatment. n represents the number of investigated mouse tissues in each group. *16 corresponds to sample 16 being an outlayer. (TIFF 1 MB)
Additional file 7: Figure S5: TQ induces the phosphorylation of pPAK1Thr212 before the phosphorylation of pERK1/2. Cells treated with 60 μM TQ were collected after 10, 15, 30, 45, 60 and 90 min. untreated cells were used as control. 30 μg of proteins were immunoblotted against pPAKThr212 and pERK1/2. GAPDH was used as loading control. Data shown are representative of two independent experiments. (TIFF 154 KB)
Additional file 8: Table S3: Hydrogen bonds formed at the PAK1 catalytic site in the absence and presence of TQ. (DOCX 18 KB)
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
About this article
Cite this article
El-Baba, C., Mahadevan, V., Fahlbusch, F.B. et al. Thymoquinone-induced conformational changes of PAK1 interrupt prosurvival MEK-ERK signaling in colorectal cancer. Mol Cancer 13, 201 (2014). https://doi.org/10.1186/1476-4598-13-201