Identification of the structural features of quinazoline derivatives as EGFR inhibitors using 3D-QSAR modeling, molecular docking, molecular dynamics simulations and free energy calculations

Abstract

Epidermal growth factor receptor (EGFR) is a promising target for the treatment of different types of malignant tumors. Therefore, a combined molecular modeling study was performed on a series of quinazoline derivatives as EGFR inhibitors. The optimum ligand-based CoMFA and CoMSIA models showed reliable and satisfactory predictability (with R2 cv=0.681, R2 ncv=0.844, R2 pred=0.8702 and R2 cv=0.643, R2 ncv=0.874, R2 pred=0.6423). The derived contour maps provide structural features to improve inhibitory activity. Furthermore, the contour maps, molecular docking, and molecular dynamics (MD) simulations have good consistency, illustrating that the derived models are reliable. In addition, MD simulations and binding free energy calculations were also carried out to understand the conformational fluctuations at the binding pocket of the receptor. The results indicate that hydrogen bond, hydrophobic and electrostatic interactions play significant roles on activity and selectivity. Furthermore, amino acids Val31, Lys50, Thr95, Leu149 and Asp160 are considered as essential residues to participate in the ligand-receptor interactions. Overall, this work would offer reliable theoretical basis for future structural modification, design and synthesis of novel EGFR inhibitors with good potency. Communicated by Ramaswamy H. Sarma