2D-QSAR driven design, molecular docking, molecular dynamics simulation and MM/GBSA studies on quinazoline derivatives for development of VEGFR-2 inhibitors

Abstract Various in silico approaches were utilized to design quinazoline derivatives as anticancer agents targeting VEGFR-2. Six 2D-QSAR models were generated via Monte Carlo optimization method of CORALSEA software. The models generated were based on hybrid optimal descriptors including Graphs and...

Full description

Saved in:
Bibliographic Details
Main Authors: Mayank Kashyap, Saurabh Gupta, Yogita Bansal, Gulshan Bansal
Format: Article
Language:English
Published: Springer 2025-05-01
Series:Discover Chemistry
Subjects:
Online Access:https://doi.org/10.1007/s44371-025-00173-4
Tags: Add Tag
No Tags, Be the first to tag this record!
_version_ 1850154586063503360
author Mayank Kashyap
Saurabh Gupta
Yogita Bansal
Gulshan Bansal
author_facet Mayank Kashyap
Saurabh Gupta
Yogita Bansal
Gulshan Bansal
author_sort Mayank Kashyap
collection DOAJ
description Abstract Various in silico approaches were utilized to design quinazoline derivatives as anticancer agents targeting VEGFR-2. Six 2D-QSAR models were generated via Monte Carlo optimization method of CORALSEA software. The models generated were based on hybrid optimal descriptors including Graphs and SMILES. All the models were observed to be robust via internal and external validation parameters. Model M4 of split 3 had the highest average correlation coefficient (r 2  = 0.7209), average $${R}_{m}^{2}$$ R m 2 value (0.6256) and lowest error (MAE = 0.465). Structural features characterised as promoters and hinderers were extracted from the M4 model and their mechanistic interpretation was employed to improve the pIC50 of least potent compound SM33 from the dataset. Removal of hinderers and addition of promoters resulted in a series of compounds, DDSL6a-DDSL6l among which DDSL6j was the most active compound with a pIC50 of 8.07. All the compounds were subsequently subjected to molecular docking against PDB: 4ASD, with DDSL6j having better LF rank score (− 14.032) and LF dG score (− 12.318). MD simulation of DDSL6j, sorafenib and SM33 carried out at 200 ns revealed sustained interactions of these compounds with all the key amino acid residues of the protein. Improved potency of compound DDSL6j was attributed to presence of additional stable H-Bond interaction. MM/GBSA calculation also indicated stable complex with binding free energy of 0.28 kcal/mol. All the designed molecules were also found to be possess good ADMET profile. These findings indicate that the quinazoline nucleus can be investigated and refined to create novel VEGFR-2 inhibitors as potential anti-angiogenic drugs in future.
format Article
id doaj-art-c359d90c0b8643bc9b19a5fb2a5035c2
institution OA Journals
issn 3005-1193
language English
publishDate 2025-05-01
publisher Springer
record_format Article
series Discover Chemistry
spelling doaj-art-c359d90c0b8643bc9b19a5fb2a5035c22025-08-20T02:25:17ZengSpringerDiscover Chemistry3005-11932025-05-012111710.1007/s44371-025-00173-42D-QSAR driven design, molecular docking, molecular dynamics simulation and MM/GBSA studies on quinazoline derivatives for development of VEGFR-2 inhibitorsMayank Kashyap0Saurabh Gupta1Yogita Bansal2Gulshan Bansal3Drug Design and Synthesis Lab, Department of Pharmaceutical Sciences and Drug Research, Punjabi UniversityDrug Design and Synthesis Lab, Department of Pharmaceutical Sciences and Drug Research, Punjabi UniversityDrug Design and Synthesis Lab, Department of Pharmaceutical Sciences and Drug Research, Punjabi UniversityDrug Design and Synthesis Lab, Department of Pharmaceutical Sciences and Drug Research, Punjabi UniversityAbstract Various in silico approaches were utilized to design quinazoline derivatives as anticancer agents targeting VEGFR-2. Six 2D-QSAR models were generated via Monte Carlo optimization method of CORALSEA software. The models generated were based on hybrid optimal descriptors including Graphs and SMILES. All the models were observed to be robust via internal and external validation parameters. Model M4 of split 3 had the highest average correlation coefficient (r 2  = 0.7209), average $${R}_{m}^{2}$$ R m 2 value (0.6256) and lowest error (MAE = 0.465). Structural features characterised as promoters and hinderers were extracted from the M4 model and their mechanistic interpretation was employed to improve the pIC50 of least potent compound SM33 from the dataset. Removal of hinderers and addition of promoters resulted in a series of compounds, DDSL6a-DDSL6l among which DDSL6j was the most active compound with a pIC50 of 8.07. All the compounds were subsequently subjected to molecular docking against PDB: 4ASD, with DDSL6j having better LF rank score (− 14.032) and LF dG score (− 12.318). MD simulation of DDSL6j, sorafenib and SM33 carried out at 200 ns revealed sustained interactions of these compounds with all the key amino acid residues of the protein. Improved potency of compound DDSL6j was attributed to presence of additional stable H-Bond interaction. MM/GBSA calculation also indicated stable complex with binding free energy of 0.28 kcal/mol. All the designed molecules were also found to be possess good ADMET profile. These findings indicate that the quinazoline nucleus can be investigated and refined to create novel VEGFR-2 inhibitors as potential anti-angiogenic drugs in future.https://doi.org/10.1007/s44371-025-00173-4QuinazolineQSARMonte CarloVEGFR-2Molecular dockingMD simulation
spellingShingle Mayank Kashyap
Saurabh Gupta
Yogita Bansal
Gulshan Bansal
2D-QSAR driven design, molecular docking, molecular dynamics simulation and MM/GBSA studies on quinazoline derivatives for development of VEGFR-2 inhibitors
Discover Chemistry
Quinazoline
QSAR
Monte Carlo
VEGFR-2
Molecular docking
MD simulation
title 2D-QSAR driven design, molecular docking, molecular dynamics simulation and MM/GBSA studies on quinazoline derivatives for development of VEGFR-2 inhibitors
title_full 2D-QSAR driven design, molecular docking, molecular dynamics simulation and MM/GBSA studies on quinazoline derivatives for development of VEGFR-2 inhibitors
title_fullStr 2D-QSAR driven design, molecular docking, molecular dynamics simulation and MM/GBSA studies on quinazoline derivatives for development of VEGFR-2 inhibitors
title_full_unstemmed 2D-QSAR driven design, molecular docking, molecular dynamics simulation and MM/GBSA studies on quinazoline derivatives for development of VEGFR-2 inhibitors
title_short 2D-QSAR driven design, molecular docking, molecular dynamics simulation and MM/GBSA studies on quinazoline derivatives for development of VEGFR-2 inhibitors
title_sort 2d qsar driven design molecular docking molecular dynamics simulation and mm gbsa studies on quinazoline derivatives for development of vegfr 2 inhibitors
topic Quinazoline
QSAR
Monte Carlo
VEGFR-2
Molecular docking
MD simulation
url https://doi.org/10.1007/s44371-025-00173-4
work_keys_str_mv AT mayankkashyap 2dqsardrivendesignmoleculardockingmoleculardynamicssimulationandmmgbsastudiesonquinazolinederivativesfordevelopmentofvegfr2inhibitors
AT saurabhgupta 2dqsardrivendesignmoleculardockingmoleculardynamicssimulationandmmgbsastudiesonquinazolinederivativesfordevelopmentofvegfr2inhibitors
AT yogitabansal 2dqsardrivendesignmoleculardockingmoleculardynamicssimulationandmmgbsastudiesonquinazolinederivativesfordevelopmentofvegfr2inhibitors
AT gulshanbansal 2dqsardrivendesignmoleculardockingmoleculardynamicssimulationandmmgbsastudiesonquinazolinederivativesfordevelopmentofvegfr2inhibitors