Synthesis, Antitumor Activities and Molecular Modelling of 4-Anilinoquinazoline Derivatives as EGFR-TK Inhibitors

New compounds of 4-anilino-6-substituted quinazoline were designed, synthesized and tested for their EGFR-TK and tumor growth inhibitory activities. The synthesized compounds were appended with amides 6 and 7, dithiocarbamate ester 8a–f or urea/thiourea 9–12 moieties at C-6 of the quinazoline ring to work as extra hydrogen bond acceptors. All the synthesized compounds were effective against EGFR-TK activity, particularly, derivatives 8a, 8f and 9 with IC 50 values of 0.14±0.003, 0.119±0.003, and 0.115±0.002 μM, respectively, showed the best activities. The three compounds were further assayed for their cytotoxicity against MCF-7, H-69, SKOV-3 and LS-174T cell lines. Multikinase enzymes inhibition activity of compound 9 was further screened including VEGFR-2, c-MER, c-MET and Her-2. Compounds 8a, 8f, and 9 were docked into the ATP binding site of EGFR-TK which also had resemblance binding pattern to erlotinib with extra binding mode with Cys-773 at the gatekeeper of the enzyme. Cell cycle analyses of MCF-7 cells treated with 8a and 9 was measured in addition to other related factors such as Bax, Bcl-2, caspase-9, and PARP-1.


Introduction
Cancer is continuing to be a major health problem in developing and developed countries [1]. WHO reported in 2020, that the incidence of population is 18.1 million cancer cases and mortality about 9.6 million deaths occurred [2]. It is regarded as the second leading cause of mortality, after cardiovascular diseases. Among women, breast cancer has stood out as the most frequently diagnosed cancer as well as the second largest cause of cancer-related deaths in women [3].
Epidermal Growth Factor Receptor (EGFR) consists of a family of four different tyrosine kinases (TKs) including EGFR (HER1), HER2, HER3 and HER4 that plays a critical role in the development of multiple cellular functions including cell growth, survival, propagation, differentiation and cell death. EGFR consists of an extracellular, a transmembrane, and an intracellular tyrosine kinase domain [4]. Upon binding of the ligands to the extracellular domain, dimerization occurs with other receptors triggering autophosphorylation of tyrosine residues with the catalytic domain thus activating downstream cell signaling pathways and subsequent signal transduction [4,5]. There are two ways to inhibit EGFR either by blocking ligand binding to the extracellular domain with monoclonal antibodies or by small-molecule inhibitors which interact with the ATP-binding region [6].
Based on the understanding of the competitive binding to the ATP catalytic active site, the rst generation was 4-anilinoquinazoline-based derivatives such as Erlotinib (Tarceva ® ) [7] and ge tinib (Iressa ® ) [8] which are reversible competitors of ATP for binding with the tyrosine kinase. Lapatinib is a reversible dual inhibitor of EGFR and HER2 and is approved for the breast cancer therapy. [9]. Dacomitinib is highly selective second-generation small-molecule inhibitor of EGFR, HER2 and HER4, that speci cally and irreversibly binds to and inhibits these receptors, resulting in inhibition of proliferation and induction of apoptosis in EGFR-expressing tumor cells [10], also Afatinib (Gilotrif ® ) [11] is an example for the second generation, which is irreversibly bound to EGFR-TKIs (Fig. 1).
In 4-anilinoquinazoline-based EGFR inhibitors, the essential parts required for the inhibition are shown in Figure 2. N-1 atom of the quinazoline ring binds with the backbone NH of Met769, in the ATP-binding pocket, via hydrogen bonding. Similarly, a water (HOH-10) molecule-induced hydrogen bonding interaction is noticed between the N-3 atom of the quinazoline ring and Thr766 side chain [4,6,12]. The aniline ring completely occupies the adjacent hydrophobic region of the EGFR kinase domain and the free NH linker of aniline is clearly optimal [6]. Electron-withdrawing lipophilic substituents on the 3-position of the aniline are favorable and only very small groups are preferred, while the substitution of the aniline side chain with larger groups will signi cantly reduce activity as limited bulk tolerance [13]. In some studies new compounds demonstrated that compounds having Br atom at position 3 of the anilino moiety showed better activity than Cl atom at the same position [14]. On the other side of the quinazoline scaffold, introducing electron-donating group with polar properties was favorable at C-6 and/or C-7 positions which lays adjacent to the hydrophilic region to improve the physicochemical properties and to provide a better pharmacokinetic feature [6]. The introduction of large groups at the 5-and 8-positions is not desirable for high potency, because larger groups at these two positions may introduce unfavorable van der Waals interaction with EGFR [15].
In this research, two series of novel 6-substituted-4-anilinoquinazolines were designed to speci cally inhibit EGFR. The rst series contains amide/dithiocarbamate moieties 6-8 and the second series contains urea/thiourea moieties 9-12, both at C-6 of the quinazoline scaffold. By this way, these new polar solubilizing moieties could target the hydrophilic regions of the protein kinase domain's ATP-binding site to develop different selectivity for the designed molecules and to improve their inhibitory effects by improving their physicochemical properties. The synthesized anilinoquinazolines were assessed for their inhibitory ability towards EGFR enzyme and MCF-7, H-69, SKOV-3 and LS-174T cancer cell lines expressing EGFR. Also, exible docking program such as MOE was used to expect the binding a nity and orientation of the target compounds 8a, 8f and 9 at the active region of the EGFR-TK ATP binding site.
Furthermore, ow cytometry was also used to evaluate their effect on the cell cycle of MCF-7 cells.

Experimental
General All chemicals and solvents were ordered from the Aldrich-Sigma and alpha company which were utilized with no more puri cation. Using precoated aluminum sheet silica gel MERCK 60 F245 (TLC technique), the reactions progressed and the spots formed were seen using a UV lamp (366 nm).
Dichloromethane/methanol (10:1), petroleum ether/ethyl acetate (1:1) and (1:3) were used as elution solvents. 1 H and 13 C NMR spectra were obtained using JEOL 500 MHz spectrometer and Bruker Avance III HD FT high resolution 400 MHz spectrometer at Faculty of Pharmacy's Unit of NMR, respectively, Mansoura University and TMS was used as an internal standard and chemical shifts were indicated as ppm and coupling constants (J) were reported in Hz. Mass spectrometry (APCI, or APCI/ASAP) m/z analyses were performed Advion compact mass spectrometer (CMS) NY |USA. General procedure for preparation of compounds (8a-8f) A solution of cyclic secondary amine (0.615 mmol, 1 eq) in DMF (3 mL) and CS 2 (3.077 mmol, 5 eq) was stirred 5 min. then added a nely divided powder of Na 3 PO 4 .H 2 O (0.615 mmol, 1 eq) and stirred at room temperature for 4 hr. Then added to intermediate N-{4-((3-bromophenyl)amino)quinazolin-6-yl}-2chloroacetamide (5) and the mixture was stirred at room temperature until the reaction was complete, then poured into ice to obtain a precipitate that was ltered, washed with water and recrystallized with petroleum ether.     General procedure for preparation of compounds (9)(10)(11)(12) A solution of N4-(3-bromophenyl)quinazoline-4,6-diamine (4) (0.1 gm, 0.317 mmol, 1 eq) in pyridine (3 mL) was reacted with adequate isocyanate or iso(thio)cyanate (0.317 mmol, 1 eq) stirred for 6 hr at room temperature. The reaction mixture was poured into cold dil. HCl and the formed precipitate was ltered, washed with water, dried and recrystallized by diethyl ether. (9) Olive green, yield = 80%, M.P = 294-296 °C. 1   The EGFR kinase assay kit was used to test the newly synthesized compounds for their ability to inhibit EGFR, which was applied in a proper 96-well format and contains a su cient amount of pure recombinant EGFR enzyme, EGFR substrate, ATP and kinase assay buffer for 100 enzyme reactions. Erlotinib was chosen as the reference EGFR inhibitor.
Add 20 μL of diluted EGFR enzyme to the wells labelled "Positive Control" and "Test Inhibitor Control" to start the reaction. Incubate for 40 min at 30 °C. Thaw Kinase-Glo Max reagent. Add 50 μL of Kinase-Glo Max reagent to every well after the 40 min reaction. Cover the plate with aluminum foil and incubate the plate at room temperature for 15 min then measure the luminescence with the microplate reader.
VEGFR-2 inhibition assay VEGFR-2 kinase activity was measured using the detection reagent Kinase-Glo® MAX. 96-well plates with puri ed recombinant VEGFR-2 enzyme, VGFR-2 substrate, ATP and kinase assay buffer were applied to the VEGFR-2 kinase Assay kit. After preparing the positive control, test inhibitor and blank wells, diluted VEGFR-2 enzyme was incorporated into the positive control and also into the test inhibitor, followed by 40 min of incubation at 30 °C. Then, each well was loaded with Kinase-Glo® MAX reagent and incubated for 45 min at room temperature. Cover the plate with aluminum foil and incubate the plate at room temperature for 15 min, then use the microplate reader to measure the luminescence.

c-MER inhibition assay
The c-MER Kinase Assay kit uses ADP-Glo® as a detection reagent to measure c-MER activity for screening and pro ling applications. The c-MER Kinase Assay Kit comes in a convenient 96-well format, with enough puri ed recombinant c-Mer enzyme, c-MER substrate peptide, ATP and kinase assay buffer for 100 enzyme reactions using erlotinib as standard.
Initiate the reaction by adding 10 μL of diluted c-MER enzyme into the "Positive Control" and "Test Inhibitor Control" wells. Incubate at 30 °C for 45 min. Then thaw the ADP-Glo reagent. After the 45 min reaction, add 25 μL of ADP-Glo reagent to every well and incubate the plate at room temperature for 45 min. Add 50 μm of kinase-Detection reagent. Cover the plate with aluminum foil and incubate it at room temperature for another 45 min. The microplate reader is used to measure luminescence. A "Blank" value is subtracted from all readings. Immediately read the sample in a luminometer or microtiter-plate capable of reading luminescence. A "Blank" value is subtracted from all readings.

Her-2 inhibition assay
The Her-2 Kinase Assay Kit uses Kinase-Glo® MAX as a detection reagent to detect Her-2 kinase activity for screening and pro ling applications. Puri ed recombinant Her-2 enzyme, Her-2 substrate, ATP and kinase assay buffer were applied to 96-well plates with Her-2 kinase test kit. After preparing the positive control, test inhibitor and blank wells, diluted Her-2 enzyme was added to both the positive control and test inhibitor, followed by 40 min incubation at 30°C. The Kinase-Glo ® MAX reagent was then applied to each well and incubated for 45 min at room temperature. Cover the plate with aluminum foil and incubate the plate at room temperature for 15 min, then use the microplate reader to measure luminescence. Reassemble each vial of MTT (M-5655) to be used with 3 mL of medium or balanced salt solution without phenol red and serum after removing cultures from the incubator into a laminar ow hood or other sterile work area. 10% of the culture medium volume was added to reconstitute MTT. Return the cultures to the incubator for 2-4 hr depending on the type of cell and the highest density of cells. A 2 hr incubation period is usually su cient, while it can be extended for low cell densities or cells with lower metabolic activity. When doing comparisons, incubation times have to be similar. Remove the cultures from the incubator once the incubation period has ended and dissolve the formazan crystals by applying the same volume of MTT Solubilization solution (M-8910) as the original culture medium. Pipetting up and down [trituration] may be required to completely dissolve the MTT formazan crystals in some cases, particularly in dense cultures. Finally, at a wavelength of 570 nm, the spectrophotometer was measured. At 690 nm, measure the background absorbance of multiwell plates and subtract it from the 570 nm reading.

Cell cycle analysis and apoptosis detection
Flow cytometry was used to scan the cell cycle analysis and investigate apoptosis. The Annexin V Apoptosis Detection Kit is based on the observation that soon after initiating apoptosis, cells translocate the membrane phosphatidylserine (PS) from the inner face of the plasma membrane to the cell surface. Once on the cell surface, PS can be easily detected by staining with a uorescent conjugate of Annexin V, a protein that has a high a nity for PS. Induce apoptosis by the desired method. 1-5 x 10 5 cells was collected by centrifugation and resuspended in 500 μL of 1X Binding Buffer. Add 5 μL of Annexin V-FITC and 5 μL of propidium iodide (PI 50 mg/mL). Incubate at room temperature for 5 min in the dark. Analyze Annexin V-FITC binding by ow cytometry using a FITC signal detector (usually FL1) and PI staining by the phycoerythrin emission signal detector.
Measurement of the effect of compounds 8a and 9 on the levels of Bax and Bcl-2

Bax assay
Except the human Bax-α standard, all reagents should be kept at room temperature for at least 30 min before use. All standards, controls and samples should be tested twice. Determine the number of wells to be utilized using the Assay Layout Sheet, then insert any remaining wells with the desiccant back into the pouch and close the ziploc. Pipette 100 μL of each sample into the appropriate wells and keep the rest of the wells at 4 °C. Tap the plate gently to mix the contents. Seal the plate and incubate at room temperature on a plate shaker for 1 hour at ~500 rpm, then empty the contents of the wells and wash each one with 400 μL of wash solution. For a total of 5 rinses, do the wash 4 more times and empty or aspirate the wells after the nal wash, then rmly tap the plate on a lint-free paper towel to eliminate any residual wash buffer.
Put 100 μL of yellow antibody into each well, except the blank one. Seal the plate and incubate at room temperature for 1 hour at ~500 rpm on a plate shaker. Evacuate the contents of the wells and wash each one with 400 μL of wash solution. For a total of 5 rinses, repeat the wash 4 more times. Fill or aspirate the wells after the last wash, then rmly tap the plate on a lint-free paper towel to eliminate any residual wash buffer. Add 100 μL of blue Conjugate to each well, except the Blank. Close the plate and incubate at room temperature for 30 min at ~500 rpm on a plate shaker. Drain the contents of the wells and wash each one by adding 400 μL of wash solution. For a total of 5 rinses, repeat the wash four more times. Fill or aspirate the wells after the last wash, then rmly tap the plate on a lint-free paper towel to eliminate any residual wash buffer. Pipet 100 μL of substrate solution into each well and on a plate shaker adjusted to 500 rpm, incubate for 30 min at room temperature. Pipet 100 μL stop solution to each well and detect the optical density at 450 nm, usually with a correction between 570 and 590 nm, with the plate reader set towards the blank wells.
If the plate reader cannot be blanked against the blank wells, calculate the result by subtracting the mean optical density of the blank wells from all the readings.

Bcl-2 assay
Wash the microwells twice with approximately 300 μL of wash buffer per well with thorough aspiration of microwell contents between washes. After washing, empty the wells and tap the microwell strips on an absorbent pad or paper towel to drain the remaining wash buffer. Utilize the microwell strips right away once they've been washed, or leave them upside down on damp absorbent paper for no more than 15 min.
To all standard wells and to the blank wells, add 100 μL of sample diluent in duplicate. Prepare doubles of the standard (1:2 dilution) with concentrations ranging from 32 ng/mL to 0.5 ng/mL. Add 100, 80, 20 μL respectively, of sample diluent, in duplicate, to the blank wells. To all wells, having the blank wells, add 50 μL of diluted biotin-conjugate. Cover with a plate cover and incubate at room temperature, on a microplate shaker at 100 rpm if available, for 2 hours. Add 100 μL of diluted Streptavidin-HRP to all wells, including the blank wells. Cover with a plate cover and incubate at room temperature, on a microplate shaker at 100 rpm if available, for 1 hour. Pipette 100 μL of mixed TMB Substrate Solution into all wells, including the blanks. Incubate the microwell strips at room temperature (18 to 25 °C) for about 15 min, if available on a rotator set at 100 rpm. Avoid direct exposure to intense light. The point, at which the substrate reaction is stopped, is often determined by the ELISA reader. Many ELISA readers record absorbance only up to 2.0 O.D. Therefore the color development within individual microwells must be watched by the person running the assay and the substrate reaction stopped before positive wells are no longer properly detectable. Stop the enzyme reaction by quickly pipetting 100 μL of stop solution into each well, including the blank wells. It is important that the stop solution is spread quickly and uniformly throughout the microwells to completely inactivate the enzyme. If the microwell strips are stored at 2-8°C in the dark, results must be read spontaneously after the stop solution is introduced, or within one hour. Read absorbance of each microwell using a spectrophotometer with a using primary wave length 450 nm (optionally 620 nm as the reference wave length; 610 nm to 650 nm) is valid and calculate the result.

Caspases-9 assay
For cell lysis, follow the cell lysate protocol. Determine the number of microwell strips required to test the desired number of samples plus the appropriate number of wells needed for running blanks and standards. Each sample, standard, blank and optional control sample should be assayed in duplicate. In duplicate, pour 100 μL of sample diluent into the blank wells. Fill the sample wells with 50 μL of Sample Diluent. Add 50 μL of each sample in duplicate to the sample wells. Add 50 μL of detection monoclonal antibody to human caspase-9 antibody to all wells. Incubate at room temperature (18 °C to 25 °C) for 2 hours, at 200 rpm, and then add 100 μL of diluted anti-rabbit-IgG-HRP to all wells, and also for the blank wells. Incubate at room temperature (18 °C to 25 °C) for 1 hour, at 200 rpm and proceed immediately to the next step. 100 μL of TMB Substrate Solution should be pipetted to all wells.
Incubate the microwell strips at room temperature (18 °C to 25 °C) for about 10 min. The color development on the plate should be monitored and the substrate reaction stopped before positive wells are no longer properly recordable. Alternatively, the color development can be monitored by the ELISA reader at 620 nm. Read the absorbance of every microwell on a spectrophotometer using 450 nm as the primary wave length (optionally 620 nm as the reference wave length; 610 nm to 650 nm is valid). Blank the plate reader according to the manufacturer's instructions with the use of the blank wells. Determine the absorbance of both the samples and the standards and calculate the result.

PARP-1 Assay
Equilibrate all reagents and materials to ambient room temperature prior to use in the procedure. Optimal results for intra-and inter-assay reproducibility will be obtained when performing incubation steps at 37°C as indicated below. Detect the required number of wells and return any residual unused wells and desiccant to the pouch. Add 100 μL of serially titrated standards, diluted samples or blank into the wells of the PARP1 Microplate. At least two replicates of each standard, sample or blank is recommended.
Cover the plate with the well plate lid and incubate at 37°C for 2 hours. Remove the plate lid and discard the liquid in the wells by rigorously icking it into an acceptable waste receptacle or aspiration. Gently blot any remaining liquid from the wells by tapping inverted on the benchtop onto paper toweling. Do not allow the wells to completely dry out at any time. Add 100 μL of prepared 1X Biotinylated PARP1 detector antibody to each well. Cover with the well-plate lid and incubate at 37°C for 60 min and discard the liquid in the wells by rigorously icking it into an acceptable waste receptacle or aspiration. Gently blot any remaining liquid from the wells by tapping inverted on the benchtop onto paper toweling. Do not allow the wells to completely dry out at any time. Wash the plate 3 times with 1X wash buffer then add 100 μL of prepared 1X Avidin-HRP Conjugate into each well and incubate at 37°C for 60 min. Discard the liquid in the wells by rigorously icking it into an acceptable waste receptacle or aspiration. Gently blot any remaining liquid from the wells by tapping inverted on the benchtop onto paper toweling. Do not allow the wells to completely dry out at any time. Wash the plate 5 times with 1X Wash Buffer as in Step 10.10. 1m17) to compare the structural similarities of our analogues to Erlotinib, which has a higher binding a nity with a lower energy pro le ranging from -8.3437 to -7.1982 kcal/mol in comparison to Erlotinib -7.9161 kcal/mol and showed the best overlay with lead compound. In terms of hydrogen bonding with the Met769 and Thr766 with the quinazoline ring, they had the same binding a nity to the receptor of concern.
Results And Discussion Chemistry In this study, the method for obtaining the desired quinazolines is illustrated in Schemes 1-3. The key  As shown in Table 1, amides 6 and 7 and dithiocarbamate ester 8a-f were evaluated for their inhibitory activities towards EGFR enzymatic assay. The majority of the compounds showed high inhibition in micro molar range (IC 50 from 0.119 to 0.46 μM) using Erlotinib as a reference of (IC 50 = 0.138±0.003 μM).
The structural modi cation of dithiocarbamte derivatives and their EGFR inhibitory action were studied and showed that compounds 8a, 8b and 8c having heterocyclic rings piperidine, morpholine and N-propyl piperazine, respectively, are the most potent when compared to other substituted piperazines which carry N-cyclohexyl and N-phenyl moieties. Interestingly, the N-pyrimidinyl piperazine derivative in 8f proved to be the most potent EGFR inhibitor (IC 50 = 0.119±0.003 μM) among this series. This may give an indication that the polarity is important for the enzyme inhibition particularly when compared to compounds 8d and 8e which contain non polar moieties with IC 50 values equal to 0.47±0.01 and 0.2±0.004, respectively. To clarify the importance of dithiocarbamte moiety, additional two compounds 6 and 7 were synthesized and the results shows that the absence of dithiocarbamte moiety reduces the inhibitory activity of both compounds with IC 50 values equal to 0.44±0.01 μM and 0.28±0.01 μM, respectively, toward EGFR enzyme (entries 8e and 8f versus entries 6 and 7, respectively). On the other hand, the urea and thiourea derivatives of series 2 were evaluated for their EGFR inhibitory activity ( Table 2). The results showed that phenyl urea derivative 12 has the half potency (IC 50 = 0.24±0.005 μM) when compared to erlotinib as a reference (IC 50 = 0.138±0.003 μM). The activity was enhanced to be equal to erlotinib by inserting the polar nitro group at para position of the phenyl moiety 9. The non-polar methyl group at para position diminishes the activity to the half (IC 50 = 0.25±0.005 μM), while the activity was retained when replacing the non-polar para methyl group in the urea compound 10 with the polar para methoxy group in the thiourea compound 11 (IC 50 = 0.16±0.003 μM).

Apoptosis assay
Apoptosis is one of the techniques used to assess the success of cancer treatments [25]. In MCF-7 cells, ow cytometry was used to determine the apoptotic activity of compounds 8a and 9 using propidium iodide (PI) staining and annexin-V-FITC [26]. The studied compound 8a produced an increase in late apoptosis from 0.23 % to 19.4 %, as well as an increase in early apoptosis from 0.36 % to 4.04 %, with necrosis percent 13.59 % versus 0.82 % in the control MCF-7 cell (Fig. 3). For compound 9, there was a rise in the late apoptosis from 0.23% to 11.84% and a rise in the early apoptosis from 0.36% to 2.35% with necrosis percent 7.15% versus 0.82% produced by the control MCF-7 cell (Fig. 3), whereas the reference drug erlotinib has an increase in the late apoptosis from 0.23% to 19.74% and an increase in the early apoptosis from 0.36% to 3.67% with necrosis percent 11.51% versus 0.82% produced by the control MCF-7 cell. It was shown that the late apoptosis percentage induced by 8a and 9 were higher than that of the early phase, implying that the apoptotic cells could be recovered to safe ones. The data collected indicated that compounds 8a and 9 inhibit cancer cell proliferation through apoptosis.

Cell cycle arrest
Compounds 8a and 9 were further assayed of their actions towards cell cycle progression at the following checkpoints G1 (start), S (metaphase), and G2/M [27] in MCF-7 cells against Erlotinib for 24 hr using propidium iodide (PI) and ow cytometer (FCM) for analysis [28]. According to the data obtained, it showed that there were cell accumulations of 37.03% and 7.26% at pre G1 and G2/M phases in MCF-7 cells treated with 8a and 21.52% and 5.87% with 9, respectively, in compared to the untreated cells (1.41% and 9.63%) ( Table 5). On the other hand, G0/G1 population was 46.51% and 57.12% after treatment with compounds 8a and 9, respectively, when compared to the controlled cells 55.42%. The results of the cell cycle analysis is consistent with the data obtained from the apoptosis assay indicate that compounds 8a and 9 are apoptotic inducers. The effect of compounds 8a and 9 on Bax and Bcl-2 Bax and Bcl-2 are main mitochondrial factors which involve in the control of the cell apoptosis [29]. So, the effects of 8a and 9 on the apoptotic factors Bax and Bcl-2 were examined in this study. Compounds 8a and 9 were incubated in the MCF-7 cancer cells for 24 hr at its IC 50 ; 7.27 and 16.55 μM, respectively. In compared to untreated cells, the investigated drugs 8a, 9 and erlotinib enhanced the levels of Bax proteins by 10, 9 and 13 folds, respectively (Table 6). While compounds 8a, 9 and erlotinib enhanced the levels of Bcl-2 proteins approximately by half compared to the untreated cells (Table 6).  (Table 6).

Molecular docking
The interaction among the tested compounds and the EGFR enzyme's ATP pocket was investigated utilizing a computational molecular docking analysis using MOE [31]. The most potent inhibitors 8a, 8f, and 9 were used as ligands in this research, with the structure of EGFR as the docking model (PDB ID: 1m17) and erlotinib's reported binding form in the X-ray crystal structure of the EGFR kinase active site (Fig. 5). The results indicate that compound 8a and 8f form two H-bonding between Met769 and Thr766 by the formation of H-bonding between the N1 and N3 of the quinazoline ring, respectively, as an illustration of dithiocarbamate derivatives. Hydrophobic region binds with Lys721 and carbon disul de linkage to form new bond with EGFR which explains the importance of the presence of the dithiocarbamte moiety (Fig. 5). While, compound 9 as an illustration of urea derivatives strongly linked to the key amino acid Met769 and Thr766 by the formation of H-bonding between the N1 and N3, respectively, and the polar NO 2 group to form covalent bond with Ser696 (Fig. 5). Also, the lead compound, erlotinib, was docked in the active site of EGFR kinase (PDB ID: 1m17), which was shown to obtain the reported binding mode of erlotinib in the 3D structure of the active site of EGFR kinase (Fig. 5), where the N1 atom of quinazoline binds via H-bond with the backbone NH of Met769.
The designed compounds 8a, 8f and 9 have higher binding a nity with a lower energy range between -8.3437 to -7.1982 kcal/mol, compared to Erlotinib -7.9161 kcal/mol. the average distance between atoms of the new pose selected and the original ligand rmsd_re ne 1.9719, 1.2455 and 1.7604, respectively, with a similar binding to the receptor of interest sharing the hydrogen bonding towards Met793 and Thr766 the quinazoline moiety showed with better overlay (Fig. 6).

Conclusion
Two new series of 4-anilino-6-substituted-quinazolines with (thio)urea or dithiocarbamate linkers were synthesized and their antitumor activity against EGFR inhibition was investigated in vitro. Among the tested compounds, 8a and 8f (with dithiocarbamate group) and 9 (with urea group) were the most potent inhibitors against EGFR in comparison to erlotinib as a reference drug. Those three compounds tested against four cell lines MCF-7, H-69, SKOV-3 and LS-174T, where 8a shows the highest cytotoxic activity, among the three tested compounds, 8f exhibits the highest activity against LS-174T and 9 exhibits the  Analysis of cell cycle caused by compounds 8a, 9 and Erlotinib in MCF-7 cells This is a list of supplementary les associated with this preprint. Click to download.