1. Introduction
The integrated stress response (ISR) pathway allows cells to adapt to stressful situations such as amino acid depletion, endoplasmic reticulum (ER) stress, viral infection and heme depletion [
1]. The ISR is comprised of four different protein kinases that respond to the different stressors and have as their only known target the eukaryotic translation initiation factor eIF2α. They are general control nonderepressible 2 (GCN2), which respond to amino acid depletion by using uncharged tRNAs as a sensor of amino acid imbalances; protein kinase R-like endoplasmic reticulum kinase (PERK), which senses ER-stress [
2]; protein kinase R (PKR), which detects double-stranded RNA [
3]; and heme-regulated inhibitor, which responds to stresses associated with hemoglobin production [
4]. In the context of this study, we focused on GCN2. Phosphorylation of eIF2α reduces overall protein biosynthesis and at the same time allows the translation of specific mRNAs with upstream open-reading frames, in particular ATF4 [
5]. ATF4 is a general activating transcription factor, which promotes the synthesis of genes involved in amino acid biosynthesis, transport, autophagy and antioxidant responses [
6,
7]. These mechanisms can restore amino acid homeostasis by improving the provision of amino acids and enhancing the recycling of proteins and cellular materials. In cancer cells rapid growth is thought to cause amino acid depletion and ER stress [
8], which has drawn attention to the role of GCN2 and PERK in cancerous cell growth [
9,
10]. In vitro, however, cancer cell lines are often unstressed and require the inhibition of other components of amino acid homeostasis to elicit the ISR [11-13]. Elements of amino acid homeostasis are frequently explored for cancer therapy. Asparaginase, for instance, is used to treat acute lymphocytic leukemia [
14]. Inhibition of amino acid transporter LAT1 [
15] is explored to treat a variety of solid tumors. Bortezomib is a proteasome inhibitor that affects protein recycling [
16], while rapamycin-derived mTORC1 inhibitors were explored in a large variety of clinical trials [
17].
Activation of GCN2 initiates a rescue mechanism for cancer cells to deal with cellular stresses elicited by other anticancer drugs, which can thwart the expected effects of targeting amino acid homeostasis [
13,
16,
17]. Thus pairing a drug, which targets amino acid homeostasis with a GCN2 inhibitor could improve its efficacy. Anticancer drug pairs can be beneficial over monotherapies as they can suppress the emergence of drug resistance and achieve sustained remission. Tumours frequently consist of a heterogenic population of malignant cells. The diversity of oncogenic profiles within a malignancy is driven by the genomic instability that is intrinsic to cancer cells and by different selection pressures arising from the immediate microenvironment of various subpopulations [
18]. One way to target diverse populations of cells is combination therapy. A substantial combination screening effort by the US National Cancer Institute, involving 5000 FDA-approved anticancer drugs in 60 human cancer cell lines, discovered two promising drug pairings: clofarabine with bortezomib, and nilotinib with paclitaxel. Both combinations are in Phase I clinical trials (NCT02211755 and NCT02379416) [
19].
Several GCN2 inhibitors have been developed and characterised for kinase selectivity and pharmacological properties. Nakamura et al. [
14] reported the discovery of GCN2iB with an IC
50 at 2.4 nM. GCN2iB also inhibits MAP2K5, STK10 and ZAK with slightly lower potencies but has a pharmacokinetic profile which makes it suitable for
in vivo studies. GZD824 (olverembatinib) is a broad-spectrum Bcr-Abl inhibitor that is orally bioavailable and recently repurposed as a GCN2 inhibitor; however, it is highly promiscuous in kinase assays [
20]. Lastly, triazolo[4,5-d]pyrimidines (TAP) constitute an entire class of GCN2 inhibitors [
21,
22]. Among these, compound 20 (TAP20) demonstrated superior potency and selectivity. In a cell-free assay, TAP20 achieved an IC
50 of 17 nM against GCN2. It was similarly potent against GSK3α/β, and, to a lesser extent, CDK9/cyclinD1.
In this study we explored synergistic relationships between 25 experimental and approved drugs and the GCN2 inhibitor TAP20 using several cell lines.
2. Materials and Methods
Cell lines: All cell lines were procured from the American Type Culture Collection (ATCC) and maintained in DMEM/F-12 with 10% heat-inactivated FBS, 2 mM glutamine and 0.1% penicillin/streptomycin. Cells were stored in a humidified atmosphere with 5% CO2 in an incubator set to 37 °C in 25 cm2 or 75 cm2 culture flasks. Once confluent, cells were trypsinised using 0.25% Trypsin-EDTA (Gibco). Experiments were conducted using cells passaged 5–20 times.
Chemicals and antibodies: Drugs and chemicals were purchased from several sources. TAP20 was synthesised and provided by Merck KGaA (Darmstadt, Germany). YH16899 was purchased from AOBIOUS (Gloucester, MA, USA); thapsigargin from Thermo Fisher Scientific (Waltham, MA, USA), Bay876 from Sigma-Aldrich (St Louis, MO, USA) and all other compounds from MedChemExpress (Monmouth Junction, NJ, USA). Bumetanide and 5-fluorouracil stocks were dissolved in water and all other compounds were dissolved in DMSO. Cetuximab was purchased from MedChemExpress and all other antibodies were from Cell Signaling Technology (Danvers, MA, USA).
Proliferation assay: Cell growth was monitored using the IncuCyte system (Essen Biosciences). To increase throughput, cell confluence was measured only at the end of the assay. Briefly, cells were seeded in 96-well flat-bottom plates at 3000 cells/well in growth medium and allowed to attach overnight. The next day, the growth medium from each well was replaced with 300 µL of growth medium with or without treatment. Cells were stored in the cell incubator (5% CO2 humidified atmosphere; 37°C) for a predefined period before being scanned by the IncuCyte system. This period was determined based on the shortest time required for the untreated control cells to reach >90% confluence. This ranged from four to ten days and depended on the cell line. Confluence was recorded as the mean from two images taken for each well. Peripheral wells were excluded from the analysis.
IC50 curves were plotted from data obtained in this assay using the logistic growth equation:
Where A1 is the maximal confluence, A2 is the minimal confluence, x is the inhibitor’s concentration, x0 is the inhibitor’s concentration at half maximal effect and p is the power value.
To evaluate synergism among drug pairings, the coefficient of drug interaction [
23] (CDI) was calculated using the following equation:
Where AB is the confluence of cells treated with the drug combination and A and B are the confluence values of cells exposed to each drug alone. CDI scores of 1 indicate an additive effect, those above 1 indicate antagonism and those below 1 indicate synergism.
Synergy analysis: Synergy diagrams were generated from data produced in the assay described above. A matrix of seven concentrations per drug was arrayed using MDA-MB-231 and HPAF-II cells seeded across three 96-well plates for each experiment. Drug concentrations were incremented by values of 1, 2 and 4 per order of magnitude. As above, peripheral wells were excluded from analysis and each plate contained untreated controls to which the confluence values of treated cells were normalised. Normalised data were then converted to synergy scores using SynergyFinder [
23] and the Loewe equation.
Western blotting: Intracellular proteins were prepared for SDS-PAGE by washing cells grown almost to full confluence in 35 mm dishes with PBS (pH 7.4) and adding 100–150 μL of RIPA lysis and extraction buffer (Thermo Fisher Scientific; Waltham, MA, USA) containing cOmplete Mini Protease Inhibitor Cocktail (Roche; Basel, Switzerland). The lysate was collected using a cell scraper and transferred to a microcentrifuge tube. After 5 min of mixing and 5 min of centrifugation at top speed, 10 μL of the supernatant was sampled and diluted at least ten-fold in water for the Bradford assay (reagent obtained from Sigma-Aldrich; St Louis, MO, USA) to measure protein concentration and to normalise the protein mass added to the wells of an SDS-PAGE gel. Lysates were then mixed with 4×Bolt LDS and 10×Bolt Sample Reducing Agent (Life Technologies) and water up to 50 µL. A volume of the mixture, equivalent to 20–30 µg protein, was then loaded into the wells of a precast NuPAGE 4-12% Bis-Tris 10- or 15-well gel fitted to a Mini Gel Tank (Life Technologies) filled with MOPS-SDS running buffer (Thermo Fisher; Waltham, MA, USA). Sample lanes were flanked with wells loaded with either SeeBlue Plus2 or Novex Sharp pre-stained protein standards (Invitrogen; Waltham, MA, USA). Proteins were separated for approximately 1 h at 150 V, after which proteins were transferred onto a 0.45-micron Amersham Protran nitrocellulose membrane (GE Healthcare; Chicago, IL, USA) in an ice-cold wet transfer system (Bio-Rad; Hercules, CA, USA) containing Towbin buffer (192 mM NaCl, 25 mM Tris base, 20% methanol; pH 8.3) for 2 h at 100 V. The blot was then incubated in 5 % (w/v) skim milk powder in PBS-Tween (0.1% Tween-20; pH 7.4) for 1 h at room temperature or overnight at 4°C. In between blocking and antibody solutions, blots were always washed thrice with PBS-Tween for 5 min at a time. Primary antibodies were always incubated with blots overnight, while secondary antibodies were incubated for 2–6 h. To visualise the proteins, either the SuperSignal West Pico PLUS or the SuperSignal West Femto Maximum Sensitivity chemiluminescent substrates (Thermo Scientific; Waltham, MA, USA) were added dropwise onto the blots. Either the ChemiDoc MP (Bio-Rad) or Fusion Solo 2 (Vilber) imaging systems were used to visualise chemiluminescence.
All blots were probed first for the protein of interest and lastly for a housekeeping protein. β-Actin was probed as loading control. To remove protein-bound antibodies from the previous immunodetection, blots were incubated at 65°C in stripping buffer (2% w/v SDS; 62.5 mM Tris; and 0.7% v/v 2-mercaptoethanol; pH 6.8) for 10 min with agitation. The blots were then washed with PBS-Tween and incubated for at least 2 h with 5% skim milk in PBS-Tween.
Matrix-dependent invasion assay: To investigate drug-pair effects on matrix-dependent invasion, 143B osteosarcoma cells were used in a Cultrex 3D spheroid invasion system (R&D Systems; Minneapolis, MN, USA) as they readily form the necessary spheroids unlike the six cell lines used in the synergy screen. Initially, 143B were seeded in 1 × spheroid formation medium (ECM) (3000 cells per well in 50 μl) into low-adhesion, round-bottom 96-well plates to induce formation of spheroids. The plates were centrifuged at 300 × g for 3 min to position cells at the bottom of the well. Working on ice, 50 μl invasion matrix was added to the spheroids after 72 h. Subsequently, plates were centrifuged at 300 × g at 4°C to position the spheroid at the bottom of their wells. Gel formation was initiated by incubation at 37°C for 1 h followed by addition of 100 μl of pre-warmed culture medium. Plates were then incubated at 37°C/5% CO2 in a humidified incubator for 3-6 days, and images were taken with a Leica M205FA stereomicroscope. Invasion perimeter area was measured using the Leica Application Suite, version 4.12.
3. Results
We reasoned that current cancer drugs and experimental anticancer compounds might cause stress to cell lines, which was counteracted by GCN2 activation. Blocking GCN2 with TAP20 would reduce the cellular response flexibility, thereby reducing cellular fitness.
We selected 24 experimental compounds with well-known targets (
Table 1). The reported target IC
50 values are all in the nanomolar range typical for modern pharmaceutical compounds. To establish a baseline for a potential synergistic effect of TAP20, we determined the IC
50 of 25 candidate drugs (
Table 2) for cell growth using a range of cell lines that have a variety of oncogenic signatures [
24]. Growth IC
50-values can be higher than the target IC
50-values due to membrane permeability, drug-resistance mechanisms and cellular reliance on the function of the target.
To determine synergy with TAP20, candidate drugs were adjusted to a concentration where inhibition was just observable (i.e. < IC
50,
Table 3). In this range, synergy is readily observed as growth rates are especially sensitive to reduced fitness. In general, we did not employ concentrations >10 µM due to insolubility or increased non-specific drug interactions. The dose-response curves for all cell lines can be found in supplementary Fig S1.
As a last preparation, we determined the response of cell lines to TAP20. TAP20 alone typically does not inhibit growth when cell lines are cultured in replete growth media because they experience little stress, which would initiate activation of GCN2. To quantify the action of TAP20 on GCN2 we starved the cell lines for 24 h in the presence of different concentrations of TAP20 (
Figure 1). In our hands, the expression of ATF4 was the most reliable indicator of GCN2 activation and we used western blotting to determine ATF4 expression at different doses of TAP20. Except for Panc 02.03 cells, ATF4 expression was suppressed by TAP20 concentrations > 1 µM. The IC
50-values were MDA-MB-231 = 0.54 µM, HPAFII = 0.56 µM, SKOV3 = 0.29 µM. For all synergy tests, TAP20 was used at 3 µM to ensure that the ISR could not be activated via GCN2.
In the synergy experiments, control cells were grown to confluency and the growth compared to that of the treated cells. A coefficient of drug interaction was calculated and used to compare the experiments (
Table 4). A score <1 (green in
Table 4) indicates a synergistic relationship, a score of 1 an additive effect (yellow) and scores >1 indicate an antagonistic relationship (red). The unprocessed growth data for each combination are shown in supplementary
Figure S2. Several trends were observed. First, instances of drug synergy were not consistently observed across all cell lines. We observed a variety of synergistic relations in MDA-MB-231, HPAFII and SKOV3 cells, while these were mostly additive in MDA-MB-468 and Panc02.03 cells. We observed only a limited number of synergistic relations in Ovcar3 cells. Second, rapamycin was antagonistic to TAP20. This was expected because inhibition of mTORC1 is known to activate GCN2 and enhance phosphorylation of eIF2α [
25].
Numerous synergistic relationships between TAP20 and the candidate drugs could be identified and grouped into several cellular processes. Some of the relationships can be readily rationalized. V-9302 for instance is an inhibitor of amino acid transporters LAT1 and SNAT2 [
26]. Both transporters play a key role in the maintenance of cytosolic amino acid levels. Inhibition of LAT1 has been shown in the past to activate GCN2 [
27]. Similarly, we expected inhibitors of amino acid homeostasis to synergize with TAP20. Halofuginone and YH16889 interfere with prolyl-tRNA synthetase and lysyl-tRNA synthetase, respectively, and are therefore mechanistically close to the tRNA-binding function of GCN2. CB-5083 inhibits p97 ATPase, a protein which is involved in unfolding of ubiquitinated proteins and membrane proteins [
28]. The unfolded proteins are then transferred to the proteasome for degradation, a process that replenishes cellular amino acid pools. Another inhibitor of proteostasis is bortezomib which blocks the proteasome itself, although synergy was observed in different cell lines compared to CB-5083. Thapsigargin, an inhibitor of the SERCA ATPase, induces ER stress which is an independent activator of GCN2 through the protein kinase PERK. However, thapsigargin only synergised with TAP20 in HPAFII and in SKOV3 cells.
A direct link between the MEK-ERK MAPK signalling pathway and the ISR has been demonstrated previously in HepG2 cells [
29]. Consistently, we found a strong synergy between TAP20 and ERK inhibitor SCH772984. The MEK inhibitor selumetinib, which inhibits the MEK-ERK pathway upstream of ERK had a similar effect, although not in HPAFII cells.
The most striking synergy we identified was between elements of cell-cycle regulation and TAP20. This was exemplified by inhibitors of cyclin-dependent kinases, such as flavopiridol, seliciclib and THZ-1 [
30] and danusertib an inhibitor of aurora kinases [
31]. While flavopiridol and seliciclib are non-specific CDK inhibitors, THZ-1 is a more selective CDK-7 inhibitor.
To further explore the scope of these synergistic relationships we performed growth experiments exploring a matrix of different drug concentrations and tested the synergistic effect on matrix-dependent cell invasion. Due to low invasiveness of the cell lines used in the growth assays, we performed the cell invasion assay with 143B osteosarcoma cells, while the synergy analysis was made using MDA-MB-231 cells and HPAF2 cells, which showed the strongest responses in this study.
Combining TAP20 with p97-ATPase inhibitor CB5083 showed a synergistic effect as shown by growth inhibition and calculation of synergy scores using Synergyfinder [
23]. In HPAFII cells the optimal effect was achieved at 1 µM CB5083 and >1 µM of TAP20 (
Figure 2).
Although TAP20 alone had little effect on growth in all cell lines, it dramatically reduced matrix-dependent invasion of 143B cells at 1–3 µM, but to evaluate synergistic effects, a lower concentration, 0.75 µM, was used (
Figure 2). In agreement with the growth experiments, spheroid size was reduced in the presence of both inhibitors (compare the growth area of the control with that of TAP20/CB5083).
The highly synergistic relationship between the CDK inhibitor flavopiridol and TAP20 was confirmed in the synergy analysis yielding scores of up to 40 in MDA-MB-231 and HPAFII cells (
Figure 3). Optimal synergy was achieved with relatively low concentrations of both compounds (MDA-MB-231 TAP20: 0.2 µM; flavopiridol: 0.04 µM (
Figure 3 a, b); HPAFII TAP20 1 µM, flavopiridol 0.03 µM (
Figure 3 c, d)). Synergy disappeared at higher drug concentrations in MDA-MB-231 cells (
Figure 3 b). As shown in
Figure 2, TAP20 (0.75 µM) had a significant effect on matrix-dependent invasion of 143B cells alone, which was similar to the effect of 50 nM flavopiridol (
Figure 3 e). This suggests that invasion generates cellular stress, which the cells ameliorate by activation of GCN2. Combination of both drugs significantly potentiated the effect. In the presence of both drugs hardly any cell protrusions from the original spheroid were observed (
Figure 3 f).
This result was confirmed using the non-specific CDK inhibitor seliciclib (
Figure 4). Inhibition of growth of MDA-MB-231 cells and synergy markedly increased when both drugs were applied at higher concentrations (
Figure 4 a, b). Synergy was confined to a narrow range of concentrations. At 10 µM, seliciclib significantly reduced matrix-dependent invasion on its own, but in combination with TAP20 (1 µM), no outgrowth was observed compared to tumor spheroids in the absence of matrix (
Figure 4 c, d).
Recently, CDK7 has been implicated in the ISR [
32]. To investigate whether this isoform might be underlying the responses observed for flavopiridol and seliciclib, we used THZ-1, a specific CDK7 inhibitor. While THZ-1 showed a strong synergistic action in MDA-MB-231 (
Figure 5 a, b) and HPAFII cells (
Figure 5 c, d) it was ineffective alone or in combination with TAP20 to suppress growth or invasion of 143B cells (
Figure 5 e, f).
We also analysed the MEK and ERK inhibitors selumetinib (
Figure 6) and SCH772984 (
Figure 7) in detail. Both compounds inhibited growth at low concentrations (
Figure 6a and
Figure 7a) and showed high synergy scores at low concentrations (SCH772984 0.4 µM, selumetinib 2 µM). However, synergy was generally low for SCH772984, which was a very effective drug independent of TAP20. At 100 nM selumetinib not only reduced invasion but also the size of the tumor spheroids (
Figure 6 c, d), suggesting that cells were dying during the treatment. A similar pattern was observed with SCH772984 (
Figure 7). Spheroids were markedly smaller than control spheroids grown in the absence of matrix (
Figure 7 c, d)
In summary, we have identified key cellular processes inhibition of which causes nutrient stress and requires metabolic flexibility, which is afforded by activation of GCN2. This can be abrogated by concurrent inhibition of GCN2. Combination of drugs is an important strategy to overcome drug resistance and to inhibit processes that require significant pathway flexibility such as metastasis.
4. Discussion
While cancer cell lines experience little nutritional stress in replete growth media, cancer cells in the microenvironment of the tumor are more likely to experience nutrient stress [
33,
34]. Evidence for this is the close association of cancer stem cells with blood vessels [
35]. the ISR is one of the cellular programs to deal with stress and to elicit responses that optimize cellular fitness [
1,
36]. Recycling of amino acids through lysosomal and proteasomal protein degradation is a key element of amino acid homeostasis [
37]. Consistently, we found synergy with the proteasomal inhibitor bortezomib, the p97-ATPase inhibitor CB-5083 and the amino acid transport inhibitor V-9302. We did not detect synergy with the PERK inhibitor AMG PERK44, suggesting that our chosen cancer cell lines do not experience significant ER stress. It is tempting to speculate that AMG PERK44 may synergize well with CB-5083 or thapsigargin, both of which cause ER stress.
A synergistic relationship between TAP20 and inhibitors of cyclin-dependent kinases may appear surprising at first sight. However, before undergoing division, cells must ensure that sufficient amino acids are available for abundant protein synthesis to duplicate cellular content. Accordingly, amino acid deprivation stops cell cycle progression by stabilizing the mRNA of endogenous cyclin-dependent kinase inhibitors p22 and p27 [
29].
CDK-7 is an unusual cyclin-dependent kinase, as it appears to activate other CDKs instead of directly regulating steps of the cell cycle. In addition, it has a role as a transcriptional activator [
38]. It allows RNA polymerase II to initiate transcription and facilitates capping of mRNAs [
39]. A direct link between CDK-7 and expression of the amino acid transporter SNAT2 was identified by Stretton et al. [
32]. Upregulation of SNAT2 mRNA in response to amino acid limitation depended on CDK-7 activity. While RNA transcription and execution of cellular programs are generally associated, activation of GCN2 silences CAP-dependent transcription allowing mRNAs with uORFs, such as SNAT2, to be translated by a CAP-independent mechanism [
40]. Sirtuin 6 was recently identified as a link between CDK-7 and the ISR [
41]. Sirtuin 6 binds to ATF4, stabilizing the latter and preventing its degradation. Low expression levels of sirtuin 6 causes cells to become sensitive to CDK7 inhibition. CDK-7 was found to enhance transcription of more than 2000 genes, including MYC [
41]. Thus, inhibition of CDK-7 will cause cellular stress, which elicits the ISR.
Related to cell cycle progression, Piecyk et al., found synergy between oxaliplatin and GCN2 inhibition in colon cancer tumoroids [
42]. This was based on GCN2 sustaining ribosomal RNA transcription in nutrient-rich conditions, which synergised with inhibition of RNA polymerase I.
Notably, inhibitors of proteostasis, such as CB-5083 and bortezomib have a significant influence on cell-cycle progression as rapid breakdown of cell-cycle proteins is required at certain stages of the cell cycle [
28,
43]. This may contribute to the synergy observed in this study.
The observed synergy with the MEK-ERK pathway is harder to explain. In the absence of stress, growth-promoting signals lead to Akt-dependent and or Erk-dependent phosphorylation of the tuberous sclerosis-2 (TSC2), which in turn causes the release of the TSC complex from the lysosome resulting in mTORC1 activation [
44]. Thus, we might have expected an antagonistic relationship similar to rapamycin. However, a direct link between the MEK-ERK MAPK signalling pathway and the ISR has been demonstrated previously in HepG2 cells [
29]. In these cells, an ERK inhibitor prevented increased transcription of ATF4 and of the amino acid response target gene SNAT2 [
45]. The study further showed that phosphorylation of eIF2α was also blocked by inhibition of ERK without ERK directly phosphorylating eIF2α. The interaction is interdependent, as lack of GCN2 prevents phosphorylation of ERK in response to amino acid limitation. Moreover, there is also a direct link between the activity of amino acid transporter SNAT2 and ERK1/2 activity [
46]. As a result, blocking the MEK-ERK axis could increase nutrient stress, which can be aggravated by simultaneously blocking GCN2.
Our results highlight the mutual relationship between nutrient stress and the ability of cancer cells to adaptively respond to stressful situations using the GCN2-ATF4 pathway. While blocking the GCN2-ATF4 axis alone may have limited efficacy to reduce cancer cell growth, its blockade in conjunction with stress-inducing drug regimens could be a promising strategy to reduce cancer cell growth and invasion.