Search for potentially biased epidermal growth factor receptor (EGFR) inhibitors through pharmacophore modelling, molecular docking, and molecular dynamics (MD) simulation approaches

Epidermal growth factor receptor (EGFR), being one of the most crucial receptor in cancer therapy, has been selected as a potential target for the present study. Ligand-based pharmacophore model (n = 30, R2=0.93 with root mean square deviation = 1.14, ΔCost = 144.27 and configuration cost = 21) was developed and validated with Fischer's randomisation (at 95% confidence), test set (n = 225, R2 pred = 0.81), external data set (n = 13, R2 pred = 0.95) and decoy set (n = 70), further the model has been used to search for novel EGFR inhibitors. The validated model was used for virtual screening of zinc database. A pool of 115,948 candidate molecules was screened through the model. Subsequently, molecules having predicted IC50<0.2 µM were selected for screening through drug-like properties filter. Based on pharmacokinetic profile (ADMET study), Lipinski's rule of five and Veber's rule, 62 molecules were shortlisted for molecular docking. Using consensus docking, five hit molecules were selected, which were further considered for molecular dynamics simulation. Additionally MM-GBSA analysis was carried which showed that affinity of hits towards the receptor of three compound mainly ZINC305, ZINC131796 and ZINC131785 were similar to the standard vanedtinib. The simulation, performed for 100 ns, revealed that two hit molecules, namely ZINC305 and ZINC131785, showing potential interactions at the ligand-binding domain of EGFR protein with good ligand-protein stability. Communicated by Ramaswamy H. Sarma.