Molecular dynamics-based insight of VEGFR-2 kinase domain: a combined study of pharmacophore modeling and molecular docking and dynamics

Background Inhibition of vascular endothelial growth factor receptor 2 (VEGFR-2) tyrosine kinase by small molecules has become a promising target in the treatment of cancer. Objective In this study, we approached pharmacophore modeling coupled with a structure-based virtual screening workflow to identify the potent inhibitors. Methods The top selected hit compounds have been rescored using the MM/GBSA approach. To understand the molecular reactivity, electronic properties, and stability of those inhibitors, we have employed density functional theory and molecular dynamics. Following that, the best 21 hit compounds have been further post-processed with a Quantum ligand partial charge-based rescoring process and further validated by implementing molecular dynamics simulation. Results The ten hit compounds have been hypothesized and considered as potent inhibitors of VEGFR-2 tyrosine kinase. This study also signifies the contribution of QM-based ligand partial charge, which is more accurate in predicting reliable free binding energy and filtering large ligand libraries to hit optimization, rather than assigning those of the force field-based method. From the binding pattern analysis of all the complexes, amino acids, such as Glu885, Cys919, Cys1045, Thr916, Thr919, and Asp1046, were found to have comprehensive interaction with the hit compounds. Conclusion Hence, this could prove to be useful as a potential inhibition site of the VEGFR-2 tyrosine kinase domain for future researchers. Moreover, this study also emphasizes the conformational changes upon ATP binding, based on either the receptor’s rigidity or flexibility.