Mostafa M Elbadawi1,2, Wagdy M Eldehna2, Amer Ali Abd El-Hafeez3,4, Warda R Somaa5, Amgad Albohy6, Sara T Al-Rashood7, Keli K Agama8, Eslam B Elkaeed9, Pradipta Ghosh4,10,11,12, Yves Pommier8, Manabu Abe1. 1. Department of Chemistry, Graduate School of Science, Hiroshima University, Hiroshima, Japan. 2. Department of Pharmaceutical Chemistry, Faculty of Pharmacy, Kafrelsheikh University, Kafrelsheikh, Egypt. 3. Pharmacology and Experimental Oncology Unit, Cancer Biology Department, National Cancer Institute, Cairo University, Cairo, Egypt. 4. Department of Cellular and Molecular Medicine, University of California San Diego, La Jolla, CA, USA. 5. Faculty of Pharmacy, Kafrelsheikh University, Kafrelsheikh, Egypt. 6. Department of Pharmaceutical Chemistry, Faculty of Pharmacy, The British University in Egypt (BUE), Cairo, Egypt. 7. Department of Pharmaceutical Chemistry, College of Pharmacy, King Saud University, Riyadh, Saudi Arabia. 8. Developmental Therapeutics Branch, Laboratory of Molecular Pharmacology, Center for Cancer Research, National Cancer Institute, NIH, Bethesda, MD, USA. 9. Department of Pharmaceutical Sciences, College of Pharmacy, AlMaarefa University, Riyadh, Saudi Arabia. 10. Department of Medicine, University of California San Diego, La Jolla, CA, USA. 11. Moores Comprehensive Cancer Center, University of California San Diego, La Jolla, CA, USA. 12. Veterans Affairs Medical Center, La Jolla, CA, USA.
Abstract
In this study, different assortments of 2-arylquinolines and 2,6-diarylquinolines have been developed. Recently, we have developed a new series of 6,7-dimethoxy-4-alkoxy-2-arylquinolines as Topoisomerase I (TOP1) inhibitors with potent anticancer activity. Utilising the SAR outputs from this study, we tried to enhance anticancer and TOP1 inhibitory activities. Though target quinolines demonstrated potent antiproliferative effect, specifically against colorectal cancer DLD-1 and HCT-116, they showed weak TOP1 inhibition which may be attributable to their non-coplanarity. Thereafter, screening against kinase panel revealed their dual inhibitory activity against EGFR and FAK. Quinolines 6f, 6h, 6i, and 20f were the most potent EGFR inhibitors (IC50s = 25.39, 20.15, 22.36, and 24.81 nM, respectively). Meanwhile, quinolines 6f, 6h, 6i, 16d, and 20f exerted the best FAK inhibition (IC50s = 22.68, 14.25, 18.36, 17.36, and 15.36 nM, respectively). Finally, molecular modelling was employed to justify the promising EGFR/FAK inhibition. The study outcomes afforded the first reported quinolines with potent EGFR/FAK dual inhibition.
In this study, different assortments of 2-arylquinolines and 2,6-diarylquinolines have been developed. Recently, we have developed a new series of 6,7-dimethoxy-4-alkoxy-2-arylquinolines as Topoisomerase I (TOP1) inhibitors with potent anticancer activity. Utilising the SAR outputs from this study, we tried to enhance anticancer and TOP1 inhibitory activities. Though target quinolines demonstrated potent antiproliferative effect, specifically against colorectal cancer DLD-1 and HCT-116, they showed weak TOP1 inhibition which may be attributable to their non-coplanarity. Thereafter, screening against kinase panel revealed their dual inhibitory activity against EGFR and FAK. Quinolines 6f, 6h, 6i, and 20f were the most potent EGFR inhibitors (IC50s = 25.39, 20.15, 22.36, and 24.81 nM, respectively). Meanwhile, quinolines 6f, 6h, 6i, 16d, and 20f exerted the best FAK inhibition (IC50s = 22.68, 14.25, 18.36, 17.36, and 15.36 nM, respectively). Finally, molecular modelling was employed to justify the promising EGFR/FAK inhibition. The study outcomes afforded the first reported quinolines with potent EGFR/FAK dual inhibition.
Cancer is a major health obstacle in the world threating the life of millions of people annually,. The universal burden of cancer has been expected to rise to 21.6 million in 2023 compared to 14.1 million in 2012 and was predicted to increase to 28.4 million in 2040 with 47% increment relative to 2020. In 2020, 19.3 million new cancer cases have been diagnosed and 10 million cancer patients passed away. As established by WHO in 2019, cancer was estimated to be the first or second dominant cause of death for the ages <70 years in 112 countries, while it was projected to be the third or fourth death cause in 23 countries. In general, the incidence and mortality of cancer are growing rapidly worldwide,. Subsequently, enormous attempts have been implemented to develop potent anticancer drugs through investigations of diverse scaffolds against numerous potential chemotherapeutic targets.Epidermal growth factor receptor (EGFR) is a member of tyrosine kinase family in which the endogenous ligand binds to the extracellular domain leading to conformational changes and dimerisation of EGFR resulting in its activation which subsequently stimulates its intrinsic intracellular protein-tyrosine kinase activity. EGFR is over-expressed in many solid tumours and is related to cancer cell proliferation, angiogenesis and metastasis, so it has a critical role in cancer growth. Therefore, EGFR has been validated as an efficient target for anticancer drug discovery. In the last two decades, different 4-anilinoquinazoline-based EFGR inhibitors, such as Gefitinib, Erlotinib, Afatinib, and Dacomitinib (Figure 1), have been FDA-approved for clinical use in treatment of non-small cell lung cancer,.
Figure 1.
Some reported EGFR and FAK inhibitors.
Some reported EGFR and FAK inhibitors.Focal adhesion kinase (FAK) is a cytoplasmic non-receptor tyrosine kinase involved in signal transductions from cell adhesions to regulate different biological cell functions including survival and cell migration,. Also, it is activated and overexpressed in diverse cancer types controlling cancer proliferation, survival and metastasis. Thus, FAK has been identified as a promising druggable target for targeted cancer therapy. Currently, several FAK inhibitors, such as 2,4-diaminopyridine derivative GSK2256098 and 2,4-diaminopyrimidine derivative Defactinib (Figure 1), are currently being evaluated in clinical trials for cancer treatment, in addition to the 2,4-diaminopyrimidine derivative TAE-226 (Figure 1) which displayed potent antitumor impact in different cancer types in vivo and in vitro and usually used as a reference drug,,. Noteworthy, it was established that the most affected colorectal cancer expressed high levels of EGFR and FAK that particularly correlated with tumour angiogenesis, cancer aggressiveness and poor prognosis,.Thus, dual EGFR/FAK inhibition mechanism is an efficient strategy to fight cancer that could be attributed to a non-overlapping downstream signalling/inhibition,. For example, the kinase inhibitor APG-2449 (Figure 1) was reported to improve the antitumor effect of Ibrutinib via EGFR/FAK inhibition mechanism in oesophageal squamous cell carcinoma. Also, combined EGFR/FAK inhibition caused higher radiosensitization than either approach alone. Interestingly, few studies have succeeded to develop dual EGFR/FAK small molecule inhibitors. In 2020, Ai et al. has exploited a fragment-based drug design approach to identify novel series of 2,4-diaminopyrimidines as potent dual EGFR/FAK inhibitors with good in vitro and in vivo antitumor effects.Quinoline is an outstanding planar heterocyclic motif playing a distinctive role in anticancer drug discovery. So far, assortments of quinoline-based small molecules have been developed and investigated against numerous biological targets for cancer treatment displaying exquisite outcomes. It is worth stressing that plenty of quinoline derivatives provoked their anticancer impact through different mechanisms of action, such as inhibition of DNA repair, tubulin polymerisation, and inhibition of various enzymes implicated in critical cancer cell proliferation prominently kinases enzymes (EGFR, VEGFR, pim-1 kinase, c-Met factor, and PI3K) which stood out as one of the most significant targets implemented in cancer therapy due to their functions in cellular signal transduction. Of special interest, Pelitinib (EKB-569, Figure 2) is a 4-anilinoquinoline derivative which is a potent irreversible inhibitor of EGFR in the clinical trials as an anticancer candidate. In addition, several 4-aminoquinoline derivatives, such as compounds I–III (Figure 2), were reported as promising EGFR inhibitors endowed with effective anticancer activities. Accordingly, quinoline stands out as a significant privileged scaffold in anticancer drug discovery to develop many efficient kinases inhibitors,.
Figure 2.
Some reported quinolines with potential anticancer activity.
Some reported quinolines with potential anticancer activity.Recently, we have developed a new series of 6,7-dimethoxy-4-alkoxy-2-arylquinolines as potential Topoisomerase I (TOP1) inhibitors. The TOP1-mediated DNA cleavage assay was utilised to assess the ability of the reported compounds to stabilise TOP1-DNA cleavage complexes (TOP1ccs). The assay outcomes revealed a moderate TOP1 inhibitory activity of compounds IVa,b (Figure 2). Interestingly, the developed quinolines showed outstanding anti-proliferative profile upon evaluation at the Developmental Therapeutics Program (DTP) of the NCI-USA. Noteworthy, the weightiness of incorporation of p-substituted phenyl at C-2, as well as propyl linker at C-4 of the quinoline scaffold, was highlighted by the SAR study.As a continuation for our previous study novel five sets of 4-propoxy-2-arylquinolines (6a–o, 8a,b, 10a,b, 12a–d, 16a–d), and 4-propoxy-2,6-diarylquinolines (20a–f) are herein designed and synthesised, exploiting the deduced SARs from the previous study, with the aim to afford more potent anticancer TOP1 inhibitors (Figure 2). Different structural modification strategies were adopted seeking to enhance both anticancer and TOP1 inhibitory activities of the lead compounds IVa,b (Figure 3).
Figure 3.
Different structural modification strategies adopted for design of target quinoline derivatives in this study.
Different structural modification strategies adopted for design of target quinoline derivatives in this study.First, diverse secondary alicyclic and aromatic amines, in addition to different primary amines were appended to the propyl linker to illuminate the influence of these moieties on the activity. Also, we introduced the electron donating methyl group as well as electron withdrawing Cl and CF3 for more elucidation of the electronic impact of p-substituent of the 2-phenyl. Likewise, the heterocycles: 2-furyl and 2-thienyl were appended to para position of the C2-phenyl moiety to explore their electronic and size impact on the desired activity. Afterward, the 6,7-dimethoxy groups were removed in some synthesised analogs to confirm their significance. Moreover, the ring system was extended via fusion of the quinoline motif with 1,3-dioxolane in attempt to enhance the planar structure requisite for DNA intercalation which may potentiate both anticancer and TOP1 poisoning effects. Finally, a structural extension approach was utilised via grafting the HBA-bearing 2-furyl and 4-methoxyphenyl moieties at C6 of quinoline scaffold, hoping to enhance the hydrophobic interactions (Figure 3). The designed 4-propoxy-2-arylquinolines were prepared employing different synthetic procedures, and then investigated for their anticancer and TOP1 inhibitory activities.Although the target quinolines demonstrated potent antiproliferative effect against different cancer cell lines, they showed no or weak TOP1 poisoning influence. Accordingly, the promising anticancer activity prompted us to search for the plausible molecular mechanism for herein reported quinolines.The diverse well-reported kinase inhibitory activities of quinoline-based small molecules, as mentioned above, motivated us to explore the potential inhibitory activity of target 4-propoxy-2-arylquinolines (6a–o, 8a,b, 10a,b, 12a–d, and 16a–d) and 4-propoxy-2,6-diarylquinolines (20a–f) against various kinases (EGFR, FAK, FRK, IGF-1R, BTK, c-Src, VEGFR-1, VEGFR-2, HER-2). Strikingly, the investigated quinolines exhibited promising dual inhibitory effect towards EGFR and FAK kinases. Moreover, the apoptotic impact of the most potent anti-proliferative agents in this study was investigated on DLD-1 cells exploiting AV/PI dual staining assay. Finally, in silico molecular modelling techniques, including docking and molecular dynamics studies, were exploited to justify and support results obtained from the biological evaluations.
Results and discussion
Chemistry
The synthetic routes adopted for the synthesis of the target 4-propoxy-2-arylquinolines (6a–o, 8a,b, 10a,b, 12a–d, 16a–d) and 4-propoxy-2,6-diarylquinolines 20a–f are illustrated in Schemes 1–3. Regarding Scheme 1, the benzamides 3a–c were prepared by reacting 4,5-dimethoxy-2′-aminoacetophenone 1 with p-substituted benzoyl derivatives 2a–c in dry THF and Et3N. The latter were cyclized in refluxing dioxane and NaOH to afford the corresponding quinolones 4a–c in excellent yields which then were subjected to O-alkylation with 1-bromo-3-chloropropane using our previously confirmed procedure to yield the respective 4-propoxy key intermediates 5a–c. Finally, these key intermediates 5a–c were converted to the target 4-prpoxy-2-arylquinolines 6a–o in good to excellent yields (70–80%) through nucleophilic substitution with the appropriate amine in dry DMF and potassium carbonate anhydrous using catalytic amount of potassium iodide at 90 °C.
Scheme 1.
Synthesis of target 4-propoxy-2-arylquinolines 6a–o; Reagents and conditions: (i) Et3N, THF, 0 °C then rt, overnight; (ii) dry dioxane, NaOH, reflux under N2, 110 °C, 4 h; (iii) KI, KOH, 1-Bromo-3-chloropropane, dry DMF, rt, 24 h; (iv) KI, K2CO3 anhydrous, appropriate amine, dry DMF, reflux, 90 °C, 12 h.
Synthesis of target 4-propoxy-2-arylquinolines 6a–o; Reagents and conditions: (i) Et3N, THF, 0 °C then rt, overnight; (ii) dry dioxane, NaOH, reflux under N2, 110 °C, 4 h; (iii) KI, KOH, 1-Bromo-3-chloropropane, dry DMF, rt, 24 h; (iv) KI, K2CO3 anhydrous, appropriate amine, dry DMF, reflux, 90 °C, 12 h.In Scheme 2(A), the bromo analog of 4-propoxy-N-methylpiperazine-2-arylquinoline 7 was obtained based on the same pathway for 6a–o using p-bromobenzoyl chloride. Then, this bromo analog 7 has been transformed into the target heterocyclic derivatives 8a,b in good yields (70–72%) under Suzuki cross coupling condition through its reaction with the appropriate boronic acid derivative in dioxane using tetrakis catalyst and 2 M sodium carbonate under N2 at 90 °C. In Scheme 2(B,C), syntheses of the target demethoxylated analogs 10a,b and 1,3-dioxolo analogs 12a–d have been accomplished in good to excellent yields (72–87%) adopting the same synthetic routes described for 6a–o using the respective starting compounds, intermediates, and amines.
Scheme 2.
Synthesis of target 4-propoxy-2-arylquinolines (A) 8a,b; (B) 10a,b; (C) 12a–d; Reagents and conditions: (i) Pd(PPh3)4, 2 M Na2CO3, Arylboronic acid, dioxane, 90 °C under N2, 16 h. (ii) KI, K2CO3 anhydrous, the respective amine, dry DMF, reflux, 90 °C, 12 h.
Synthesis of target 4-propoxy-2-arylquinolines (A) 8a,b; (B) 10a,b; (C) 12a–d; Reagents and conditions: (i) Pd(PPh3)4, 2 M Na2CO3, Arylboronic acid, dioxane, 90 °C under N2, 16 h. (ii) KI, K2CO3 anhydrous, the respective amine, dry DMF, reflux, 90 °C, 12 h.Concerning Scheme 3, initiated from 5-bromo-2′-aminoacetophenone 13, the target 6-bromo analogs 16a–d have been afforded in good to excellent yields (78–87%) applying the synthetic routes utilised for 6a–o. Moreover, 5-bromo-2′-aminoacetophenone 13 was transferred to the 6-aryl derivatives 17a,b
via Suzuki coupling reaction by its reaction with the appropriate boronic acid derivative in dioxane/water in the presence of tetrakis catalyst and potassium carbonate under N2 at 100 °C. Subsequently, these 6-aryl derivatives 17a,b have been similarly converted to the target 4-propoxy-2,6-diarylquinolnes 20a–f in good to excellent yields (74–89%) exploiting the same experimental pathways adopted for 6a–o.
Scheme 3.
Synthesis of the target 6-bromo-4-propoxy-2-arylquinolines 16a–d and 4-propoxy-2,6-diarylquinolines 20a–f; Reagents and conditions: (i) the respective benzoyl chloride, Et3N, THF, 0 °C then rt, overnight; (ii) dry dioxane, NaOH, reflux under N2, 110 °C, 4 h; (iii) KI, KOH, 1-Bromo-3-chloropropane, dry DMF, rt, 24 h; (iv) KI, K2CO3 anhydrous, the respective amine, dry DMF, reflux, 90 °C, 12 h; (v) K2CO3, Pd(PPh3)4, Arylboronic acid, dioxane/H2O (1:1), 100 °C under N2, 4 h.
Synthesis of the target 6-bromo-4-propoxy-2-arylquinolines 16a–d and 4-propoxy-2,6-diarylquinolines 20a–f; Reagents and conditions: (i) the respective benzoyl chloride, Et3N, THF, 0 °C then rt, overnight; (ii) dry dioxane, NaOH, reflux under N2, 110 °C, 4 h; (iii) KI, KOH, 1-Bromo-3-chloropropane, dry DMF, rt, 24 h; (iv) KI, K2CO3 anhydrous, the respective amine, dry DMF, reflux, 90 °C, 12 h; (v) K2CO3, Pd(PPh3)4, Arylboronic acid, dioxane/H2O (1:1), 100 °C under N2, 4 h.The structures of all synthesised compounds have been authenticated employing one- dimensional (1D) and two-dimensional (2D) NMR (for 2D NMR, see Supplementary Data), in addition to high resolution mass spectrometry (HRMS).
2.2. In vitro anticancer activity
Antiproliferative activity against different cancer cell lines
The in vitro antiproliferative activity of the synthesised 4-propoxy-2-arylquinolines (6a–o, 8a,b, 10a,b, 12a–d, and 16a–d) and 4-propoxy-2,6-diarylquinolines (20a–f) was preliminary investigated by their screening at one dose level (10 µM) against five cancer cell lines representing three cancer types; colorectal cancer (DLD-1 and HCT-116), breast cancer (MDMBA-231 and MCF-7) and cervical cancer (HeLa). After the incubation of the tested compounds for 24 h, the percent growth inhibition (GI%) was calculated. For DLD-1 cell line, twelve compounds exerted good to excellent growth inhibition ranging from 71.94 to 95.36%, from them compounds 6f, 6h, 6i, 16d, 20e, and 20f displayed the highest GI%. Similarly, HCT-116 cell line established good to excellent sensitivity to the tested compounds and thirteen compounds exhibited GI% ranging from 64.26 to 97.48%, from them compounds 6d, 6f, 6h, 6i, 12c, 16b, 16d, 20e, and 20f demonstrated the best sensitivity (Table 1).
Table 1.
% Growth inhibition (GI%) of all target compounds 6a–o, 8a,b, 10a,b, 12a–d, 16a–d, and 20a–f against different cancer cell lines at 10 µM dose level.
Compounds
GI%
DLD1
HCT-116
MDMBA-231
MCF-7
Hela
6a
22.56
31.34
9.38
17.34
20.15
6b
36.25
42.81
19.67
23.49
12.81
6c
28.22
41.36
19.38
15.76
4.96
6d
75.64
82.61
63.39
71.97
26.98
6e
15.79
22.03
4.67
5.34
22.36
6f
95.36
88.46
75.34
94.36
69.17
6g
25.36
22.16
48.15
76.38
46.28
6h
92.36
89.34
84.76
90.36
90.78
6i
88.76
82.11
81.60
89.49
78.64
6j
6.31
2.46
8.42
14.67
5.67
6k
18.69
22.15
8.64
17.49
11.32
6l
75.34
69.15
38.36
44.71
34.29
6m
22.16
24.31
8.16
22.09
30.19
6n
81.34
64.26
14.39
58.15
71.12
6o
12.36
1.08
8.64
5.36
9.16
8a
5.15
1.02
4.40
7.22
11.78
8b
74.25
52.10
36.87
74.90
50.81
10a
22.51
18.64
37.15
14.26
25.64
10b
35.72
24.10
45.35
64.81
55.49
12a
21.25
16.25
4.98
11.72
22.79
12b
72.95
31.28
14.28
64.38
25.39
12c
71.94
80.49
61.94
43.27
68.22
12d
22.58
54.37
9.36
26.54
36.67
16a
60.08
76.17
51.24
22.97
47.35
16b
55.82
87.61
45.31
19.38
35.94
16c
12.25
18.14
20.11
4.36
8.67
16d
90.11
97.48
80.24
86.38
85.67
20a
44.98
36.71
29.08
15.23
22.06
20b
53.15
71.05
48.31
29.37
34.02
20c
41.34
38.64
32.11
76.82
22.08
20d
24.65
36.91
27.05
15.34
11.97
20e
90.65
89.64
56.29
73.61
65.73
20f
94.18
91.37
82.58
87.36
79.64
% Growth inhibition (GI%) of all target compounds 6a–o, 8a,b, 10a,b, 12a–d, 16a–d, and 20a–f against different cancer cell lines at 10 µM dose level.Regarding breast cancer cell lines, GI% ranged from 75.34 to 84.76% and 71.97 to 87.36% for MDMBA-231 and MCF-7 cell lines, respectively. While five compounds showed good inhibitory activity towards MDMBA-231, ten compounds possessed good inhibition to the growth of MCF-7 cancer cell line. For HeLa cell line, seven compounds had good GI% from 65.73 to 90.78%. Based on the preceding screening results, it was revealed that colorectal cancer cell lines (DLD-1 and HCT-116) were the most sensitive to the tested compounds, therefore they were selected for further antiproliferative assay at six doses levels (Table 1).The growth inhibitory activity of the tested 4-propoxy-2-arylquinolines (6a–o, 8a,b, 10a,b, 12a–d, and 16a–d) and 4-propoxy-2,6-diarylquinolines 20a–f exerted on the colorectal cancer cell lines (DLD-1 and HCT-116) was assessed using MTT assay. DLD-1 and HCT-116 cell lines were incubated for 24 h with increasing concentrations (0.5, 1, 10, 30, 50, and 100 µM) of the tested compounds. Gefitinib (EGFR inhibitor) and TAE226 (FAK inhibitor) were used as reference drugs. The results were presented as half maximal growth inhibitory concentration (IC50) which represents the concentration of a drug exhibiting 50% growth inhibition of the cell line compared to the negative control (Table 2).
Table 2.
The half maximal growth inhibitory concentration (IC50) of all target compounds 6a–o, 8a,b, 10a,b, 12a–d, 16a–d, and 20a–f against two colorectal cancer cell lines (DLD1 and HCT-116) compared to Gefitinib and TAE226.
Compounds
IC50 (µM)a
Compounds
IC50 (µM)
DLD1
HCT-116
DLD1
HCT-116
6a
>100
65.36±4.37
10b
40.28±2.09
>100
6b
19.37±2.15
16.79±1.56
12a
>100
>100
6c
46.95±2.46
13.70±1.29
12b
8.36±1.08
14.36±3.69
6d
8.15±1.05
4.67±0.85
12c
9.12±2.11
6.37±1.09
6e
>100
>100
12d
74.94±5.46
9.20±2.49
6f
2.25±0.96
3.09±1.05
16a
8.22±1.12
6.15±1.30
6g
75.36±5.46
46.82±3.16
16b
8.69±2.64
4.22±1.05
6h
1.79±0.21
3.28±0.67
16c
>100
88.05±5.24
6i
2.48±0.86
5.68±1.42
16d
2.18±0.52
2.43±0.71
6j
>100
>100
20a
13.05±3.69
24.11±2.46
6k
>100
>100
20b
9.65±2.15
7.34±1.82
6l
6.34±0.52
8.11±1.04
20c
15.97±1.10
12.49±3.05
6m
54.21±2.11
43.08±1.05
20d
65.30±5.36
29.47±4.15
6n
6.11±1.80
8.69±0.94
20e
4.46±0.65
4.75±1.02
6o
>100
>100
20f
2.09±0.14
2.96±0.12
8a
>100
>100
Gefitinib
10.24±2.10
6.94±1.24
8b
7.34±1.22
9.84±0.68
TAE226
0.12±0.05
0.17±0.04
10a
>100
>100
aIC50 values are the mean of three separate experiments ± SD.
The half maximal growth inhibitory concentration (IC50) of all target compounds 6a–o, 8a,b, 10a,b, 12a–d, 16a–d, and 20a–f against two colorectal cancer cell lines (DLD1 and HCT-116) compared to Gefitinib and TAE226.aIC50 values are the mean of three separate experiments ± SD.The antiproliferative investigations against DLD-1 revealed that compound 6h emerged as the most potent counterpart showing IC50 = 1.79 µM surpassing the activity of Gefitinib by 5-folds which possessed IC50 = 10.24 µM. Thereafter, compounds 6f, 6i, 16d, and 20f displayed 4-folds superior activity compared to Gefitinib with IC50 values = 2.25, 2.48, 2.18, and 2.09 µM, respectively. In addition, compounds 6d, 6l, 6n, 8b, 12b, 12c, 16a, 16b, 20b, and 20e exerted better inhibitory activity than Gefitinib demonstrating IC50 values = 8.15, 6.34, 6.11, 7.34, 8.36, 9.12, 8.22, 8.69, 9.65, and 4.46 µM, respectively. The rest of compounds had week to moderate or no activity relative to Gefitinib.Concerning HCT-116 cell line, it was found that compounds 6f, 6h, 16d, and 20f established twice inhibitory activity compared to Gefitinib (IC50 = 6.94 µM) displaying IC50 values = 3.09, 3.28, 2.43, and 2.96 µM, respectively. Furthermore, compounds 6d, 6i, 12c, 16a, 16b, and 20e exhibited higher growth inhibitory activities than Gefitinib with IC50 values = 4.67, 5.68, 6.37, 6.15, 4.22, and 4.75 µM, respectively. The remaining compounds possessed week to moderate or no activity compared to Gefitinib.It is noteworthy that appending p-CF3 to the 2-phenyl enhanced the antiproliferative activity compared to p-Cl and p-CH3, except for the cyclohexylamine analogs 6b, 6g, and 6l in which the p-CH3 substituent is preferred for antiproliferative activity. Remarkably, the incorporation of p-(2-furyl) to the 2-phenyl abolished the growth inhibitory activity, while the introduction of p-(2-thienyl) increased the activity compared to p-Cl and p-CH3 analogs. In the context of impact of amine substituents on the antiproliferative activity, it was proved that 4-amino-N-methylpiperidine is preferred for anticancer activity, then piperazine, tetrahydrofurfurylamine and cyclohexylamine. Notably, grafting imidazole along with 6,7-dimethoxy substituents abolished the antiproliferative activity for all derivatives.Besides, the removal of 6,7-dimethoxy groups from the quinoline scaffold decreased or abolished activity. While the fusion of 1,3-dioxolo to the quinoline scaffold along with p-CF3 substitution on phenyl dramatically potentiated the anticancer activity of the imidazole derivatives, the replacement of imidazole with morpholine decreased the activity of p-CF3 analogs and markedly increased the activity of p-Cl counterparts. On the other hand, appending of electron withdrawing (Br) to position 6 of quinoline parallel with p-Cl substitution on 2-phenyl greatly enhanced the activity of the imidazole derivatives while, p-F substitution on the 2-phenyl extremely decreased the anticancer activity of imidazole derivatives. Also, the replacement of imidazole with morpholine almost had no impact on the p-Cl analogs, but tremendously elevated the anticancer activity of p-F derivatives.Furthermore, the incorporation of 4-methoxyphenyl or 2-furyl to position 6 of quinoline along with p-Cl or F on the 2-phenyl enhanced the activity of the imidazole derivatives compared to the dimethoxy analogs, but the 2-furyl derivatives exhibited better activity. Moreover, the replacement of imidazole with morpholine elevated the antiproliferative activity of p-Cl analogs while, diminished the activity of p-F derivatives.Finally, the deduced structure activity relationships indicated that the substitution pattern on positions 6 and 7 of quinoline and position 4 of the 2-phenyl, in addition to the amine substitution on the 4-propoxy linker are crucial elements for the anticancer activity. In general, incorporation of dimethoxy groups at positions 6 and 7 of quinoline along with p-CF3 at 2-phenyl and 4-amino-N-methylpiperidine, piperazine or tetrahydrofurfurylamine on the 4-propoxy linker, in addition to grafting of electron withdrawing (Br) or 2-furyl at position 6 of quinoline along with p-F or Cl on the 2-phenyl and morpholine on the propoxy linker resulted in the most potent antiproliferative agents in this study.
The apoptotic impact of the most potent antiproliferative agents 6f, 6h, 6i, 16d, and 20f on DLD-1 colorectal cancer cell line was investigated exploiting AV/PI dual staining assay. The assay outcomes proved that the tested compounds elicited apoptosis of such cell line as indicated by significant rise in the total percentage of AV positive apoptotic DLD-1 cells compared to the control (Figure 4). Compound 6h increased the total percentage of apoptotic cells from 7% for the control to 90.33%. Also, compounds 6f, 6i, 16d, and 20f exerted potential apoptotic effect elevating the total percentage of apoptotic cells to 84.33, 71.66, 51.66, and 80.66%, respectively.
Figure 4.
Influence of the promising compounds on the total percentage of AV-FITC positive staining in DLD1 cancer cell line.
Influence of the promising compounds on the total percentage of AV-FITC positive staining in DLD1 cancer cell line.
Topoisomerase I-mediated DNA cleavage assay
TOP1 poisoning activity of all target compounds has been estimated utilising TOP1-mediated DNA cleavage assay that determines the TOP1 poisoning activity relative to 1 µM Camptothecin (CPT). The tested compounds were incubated at 0.1, 1, 10, and 100 µM with recombinant human TOP1 enzyme and a 3′-[32P]-labeled 117-bp DNA oligonucleotide. TOP1 poisoning agents specifically trap TOP1ccs leading to their stabilisation and DNA cleavage. The drug-induced stabilised TOP1ccs are visualised by gel electrophoresis demonstrating specific DNA cleavage patterns. Then a semiquantitative scoring system by visual comparison between lanes induced by the target compounds and 1 µM CPT was used to score the compounds (Table 3, see Supplementary Data for the results of gel electrophoresis).
Table 3.
The TOP1 inhibitory activity of all target compounds 6a–o, 8a,b, 10a,b, 12a–d, 16a–d, and 20a–f compared to Camptothecin (CPT).
The TOP1 inhibitory activity of all target compounds 6a–o, 8a,b, 10a,b, 12a–d, 16a–d, and 20a–f compared to Camptothecin (CPT).aScoring: 0: no activity; −/+: 0–25% 1 µM CPT; +: 25–50% 1 µM CPT.Compounds 6a–e, 6h, 6l, 6m, 12a, and 12b exhibited weak TOP1 poisoning effect (-/+) displaying DNA cleavage activity equal to 0–25% of the activity of 1 µM CPT, while compound 16c demonstrating activity (+) equals to 25–50% of the activity of 1 µM CPT. The rest of compounds possessed no cleavage activity. Despite the target compounds showing promising antiproliferative activity against cancer cell lines, their TOP1 poisoning activities were not encouraging for further development as potent TOP1 inhibitors. Accordingly, the synthesised compounds have been evaluated for the plausible mechanism by which they provoked the antiproliferative activity.
Kinase inhibitory activities of target quinolines
Kinase profiling
The non-significant TOP1 poisoning activities of target compounds obtained from Topoisomerase I-mediated DNA cleavage assay motivated us to search for the plausible molecular mechanism for herein reported 4-propoxy-2-arylquinolines.The potential inhibitory activity of the target 4-propoxy-2-arylquinolines (6a–o, 8a,b, 10a,b, 12a–d, 16a–d) and 4-propoxy-2,6-diarylquinolines (20a–f) was explored against a panel of nine kinases representing different signalling pathways; EGFR, FAK, FRK, IGF-1R, BTK, c-Src, VEGFR-1, VEGFR-2 and HER-2 (see Supplementary Data, Table S1). The half maximal inhibitory concentration (IC50) values were calculated for each kinase and presented in Table S1 and Table 4. Strikingly, the screening outcomes revealed that the investigated quinolines exhibited promising dual inhibitory effect towards EGFR and FAK kinases.
Table 4.
The half maximal inhibitory concentration (IC50) of all target compounds 6a–o, 8a,b, 10a,b, 12a–d, 16a–d, and 20a–f against EGFR and FAK kinase activity compared to Gefitinib and TAE226.
Compounds
IC50 (nM)a
Compounds
IC50 (nM)
EGFR
FAK
EGFR
FAK
6a
142.64±2.54
214.36±1.09
10b
179.64±9.12
164.74±5.37
6b
45.26±5.36
98.16±4.67
12a
124.97±7.94
225.46±14.02
6c
85.67±6.46
45.70±3.40
12b
95.36±2.05
70.85±3.16
6d
46.37±4.09
36.97±2.34
12c
246.70±12.29
111.06±8.94
6e
156.72±11.36
211.08±8.96
12d
450.16±4.25
273.16±2.84
6f
25.39±3.49
22.68±2.38
16a
222.15±8.25
44.15±3.26
6g
365.49±14.82
145.71±10.54
16b
313.34±15.34
50.36±4.81
6h
20.15±1.07
14.25±2.72
16c
485.46±11.37
224±10.46
6i
22.36±2.05
18.36±3.17
16d
35.03±2.64
17.36±2.15
6j
258.34±11.94
186.46±6.22
20a
121.74±9.40
77.25±4.37
6k
410.38±12.73
157.84±8.73
20b
245.11±12.34
63.25±3.25
6l
34.91±3.76
26.37±2.81
20c
362.30±5.26
91.03±5.85
6m
154.29±12.80
172.49±13.67
20d
146.95±8.37
125.38±3.15
6n
41.82±2.34
44.36±2.94
20e
33.65±1.02
25.36±3.48
6o
256.19±6.94
204.84±8.04
20f
24.81±2.71
15.36±0.98
8a
349.37±14.05
198.32±12.32
Gefitinib
48.52±3.64
–
8b
35.48±1.50
29.79±2.37
TAE226
–
4.60±0.94
10a
244.30±8.41
298.74±1.94
aIC50 values are the mean of three separate experiments ± SD.
The half maximal inhibitory concentration (IC50) of all target compounds 6a–o, 8a,b, 10a,b, 12a–d, 16a–d, and 20a–f against EGFR and FAK kinase activity compared to Gefitinib and TAE226.aIC50 values are the mean of three separate experiments ± SD.
EGFR and FAK kinase inhibitory activity
All the newly prepared 4-propoxy-2-arylquinolines (6a–o, 8a,b, 10a,b, 12a–d, and 16a–d) and 4-propoxy-2,6-diarylquinolines (20a–f) were examined for their potential EGFR and FAK inhibitory activities. Gefitinib and TAE-226 were used as reference EGFR and FAK inhibitors, respectively. The results are reported as half maximal inhibitory concentration values (IC50), as determined from triplicate measurements and are presented in Table 4.Results in Table 4 revealed that the examined quinolines displayed moderate to potent inhibitory activity towards EGFR (IC50 values ranging between 20.15 ± 1.07 and 485.46 ± 11.37 nM, Table 4). In particular, trifluoromethyl phenyl-bearing 6,7-dimethoxy-2-arylquinolines 6f, 6h, and 6i, as well as 6-furanyl-2-arylquinoline 20f emerged as the most efficient EGFR inhibitors with two-digits nanomolar IC50s (IC50 = 25.39 ± 3.49, 20.15 ± 1.07, 22.36 ± 2.05, and 24.81 ± 2.71 nM, respectively). Notably, these four derivatives displayed 2-fold higher activity than the reference EGFR inhibitor Gefitinib (IC50 = 48.52 ± 3.64 nM). In addition, compounds 6b, 6d, 6l, 6n, 8b, 16d, and 20e exhibited potent EGFR inhibitory activity, as the measured IC50 values ranged between 33.65 ± 1.02 and 46.37 ± 4.09 nM, which are slightly improved or comparable to that of the reference drug Gefitinib (Table 4). Moreover, compounds 6c and 12b showed 2-fold decreased activity (IC50 = 85.67 ± 6.46 and 95.36 ± 2.05 nM, respectively) than Gefitinib against EGFR. The remaining examined quinolines possessed moderate EGFR inhibitory activity (IC50 range: 121.74 ± 9.40–485.46 ± 11.37 nM) compared to Gefitinib (Table 4). Strikingly, the inclusion of 4-amino-N-methylpiperidine, tetrahydrofurfurylamine and N-methylpiperazine along with 2-(p-CF3 phenyl) and 6,7-dimethoxy substituents (6f, 6h, and 6i), in addition to the grafting of morpholine together with 2-(p-Cl phenyl) and 6–(2-furyl) 20f afforded the most potent EGFR inhibitors in this study displaying IC50 range from 20.15 ± 1.07 to 25.39 ± 3.49 nM.On the other hand, as depicted in Table 4, FAK kinase was efficiently inhibited by all 4-propoxy-2-arylquinolines (6a–o, 8a,b, 10a,b, 12a–d, and 16a–d) and 4-propoxy-2,6-diarylquinolines (20a–f) herein reported in the nanomolar range (IC50 range: 14.25 ± 2.72–298.74 ± 1.94 nM). Superiorly, p-CF3-phenyl-bearing 6,7-dimethoxy-2-arylquinolines 6f, 6h, and 6i, as well as, morpholine-bearing 6-bromo-2-arylquinoline 16d and 6-furanyl-2-arylquinoline 20f were the most potent FAK inhibitors in this study with IC50 values equal 22.68 ± 2.38, 14.25 ± 2.72, 18.36 ± 3.17, 17.36 ± 2.15, and 15.36 ± 0.98 nM, respectively (Table 4). Moreover, compounds 6c, 6d, 6l, 6n, 8b, 16a, 16b, and 20e exerted potent FAK inhibitory activity with IC50 spanning in the range 25.36 ± 3.48–50.36 ± 4.81 nM.Further analysis of the obtained results in Table 4 revealed that compounds 6b, 12b, and 20a–c exhibited two-digit nanomolar IC50s; 98.16 ± 4.67, 70.85 ± 3.16, 77.25 ± 4.37, 63.25 ± 3.25, and 91.03 ± 5.85 nM, respectively, whereas the remaining derivatives displayed moderate inhibitory activity against FAK kinase (IC50 range: 111.06 ± 8.94–298.74 ± 1.94 nM) (Table 4). Interestingly, the incorporation of 4-amino-N-methylpiperidine, tetrahydrofurfurylamine and N-methylpiperazine along with 2-(p-CF3 phenyl) and 6,7-dimethoxy substituents (6f, 6h, and 6i), besides the appending of morpholine with 2-(p-F phenyl) and 6-Br 16d, as well as the addition of morpholine in conjunction with 2-(p-Cl phenyl) and 6–(2-furyl) 20f provided the most potent FAK inhibitors in this study demonstrating IC50 values ranging from 14.25 ± 2.72 to 22.68 ± 2.38 nM.It is worth stressing that 4-propoxy-2-arylquinolines 6f, 6h, 6i, and 20f emerged not only as the most potent dual EGFR/FAK inhibitors in this study, but also as the most efficient anti-proliferative agents towards the examined colorectal cancer (DLD-1 and HCT-116) cell lines.
In silico molecular docking
Docking into EGFR binding site
The molecular docking approach was utilised to investigate the potential binding of herein reported 4-propoxy-2-arylquinolines to EGFR binding site (PDB: 1M17). The docking procedure was validated through the redocking of the co-crystalised ligand. The correct pose was predicted accurately with RMSD of 1.498 between the docked and co-crystalised ligand using DockRMSD server (Figure 5(a)). In addition, docking was able to maintain hydrogen bonding seen in the co-crystalised ligand with NH of M769 (2.7 Å) and with the NH of G772 (3.2 Å). Furthermore, hydrophobic interactions with residues in the active site were also maintained. These include interactions between K721 and ethyne benzene moiety, and L694, L768, and L820 with hydrophobic part of the quinazoline ring (Figure 5(b)).
Figure 5.
Docking of target quinolines in the active site EGFR. (a) validation of docking procedure showing overlapping of crystalised (blue) and docked (pink) poses; (b) interactions of Erlotinib with EGFR; (c) docking pose of 20c; (d) docking pose of 20e.
Docking of target quinolines in the active site EGFR. (a) validation of docking procedure showing overlapping of crystalised (blue) and docked (pink) poses; (b) interactions of Erlotinib with EGFR; (c) docking pose of 20c; (d) docking pose of 20e.Docking scores of the tested compounds with EGFR are shown in Table 5. The docking score of the co-crystalised ligand was −7.2 kcal/mol, whereas all the tested quinolines have shown better docking scores than that of the co-crystalised ligand (−7.9 to −9.7 kcal/mol). Best docking scores were seen with quinolines form the series 20. Compound 20c showed the best binding energy to EGFR with docking score of −9.7 kcal/mol and its binding pose is shown in Figure 5(c). The compound has formed 2 hydrogen bonds between imidazole ring and both K721 (3.5 Å) and T766 (3.1 Å). Also, several hydrophobic interactions have been seen, which included interactions between quinoline ring and side chains of V702 and T830. The phenyl ring at position 2 formed hydrophobic interactions with L694 and L820 which are also common with the co-crystalised ligand. In addition, π- π stacking between the ring at quinoline position 6 and F699.
Table 5.
Docking results of target compounds with EGFR and FAK.
Compound
Docking Score (kcal/mol)
EGFR (PDB: 1M17)
FAK (PDB: 2JKM)
6a
−8.1
−7.7
6b
−8.5
−7.3
6c
−8.2
−7.5
6d
−7.9
−7
6e
−8.4
−7.6
6f
−8.2
−7.4
6g
−8.5
−7.6
6h
−8.2
−7.6
6i
−8.1
−7.3
6j
−8.6
−8.0
6k
−8.1
−7.7
6l
−8.7
−7.4
6m
−8.3
−7.5
6n
−8.2
−7.1
6o
−8.4
−7.7
8a
−8.3
−7.6
8b
−8.4
−7.5
10a
−8.8
−7.9
10b
−8.3
−8.0
12a
−8.5
−8.0
12b
−8.5
−8.3
12c
−8.8
−8.4
12d
−8.7
−8.3
16a
−8.1
−7.7
16b
−8.2
−8.0
16c
−8.3
−7.7
16d
−8.3
−7.9
20a
−9.3
−8.3
20b
−9.1
−8.2
20c
−9.7
−8.3
20d
−9.2
−8.2
20e
−9.3
−8.2
20f
−8.8
−7.7
Co-crystalised ligand
−7.2 (Erlotinib)
−8.5 (AZW592)
Docking results of target compounds with EGFR and FAK.Another compound from this series is compound 20e which was selected as a representative example because it has shown good biological results with both EGFR and FAK. The docking pose of this compound is shown in Figure 5(d) showing a similar docking pose to 20c. The same hydrogen bonds with K721 (3.5 Å) and T766 (3.1 Å) were maintained as well as hydrophobic interactions with V702 and T830 as well as with L694 and L820. The π–π stacking was also seen between F699 and the furan ring of 20e. This compound, in complex with EGFR, was subjected to further investigation using molecular dynamics to study the stability of its complex with EGFR as will be discussed later.
Docking into FAK binding site
Potential binding of target quinolines to FAK was also investigated using docking studies (PDB: 2JKM). Initially, the docking procedure was validated through the redocking of the co-crystalised ligand (AZW592). The docking searching algorithm was able to correctly predict the binding pose with acceptable accuracy with RMSD of 1.318 between the docked and co-crystalised ligands as predicted by DockRMSD server (Figure 6(a)). The docked structure was able to maintain same hydrogen bonds that are in the crystal structure including those between the sulphamoyl moiety oxygen and the terminal amino group of K454 (Å) and the hydrogen bond with the α-carbonyl group of C502. In addition, several hydrophobic interactions have been also seen with residues in the active site including I428, V436, V484, L501, and L553 (Figure 6(b)).
Figure 6.
Docking of target compounds in the active site FAK. (a) validation of docking procedure showing overlapping of crystalised (blue) and docked (pink) poses; (b) interactions of AZW592 with FAK; (c) general binding of 16 and 20 compound series; (d) docking pose of 20e.
Docking of target compounds in the active site FAK. (a) validation of docking procedure showing overlapping of crystalised (blue) and docked (pink) poses; (b) interactions of AZW592 with FAK; (c) general binding of 16 and 20 compound series; (d) docking pose of 20e.Next, target synthesised quinolines were docked in the active site of the FAK after their preparation. The docking scores of tested compounds are shown in Table 5. The docking score of co-crystalised ligand (AZW592) was found to be −8.5 kcal/mol. Some of tested compounds have shown scores that are comparable to the co-crystalised ligand. Best results were seen with 12 and 20 and some of the 16 series. Docking poses of these compounds were similar with most of the docked compounds as can be seen in Figure 6(c) which shows docking pose of some compounds from these series.Docking pose of compound 20e which was chosen as representative example is shown in Figure 6(d). The compound was able to form 3 hydrogen bonds with N551 (3.1 Å), D564 (3.1 Å), and E506 (3.3 Å). In addition, several hydrophobic interactions were also seen, such as the hydrophobic interaction between quinoline ring and L553 and between phenyl ring at position 2 and I428 which are common with the co-crystallized ligand. This pose was selected for further investigation of the complex stability using molecular dynamics study.
Molecular dynamics (MD) simulation
The stability of compound 20e complexes with both EGFR and FAK was investigated using 100 ns molecular dynamics studies. With each complex, the results were compared with the co-crystalised ligand complex as a control and with the apoprotein (the protein alone with no ligands). The missing loops in both targets were built using Swiss-Model server before starting the dynamics to ensure correct results. All complexes were equilibrated under NVT then NPT conditions for 1 ns each and the analysis was done on the production run.Analysis of the production runs trajectories for 20e in the active site of EGFR demonstrated stability comparable to the co-crystalised ligand. Radius of gyration (R) is a measure of the compactness of the complexes. Stable R suggested the stability of the protein or complex under investigation. Figure 7(a) shows a plot of R of 20e, co-crystalised ligand and apoprotein. The average Rg was found to be 2.01 ± 0.02, 2.03 ± 0.02, and 2.02 ± 0.01 nm for apoprotein, control, and 20e, respectively. In addition, Root mean square fluctuation (RMSF) of protein residue (Figure 7(b)) for all the three complexes showing similar patterns. The average RMSF for co-crystalised ligand and 20e was found to be 0.19 and 0.18 nm, respectively which is slightly higher than that of the apoprotein (0.15 nm). Although root mean square deviation (RMSD) of ligand heavy atoms for 20e is slightly higher than that of the co-crystalised ligand (Figure 7(c)), the value is <1 nm for most of the trajectory. This value cannot be calculated for the apoprotein as it has no ligand in the system. Finally, the number of hydrogen bonds between ligands and protein are shown in Figure 7(d), which showed that 20e formed extra hydrogen bonds during at least 50% of the production run time. These results suggested that 20e complex with EGFR is at least of comparable stability when compared to the complex with the co-crystalised ligand; Erlotinib.
Figure 7.
Molecular dynamics analysis of the production run trajectory of 20e in the active site of EGFR compared to control and apoprotein. (a) Radius of gyration; (b) Root mean square fluctuation of residues; (c) Root mean square deviation of ligand heavy atoms; (d) Number of hydrogen bonds between ligand and protein.
Molecular dynamics analysis of the production run trajectory of 20e in the active site of EGFR compared to control and apoprotein. (a) Radius of gyration; (b) Root mean square fluctuation of residues; (c) Root mean square deviation of ligand heavy atoms; (d) Number of hydrogen bonds between ligand and protein.Complexes of FAK with 20e, its co-crystalised ligand and the apoprotein showed similar pattern to that of the EGFR (Figure 8). This includes the radius of gyration (Rg), which showed an average of 2.00 ± 0.01, 1.99 ± 0.01, and 1.96 ± 0.01 nm for apoprotein, co-crystalised ligand, and 20e, respectively (Figure 8(a)). Also, RMSF of protein residues was found to follow similar patterns for all the three studied systems (Figure 8(b)). The average RMSF for the three systems was found to be 0.12 ± 0.07, 0.13 ± 0.09, and 0.12 ± 0.07 nm for apoprotein, co-crystalised ligand, and 20e, respectively. In addition, plotting of RMSD of ligand heavy atoms (Figure 8(c)) showed minimal fluctuation for both with and average RMSD of 0.19 ± 0.05 and 0.57 ± 0.16 nm for co-crystalised ligand and 20e, respectively. Although, the value for 20e is higher but it is within acceptable range (<1 nm). Finally, plotting of hydrogen bonds between ligands and target protein (Figure 8(d)) showed that the number of hydrogen bonds is higher in case of the co-crystalised ligand compared to 20e. Being said, 20e was still able to maintain an average of 1.93 ± 0.45 hydrogen bonds during the 100 ns production run.
Figure 8.
Molecular dynamics analysis of the production run trajectory of 20e in the active site of FAK compared to control and apoprotein. (a) Radius of gyration; (b) Root mean square fluctuation of residues; (c) Root mean square deviation of ligand heavy atoms; (d) Number of hydrogen bonds between ligand and protein.
Molecular dynamics analysis of the production run trajectory of 20e in the active site of FAK compared to control and apoprotein. (a) Radius of gyration; (b) Root mean square fluctuation of residues; (c) Root mean square deviation of ligand heavy atoms; (d) Number of hydrogen bonds between ligand and protein.These results collectively suggested the stability of compound 20e complexes with both EGFR and FAK compared to the corresponding co-crystalised ligands in each target. This in general supported the dual mechanism similar to the enzymatic inhibition data.
Conclusion
Different series of 4-propoxy-2-arylquinolines (6a–o, 8a,b, 10a,b, 12a–d, and 16a–d) and 4-propoxy-2,6-diarylquinolines (20a–f) have been designed and synthesised as potential anticancer agents. The quinolines 6f, 6h, 6i, 16d, and 20f demonstrated the most potent antiproliferative effect against DLD-1 colorectal cancer with respective IC50 values = 2.25, 1.79, 2.48, 2.18, and 2.09 µM with 4- to 5-folds potency compared to Gefitinib (IC50 = 10.24 µM). Additionally, compounds 6f, 6h, 16d, and 20f possessed twice growth inhibitory impact as Gefitinib (IC50 = 6.94 µM) displaying IC50 values = 3.09, 3.28, 2.43, and 2.96 µM against HCT-116 cell line, respectively. Moreover, compounds 6f, 6h, 6i, 16d, and 20f significantly elevated the total percentage of DLD-1 apoptotic cells. Furthermore, the quinolines 6f, 6h, 6i, and 20f exerted potent EGFR inhibitory effects with IC50 values = 25.39 ± 3.49, 20.15 ± 1.07, 22.36 ± 2.05, and 24.81 ± 2.71 nM, respectively compared to Gefitinib (IC50 = 48.52 ± 3.64 nM). In a similar fashion, the quinolines 6f, 6h, 6i, 16d, and 20f displayed the best FAK inhibitory actions with IC50 values = 22.68 ± 2.38, 14.25 ± 2.72, 18.36 ± 3.17, 17.36 ± 2.15, and 15.36 ± 0.98 nM, respectively. Molecular docking and molecular dynamics simulation rationalised EGFR/FAK dual inhibition providing different quinolines being as the first reported quinolines possessing potential EGFR/FAK dual inhibition. The latter compounds can be used as lead compounds for the development of more potent EGFR/FAK dual inhibitors as potential anticancer agents.
Experimental
General
Melting points have been measured by Yanaco melting point device and were uncorrected. NMR spectra were measured using Bruker Advance III HD at 400 MHz for 1H NMR, 100 MHz for 13C NMR and 376 MHz for 19F NMR in deuterated CDCl3 or DMSO-d6 using tetramethyl silane (TMS) as an internal standard. Coupling constant values (J) were determined in Hertz (Hz) and chemical shifts (δ) were expressed in ppm. High resolution mass spectrometry (HRMS) have been measured by Thermo Fisher Scientific LTQ Orbitrap XL spectrophotometer using electrospray ionisation (ESI) and the results were expressed as [M + H]+ or [M + Na]+ at Natural Science Research and Development Centre, Hiroshima University, Japan. The purities of all biologically tested compounds were determined by HPLC and were found to be ≥95%. HPLC analysis was performed utilising JASCO 880-PU HPLC system (Japan spectroscopic Co. Ltd) connected to a diode array detector with detection at a wavelength of 254 nm. The column exploited in the HPLC analysis was Inertsil ODS-3 column with dimensions of 250 × 4.6 mm and 5 µm particle size (GL SCIENCES INC., Japan). The mobile phase employed for HPLC analysis was acetonitrile/water/TFA (29.9/70/0.1, v/v) at a flow rate of 0.5 ml/min. The reactions have been monitored by thin layer chromatography (TLC) using Merck silica gel 60F254 aluminium sheets. Column chromatography has been performed utilising silica gel 60 N, 63–210 µm that was purchased from Kanto Chemical Co. Inc., Japan using dichloromethane/methanol (100/0 to 90/10, v/v) and hexane/ethylacetate (100/0 to 90/10, v/v). Unless otherwise stated, all chemicals and solvents were available commercially and have been used without further purification.
Synthesis of 1-(2-amino-4,5-dimethoxyphenyl)ethan-1-one (1)
General procedure for the synthesis of benzamides (3a–c)
In dry THF (8 ml) and Et3N (2 ml), 1–(2-amino-4,5-dimethoxyphenyl)ethan-1-one 1 (0.976 g, 5 mmol) was dissolved and cooled in ice bath. Then, a solution of the respective p-substituted benzoyl chloride 2a–c (5.1 mmol) in dry THF (2 ml) was added dropwise while cooling in ice bath. The reaction mixture was stirred in ice bath for 30 min and then overnight at room temperature. After that, reaction mixture was poured into ice/water and the resulting solid was filtered off and washed excessively with water and methanol to afford the corresponding 3a–c.The spectral characterisation of the benzamides 3a,b were reported in our previous study.
General procedures for synthesis of the quinolones (4a–c)
Under N2 atmosphere, the benzamides 3a–c and three equivalents NaOH were refluxed in dry dioxane at 110 °C for 4 h and then cooled to room temperature. Small amount of water and excess amount of hexane were added to the reaction mixture. The resulting mixture was subjected to sonication for 2 min and then neutralised using 1 M HCl. The separated solid was filtered off and washed excessively with water to give the corresponding quinolones 4a–c.The spectral data of the quinolones 4a,b have been reported in our previous study.
General procedure for synthesis of the key intermediates (5a–c)
The quinolones 4a–c (3 mmol), KI (0.498 g, 3 mmol) and KOH (1.009 g, 18 mmol) were stirred for 2 h in dry DMF (30 ml) and then 1-bromo-3-chloropropane (2.361 g, 15 mmol) was added to the mixture and stirred at room temperature for 24 h. Thereafter, the mixture was poured into ice/water and the separated solid was filtered off and washed with water then hexane to furnish the key intermediates 5a–c which were used without further purification.The spectral data of 5a,b have been reported in our previous study.
General procedure for synthesis of the target 4-propoxy-2-arylquinolines (6a–o)
To a stirred mixture of 5a–c (1 mmol), anhydrous K2CO3 (1.38 g, 10 mmol) and KI (0.83 g, 5 mmol) in dry DMF (20 ml), the respective amine (10 mmol) was added. Then, the mixture was refluxed at 90 °C for 12 h and poured into ice/water (50 ml). The separated solid was filtered off then washed with water and hexane. The products were purified by silica gel column chromatography using DCM/MeOH to furnish the pure target compounds 6a–o.
General procedure for synthesis of heterocyclic analogs of the bromo derivative (8a,b)
The bromo derivative 7 (125 mg, 0.25 mmol) was taken with 2-furylboronic acid or 2-thienylboronic acid (0.5 mmol) and Pd(PPh3)4 (0.05 equivalent, 15 mg), then dioxane (5 ml) and 2 M Na2CO3 (0.3 ml) were added to the mixture under N2 atmosphere. The reaction mixture was refluxed under N2 at 90 °C for 16 h. Then, the reaction mixture was poured into ice/water (50 ml), the aqueous layer was extracted with ethyl acetate (50 × 3) and the organic layers were washed with water and brine. After evaporation of the organic solvent under vacuum, the residue was purified by silica gel column chromatography using DCM/MeOH to furnish 8a,b in pure form.
Synthesis of the target demethoxylated 4-propoxy-2-arylquinolines (10a,b)
The demethoxylated 4-propoxy-2-arylquinolines 10a,b were prepared according to the general synthetic procedure for 6a–o using the respective key intermediate 9.
Synthesis of 1,3-dioxoloarylquinolines key intermediates (11a,b)
Starting from 6′-amino-3′,4′-(methylenedioxy)acetophenone and the respective p-substituted benzoyl chloride, the key intermediates 11a,b were prepared according to the synthetic route for 5a–c.
Synthesis of the target propoxy derivatives of 1,3-dioxoloarylquinolines (12a–d)
The target dioxolo derivatives 12a–d have been synthesised utilising the synthetic procedures used for 6a–o using imidazole or morpholine with the appropriate key intermediate 11a,b.
Synthesis of 1-(2-amino-5-bromophenyl)ethan-1-one (13)
The 5-bromo derivative of 2-aminoacetophenone 13 has been synthesised based on the reported procedure,.Yellow solid, Yield 96%, m.p. 77–79 °C; 1H NMR (400 MHz, CDCl3) δ (ppm): 2.55 (s, 3H, CH3), 6.29 (s, 2H, NH2), 6.54 (d, 1H, phenyl CH, J = 8.8 Hz), 7.31 (dd, 1H, phenyl CH, J = 2.3, 8.8 Hz), 7.78 (d, 1H, phenyl CH, J = 2.3 Hz); 13 C NMR (100 MHz, CDCl3) δ (ppm): 27.83, 106.63, 118.99, 119.40, 134.11, 136.99, 149.08, 199.64; HRESIMS (m/z): [M + H]+ Calcd for C8H9CBrNO, 213.98620; found, 213.98599.
Synthesis of 6-bromo-2-arylquinolones (14a,b)
1-(2-amino-5-bromophenyl)ethan-1-one 13 was benzoylated with p-chloro or fluorobenzoyl chloride using the same procedure for 3a–c. Then, the resulted benzoyl derivatives have been subjected to ring closure reaction according to the synthetic route for 4a–c to afford the corresponding 6-bromo-2-arylquinolones 14a,b.
Synthesis of the target 6-bromo-4-propoxy-2-arylquinolines (16a–d)
The synthesis of the target 6-bromo-4-propoxy-2-arylquinolines 16a–d has been accomplished based on the general synthetic route for 6a–o utilising imidazole or morpholine with the respective key intermediate 15a,b.
To a mixture of 1–(2-amino-5-bromophenyl)ethan-1-one 13 (1.07 g, 5 mmol), 4-methoxyphenylboronic acid or 2-furylboronic acid (5.5 mmol), K2CO3 (2.28 g, 16.5 mmol) and Pd(PPh3)4 (0.02 equivalent, 116 mg), dioxane (14 ml) and H2O (14 ml) were added under N2 atmosphere. The reaction mixture was refluxed under N2 at 100 °C for 4 h. Then, the reaction mixture was poured into ethylacetate (50 ml) and the organic layer was separated. The aqueous layer was extracted with ethylacetate (30 × 3) and the organic layers were collected and washed with 1 M HCl (100 × 3) then brine ((100 × 3). After evaporation of the organic solvent under vacuum, the residue was purified by silica gel column chromatography using hexane/ethylacetate to afford the corresponding compounds 17a,b.
The synthesis of 2,6-diarylquinolones 18a–c has been accomplished using 5-aryl-2-aminoacetophenones 17a,b and the respective p-chloro or fluorobenzoyl chloride according to the synthetic route used for 6-bromo-2-arylquinolones 14a,b. The 2,6-diarylquinolones 18a,b have poor solubility for NMR spectral analysis, so the NMR spectral analysis of the soluble 6-furyl analog 18c was used for their structural authentication in addition to HRMS for 18a,b. Moreover, the corresponding 4-propoxy analogs 19a–c exhibited good solubility and their spectral characterisation was enough for structural confirmation of 2,6-diarylquinolones 18a–c.
Synthesis of 4-propoxy-2,6-diarylquinolines key intermediates (19a–c)
The key intermediates 4-propoxy-2,6-diarylquinolines 19a–c were synthesised from 2,6-diarylquinolones 18a–c utilising the synthetic route used for 5a–c.
Synthesis of the target 4-propoxy-2,6-diarylquinolines (20a–f)
The synthesis of the target 4-popoxy-2,6-diarylquinolines 20a–f has been conducted by reaction of imidazole or morpholine with the corresponding key intermediate 19a–c under the same conditions used for synthesis of 6a–o.
The antiproliferative assay against five cancer cell lines representing three different tumour subpanels, including colorectal (DLD-1, HCT-116), breast (MDMBA-231, MCF-7), and cervical (HeLa) cell lines was conducted using MTT assay. Cells were seeded at 1 × 104 cells/well and cultured overnight in a 96-well plate. The cells were treated with either 10 µМ of tested quinolines 6a–o, 8a,b, 10a,b, 12a–d, 16a–d, and 20a–f, or DMSO as a negative control. After 24 h, the cells washed with phosphate buffered saline (PBS, Invitrogen Gibco) and incubated with 20 µl of MTT solution (2 mg/ml) for 4 h at 37 °C. Then, 150 µl DMSO was used to solubilise MTT formazan crystals. Finally, the plates were shaken, and the optical density was determined at 570 nm using ELISA plate reader. At least, three independent experiments were performed. Percentage of growth inhibition was determined as (1-[OD of treated cells/OD of control cells]). On the other hand, using the MTT assay, we tested the effect of different concentrations (0.5, 1, 10, 30, 50, and 100 µM)) of the synthesised compounds on colorectal cancer cell lines (DLD-1 and HCT-116), using DMSO as a negative control, whereas Gefitinib and TAE226 were used as positive controls. The IC50 values were calculated using Prism v0.8 software (GraphPad Software Inc., La Jolla, CA).
Apoptosis assay
The apoptotic effect of the most potent antiproliferative agents 6f, 6h, 6i, 16d, and 20f on DLD-1 colorectal cancer cell line was investigated using the annexin V/propidium iodide (AV/PI) staining kit (BioLegend, San Diego, CA, USA) according to the manufacturer's instructions. DLD-1 cells were treated with 3 µM of the most potent antiproliferative compounds 6f, 6h, 6i, 16d, and 20f or DMSO as a negative control then incubated for 24 h. Flow cytometric analysis was conducted using FACS Calibur flow cytometer (BD Biosciences, San Jose, CA, USA).A 3′-[32P]-labeled 117-bp DNA substrate oligonucleotide was prepared as described previously. Radiolabeled DNA was incubated with recombinant human TOP1 in 20 µL reaction buffer (10 mmol/L Tris-HCl, pH 7.5, 50 mmol/L KCl, 5 mmol/L MgCl2, 0.1 mmol/L EDTA, and 15 µg/mL BSA) at 30 °C for 20 min in the presence of the indicated drug concentrations. Reactions were terminated by adding SDS (0.5% final concentration) followed by the addition of two volumes of loading dye (80% formamide, 10 mmol/L sodium hydroxide, 1 mmol/L sodium EDTA, 0.1% xylene cyanol, and 0.1% bromophenol blue). Aliquots of reaction mixtures were subjected to 20% denaturing PAGE. Gels were dried and visualised by using PhosphorImager and Image Quant software (Molecular Dynamics).
Kinases inhibitory assay
The IC50 values of the tested compounds, Gefitinib and TAE226 on different nine kinases (EGFR, FAK, FRK, IGF-1R, BTK, c-Src, VEGFR-1, VEGFR-2, HER-2) were estimated utilising Z`-LYTE® technology, which is based on FRET (Invitrogen/Life Technologies).Ligands were converted to 3D structures and minimised using Avogadro. Ligands were prepared and converted to pdbqt files using PyRx. Protein targets were downloaded from the protein data bank under the codes 1M17 for EGFR and 2JKM for FAK. Co-crystalised ligands were extracted form pdb files and prepared similar to the tested ligands. Docking was done using Autodock Vina in a grid box of 253 Å3 centred on the co-crystalised ligand with exhaustiveness of 16. Visualisation and 3D images were prepared using PyMol.Missing loops in the 3D structures of the protein were constructed using Swiss-Model before starting molecular dynamics steps. All atom molecular dynamics simulations were performed using GROMACS 2020.3 for the selected protein-ligand complexes as reported earlier. In brief, SwissParam server was used for ligands parameterisation while Charmm36 all-atom force field was used to generate topology files for the protein. Ligand coordinates obtained from docking studies were used to build complexes which were boxed in a dodecahedron box and then solvated with TIP3P explicit water. Systems were neutralised by the addition of required number of Na+ or Cl− ions. Systems energy was minimised with a maximum force of 1000 kJ mol−1 nm−1 using steepest descent algorithm. Equilibration using NVT and NPT ensembles for 1 ns each was done afterward then production run was done for 50 ns. Temperature was kept at 300 K using the V-rescale algorithm while pressure was controlled using the Parrinello-Rahman barostat as required. The LINear Constraint Solver (LINCS) algorithm and Particle mesh Ewald (PME) method were used for bond’s length constraints and long-range electrostatics calculations, respectively. Two femtosecond timestep was used for all simulations. Van der Waals distance cut-off (rvdw) was set to 1.2 nm. Trajectories from the production run were used for analysis using trajconv after correction of periodic boundary condition (PBC).Click here for additional data file.
Authors: Marcus D Hanwell; Donald E Curtis; David C Lonie; Tim Vandermeersch; Eva Zurek; Geoffrey R Hutchison Journal: J Cheminform Date: 2012-08-13 Impact factor: 5.514