Computational Investigation of 1, 3, 4 Oxadiazole Derivatives as Lead Inhibitors of VEGFR 2 in Comparison with EGFR: Density Functional Theory, Molecular Docking and Molecular Dynamics Simulation Studies

Muhammad Sajjad Bilal,Syeda Abida Ejaz,Seema Zargar,Naveed Akhtar,Tanveer A Wani,Naheed Riaz,Adullahi Tunde Aborode,Farhan Siddique,Nojood Altwaijry,Hamad M Alkahtani,Haruna Isiyaku Umar
DOI: https://doi.org/10.3390/biom12111612
IF: 6.064
2022-11-01
Biomolecules
Abstract:Vascular endothelial growth factor (VEGF) is an angiogenic factor involved in tumor growth and metastasis. Gremlin has been proposed as a novel therapeutic pathway for the treatment of renal inflammatory diseases, acting via VEGFR 2 receptor. To date, most FDA-approved tyrosine kinase (TK) inhibitors have been reported as dual inhibitors of EGFR and VEGFR 2. The aim of the present study was to find the potent and selective inhibitor of VEGFR 2 specifically for the treatment of renal cancer. Fourteen previously identified anti-inflammatory compounds i.e., 1, 3, 4 oxadiazoles derivatives by our own group were selected for their anti-cancer potential, targeting the tyrosine kinase (TK) domain of VEGFR2 and EGFR. A detailed virtual screening-based study was designed viz density functional theory (DFT) study to find the compounds' stability and reactivity, molecular docking for estimating binding affinity, SeeSAR analysis and molecular dynamic simulations to confirm protein ligand complex stability and ADMET properties to find the pharmacokinetic profile of all compounds. The DFT results suggested that among all the derivatives, the 7g, 7j, and 7l were chemically reactive and stable derivatives. The optimized structures obtained from the DFTs were further selected for molecular docking, and the results suggested that 7g, 7j and 7l derivatives as the best inhibitors of VEGFR 2 with binding energy values -46.32, -48.89 and -45.01 kJ/mol. The Estimated inhibition constant (IC50) of hit compound 7j (0.009 µM) and simulation studies of its complexes confirms its high potency and best inhibitor of VEGFR2. All the derivatives were also docked with EGFR, where they showed weak binding energies and poor interactions, important compound 7g, 7j and 7i exhibited binding energy of -31.01, -33.23 and -34.19 kJ/mol respectively. Furthermore, the anticancer potential of the derivatives was confirmed by cell viability (MTT) assay using breast cancer and cervical cancer cell lines. At the end, the results of ADMET studies confirmed these derivatives as drug like candidates. Conclusively, the current study suggested substituted oxadiazoles as the potential anticancer compounds which exhibited more selectivity towards VEGFR2 in comparison to EGFR. Therefore, the identified lead molecules can be used for the synthesis of more potent derivatives of VEGFR2, along with extensive in vitro and in vivo experiments, that can be used to treat various cancers, especially renal cancers, and to prevent angiogenesis due to aberrant expression of VEGFR2.
What problem does this paper attempt to address?