Exploring binding mechanisms of VEGFR2 with three drugs lenvatinib, sorafenib, and sunitinib by molecular dynamics simulation and free energy calculation

From Deep Learning to the Discovery of Promising VEGFR-2 Inhibitors

Vascular endothelial growth factor receptor 2 (VEGFR-2) stands as a prominent therapeutic target in oncology, playing a critical role in angiogenesis, tumor growth, and metastasis. FDA-approved VEGFR-2 inhibitors are associated with diverse side effects. Thus, finding novel and more effective inhibitors is of utmost importance. In this study, a deep learning (DL) classification model was first developed and then employed to select putative active VEGFR-2 inhibitors from an in-house chemical library including 187 druglike compounds. A pool of 18 promising candidates was shortlisted and screened against VEGFR-2 by using molecular docking. Finally, two compounds, RHE-334 and EA-11, were prioritized as promising VEGFR-2 inhibitors by employing PLATO, our target fishing and bioactivity prediction platform. Based on this rationale, we prepared RHE-334 and EA-11 and successfully tested their anti-proliferative potential against MCF-7 human breast cancer cells with IC50 values of 26.78±4.02 and 38.73±3.84 μM, respectively. Their toxicities were instead challenged against the WI-38. Interestingly, expression studies indicated that, in the presence of RHE-334, VEGFR-2 was equal to 0.52±0.03, thus comparable to imatinib equal to 0.63±0.03. In conclusion, this workflow based on theoretical and experimental approaches demonstrates effective in identifying VEGFR-2 inhibitors and can be easily adapted to other medicinal chemistry goals.

Identification of indole-based natural compounds as inhibitors of PARP-1 against triple-negative breast cancer: a computational study

Triple-negative breast cancer (TNBC) is the most aggressive kind of breast cancer known to mankind. It is a heterogeneous disease that is formed due to the missing estrogen, progesterone and human epidermal growth factor 2 receptors. Poly(ADP-ribose) polymerase-1 (PARP-1) protein helps in the development of TNBC by repairing the cancer cells, which proliferate and spread metastatically. To determine the potential PARP-1 inhibitors (PARPi), 0.2 million natural products from Universal Natural Product Database were screened using molecular docking and six hit compounds were selected based on their binding affinity towards PARP-1. The bio-availability and drug-like properties of these natural products were evaluated using ADMET analysis. Molecular dynamics simulations were conducted for these complexes for 200 ns to examine their structural stability and dynamic behaviour and further compared with the complex of talazoparib (TALA), an FDA-approved PARPi. Using MM/PBSA calculations, we conclude that the complexes HIT-3 and HIT-5 (-25.64 and -23.14 kcal/mol, respectively) show stronger binding energies with PARP-1 than TALA with PARP-1 (-10.74 kcal/mol). Strong interactions were observed between the compounds and hotspot residues, Asp770, Ala880, Tyr889, Tyr896, Ala898, Asp899 and Tyr907, of PARP-1 due to the existence of various types of non-covalent interactions between the compounds and PARP-1. This research offers critical information about PARPi, which could potentially be incorporated into the treatment of TNBC. Moreover, these findings were validated by comparing them with an FDA-approved PARPi.

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

Considering the pivotal role of angiogenesis in solid tumor progression, we developed a novel series of quinazoline-thiazole hybrids (SA01-SA07) as antiproliferative and anti-angiogenic agents. Four out of the seven compounds displayed superior antiproliferative activity (IC50 =1.83-4.24 µM) on HepG2 cells compared to sorafenib (IC50 = 6.28 µM). The affinity towards the VEGFR2 kinase domain was assessed through in silico prediction by molecular docking, molecular dynamics studies, and MM-PBSA. The series displayed a high degree of similarity to sorafenib regarding the binding pose within the active site of VEGFR2, with a different orientation of the 4-substituted-thiazole moieties in the allosteric pocket. Molecular dynamics and MM-PBSA evaluations identified SA05 as the hybrid forming the most stable complex with VEGFR2 compared to sorafenib. The impact of the compounds on vascular cell proliferation was assessed on EA.hy926 cells. Six compounds (SA01-SA05, SA07) displayed superior anti-proliferative activity (IC50 = 0.79-5.85 µM) compared to sorafenib (IC50 = 6.62 µM). The toxicity was evaluated on BJ cells. Further studies of the anti-angiogenic effect of the most promising compounds, SA04 and SA05, through the assessment of impact on EA.hy296 motility using a wound healing assay and in ovo potential in a CAM assay compared to sorafenib, led to the confirmation of the anti-angiogenic potential.

New benzothiazole hybrids as potential VEGFR-2 inhibitors: design, synthesis, anticancer evaluation, and in silico study

A new series of 2-aminobenzothiazole hybrids linked to thiazolidine-2,4-dione 4a-e, 1,3,4-thiadiazole aryl urea 6a-d, and cyanothiouracil moieties 8a-d was synthesised. The in vitro antitumor effect of the new hybrids was assessed against three cancer cell lines, namely, HCT-116, HEPG-2, and MCF-7 using Sorafenib (SOR) as a standard drug. Among the tested compounds, 4a was the most potent showing IC50 of 5.61, 7.92, and 3.84 µM, respectively. Furthermore, compounds 4e and 8a proved to have strong impact on breast cancer cell line with IC50 of 6.11 and 10.86 µM, respectively. The three compounds showed a good safety profile towards normal WI-38 cells. Flow cytometric analysis of the three compounds in MCF-7 cells revealed that compounds 4a and 4c inhibited cell population in the S phase, whereas 8a inhibited the population in the G1/S phase. The most promising compounds were subjected to a VEGFR-2 inhibitory assay where 4a emerged as the best active inhibitor of VEGFR-2 with IC50 91 nM, compared to 53 nM for SOR. In silico analysis showed that the three new hybrids succeeded to link to the active site like the co-crystallized inhibitor SOR.

Novel VEGFR-2 inhibitors as antiangiogenic and apoptotic agents via paracrine and autocrine cascades: Design, synthesis, and biological evaluation

Appertaining to its paracrine and autocrine signaling loops, VEGFR-2 succeeded in grabbing attention as one of the leading targets in cancer treatment. Based on the foregoing and our comprehensive studies regarding pharmacophoric features and activity of sorafenib, novel phenylpyridazinone based VEGFR-2 inhibitors 4, 6a-e, 7a,b, 9a,b, 12a-c, 13a,b, 14a,b, 15a,b, and 17a-d were optimized. An assortment of biological assays was conducted to assess the antiangiogenic and apoptotic activities of the synthesized derivatives. In vitro VEGFR-2 kinase assay verified the inhibitory activity of the synthesized derivatives with IC50 values from 49.1 to 418.0 nM relative to the reference drug sorafenib (IC50 = 81.8 nM). Antiproliferative activity against HUVECs revealed that compounds 2-{2-[2-(6-oxo-3-phenylpyridazin-1(6H)-yl)acetyl]hydrazineyl}-N-(p-tolyl)acetamide (12c) and 2-[(5-mercapto-4-methyl-4H-1,2,4-triazol-3-yl)methyl]-6-phenylpyridazin-3(2H)-one (13a) possessed superior activity (IC50 values = 11.5 and 12.3 nM, respectively) in comparison to sorafenib (IC50 = 23.2 nM). For the purpose of appraising their antiproliferative effect, derivatives 12c and 13a were exposed to cell cycle analysis, apoptotic, cell invasion and migration assays in addition to determination of VEGFR-2 in protein level. Moreover, cytotoxicity as well as selectivity index against WI-38 cell line was measured to examine safety of derivatives 12c and 13a. After that, molecular docking study was executed on the top five compounds in the in vitro VEGFR-2 kinase assay 6d, 12c, 13a, 14a and 17c to get a deep perception on binding mode of the synthesized compounds and correlate the design strategy with biological results. Finally, physicochemical, pharmacokinetic properties, and drug-likeness studies were performed on the top five derivative in in vitro VEGFR-2 kinase assay.

The use of machine learning modeling, virtual screening, molecular docking, and molecular dynamics simulations to identify potential VEGFR2 kinase inhibitors

Targeting the signaling pathway of the Vascular endothelial growth factor receptor-2 is a promising approach that has drawn attention in the quest to develop novel anti-cancer drugs and cardiovascular disease treatments. We construct a screening pipeline using machine learning classification integrated with similarity checks of approved drugs to find new inhibitors. The statistical metrics reveal that the random forest approach has slightly better performance. By further similarity screening against several approved drugs, two candidates are selected. Analysis of absorption, distribution, metabolism, excretion, and toxicity, along with molecular docking and dynamics are performed for the two candidates with regorafenib as a reference. The binding energies of molecule1, molecule2, and regorafenib are - 89.1, - 95.3, and - 87.4 (kJ/mol), respectively which suggest candidate compounds have strong binding to the target. Meanwhile, the median lethal dose and maximum tolerated dose for regorafenib, molecule1, and molecule2 are predicted to be 800, 1600, and 393 mg/kg, and 0.257, 0.527, and 0.428 log mg/kg/day, respectively. Also, the inhibitory activity of these compounds is predicted to be 7.23 and 7.31, which is comparable with the activity of pazopanib and sorafenib drugs. In light of these findings, the two compounds could be further investigated as potential candidates for anti-angiogenesis therapy.

An oleanolic acid derivative, K73-03, inhibits pancreatic cancer cells proliferation in vitro and in vivo via blocking EGFR/Akt pathway

Oleanolic acid (OA) and its derivatives show potent anticancer function. Pancreatic cancer (PC) is the fourth core motive of cancer-related deaths worldwide. Epidermal growth factor receptor (EGFR) has been implicated in PC and has been validated as a therapeutic target. Our study demonstrated that K73-03, an OA derivative, was identified as a potent inhibitor of EGFR by using reverse pharmacophore screening and molecular dynamics simulation assays. Moreover, Western blot analysis showed that K73-03 markedly suppressed the levels of phosphorylated-EGFR (p-EGFR) and phosphorylated-Akt (p-Akt). The inhibitory effect of K73-03 on PC cells was assessed in vitro and in vivo. Mechanistically, K73-03 effectively inhibited the cell proliferation of PC cells, and induced apoptosis and autophagy of ASPC-1 cells in a dose-dependent manner. Additionally, pretreatment with chloroquine, an autophagy inhibitor, significantly inhibited K73-03-induced autophagy and enhanced K73-03-induced apoptotic cell death. K73-03 also strongly repressed ASPC-1 cells xenograft growth in vivo. Thus, all these findings provided new clues about OA analog K73-03 as an effective anticancer agent targeted EGFR against ASPC-1 cells, it is worth further evaluation in the future.

Optimizing the Calculation of Free Energy Differences in Nonequilibrium Work SQM/MM Switching Simulations

A key step during indirect alchemical free energy simulations using quantum mechanical/molecular mechanical (QM/MM) hybrid potential energy functions is the calculation of the free energy difference ΔAlow→high between the low level (e.g., pure MM) and the high level of theory (QM/MM). A reliable approach uses nonequilibrium work (NEW) switching simulations in combination with Jarzynski's equation; however, it is computationally expensive. In this study, we investigate whether it is more efficient to use more shorter switches or fewer but longer switches. We compare results obtained with various protocols to reference free energy differences calculated with Crooks' equation. The central finding is that fewer longer switches give better converged results. As few as 200 sufficiently long switches lead to ΔAlow→high values in good agreement with the reference results. This optimized protocol reduces the computational cost by a factor of 40 compared to earlier work. We also describe two tools/ways of analyzing the raw data to detect sources of poor convergence. Specifically, we find it helpful to analyze the raw data (work values from the NEW switching simulations) in a quasi-time series-like manner. Principal component analysis helps to detect cases where one or more conformational degrees of freedom are different at the low and high level of theory.

Theoretical investigations on the effects of mutations in important residues of NS1B on its RNA-binding using molecular dynamics simulations

NS1B protein plays an important role in countering host antiviral defense and virulence of influenza virus B, considered as the promising target. The first experimental structure of the NS1B protein has recently been determined, was able to bind to double-stranded RNA (dsRNA). However, few studies attempt to investigate the RNA-binding mechanism of the NS1B. In this study, we provide our understanding of the structure-function relationship, dynamics and RNA-binding mechanism of the NS1B protein by performing molecular dynamics simulations combined and MM-GBSA calculations on the NS1B-dsRNA complex. 12 key residues are identified for RNA-binding by forming hydrogen bonds with the. Our results also demonstrate that mutations (R156A, K160A, R208A and K221A) can cause the local structure changes of NS1B CTD and the hydrogen bonds between NS1B CTD and RNA disappearance, which may be the main reasons for the decrease in RNA-binding affinity. These results mentioned will help us understanding the RNA-binding mechanism and could provide some medicinal chemistry insights chances for rational drug design targeting NS1B protein.

Synthesis and Characterization of Novel Copper(II)-Sunitinib Complex: Molecular Docking, DFT Studies, Hirshfeld Analysis and Cytotoxicity Studies

The main goal of this work was to report the synthesis, characterization, and cytotoxicity study of a novel copper(II)-sunitinib complex, CuSun. It has been synthesized and characterized in solid state and in solution by different methods (such as DFT, FTIR, Raman, UV-vis, EPR, NMR, etc.). The solid-state molecular structure of trichlorosunitinibcopper(II), where sunitinib: N-[2-(diethylamino)ethyl]-5-[(Z)-(5-fluoro-2-oxo-1H-indol-3-ylidene)methyl]-2,4-dimethyl-1H-pyrrole-3-carboxamide, for short Cu(Sun)Cl3, was determined by X-ray diffraction. It crystallizes in the triclinic space group P-1 with a = 7.9061(5) Å, b = 12.412(1) Å, c = 13.7005(8) Å, α = 105.021(6)°, β = 106.744(5)°, γ = 91.749(5)°, and Z = 2 molecules per unit cell. Also, we have found π-π interactions and classic and non-classic H-bonds in the crystal structure by using Hirshfeld surface analysis. In the speciation studies, the complex has dissociated in protonated sunitinib and chlorocomplex of copper(II), according to 1HNMR, EPR, UV-vis and conductimetric analysis. Molecular docking of the complex in both, ATP binding site and allosteric site of VEGFR2 have shown no improvement in comparison to the free ligand. Besides, cytotoxicity assay on HepG2 cell line shows similar activity for complex and ligand in the range between 1–25 μM supporting the data obtained from studies in solution.

Synthesis, in vitro anticancer activity and in silico studies of certain isoxazole-based carboxamides, ureates, and hydrazones as potential inhibitors of VEGFR2

The ensuing research presents the results of in vitro anticancer activity of novel 28 compounds of isoxazole-based carboxamides 3(a-d); ureates 4(a-g), 5, 6, 7a,b, 8; and hydrazones 9(a-f), 10(a-d), 11a,b as potential inhibitors of VEGFR2. The carboxamides and ureates were synthesized by converting 5-(aryl)-isoxzaole-3-carbohydrazides 1a,b to the corresponding carbonylazides 2a,b followed by treatment with the appropriate amines. The hydrazones were directly obtained through condensation of the carbohydrazide 1a,b with aldehydes and/or ketones. The structures of the target compounds were confirmed by elemental and spectral analyses. A preliminary in vitro anticancer screening of solutions (10^-5M) on 60 cancer cell lines (NCI, USA) revealed that the carboxamide 3c is the most promising growth inhibitor. Explicitly, 3c showed potent anticancer activity at 10µ M against leukemia (HL-60(TB), K-562 and MOLT-4), colon cancer (KM12) and melanoma (LOX IMVI) cell lines with %GI range = 70.79-92.21. Evaluation of growth inhibitory activity of the synthesized compounds against hepatocellular carcinoma (HepG2), that overexpresses VEGFR2, showed superior activity of compounds 8, 10a and 10c with IC₅₀ in sub micromolar concentrations of 0.84, 0.79 and 0.69 μM, respectively, which is better than that of the reference drug, Sorafenib (IC₅₀ = 3.99 µM). Moreover, these compounds displayed high selective cytotoxicity for HepG2 cancer cells over the nontumorigenic THLE2 liver cells (SI range = 26.37-38.60) which reflect their safety. The results of VEGFR2 kinase inhibition assay demonstrate that, compounds 8 and 10a are the most active inhibitors with IC₅₀ = 25.7 and 28.2 nM, respectively, (Sorafenib IC₅₀ = 28.1 nM). Molecular docking of the synthesized derivatives to VEGFR2 (PDB: 3WZE) showed similar binding modes to that of the co-crystallized ligand, sorafenib. Moreover, the results of computational assessment of ADME and drug-likeness characteristics inspire further investigations of the new isoxazole-based derivatives to afford more potent, safe and orally active VEGFR2 inhibitors as potential anticancer drug candidates.

Structural Basis for Vascular Endothelial Growth Factor Receptor Activation and Implications for Disease Therapy

Vascular endothelial growth factors (VEGFs) bind to membrane receptors on a wide variety of cells to regulate diverse biological responses. The VEGF-A family member promotes vasculogenesis and angiogenesis, processes which are essential for vascular development and physiology. As angiogenesis can be subverted in many disease states, including tumour development and progression, there is much interest in understanding the mechanistic basis for how VEGF-A regulates cell and tissue function. VEGF-A binds with high affinity to two VEGF receptor tyrosine kinases (VEGFR1, VEGFR2) and with lower affinity to co-receptors called neuropilin-1 and neuropilin-2 (NRP1, NRP2). Here, we use a structural viewpoint to summarise our current knowledge of VEGF-VEGFR activation and signal transduction. As targeting VEGF-VEGFR activation holds much therapeutic promise, we examine the structural basis for anti-angiogenic therapy using small-molecule compounds such as tyrosine kinase inhibitors that block VEGFR activation and downstream signalling. This review provides a rational basis towards reconciling VEGF and VEGFR structure and function in developing new therapeutics for a diverse range of ailments.

Design, synthesis, biological evaluation, and modeling studies of novel conformationally-restricted analogues of sorafenib as selective kinase-inhibitory antiproliferative agents against hepatocellular carcinoma cells

Sorafenib is one of the clinically used anticancer agents that inhibits several kinases. In this study, novel indole-based rigid analogues of sorafenib were designed and synthesized in order to enhance kinase selectivity and hence minimize the side effects associated with its use. The target compounds possess different linkers; urea, amide, sulfonamide, or thiourea, in addition to different terminal aryl moieties attached to the linker in order to investigate their impact on biological activity. They were tested against Hep3B, Huh7, and Hep-G2 hepatocellular carcinoma (HCC) cell lines to study their potency. Among all the tested target derivatives, compound 1h exerted superior antiproliferative potency against all the three tested HCC cell lines compared to sorafenib. Based on these preliminary results, compound 1h was selected for further biological and in silico investigations. Up to 30 μM, compound 1h did not inhibit 50% of the proliferation of WI-38 normal cells, which indicated promising selectivity against HCC cells than normal cells. In addition, compound 1h exerted superior kinase selectivity than sorafenib. It is selective for VEGFR2 and VEGFR3 angiogenesis-related kinases, while sorafenib is a multikinase inhibitor. Superior kinase selectivity of compound 1h to sorafenib can be attributed to its conformationally-restricted indole nucleus and the bulky N-methylpiperazinyl moiety. Western blotting was carried out and confirmed the ability of compound 1h to inhibit VEGFR2 kinase inside Hep-G2 HCC cells in a dose-dependent pattern. Compound 1h induces apoptosis and necrosis in Hep-G2 cell line, as shown by caspase-3/7 and lactate dehydrogenase (LDH) release assays, respectively. Moreover, compound 1h is rather safe against hERG. Thus, we could achieve a more selective kinase inhibitor than sorafenib with retained or even better antiproliferative potency against HCC cell lines. Furthermore, molecular docking and dynamic simulation studies were carried out to investigate its binding mode with VEGFR2 kinase. The molecule has a unique orientation upon binding with the kinase.

Design, synthesis, biological activity evaluation of 3-(4-phenyl-1H-imidazol-2-yl)-1H-pyrazole derivatives as potent JAK 2/3 and aurora A/B kinases multi-targeted inhibitors

In this study, a series of 3-(4-phenyl-1H-imidazol-2-yl)-1H-pyrazole derivatives were designed, synthesized, and evaluated for their biological activities. Upon performing kinase assays, most of the compounds exhibited potent inhibition against JAK2/3 and Aurora A/B with the IC₅₀ values ranging from 0.008 to 2.52 μM. Among these derivatives, compound 10e expressed the most moderate inhibiting activities against all the four kinases with the IC₅₀ values of 0.166 μM (JAK2), 0.057 μM (JAK3), 0.939 μM (Aurora A), and 0.583 μM (Aurora B), respectively. Moreover, most of the derived compounds exhibited potent cytotoxicity against human chronic myeloid leukemia cells K562 and human colon cancer cells HCT116, while compound 10e expressed antiproliferative activities against K562 (IC₅₀₌6.726  μM). According to western blot analysis, compound 10e down-regulated the phosphorylation of STAT3, STAT5, Aurora A, and Aurora B in a dose-dependent manner in K562 and HCT116 cells. Cell cycle analysis revealed that compound 10e inhibited the proliferation of cells by inducing cell cycle arrest in the G2 phase. The molecular modeling suggested that compound 10e could maintain a binding mode similar to the binding mode of AT9832, a common JAK 2/3 and Aurora A/B kinases multi-target kinase inhibitor. Therefore, compound 10e might be a potential agent for cancer therapy deserving further research.