Unleashing Nature’s potential: a computational approach to discovering novel VEGFR-2 inhibitors from African natural compound using virtual screening, ADMET analysis, molecular dynamics, and MMPBSA calculations

Abstract

One of the characteristic features of cancer is angiogenesis, the process by which new, aberrant blood vessels are formed from pre-existing blood vessels. The process of angiogenesis begins when VEGF binds to its receptor, the VEGF receptor (VEGFR). The formation of new blood vessels provides nutrients that can promote the growth of cancer cells. When it comes to new blood vessel formation, VEGFR2 is a critical player. Therefore, inhibiting VEGFR2 is an effective way to target angiogenesis in cancer treatment. The aim of our research was to find new VEGFR-2 inhibitors by performing a virtual screening of 13313 from African natural compounds using different in silico techniques. Using molecular docking calculations and ADMET properties, we identified four compounds that exhibited a binding affinity ranging from −11.0 kcal/mol to −11.5 Kcal/mol when bound to VEGFR-2. These four compounds were further analyzed with 100 ns simulations to determine their stability and binding energy using the MM-PBSA method. After comparing the compounds with Regorafenib, a drug approved for anti-angiogenesis treatment, it was found that all the candidates (EANPDB 252, NANPDB 4577, and NANPDB 4580), with the exception of EANPDB 76, could target VEGFR-2 similarly effectively to Regorafenib. Therefore, we recommend three of these agents for anti-angiogenesis treatment because they are likely to deactivate VEGFR-2 and thus inhibit angiogenesis. However, it should be noted that the safety and suitability of these agents for clinical use needs further investigation, as the computer-assisted study did not include in vitro or in vivo experiments.

QSAR, ADMET, molecular docking, and dynamics studies of 1,2,4-triazine-3(2H)-one derivatives as tubulin inhibitors for breast cancer therapy

“…In the simulation, this approach assesses protein flexibility and can identify residues with high or low flexibility. Root-mean-square-fluctuations (RMSF) were calculated for each complex to reveal which residues changed the most during the MD simulation 47 , 53 . Higher RMSF values for biomolecular system corresponds to a lower residual stability and vice versa.…”

Discovery of A Novel Series of Quinazoline–Thiazole Hybrids as Potential Antiproliferative and Anti-Angiogenic Agents

“…The search space was defined as a cube with the coordinates of the center x = −21.966, y = −0.473, and z = −11.442 and sides equal to 20, according to an adaptation of a previously reported protocol [42]. The parameters' exhaustiveness and num_modes were set to 50 and 20, respectively.…”

Potent VEGFR-2 inhibitors for resistant breast cancer: a comprehensive 3D-QSAR, ADMET, molecular docking and MMPBSA calculation on triazolopyrazine derivatives

“…Water molecules and co-crystallized small molecules were removed from the protein structure, the polar hydrogen and Kollmann charges were added to the structure (Shivanika et al, 2020). The grid box spacing was set to 0.375 Å, the center to (−24.611 Å, −0.388 Å, −10.929 Å), and the lattice size to 20 Å × 20 Å × 20 Å where the co-crystallized ligand interacts with the active residues (Baammi et al, 2023a). Nine poses were constructed for each protein-ligand complex based on docking affinity.…”