Exploring the interaction between epidermal growth factor receptor tyrosine kinase and some of the synthesized inhibitors using combination of in-silico and in-vitro cytotoxicity methods

Rezvan Rezaee Nasab , Mahboubeh Mansourian , Farshid Hassanzadeh, Mohsen Shahlaei


Quinazoline derivatives are potent inhibitors of human epidermal growth factor receptor (EGFR) as anticancer agents. In this study, the cytotoxic effects of a new series of synthesized quinazoline derivatives were evaluated using MTT assay against MCF-7 and HT-29 cell lines. Using molecular docking, the binding modes of all compounds were analyzed at the binding site of EGFR. Based on the results, the compounds L1, L2, L4, L5, L6, L7, L10, L15, and L18 may be promising EGFR inhibitors based on docking score and hydrogen bonds. Consistent with the experimental data, Met769 is recognized as a key residue in the binding of potential inhibitors. According to the MTT cytotoxicity assays, Lipinski’s rule of five (RO5), absorption, distribution, metabolism, excretion, and toxicity (ADMET) parameters, and docking studies, three compounds L4, L15, and L10 with IC50 values of 80, 60, and 1 µM against the MCF-7 were selected for further comparative assessments. The dynamics of free EGFR, and selected ligand-EGFR complexes were investigated using molecular dynamics (MD) simulation studies. The results indicated that the three compounds bound to EGFR active site in a stable manner during the simulation through the formation of new hydrogen bonds with Phe699, Leu694, Gly700, Lys721, Met769, Arg817, and Asp831 with the superiority of compound L15. These features can promote future drug candidate designing to produce better derivatives in the search for the anticancer agents.


Keywords: 4-Anilinoquinazoline; EGFR; Molecular docking; Molecular dynamics simulation; 4(3H)-quinazolinones.

Full Text:



Xu H, Yu Y, Marciniak D, Rishi AK, Sarkar FH, Kucuk O, et al. Epidermal growth factor receptor (EGFR)-related protein inhibits multiple members of the EGFR family in colon and breast cancer cells. Mol Cancer Ther. 2005;4(3):435-342.

Nasab RR, Hassanzadeh F, Khodarahmi GA, Mirzaei M, Rostami M. Synthesis, characterization, cytotoxic screening, and density functional theory studies of new derivatives of quinazolin-4(3H)-one Schiff bases. Res Pharm Sci. 2017;12(6):444-455.

Sun J, Wang XY, Lv PC, Zhu HL. Discovery of a series of novel phenylpiperazine derivatives as EGFR TK inhibitors. Sci Rep. 2015;5: 13934-13945.

Abuelizz HA, Marzouk M, Ghabbour H, Al-Salahi R. Synthesis and anticancer activity of new quinazoline derivatives. Saudi Pharm J. 2017;25(7):1047-1054.

Ismail RS, Ismail NS, Abuserii S, El Ella DAA. Recent advances in 4-aminoquinazoline based scaffold derivatives targeting EGFR kinases as anticancer agents. Future Journal of Pharmaceutical Sciences. 2016;2(1):9-19.

Jafari E, Khajouei MR, Hassanzadeh F, Hakimelahi GH, Khodarahmi GA. Quinazolinone and quinazoline derivatives: recent structures with potent antimicrobial and cytotoxic activities. Res Pharm Sci. 2016;11(1):1-14.

Al-Suwaidan IA, Abdel-Aziz AA, Shawer TZ, Ayyad RR, Alanazi AM, El-Morsy AM, et al. Synthesis, antitumor activity and molecular docking study of some novel 3-benzyl-4 (3H) quinazolinone analogues. J Enzyme Inhib Med Chem. 2016;31(1):78-89.

Babu YR, Bhagavanraju M, Reddy GD, Peters GJ, Prasad VV. Design and synthesis of quinazolinone tagged acridones as cytotoxic agents and their effects on EGFR tyrosine kinase. Arch Pharm (Weinheim). 2014;347(9):624-634.

Ahmed MF, Belal A, Youns M. Design, synthesis, molecular modeling and anti-breast cancer activity of novel quinazolin-4-one derivatives linked to thiazolidinone, oxadiazole or pyrazole moieties. Med Chem Res. 2015;24(7):2993-3007.

Hosseinzadeh L, Aliabadi A, Rahnama M, Sadeghi HMM, Khajouei MR. Synthesis and cytotoxic evaluation of some new 3-(2-(2-phenylthiazol-4-yl) ethyl)-quinazolin-4 (3H) one derivatives with potential anticancer effects. Res Pharm Sci. 2017;12(4):290-298.

Nasab RR, Hassanzadeh F, Khodarahmi GA, Rostami M, Mirzaei M, Jahanian-Najafabadi A, et al. Docking study, synthesis and antimicrobial evaluation of some novel 4-anilinoquinazoline derivatives. Res Pharm Sci. 2017;12(5):425-433.

Nasab RR, Mansourian M, Hassanzadeh F. Synthesis, antimicrobial evaluation and docking studies of some novel quinazolinone Schiff base derivatives. Res Pharm Sci. 2018;13(3):213-221.

Jayakanthan M, Jubendradass R, D’Cruz SC, Mathur PP. A use of homology modeling and molecular docking methods: to explore binding mechanisms of nonylphenol and bisphenol A with antioxidant enzymes. Methods Mol Biol. 2015;1268:273-289.

Ali AA, Lee YR, Chen TC, Chen CL, Lee CC, Shiau CY, et al. Novel anthra [1,2-c][1,2,5] thiadiazole-6,11-diones as promising anticancer lead compounds: Biological evaluation, characterization & molecular targets determination. PLoS One. 2016;11(4):e0154278.

Lipinski CA. Lead-and drug-like compounds: the rule-of-five revolution. Drug Discov Today Technol. 2004;1(4):337-341.

Pires DE, Blundell TL, Ascher DB. pkCSM: predicting small-molecule pharmacokinetic and toxicity properties using graph-based signatures. J Med Chem. 2015;58(9):4066-4072.

Daina A, Michielin O, Zoete V. SwissADME: a free web tool to evaluate pharmacokinetics, drug-likeness and medicinal chemistry friendliness of small molecules. Sci Rep. 2017;7:42717-42729.

Drwal MN, Banerjee P, Dunkel M, Wettig M, Preissner R. ProTox: a web server for the in silico prediction of rodent oral toxicity. Nucleic Acids Res. 2014;42(Web Server issue):W53-W58.

Morris GM, Huey R, Lindstrom W, Sanner MF, Belew RK, Goodsell DS, et al. AutoDock4 and AutoDockTools4: Automated docking with selective

receptor flexibility. J Comput Chem. 2009;30(16):2785-2791.

Berman HM, Battistuz T, Bhat TN, Bluhm WF, Bourne PE, Burkhardt K, et al. The protein data bank. Acta Crystallogr D Biol Crystallogr. 2002;58(1):899-907.

Stamos J, Sliwkowski MX, Eigenbrot C. Structure of the epidermal growth factor receptor kinase domain alone and in complex with a 4-anilinoquinazoline inhibitor. J Biol Chem. 2002;277(48):46265-4672.

Mansourian M, Fassihi A, Saghaie L, Madadkar-Sobhani A, Mahnam K, Abbasi M. QSAR and docking analysis of A2B adenosine receptor antagonists based on non-xanthine scaffold. Med Chem Res. 2015;24(1):394-407.

Coleman WF, Arumainayagam CR. HyperChem 5 (by Hypercube, Inc.). J Chem Educ. 1998;75(4):416.

Mansourian M, Madadkar-Sobhani A, Mahnam K, Fassihi A, Saghaie L. Characterization of adenosine receptor in its native environment: insights from molecular dynamics simulations of palmitoylated/glycosylated, membrane-integrated human A(2B) adenosine receptor. J Mol Model. 2012;18(9):4309-4324.

SchuÈttelkopf AW, Van Aalten DM. PRODRG: a tool for high-throughput crystallography of protein-ligand complexes. Acta Crystallogr D Biol Crystallogr. 2004;60(Pt 8):1355-1363.

Hess B, Kutzner C, van der Spoel D, Lindahl E. GROMACS 4: algorithms for highly efficient, load-balanced, and scalable molecular simulation. J Chem Theory Comput. 2008;4(3):435-447.

Sadeghian-Rizi S, Khodarahmi GA, Sakhteman A, Jahanian-Najafabadi A, Rostami M, Mirzaei M, et al. Biological evaluation, docking and molecular dynamic simulation of some novel diaryl urea derivatives bearing quinoxalindione moiety. Res Pharm Sci. 2017;12(6):500-509

Berendsen HJC, Postma JPM, van Gunsteren WF, Hermans J. Interaction Models for Water in Relation to Protein Hydration. in: Intermolecular Forces. Pullman B, editor. Vol 14. Dordrecht: Springer; 1981. pp. 331-342.

Berendsen HJC, Postma JPM, van Gunsteren WF, DiNola A, Haak JR. Molecular dynamics with coupling to an external bath. J Chem Phys. 1984;81(8):3684-3690.

Hess B, Bekker H, Berendsen HJC, Fraaije JGEM. LINCS: A linear constraint solver for molecular simulations. J Comput Chem. 1998;18(12):1463-1472.

Darden TA, York DM, Pedersen LG. Particle mesh Ewald: An N log (N) method for Ewald sums in large systems. J Chem Phys. 1993;98(12):10089-10092.

Mansourian M, Mahnam K, Madadkar-Sobhani A, Fassihi A, Saghaie L. Insights into the human A1 adenosine receptor from molecular dynamics simulation: structural study in the presence of lipid membrane. Med Chem Res. 2015;24(10):3645-3659.

Humphrey W, Dalke A, Schulten K. VMD: visual molecular dynamics. J Mol Graph. 1996;14(1):33-38.

Wallace AC, Laskowski RA, Thornton JM. LIGPLOT: a program to generate schematic diagrams of protein-ligand interactions. Protein Eng. 1995;8(2):127-134.

Songtawee N, Gleeson MP, Choowongkomon K. Computational study of EGFR inhibition: molecular dynamics studies on the active and inactive protein conformations. J Mol Model. 2013;19(2):497-509.


  • There are currently no refbacks.

Creative Commons Attribution-NonCommercial 3.0

This work is licensed under a Creative Commons Attribution-NonCommercial 3.0 Unported License which allows users to read, copy, distribute and make derivative works for non-commercial purposes from the material, as long as the author of the original work is cited properly.