Manijeh Nematpour1, Elham Rezaee1, Maryam Nazari1, Omid Hosseini2, Sayyed Abbas Tabatabai1. 1. Department of Pharmaceutical Chemistry, School of Pharmacy, Shahid Beheshti University of Medical Sciences, Tehran, Iran. 2. Central Research Labretories, Shahid Beheshti University of Medical Sciences, Tehran, Iran.
Abstract
Impaired cell cycle regulation and disturbance in signal transduction pathway are two major causes of a condition defined as cancer, one of the significant reasons for mortality worldwide. Epidermal growth factor receptor tyrosine kinase inhibitors (EGFR-TKIs) have been commonly used as anticancer agents, and the majority of this medications possess quinazoline moiety as a heteroaromatic core. In this study, two novel series of EGFR-TKIs containing quinazolinone core were designed and synthesized. Most compounds showed reasonable inhibitory activity against EGFR-TK compared to that of erlotinib, a reversible inhibitor of this enzyme. Compound 8b, 2-((2-chlorobenzyl)amino)-6-phenoxyquinazolin-4(1H)-one, with an IC50 value of 1.37 nM exhibited the highest potency. Molecular docking study of compound 8b showed that it had the same direction of erlotinib and formed proper hydrogen bonds and hydrophobic interactions with the important amino acid residues of the active site. Based on in-silico calculations of ADME properties, our novel compounds have the potential to be orally active agents.
Impaired cell cycle regulation and disturbance in signal transduction pathway are two major causes of a condition defined as cancer, one of the significant reasons for mortality worldwide. Epidermal growth factor receptor tyrosine kinase inhibitors (EGFR-TKIs) have been commonly used as anticancer agents, and the majority of this medications possess quinazoline moiety as a heteroaromatic core. In this study, two novel series of EGFR-TKIs containing quinazolinone core were designed and synthesized. Most compounds showed reasonable inhibitory activity against EGFR-TK compared to that of erlotinib, a reversible inhibitor of this enzyme. Compound 8b, 2-((2-chlorobenzyl)amino)-6-phenoxyquinazolin-4(1H)-one, with an IC50 value of 1.37 nM exhibited the highest potency. Molecular docking study of compound 8b showed that it had the same direction of erlotinib and formed proper hydrogen bonds and hydrophobic interactions with the important amino acid residues of the active site. Based on in-silico calculations of ADME properties, our novel compounds have the potential to be orally active agents.
Nowadays, cancer is one of the leading causes of death worldwide, and discovering new anticancer agents with high selectivity and low adverse effects is a significant challenge encountered by researchers (1, 2). Receptor protein tyrosine kinases (RPTKs), also called phosphotyrosine kinases, have a considerable role in the signal transduction pathway and are responsible for regulating various functions in normal cells, including cell growth, proliferation, metabolism, differentiation, and apoptosis (3, 4). The alteration of this enzyme via mutations would lead to dysregulated cell cycle and uncontrolled cell duplication, defined as cancer (5-7). The FDA approved 47 kinase inhibitors over the last 30 years that have developed as anticancer agents (8).The human epidermal growth factor receptors (EGFRs) belong to a class of RPTKs, including four members: EGFR/ERBB1/HER1, ERBB2/HER2, ERBB3/HER3 and ERBB4/HER4 (9). This type of receptors can perform a key role in signal transduction for adjusting cell-growth and developing many types of solid tumors, including head and neck, lung, breast, bladder, prostate, and kidney cancers (10-12).As a relatively lethal disease, lung cancer’s manifestations are categorized into two main groups: small-cell lung carcinoma (SCLC) and non-small-cell lung carcinoma (NSCLC) (13, 14). This invasive and metastatic cancer type and stage would guide the specialist to select proper treatment, including surgery, radiation, and chemotherapy. Erlotinib (Tarceva®), an EGFR tyrosine kinase inhibitor (EGFR-TKI), has been approved as the first-line therapy for EGFR mutated progressive or metastatic NSCLC (15). Quinazolines are ubiquitous heterocyclic scaffolds found in various pharmaceuticals and natural products (16). The quinazoline-containing derivatives have been contributed to the variety of EGFR-TKIs. There have been many novel quinazoline derivatives possessing high inhibitory activity compared to FDA approved drugs in the last few years. Some novel and potent EGFR-TKIs with quinazoline moiety were illustrated in Figure 1. Compound 1 is a semicarbazone fused quinazoline group derivative with an IC50 value of 0.05 nM (17). The concept behind the designing compound 2 (IC50 = 0.76 nM) was changing the size of alkoxyalkane chains at 7-position of quinazoline core (18). Compound 3 (IC50 = 3 nM) bearing 2,3-dihydro-indolequinoline core demonstrated strong potency as well as excellent antiproliferative activity (19). Compound 4 (IC50 = 1.79 nM) was designed as dual EGFR and α-glucosidase inhibitor possessing quinazoline moiety hybridized with 1-deoxynojirimycin segment (20).
Figure 1.
The structures of some potent EGFR-TKIs containing quinazoline moiety.
Based on studies, drug resistance in cancer and ineffective chemotherapy mainly stem from EGFR mutation (replacement of the threonine 790 with methionine 790), leading to the progression of cancer stages. Although this factor could be conquered through developing inhibitors binding covalently to the ATP binding site of EGFR (21), these inhibitors suffer from weak selectivity due to their affinity to other kinases (22). The alternative strategy introduces reversible EGFR-TKIs, which have been desirable recently, and researchers are working on them.Considering the necessary features of erlotinib as a reversible inhibitor and related structures, we designed two series of novel quinazolinone derivatives as EGFR-TKIs illustrated in Figure 2. According to the essential pharmacophoric features of erlotinib (3), a hydrophobic head and tail on a heteroaromatic system and a hydrogen acceptor spacer that bridges between head and heteroaromatic system are required to locate a compound into the active site of EGFR tyrosine kinase (EGFR-TK). Therefore, the modifications of our novel series of quinazolinone derivatives (compounds 8 and 9) included the introduction of various N-aryl and -alkyl groups (hydrophobic head) with different substituents at the 2-position of quinazolinone moiety (bioisoster of quinazoline group) and the amino group served as a spacer. Besides, hydrogen at the 6-position was replaced with phenoxy group in 8a-o derivatives to resemble a hydrophobic tail (Figure 2). Finally, the activities of the novel compounds toward recombinant EGFR-TK were evaluated.
Figure 2.
The design strategy of our novel compounds as EGFR-TKIs based on pharmacophoric features of erlotinib
2. Methods
2.1. Chemistry
All reagents and solvents were purchased from Aldrich or Merck Company and utilized without any purification. IR spectra were measured on a Perkin Elmer 834 spectrometer. The absorptions were reported on the wavenumber (cm-1) scale, in the range 400 - 4000 cm-1. All the 1H-NMR and 13C-NMR spectra were recorded on a Bruker Avance II spectrophotometer using CDCl3 as a solvent and tetramethylsilane as an internal standard. Finnigan-MAT-8430EI-MS was applied to obtain mass spectra; at 70 eV; in m/z (rel.%). Melting points were obtained using Electrothermal 9100 apparatus and were uncorrected.
2.2. General Procedure for the Synthesis of Compounds 8a-o
A mixture of 2-amino-5-fluorobenzoic acid (5.0 mmol), phenol (5.2 mmol), CuI (0.1 mmol), NaH (2.0 mmol) in DMF (3 mL) was refluxed at 110°C for 2 h. After the addition of ethyl chloroformate (6.0 mmol) to the medium, Et3N (2.0 mmol) was added dropwise over 10 min at room temperature and refluxed at 110°C for 3 h. Then, the mixture was cooled, diluted with CH2Cl2 (2 mL) and aqueous NH4Cl solution (3 mL), and stirred for 30 min. The aqueous layer was extracted with CH2Cl2, dried over sodium sulfate anhydrous, and concentrated under reduced pressure. The precipitate was recrystallized from EtOH (30 mL) to afford 6-phenoxyisatoic anhydride.To a mixture of appropriate aniline or benzylamine (3.0 mmol) and trichloroacetonitrile (3.5 mmol), a mixture of synthesized 6-phenoxyisatoic anhydride (3.5 mmol) and CuO (0.1 mmol) in DMF (3 mL) was slowly added and refluxed at 110°C for 4 h. The end of the reaction was controlled by TLC (eluent: AcOEt/ hexane 1:3). Then, the mixture was diluted with CH2Cl2 (2 mL) and aqueous NH4Cl solution (3 mL) and stirred for 30 min. The aqueous layer was extracted with CH2Cl2 dried over sodium sulfate anhydrous and concentrated under reduced pressure. The precipitate was washed with diethyl ether to afford final product 8.
2.3. General Procedure for the Synthesis of Compounds 9a-o
A mixture of 2-amino-benzoic acid (5.0 mmol) and ethyl chloroformate (6.0 mmol) in DMF (3 ml) was refluxed at 110°C for 3 h. Et3N (2.0 mmol) was added to the reaction dropwise over 10 min at room temperature and refluxed at 110 °C for 3 h. Then, the mixture was cooled, diluted with CH2Cl2 (2 mL) and aqueous NH4Cl solution (3 mL), and stirred for 30 min. The aqueous layer was extracted with CH2Cl2, dried over sodium sulfate anhydrous, and concentrated under reduced pressure. The precipitate was recrystallized from EtOH (30 mL) to afford isatoic anhydride.To a mixture of appropriate aniline or benzylamine (3.0 mmol) and trichloroacetonitrile (3.5 mmol), a mixture of synthesized isatoic anhydride (3.5 mmol) and CuO (0.1 mmol) in DMF (3 mL) was slowly added and refluxed at 110°C for 4 h. TLC (eluent: AcOEt/ hexane 1:3) was applied to determine the end of the reaction. The mixture was diluted with CH2Cl2 (2 mL) and aqueous NH4Cl solution (3 mL) and stirred for 30 min. The aqueous layer was extracted with CH2Cl2 dried over sodium sulfate anhydrous and concentrated under reduced pressure. The precipitate was washed with diethyl ether to afford final product 9.
The biological activities of the synthesized compounds, 8a-o and 9a-o, were carried out using EGFR-TK assay kit (BPS biosciences) and kinase-Glo luminescent assay (Promega). In order to start the enzymatic reaction, EGFRT790M and their substrates were incubated with the synthesized compounds in kinase assay buffer for 60 min at room temperature. Then detection reagent (ADP-GloTM reagent) was added to each well, followed by incubation at room temperature for 40 min. With the addition of kinase detection reagent, the mixture was incubated for further 30 min. The IC50 values were determined after the luminescence record, and all samples and controls were tested in triplicate.
2.5. Docking Study
Crystallographic structure of EGFR tyrosine kinase was retrieved from Protein Data Bank (PDB code: 4G5P, resolution 3.17 Å). The docking analysis was performed using AutoDock Tools version 1.5.6rc3 (http://mgltools.scripps.edu) software to evaluate the binding mode of the compound 8b against EGFRT790M. At first, the original ligand and water molecules were removed, and polar hydrogen atoms and Kollman united partial atom charges were added onto the protein. The 2D structure of the synthesized compound and erlotinib were sketched using ChemBioDraw Ultra 12.0 and optimized by MM+ force field using HyperChem8 (http://www.hyper.com). Each docking system was performed by 100 runs of the Autodock using the Lamarckian genetic algorithm (LGA). The lowest docking-energy conformation of the highest populated cluster was considered the most stable orientation and selected for analysis. Finally, a cluster analysis was performed on the docking results using a root mean square deviation (RMSD) of 0.5 Å. Graphic manipulations and visualizations were done by Pymol software version 1.5.0.1 (http://pymol.findmysoft.com).
2.6. ADME Properties
The ADME properties of the synthesized compounds were predicted using the Molinspiration online property calculation toolkit (Molinspiration, 2018). Topological polar surface area (TPSA), number of rotatable bonds (n-ROTB), molecular volume (MV), molecular weight (MW), the logarithm of the partition coefficient (miLog P), number of hydrogen bond acceptors (n-ON), number of hydrogen bond donors (n-OHNH), and Lipinski's rule of five were calculated (23). Additionally, intestinal absorption (% ABS) was calculated according to the following equation: ABS% = 109 − (0.345 × TPSA) (24).
3. Results and Discussion
3.1. Chemistry
Figure 3 displayed the general synthetic pathway of the novel compounds. 2-Amino-5-fluorobenzoic acid 1 was reacted with phenol 2 through a nucleophilic aromatic substitution using sodium hydride as the base under reflux conditions to produce compound 3. Isatoic anhydride derivatives (5, 6) were achieved by reacting the compounds 3 and 2-aminobenzoic acid 4 with ethyl chloroformate and triethylamine. The final products 8(a-o) and 9(a-o) were afforded by the treatment of the above isatoic anhydrides with corresponding aniline or benzylamine in the presence of trichloroacetonitrile and catalytic amount of CuO under reflux condition (25).
The inhibitory activities of the novel compounds were assessed using EGFRT790M/ADP-GloTM kinase assay system. Erlotinib, a reversible EGFR-TKI, was used as a standard compound. As shown in Figure 4, most synthesized compounds provided nanomolar range inhibition against EGFR-TK. Some derivatives of compound 8 (i.e., compounds 8a, 8b, 8c, 8d and 8f) exhibited higher inhibitory activities than the other compounds, compared to erlotinib (IC50 = 1.5 nM). The most potent compound was 8b, with an IC50 value of 1.37 nM. By comparing the IC50 values of compounds 8b and 9b, it can be concluded that introducing a phenoxy group at the 6-position could be beneficial for the activity owing to the proper fitting of the compound into a hydrophobic pocket. According to compounds 8c, 8f, and 9c, the presence of the electron-donating substituents in para position of benzyl ring in phenoxy series were more desirable than electron-withdrawing substituents; however, this finding cannot be extended to ortho position as compound 8b exhibited the highest activity with an electron withdrawing substituent. Substitution of fluorine in para position of either benzyl or phenyl ring (i.e., compounds 8m and 9m) was not favorable for the activity since fluorine, a hydrophilic atom could be surrounded by water molecules, leading to a difficult location into the hydrophobic pocket of the enzyme.
Figure 4.
Inhibitory activities of novel quinazolin-4(1H)-one derivatives against EGFR tyrosine kinase (NA: no activity).
Moreover, compounds 8j and 9j possessing 3-nitro substituent on phenyl ring were potent analogs in both series. It could be explained that the nitro group could afford a hydrogen bond with polar amino acids in the active site. Generally, due to close variations of IC50s, the judgment about the activity trend was a challenge we were encountered, and we did not examine electron-donating substituents in phenyl bearing series, which is a limitation for our study.
3.3. Docking Study
A docking study was performed to identify the binding mode of compound 8b and EGFR-TK active site. As shown in Figure 5, docking study revealed that this structure had a good affinity to the active site of the enzyme and fitted properly. The amide group of quinazolinone moiety and nitrogen in the hydrophobic head had suitable distances from Met793 for hydrogen bonding interaction. Additional hydrogen bonds formed between the oxygen of the phenoxy group and residues of Met790 and Thr854. As compound 8b was surrounded by lipophilic amino acids, hydrophobic interactions were observed between 8b and Lys745, Val726, Leu718, Ile744. Furthermore, the superimposition of compound 8b and erlotinib demonstrated that these two structures were located in the same direction in the active site as the nitrogen atoms at the 1-position of erlotinib and 3-position of 8b were completely aligned and formed hydrogen bonds with Met793.
Figure 5.
(A) Compound 8b (green sticks) in the active site of EGFRT790M receptor (amino acid residues in blue PDB: 4G5P). The amide group of quinazolinone moiety and nitrogen in the hydrophobic head had suitable distances from Met793 for hydrogen bonding interaction. Additional hydrogen bond formed between oxygen of phenoxy group and Met790 and Thr854. The hydrophobic interactions observed between 8b and Lys745, Val726, Leu718, and Ile744. (B) The superimposition of compound 8b (green sticks) and erlotinib (salmon pink sticks) in the active site of EGFRT790M receptor.
3.4. ADME Properties
To predict oral bioavailability and druglikeness of our novel compounds, ADME calculations were performed using molinspiration calculator (Table 1). According to the Lipinski's rule of five (miLog P (octanol–water partition coefficient) ≤ 5, molecular weight ≤ 500, number of hydrogen bond acceptors ≤ 10, and number of hydrogen bond donors ≤ 5), an orally active drug candidate has no more than one violation of the criteria (26). All the synthesized compounds exhibited an excellent ABS% ranging from 70.07% to 89.06%. They could be orally active agents through fulfilling the criteria.
Table 1.
Pharmacokinetic Parameters for Appropriate Oral Bioavailability of the Synthesized Compounds 8a-o and 9a-o
Number
%ABS
TPSA (A2)
n ROTB
MV
MW
miLogP
n-ON Acceptors
n-OHNH Donors
Lipinski's Violations
8a
85.88
67.02
5
308.85
343.39
4.71
5
2
0
8b
85.88
67.02
5
322.39
377.83
5.34
5
2
1
8c
82.69
76.25
6
334.40
373.41
4.77
6
2
0
8d
85.88
67.02
5
322.39
377.83
5.39
5
2
1
8e
85.88
67.02
5
313.79
361.38
4.88
5
2
0
8f
85.88
67.02
5
325.42
357.41
5.16
5
2
1
8g
85.88
67.02
4
292.05
329.36
5.03
5
2
1
8h
70.07
112.84
5
315.39
374.36
4.94
8
2
0
8i
70.07
112.84
5
315.39
374.36
4.99
8
2
0
8j
70.07
112.84
5
315.39
374.36
4.96
8
2
0
8k
85.88
67.02
4
309.94
408.25
5.84
5
2
1
8l
85.88
67.02
4
319.12
398.25
6.31
5
2
1
8m
85.88
67.02
4
296.98
347.35
5.19
5
2
1
8n
85.88
67.02
4
305.59
363.80
5.71
5
2
1
8o
85.88
67.02
4
305.59
363.80
5.68
5
2
1
9a
89.06
57.78
3
228.46
251.29
2.98
4
2
0
9b
89.06
57.78
3
242.00
285.73
3.61
4
2
0
9c
85.88
67.02
4
254.01
281.31
3.04
5
2
0
9d
89.06
57.78
3
242.00
285.73
3.66
4
2
0
9e
89.06
57.78
3
233.39
269.28
3.15
4
2
0
9f
89.06
57.78
3
245.02
265.32
3.43
4
2
0
9g
89.06
57.78
2
211.66
237.26
3.30
4
2
0
9h
73.25
103.61
3
234.99
282.26
3.21
7
2
0
9i
73.25
103.61
3
234.99
282.26
3.26
7
2
0
9j
73.25
103.61
3
234.99
282.26
3.23
7
2
0
9k
89.06
57.78
2
229.54
316.16
4.11
4
2
0
9l
89.06
57.78
2
238.73
306.15
4.58
4
2
0
9m
89.06
57.78
2
216.59
255.25
3.46
4
2
0
9n
89.06
57.78
2
225.19
271.71
3.98
4
2
0
9o
89.06
57.78
2
225.19
271.71
3.95
4
2
0
Erlotinib
83.31
74.45
10
362.06
393.44
2.79
7
1
0
Abbreviations: ABS%, percentage absorption; TPSA, topological polar surface area; n-ROTB, number of rotatable bonds; MV, molecular volume; MW, molecular weight; miLog P, logarithm of partition coefficient of the compound between n-octanol and water; n-ON acceptors, number of hydrogen bond acceptors; n-OHNH donors, number of hydrogen bonds donors.
Abbreviations: ABS%, percentage absorption; TPSA, topological polar surface area; n-ROTB, number of rotatable bonds; MV, molecular volume; MW, molecular weight; miLog P, logarithm of partition coefficient of the compound between n-octanol and water; n-ON acceptors, number of hydrogen bond acceptors; n-OHNH donors, number of hydrogen bonds donors.
4. Conclusions
In summary, we designed and synthesized thirty qinazolinone derivatives as EGFR-TKIs. All the compounds were evaluated for their inhibitory activities against EGFRT790M. Most of our novel compounds provided nanomolar range inhibition. Compounds 8a, 8b, 8c, 8d, and 8f exhibited more potent inhibitory activities than erlotinib. Remarkably, Compound 8b with 2-chlorobenzyl group in phenoxy derivatives showed the exemplary inhibitory activity with an IC50 value of 0.37 nM. The docking study of 8b demonstrated that it had a good affinity to the active site of the enzyme and fitted properly. According to our findings, attaching a phenoxy group at the 6-position could be beneficial for the activity due to the proper disposition of the compound into the hydrophobic pocket of the enzyme. Based on ADME prediction analysis, the synthesized compounds possessed excellent druglikeness properties to be formulated for the oral route of administration. The extension of this work would provide valuable facts for the further design of EGFR-TKIs.
Authors: Yuan H Zhao; Michael H Abraham; Joelle Le; Anne Hersey; Chris N Luscombe; Gordon Beck; Brad Sherborne; Ian Cooper Journal: Pharm Res Date: 2002-10 Impact factor: 4.200
Authors: Wagdy M Eldehna; Mahmoud F Abo-Ashour; Alessio Nocentini; Paola Gratteri; Ibrahim H Eissa; Mohamed Fares; Omnia E Ismael; Hazem A Ghabbour; Mahmoud M Elaasser; Hatem A Abdel-Aziz; Claudiu T Supuran Journal: Eur J Med Chem Date: 2017-08-01 Impact factor: 6.514
Authors: Eman A Fayed; Sally I Eissa; Ashraf H Bayoumi; Nirvana A Gohar; Ahmed B M Mehany; Yousry A Ammar Journal: Mol Divers Date: 2018-08-11 Impact factor: 2.943
Figure 1.
Figure 2.
Figure 3.
Figure 4.
Figure 5.