Quantitative structure activity relationship and molecular simulations for the exploration of natural potent VEGFR-2 inhibitors: an in silico anti-angiogenic study

VEGFR-2 has recently become an eye-catching molecular target for the novel therapeutic designs against cancer for its well known role in persuading angiogenesis in tumor cells. The current study set sights on the exploration of novel potent natural compound targeting VEGFR-2 via computational ligand-based modeling and database screening followed by binding pattern analysis, reactivity site prediction and MD simulation studies. The known 53 VEGFR-2 inhibitors (with IC50 ranging from 0.7nM to 9700nM) were headed for development of Ligand based pharmacophore model using 3D QSAR pharmacophore generation module of DS Client. Training set inhibitors (23 compounds) were exploited to create pharmacophore model based on their chemical features. The model was validated through 30 test set inhibitors and exploited further for screening of 62,082 natural compounds from InterBioscreen natural compound database. Screened compounds further went through Drug-Likeliness study, ADMET prediction, Binding pattern analysis, In silico prediction of reactivity sites, Biological activity spectra prediction, pan assay interference compound identification and MD simulation analysis. Out of 5 screened compounds, Compound A and Compound B exhibited highest binding energy judged against the standard drug “Sorafenib”. On further conducting reactivity site prediction, BAS prediction, and pan assay interference compound identification, Compound B exhibited better result which was carried forward for MD simulation study for 50ns. MD simulation results suggested that Compound B exhibited more stable binding to the active site of VEGFR-2 without causing any conformational changes in protein-ligand complex. Thereby, the investigation proposes Compound B to hold potent antiangiogenic potential targeting VEGFR-2. (Figure presented.) Communicated by Ramaswamy H. Sarma.