In silico discovery of potential VEGFR-2 inhibitors from natural derivatives for anti-angiogenesis therapy

Interaction of plant phenol vanillin with human hemoglobin: A spectroscopic and molecular docking study

Vanillin is a phenolic aldehyde widely used as a flavouring agent in the food industry. Vanillin has many health benefits and has gained attention in pharmacological industries also, due to its antioxidant properties and non-toxic nature. The interaction of vanillin with human hemoglobin (hHb), an abundant tetrameric heme protein, was investigated by several spectroscopic techniques and molecular modeling methods. UV-visible spectra showed that the binding of vanillin to hHb induces structural changes due to alterations in the micro-environment of hHb. Vanillin quenches the intrinsic fluorescence of hHb by the dynamic mechanism, which was confirmed by both temperature dependent and time resolved fluorescence studies. Vanillin binds spontaneously to hHb at a single site and the binding is stabilized by hydrogen bonds and hydrophobic interactions. The circular dichroism spectra showed that the binding of vanillin altered the secondary structure of hHb due to change in its alpha-helical content. Molecular docking identified the amino acids of hHb involved in binding to vanillin and also that the free energy change of the binding reaction is -5.5 kcal/mol. Thus, our results indicate that vanillin binds spontaneously to hHb at a single site and alters its secondary structure. This will help in understanding the potential use of vanillin and related antioxidants as therapeutic agents in various hematological disorders.

Biochemical and molecular docking-based strategies of Acalypha indica and Boerhavia diffusa extract by targeting bacterial strains and cancer proteins

Antibiotic-resistant microbes have emerged around the world, presenting a risk to health. Plant-derived drugs have become a potential source for the production of antibiotic-resistant drugs and cancer therapies. In this study, we investigated the antibacterial, cytotoxic and antioxidant properties of Acalypha indica and Boerhavia diffusa, and conducted in silico molecular docking experiments against EGFR and VEGFR-2 proteins. The metabolic extract of A. indica inhibited Streptococcus iniae and Staphylococcus sciuri with inhibition zones of 21.66 ± 0.57 mm and 20.33 ± 0.57 mm, respectively. The B. diffusa leaf extract produced inhibition zones of 20.3333 ± 0.5773 mm and 20.33 ± 0.57 mm against Streptococcus iniae and Edwardsiella anguillarum, respectively. A. indica and B. diffusa extracts had toxicities of 162.01 μg/ml and 175.6 μg/ml, respectively. Moreover, B. diffusa (IC50 =154.42 µg/ml) leaf extract exhibited moderately higher antioxidant activity compared with the A. indica (IC50 = 218.97 µg/ml) leaf extract. Multiple interactions were observed at Leu694, Met769 and Leu820 sites for EGFR and at Asp1046 and Cys1045 sites for VEGFR during the molecular docking study. CID-235030, CID-70825 and CID-156619353 had binding energies of -7.6 kJ/mol, -7.5 kJ/mol and -7.6 kJ/mol, respectively, with EGFR protein. VEGFR-2 protein had docking energies of -7.5 kJ/mol, -7.6 kJ/mol and -7.3 kJ/mol, respectively, for CID-6420353, CID-156619353 and CID-70825 compounds. The MD simulation trajectories revealed the hit compound; CID-235030 and EGFR complex, CID-6420353 and VEGFR-2 exhibit stable profile in the root mean square deviation (RMSD), radius of gyration (Rg), solvent accessible surface area (SASA), hydrogen bond and root mean square fluctuation (RMSF) and the binding free energy by MM-PBSA method. This study indicates that methanol extracts of A. indica and B. diffusa may play a crucial role in developing antibiotic-resistant and cancer drugs.Communicated by Ramaswamy H. Sarma.

Novel quinoxaline-based VEGFR-2 inhibitors to halt angiogenesis

Vascular endothelial growth factor receptor-2 is a dynamic target for therapeutic intervention in various types of cancer. This study was aimed at exploring the VEGFR-2 inhibitory activity of a novel library of quinoxalin-2-one derivatives such as 3-furoquinoxaline carboxamides, 3-pyrazolylquinoxalines, and 3-pyridopyrimidyl-quinoxalines. Among them, 6c, 7a, and 7d-f produced remarkable cytotoxicity against HCT-116 (IC50's 4.28-9.31 µM) and MCF-7 (IC50's 3.57-7.57 µM) cell lines using the MTT assay and doxorubicin (DOX) as a reference standard. Interestingly, results of cytotoxicity towards the human fibroblast cell line WI38 revealed that these hits demonstrated higher selectivity indices towards both HCT-116 (SI 8.69-23.19) and MCF-7 (SI 9.48-27.80) than DOX, SI 0.72 and 0.90, respectively. Then, these hits were subjected to a mechanistic study; they showed direct inhibition of VEGFR-2. Impressively, compound 7f displayed 1.2 times the VEGFR-2 inhibitory activity of sorafenib. The antiangiogenic potential of 7f was proved via lowering the level of VEGF-A, than that of control. It as well, exhibited scratch closure percent of 61.8%, compared with 74.5% of control at 48 hrs, indicating the potential anti-migratory effect of the compound 7f. It significantly increased the expression of tumor suppressor gene (p53) on MCF-7 cells by almost 18 folds and upregulated the caspase-3 level by 10.7 folds, compared to the control. Cell cycle analysis revealed cell cycle arrest at G2/M together with a PreG increase which indicated apoptosis induction potential. Annexin V-FITC apoptosis results proposed the two modes of cell death (apoptosis and necrosis) as an inherent mechanism of cytotoxicity of compound 7f. Molecular docking further supported the mechanism showing the affinity of target compounds for VEGFR-2 active site. Moreover, physicochemical and drug-like properties were assessed from the ADME properties.

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.

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.

Inhibitory Effects of Myriocin on Non-Enzymatic Glycation of Bovine Serum Albumin

Advanced glycation end products (AGEs) are the compounds produced by non-enzymatic glycation of proteins, which are involved in diabetic-related complications. To investigate the potential anti-glycation activity of Myriocin (Myr), a fungal metabolite of Cordyceps, the effect of Myr on the formation of AGEs resulted from the glycation of bovine serum albumin (BSA) and the interaction between Myr and BSA were studied by multiple spectroscopic techniques and computational simulations. We found that Myr inhibited the formation of AGEs at the end stage of glycation reaction and exhibited strong anti-fibrillation activity. Spectroscopic analysis revealed that Myr quenched the fluorescence of BSA in a static process, with the possible formation of a complex (approximate molar ratio of 1:1). The binding between BSA and Myr mainly depended on van der Waals interaction, hydrophobic interactions and hydrogen bond. The synchronous fluorescence and UV-visible (UV-vis) spectra results indicated that the conformation of BSA altered in the presence of Myr. The fluorescent probe displacement experiments and molecular docking suggested that Myr primarily bound to binding site 1 (subdomain IIA) of BSA. These findings demonstrate that Myr is a potential anti-glycation agent and provide a theoretical basis for the further functional research of Myr in the prevention and treatment of AGEs-related diseases.

Interaction of the insecticide bioallethrin with human hemoglobin: biophysical, in silico and enzymatic studies

Bioallethrin is an insecticide that is widely used in households resulting in human exposure. Bioallethrin is cytotoxic to human erythrocytes. Here we have studied the interaction of bioallethrin with human hemoglobin (Hb) using in silico and biophysical approaches. Incubation of Hb (5 μM) with bioallethrin (1-50 µM) led to increase in absorbance at 280 nm while the Soret band at 406 nm was slightly reduced. The intrinsic fluorescence of Hb was enhanced with the appearance of a new peak around 305 nm. Synchronous fluorescence showed that the binding of bioallethrin to Hb mainly affects the tyrosine microenvironment. The structural changes in Hb were confirmed with a significant shift in CD spectra and about 25% loss of α-helix. Molecular docking and visualisation through Discovery studio confirmed the formation of Hb-bioallethrin complex with a binding energy of -7.3 kcal/mol. Molecular simulation showed the stability and energy dynamics of the binding reaction between bioallethrin and Hb. The structural changes induced by bioallethrin led to inhibition of the esterase activity of Hb. In conclusion, this study shows that bioallethrin forms a stable complex with human Hb which may lead to loss of Hb function in the body.Communicated by Ramaswamy H. Sarma.

Target-Based Small Molecule Drug Discovery for Colorectal Cancer: A Review of Molecular Pathways and In Silico Studies

Colorectal cancer is one of the most prevalent cancer types. Although there have been breakthroughs in its treatments, a better understanding of the molecular mechanisms and genetic involvement in colorectal cancer will have a substantial role in producing novel and targeted treatments with better safety profiles. In this review, the main molecular pathways and driver genes that are responsible for initiating and propagating the cascade of signaling molecules reaching carcinoma and the aggressive metastatic stages of colorectal cancer were presented. Protein kinases involved in colorectal cancer, as much as other cancers, have seen much focus and committed efforts due to their crucial role in subsidizing, inhibiting, or changing the disease course. Moreover, notable improvements in colorectal cancer treatments with in silico studies and the enhanced selectivity on specific macromolecular targets were discussed. Besides, the selective multi-target agents have been made easier by employing in silico methods in molecular de novo synthesis or target identification and drug repurposing.

Design, synthesis and evaluation of novel tetrahydropyridothienopyrimidin-ureas as cytotoxic and anti-angiogenic agents

The novel derivatives of tetrahydropyridothienopyrimidine-based compounds have been designed and efficiently synthesized with good yields through seven steps reaction. The anticancer activity of compounds 11a-y has been evaluated against MCF-7, PC-3, HEPG-2, SW-480, and HUVEC cell lines by MTT assay. The target compounds showed IC₅₀ values between 2.81–29.6 μg/mL and were compared with sorafenib as a reference drug. Among them, compound 11n showed high cytotoxic activity against four out of five examined cell lines and was 14 times more selective against MRC5. The flow cytometric analysis confirmed the induction of apoptotic cell death by this compound against HUVEC and MCF-7 cells. In addition, 11n caused sub-G1 phase arrest in the cell cycle arrest. Besides, this compound induced anti-angiogenesis in CAM assay and increased the level of caspase-3 by 5.2 fold. The western-blot analysis of the most active compound, 11n, revealed the inhibition of VEGFR-2 phosphorylation. Molecular docking study also showed the important interactions for compound 11n.

In Silico and In Vitro Studies for Benzimidazole Anthelmintics Repurposing as VEGFR-2 Antagonists: Novel Mebendazole-Loaded Mixed Micelles with Enhanced Dissolution and Anticancer Activity

Cancer is a leading cause of death worldwide and its incidence is unfortunately anticipated to rise in the next years. On the other hand, vascular endothelial growth factor receptor 2 (VEGFR-2) is highly expressed in tumor-associated endothelial cells, where it affects tumor-promoting angiogenesis. Therefore, VEGFR-2 is considered one of the most promising therapeutic targets for cancer treatment. Furthermore, some FDA-approved benzimidazole anthelmintics have already shown potential anticancer activities. Therefore, repurposing them against VEGFR-2 can provide a rapid and effective alternative that can be implicated safely for cancer treatment. Hence, 13 benzimidazole anthelmintic drugs were subjected to molecular docking against the VEGFR-2 receptor. Among the tested compounds, fenbendazole (FBZ, 1), mebendazole (MBZ, 2), and albendazole (ABZ, 3) were proposed as potential VEGFR-2 antagonists. Furthermore, molecular dynamics simulations were carried out at 200 ns, giving more information on their thermodynamic and dynamic properties. Besides, the anticancer activity of the aforementioned drugs was tested in vitro against three different cancer cell lines, including liver cancer (HUH7), lung cancer (A549), and breast cancer (MCF7) cell lines. The results depicted potential cytotoxic activity especially against both HUH7 and A549 cell lines. Furthermore, to improve the aqueous solubility of MBZ, it was formulated in the form of mixed micelles (MMs) which showed an enhanced drug release with better promising cytotoxicity results compared to the crude MBZ. Finally, an in vitro quantification for VEGFR-2 concentration in treated HUH7 cells has been conducted based on the enzyme-linked immunosorbent assay. The results disclosed that FBZ, MBZ, and ABZ significantly (p < 0.001) reduced the concentration of VEGFR-2, while the lowest inhibition was achieved in MBZ-loaded MMs, which was even much better than the reference drug sorafenib. Collectively, the investigated benzimidazole anthelmintics could be encountered as lead compounds for further structural modifications and thus better anticancer activity, and that was accomplished through studying their structure-activity relationships.

Glioma-Targeted Therapeutics: Computer-Aided Drug Design Prospective

Amongst the several types of brain cancers known to humankind, glioma is one of the most severe and life-threatening types of cancer, comprising 40% of all primary brain tumors. Recent reports have shown the incident rate of gliomas to be 6 per 100,000 individuals per year globally. Despite the various therapeutics used in the treatment of glioma, patient survival rate remains at a median of 15 months after undergoing first-line treatment including surgery, radiation, and chemotherapy with Temozolomide. As such, the discovery of newer and more effective therapeutic agents is imperative for patient survival rate. The advent of computer-aided drug design in the development of drug discovery has emerged as a powerful means to ascertain potential hit compounds with distinctively high therapeutic effectiveness against glioma. This review encompasses the recent advances of bio-computational in-silico modeling that have elicited the discovery of small molecule inhibitors and/or drugs against various therapeutic targets in glioma. The relevant information provided in this report will assist researchers, especially in the drug design domains, to develop more effective therapeutics against this global disease.

Phthalazine-based VEGFR-2 inhibitors: Rationale, design, synthesis, in silico, ADMET profile, docking, and anticancer evaluations

In the designed compounds, a new linker was inserted in the form of fragments with verified VEGFR-2 inhibitory potential, including an α,β-unsaturated ketonic fragment, pyrazole, and pyrimidine. Also, new distal hydrophobic moieties were attached to these linkers that are expected to increase the hydrophobic interaction with VEGFR-2 and, consequently, the affinity. These structural optimizations have led us to identify the novel dihydropyrazole derivative 6_e as a promising hit molecule. All the new derivatives were evaluated to assess their anticancer activity against three human cancer cell lines, including HepG2, HCT-116, and MCF-7. The results of the in vitro anticancer evaluation study revealed the moderate to excellent cytotoxicity of 6_c , 6_e , 6_g , and 7_b , with IC₅₀ values in the low micromolar range. The inhibitory activity of VEGFR-2 was investigated for 16 of the designed compounds. The enzyme assay results of the new compounds were compared with those of sorafenib as a reference VEGFR-2 inhibitor. The obtained results demonstrated that our derivatives are potent VEGFR-2 inhibitors. The most potent derivatives 6_c , 6_e , 6_g , and 7_b showed IC₅₀ values in the range of 0.11-0.22 µM. Molecular docking and pharmacokinetic studies were also conducted to rationalize the VEGFR-2 inhibitory activity and to evaluate the ability of the most potent derivatives to be developed as good drug candidates.

Synthesis, SAR study, and bioactivity evaluation of a series of Quinoline-Indole-Schiff base derivatives: Compound 10E as a new Nur77 exporter and autophagic death inducer

We previously reported 5-((8-methoxy-2-methylquinolin-4-yl)amino)-1H-indole- 2-carbohydrazide derivatives as new Nur77 modulators. In this study, we explored whether the 8-methoxy-2-methylquinoline moiety and bicyclic aromatic rings at the N'-methylene position were critical for their antitumor activity against hepatocellular carcinoma (HCC). For this purpose, a small library of 5-substituted 1H-indole-2-carbohydrazide derivatives was designed and synthesized. We found that the 8-methoxy-2-methylquinoline moiety was a fundamental structure for its biological function, while the introduction of the bicyclic aromatic ring into the N'-methylene greatly improved its anti-tumor effect. We found that the representative compound 10E had a high affinity to Nur77. The K_D values were in the low micromolar (2.25-4.10 μM), which were coincident with its IC₅₀ values against the tumor cell lines (IC₅₀ < 3.78 μM). Compound 10E could induce autophagic cell death of liver cancer cells by targeting Nur77 to mitochondria while knocking down Nur77 greatly impaired anti-tumor effect. These findings provide an insight into the structure-activity relation of Quinoline-Indole-Schiff base derivatives and further demonstrate that antitumor agents targeting Nur77 may be considered as a promising strategy for HCC therapy.

Synthesis and Rational Design of New Appended 1,2,3-Triazole-uracil Ensembles as Promising Anti-Tumor Agents via In Silico VEGFR-2 Transferase Inhibition

Angiogenesis inhibition is a key step towards the designing of new chemotherapeutic agents. In a view to preparing new molecular entities for cancer treatment, eighteen 1,2,3-triazole-uracil ensembles 5a-r were designed and synthesized via the click reaction. The ligands were well characterized using 1H-, 13C-NMR, elemental analysis and ESI-mass spectrometry. The in silico binding propinquities of the ligands were studied sequentially in the active region of VEGFR-2 using the Molegro virtual docker. All the compounds produced remarkable interactions and potentially inhibitory ligands against VEGFR-2 were obtained with high negative binding energies. Drug-likeness was assessed from the ADME properties. Cytotoxicity of the test compounds was measured against HeLa and HUH-7 tumor cells and NIH/3T3 normal cells by MTT assay. Compound 5h showed higher growth inhibition activity than the positive control, 5-fluorouracil (5-FU), against both HeLa and HUH-7 cells with IC50 values of 4.5 and 7.7 μM respectively. Interestingly, the compounds 5a-r did not show any cytotoxicity towards the normal cell lines. The results advance the position of substituted triazoles in the area of drug design with no ambiguity.

Targeting the Autophagy Specific Lipid Kinase VPS34 for Cancer Treatment: An Integrative Repurposing Strategy

The impact of autophagy on cancer treatment and its corresponding responsiveness has galvanized the scientific community to develop novel inhibitors for cancer treatment. Importantly, the discovery of inhibitors that targets the early phase of autophagy was identified as a beneficial choice. Despite the number of research in recent years, screening of the DrugBank repository (9591 molecules) for the Vacuolar protein sorting 34 (VPS34) has not been reported earlier. Therefore, the present study was designed to identify potential VPS34 antagonists using integrated pharmacophore strategies. Primarily, an energy-based pharmacophore and receptor cavity-based analysis yielded five (DHRRR) and seven featured (AADDHRR) pharmacophore hypotheses respectively, which were utilized for the database screening process. The glide score, the binding free energy, pharmacokinetics and pharmacodynamics properties were examined to narrow down the screened compounds. This analysis yielded a hit molecule, DB03916 that exhibited a better docking score, higher binding affinity and better drug-like properties in contrast to the reference compound that suffers from a toxicity property. Importantly, the result was validated using a 50 ns molecular dynamics simulation study. Overall, we conclude that the identified hit molecule DB03916 is believed to serve as a prospective antagonist against VPS34 for cancer treatment.

Quinazolin-4(3H)-one based agents bearing thiadiazole-urea: Synthesis and evaluation of anti-proliferative and antiangiogenic activity

A series of quinazolin-4(3H)-one based agents containing thiadiazole-urea were designed, synthesized, and biologically evaluated. The proliferation rate of PC3 cells was moderately reduced by compound 9f (IC₅₀ = 17.7 μM)which was comparable with sorafenib (IC₅₀ = 17.3 μM). There was also a significant reduction in the number of HUVEC cells, when they were exposed to compound 9y (IC₅₀ = 6.1 μM). To test the potential of compounds in inducing apoptosis, Annexin V-FITC/propidium iodide double staining assay was used. After the treatment of HUVEC cells with 9f, they underwent apoptotic effects. A substantial effort was dedicated to gathering comprehensive data across CAM assay. These data showed that 9f moderately inhibits the growth of corresponding blood vessels. Finally, the outcomes of Western blotting proposed a mechanism of action, by which the phosphorylation of VEGFR-2 is inhibited by compounds 9f and 9y.

Ligand-Based Pharmacophore Modeling, Molecular Docking, and Molecular Dynamic Studies of Dual Tyrosine Kinase Inhibitor of EGFR and VEGFR2

Epidermal growth factor receptor (EGFR) and vascular endothelial growth factor receptor 2 (VEGFR2) play an important role in cancer growth. Both of them have close relationships. Expression of EGFR will induce an angiogenic factor (VEGF) release for binding with VEGFR2. However, the existence of VEGF up-regulation independent of EGFR leads to cancer cell resistance to anti-EGFR. Therefore, a therapeutic approach targeting EGFR and VEGFR2 simultaneously may improve the outcome of cancer treatment. The present study was designed to identify potential compounds as a dual inhibitor of EGFR and VEGFR2 by the computational method. Firstly, the ligand-based pharmacophore model for each target was setup to screen of ZINC database of purchasable compounds. The hit compounds obtained by pharmacophore screening were then further screened by molecular docking studies. Taking erlotinib (EGFR inhibitor) and axitinib (VEGFR2 inhibitor) as reference drugs, six potential compounds (ZINC08398597, ZINC12047553, ZINC16525481, ZINC17418102, ZINC21942954, and ZINC38484632) were selected based on their docking scores and binding interaction. However, molecular dynamics simulations demonstrated that only ZINC16525481 and ZINC38484632 which have good binding free energy and stable hydrogen bonding interactions with EGFR and VEGFR2. The result represents a promising starting point for developing potent dual tyrosine kinases inhibitor of EGFR and VEGFR2.

Design, synthesis and evaluation of novel thienopyrimidine-based agents bearing diaryl urea functionality as potential inhibitors of angiogenesis

Inhibition of angiogenesis is a promising strategy for the treatment of cancer. Herein, we describe the design and synthesis of thieno[2,3-d]pyrimidine-1,3,4-thiadiazole-aryl urea derivatives 11a-m to evaluate their efficacy in the chick chorioallantoic membrane (CAM) assay. Among target agents, 11i had a considerable activity against prostate cancer cell line, PC3 (IC₅₀ = 3.6 μM). Moreover, induction of apoptosis, good inhibitory activity against the growth of capillary blood vessels, and inhibition of VEGFR-2 phosphorylation were noticeable parameters which convinced us that 11i could be considered as a promising candidate for the discovery of novel drugs to treat tumors, particularly prostate cancer.

Molecular simulation strategies for the discovery of selective inhibitors of β-catenin

Tumor dissemination and relapse in lung cancer were found to be due to the existence of cancer stem cells. In particular, the [Formula: see text]-catenin pathway is found to be one of the crucial pathways in maintaining the stem-like properties of the cells. Thus, targeting the [Formula: see text]-catenin family of proteins is a significant therapeutic route in the treatment of lung cancer. Therefore, in the present study, a pharmacophore-based drug repurposing approach was accomplished to pinpoint potent [Formula: see text]-catenin inhibitors from the DrugBank database. Primarily, ligand-based pharmacophore hypothesis (AAHHR) was generated using existing [Formula: see text]-catenin inhibitors available in the literature and utilized for library screening. Subsequently, the inhibitory activity of the screened compounds was examined by the hierarchical docking process and the Prime MM-GBSA algorithm. Moreover, quantum chemical calculations and molecular dynamics simulations were executed to analyze the inhibitory effects of the screened hit molecule. The results indicate that hit molecule, DB08047 was found to possess better binding free energy, favorable ligand strain energy, satisfactory pharmacokinetic properties and superior free energy landscape profile. Eventually, the pIC[Formula: see text] values of the lead compounds were predicted by the AutoQSAR algorithm. It is noteworthy to mention that DB08047 was found to possess pyrazole scaffolds which could downregulate [Formula: see text]-catenin pathway and thus facilitate the controlled cell growth/inhibit tumor growth.

In Silico Investigation of the Anti-Tumor Mechanisms of Epigallocatechin-3-Gallate

The EGCG, an important component of polyphenol in green tea, is well known due to its numerous health benefits. We employed the reverse docking method for the identification of the putative targets of EGCG in the anti-tumor target protein database and these targets were further uploaded to public databases in order to understand the underlying pharmacological mechanisms and search for novel EGCG-associated targets. Similarly, the pharmacological linkage between tumor-related proteins and EGCG was manually constructed in order to provide greater insight into the molecular mechanisms through a systematic integration with applicable bioinformatics. The results indicated that the anti-tumor mechanisms of EGCG may involve 12 signaling transduction pathways and 33 vital target proteins. Moreover, we also discovered four novel putative target proteins of EGCG, including IKBKB, KRAS, WEE1 and NTRK1, which are significantly related to tumorigenesis. In conclusion, this work may provide a useful perspective that will improve our understanding of the pharmacological mechanism of EGCG and identify novel potential therapeutic targets.

Discovery of Small Molecules that Target Vascular Endothelial Growth Factor Receptor-2 Signalling Pathway Employing Molecular Modelling Studies

Angiogenesis is defined as the formation of new blood vessels and is a key phenomenon manifested in a host of cancers during which tyrosine kinases play a crucial role. Vascular endothelial growth factor receptor-2 (VEGFR-2) is pivotal in cancer angiogenesis, which warrants the urgency of discovering new anti-angiogenic inhibitors that target the signalling pathways. To obtain this objective, a structure-based pharmacophore model was built from the drug target VEGFR-2 (PDB code: 4AG8), complexed with axitinib and was subsequently validated and employed as a 3D query to retrieve the candidate compounds with the key inhibitory features. The model was escalated to molecular docking studies resulting in seven candidate compounds. The molecular docking studies revealed that the seven compounds displayed a higher dock score than the reference-cocrystallised compound. The GROningen MAchine for Chemical Simulations (GROMACS) package guided molecular dynamics (MD) results determined their binding mode and affirmed stable root mean square deviation. Furthermore, these compounds have preserved their key interactions with the residues Glu885, Glu917, Cys919 and Asp1046. The obtained findings deem that the seven compounds could act as novel anti-angiogenic inhibitors and may further assist as the prototype in designing and developing new inhibitors.

S, Balakrishnan A, Katakia YT, Chatterjee S, Nagasundaram N. A marine sponge associated fungal metabolite monacolin X suppresses angiogenesis by down regulating VEGFR2 signaling

Cancer is one of the leading causes of global death and there is an urgent need for the development of cancer treatment; targeting VEGFR2 could be one of the promising therapies. In the present study, previously isolated marine fungal metabolite monacolin X, suppresses in vitro angiogenic characteristics such as proliferation, migration, adhesion, invasion and tube formation of HUVECs when stimulated by VEGF, at a non-toxic concentration. Monacolin X downregulated VEGFR2, PKCα and PKCη mRNA expression. Further, monacolin X inhibited in vivo angiogenesis in CAM assay, vascular sprouting in aortic ring, decreased ISV and SIV length and diameter in Tg (Kdr:EGFP)/ko1 zebrafish embryos. Monacolin X showed reduced protein expression of pVEGFR2, pAKT1, pMAPKAPK2, pFAK and pERK1 in breast cancer lines and in DMBA induced mammary carcinoma in SD rats showed tumor regression and anti-angiogenesis ability via decrease pVEGFR2 and pAKT1 protein expression. In silico studies also revealed monacolin X ability to bind to crucial amino acid Cys 919 in the active site of VEGFR2 suggesting it to be a potent VEGFR2 inhibitor.

Cancer is one of the leading causes of global deaths and there is an urgent need for the development cancer treatment; targeting VEGFR2 could be one of the promising therapies.

Preclinical characterization of anlotinib, a highly potent and selective vascular endothelial growth factor receptor-2 inhibitor

Abrogating tumor angiogenesis by inhibiting vascular endothelial growth factor receptor‐2 (VEGFR2) has been established as a therapeutic strategy for treating cancer. However, because of their low selectivity, most small molecule inhibitors of VEGFR2 tyrosine kinase show unexpected adverse effects and limited anticancer efficacy. In the present study, we detailed the pharmacological properties of anlotinib, a highly potent and selective VEGFR2 inhibitor, in preclinical models. Anlotinib occupied the ATP‐binding pocket of VEGFR2 tyrosine kinase and showed high selectivity and inhibitory potency (IC₅₀ <1 nmol/L) for VEGFR2 relative to other tyrosine kinases. Concordant with this activity, anlotinib inhibited VEGF‐induced signaling and cell proliferation in HUVEC with picomolar IC₅₀ values. However, micromolar concentrations of anlotinib were required to inhibit tumor cell proliferation directly in vitro. Anlotinib significantly inhibited HUVEC migration and tube formation; it also inhibited microvessel growth from explants of rat aorta in vitro and decreased vascular density in tumor tissue in vivo. Compared with the well‐known tyrosine kinase inhibitor sunitinib, once‐daily oral dose of anlotinib showed broader and stronger in vivo antitumor efficacy and, in some models, caused tumor regression in nude mice. Collectively, these results indicate that anlotinib is a well‐tolerated, orally active VEGFR2 inhibitor that targets angiogenesis in tumor growth, and support ongoing clinical evaluation of anlotinib for a variety of malignancies.

Exploration for novel inhibitors showing back-to-front approach against VEGFR-2 kinase domain (4AG8) employing molecular docking mechanism and molecular dynamics simulations

Background

Angiogenesis is a process of formation of new blood vessels and is an important criteria demonstrated by cancer cells. Over a period of time, these cancer cells infect the other parts of the healthy body by a process called progression. The objective of the present article is to identify a drug molecule that inhibits angiogenesis and progression.

Methods

In this pursuit, ligand based pharmacophore virtual screening was employed, generating a pharmacophore model, Hypo1 consisting of four features. Furthermore, this Hypo1 was validated recruiting, Fischer’s randomization, test set method and decoy set method. Later, Hypo1 was allowed to screen databases such as Maybridge, Chembridge, Asinex and NCI and were further filtered by ADMET filters and Lipinski’s Rule of Five. A total of 699 molecules that passed the above criteria, were challenged against 4AG8, an angiogenic drug target employing GOLD v5.2.2.

Results

The results rendered by molecular docking, DFT and the MD simulations showed only one molecule (Hit) obeyed the back-to-front approach. This molecule displayed a dock score of 89.77, involving the amino acids, Glu885 and Cys919, Asp1046, respectively and additionally formed several important hydrophobic interactions. Furthermore, the identified lead molecule showed interactions with key residues when challenged with CDK2 protein, 1URW.

Conclusion

The lead candidate showed several interactions with the crucial residues of both the targets. Furthermore, we speculate that the residues Cys919 and Leu83 are important in the development of dual inhibitor. Therefore, the identified lead molecule can act as a potential inhibitor for angiogenesis and progression.

Electronic supplementary material

The online version of this article (10.1186/s12885-018-4050-1) contains supplementary material, which is available to authorized users.

Inherent formulation issues of kinase inhibitors

The small molecular Kinase Inhibitor (smKI) drug class is very promising and rapidly expanding. All of these drugs are administered orally. The clear relationship between structure and function has led to drugs with a general low intrinsic solubility. The majority of the commercial pharmaceutical formulations of the smKIs are physical mixtures that are limited by the low drug solubility of a salt form. This class of drugs is therefore characterized by an impaired and variable bioavailability rendering them costly and their therapies suboptimal. New formulations are sparingly being reported in literature and patents. The presented data suggests that continued research into formulation design can help to develop more efficient and cost-effective smKI formulation. Moreover, it may also be of help in the future design of the formulations of new smKIs.

Discovery of a Potential HER2 Inhibitor from Natural Products for the Treatment of HER2-Positive Breast Cancer

Breast cancer is one of the most lethal types of cancer in women worldwide due to the late stage detection and resistance to traditional chemotherapy. The human epidermal growth factor receptor 2 (HER2) is considered as a validated target in breast cancer therapy. Even though a substantial effort has been made to develop HER2 inhibitors, only lapatinib has been approved by the U.S. Food and Drug Administration (FDA). Side effects were observed in a majority of the patients within one year of treatment initiation. Here, we took advantage of bioinformatics tools to identify novel effective HER2 inhibitors. The structure-based virtual screening combined with ADMET (absorption, distribution, metabolism, excretion and toxicity) prediction was explored. In total, 11,247 natural compounds were screened. The top hits were evaluated by an in vitro HER2 kinase inhibition assay. The cell proliferation inhibition effect of identified inhibitors was evaluated in HER2-overexpressing SKBR3 and BT474 cell lines. We found that ZINC15122021 showed favorable ADMET properties and attained high binding affinity against HER2. Moreover, ZINC15122021 showed high kinase inhibition activity against HER2 and presented outstanding cell proliferation inhibition activity against both SKBR3 and BT474 cell lines. Results reveal that ZINC15122021 can be a potential HER2 inhibitor.

In Silico Screening Identifies a Novel Potential PARP1 Inhibitor Targeting Synthetic Lethality in Cancer Treatment

Synthetic lethality describes situations in which defects in two different genes or pathways together result in cell death. This concept has been applied to drug development for cancer treatment, as represented by Poly (ADP-ribose) polymerase (PARPs) inhibitors. In the current study, we performed a computational screening to discover new PARP inhibitors. Among the 11,247 compounds analyzed, one natural product, ZINC67913374, stood out by its superior performance in the simulation analyses. Compared with the FDA approved PARP1 inhibitor, olaparib, our results demonstrated that the ZINC67913374 compound achieved a better grid score (-86.8) and amber score (-51.42). Molecular dynamics simulations suggested that the PARP1-ZINC67913374 complex was more stable than olaparib. The binding free energy for ZINC67913374 was -177.28 kJ/mol while that of olaparib was -159.16 kJ/mol. These results indicated ZINC67913374 bound to PARP1 with a higher affinity, which suggest ZINC67913374 has promising potential for cancer drug development.

Discover natural compounds as potential phosphodiesterase-4B inhibitors via computational approaches

cAMP, intracellular cyclic adenosine monophosphate, is a ubiquitous second messenger that plays a key role in many physiological processes. PDE4B which can reduce the cAMP level by hydrolyzing cAMP to 5'-AMP has become a therapeutic target for the treatment of human diseases such as respiratory disorders, inflammation diseases, neurological and psychiatric disorders. However, the use of currently available PDE4B inhibitors is restricted due to serious side effects caused by targeting PDE4D. Hence, we are attempting to find out subfamily-selective PDE4B inhibitors from natural products, using computer-aided approaches such as virtual screening, docking, and molecular dynamics simulation. Finally, four potential PDE4B-selective inhibitors (ZINC67912770, ZINC67912780, ZINC72320169, and ZINC28882432) were found. Compared to the reference drug (roflumilast), they scored better during the virtual screening process. Binding free energy for them was -317.51, -239.44, -215.52, and -165.77 kJ/mol, better than -129.05 kJ/mol of roflumilast. The pharmacophore model of the four candidate inhibitors comprised six features, including one hydrogen bond donor, four hydrogen bond acceptors, and one aromatic ring feature. It is expected that our study will pave the way for the design of potent PDE4B-selective inhibitors of new drugs to treat a wide variety of diseases such as asthma, COPD, psoriasis, depression, etc.

Molecular Dynamics: New Frontier in Personalized Medicine

The field of drug discovery has witnessed infinite development over the last decade with the demand for discovery of novel efficient lead compounds. Although the development of novel compounds in this field has seen large failure, a breakthrough in this area might be the establishment of personalized medicine. The trend of personalized medicine has shown stupendous growth being a hot topic after the successful completion of Human Genome Project and 1000 genomes pilot project. Genomic variant such as SNPs play a vital role with respect to inter individual's disease susceptibility and drug response. Hence, identification of such genetic variants has to be performed before administration of a drug. This process requires high-end techniques to understand the complexity of the molecules which might bring an insight to understand the compounds at their molecular level. To sustenance this, field of bioinformatics plays a crucial role in revealing the molecular mechanism of the mutation and thereby designing a drug for an individual in fast and affordable manner. High-end computational methods, such as molecular dynamics (MD) simulation has proved to be a constitutive approach to detecting the minor changes associated with an SNP for better understanding of the structural and functional relationship. The parameters used in molecular dynamic simulation elucidate different properties of a macromolecule, such as protein stability and flexibility. MD along with docking analysis can reveal the synergetic effect of an SNP in protein-ligand interaction and provides a foundation for designing a particular drug molecule for an individual. This compelling application of computational power and the advent of other technologies have paved a promising way toward personalized medicine. In this in-depth review, we tried to highlight the different wings of MD toward personalized medicine.

Screening of multi-targeted natural compounds for receptor tyrosine kinases inhibitors and biological evaluation on cancer cell lines, in silico and in vitro

Receptors for growth factors encompass within the superfamily of receptor tyrosine kinases and are known to regulate numerous biological processes including cellular growth, proliferation, metabolism, survival, cell differentiation and apoptosis. These receptors have recently caught the attention of the researchers as an attractive target to combat cancer owing to the evidence suggesting their over-expression in cancer cells. Therefore, we studied receptor-based molecular docking of IR (PDB; 3ETA), IGF1R (PDB; 1K3A), EGFR (PDB; 1M17), VEGFIR (PDB; 3HNG), and VEGFIIR (PDB; 2OH4) against natural compounds. Further, in vitro investigation of the biological effect of lead molecules in an array of cancer cell lines was done. All selected natural compounds were docked with the X-ray crystal structure of selected protein by employing GLIDE (Grid-based Ligand Docking with Energetics) Maestro 9.6. InterBioScreen natural compounds docked with each selected protein molecules by using GLIDE high throughput virtual screening. On the basis of Gscore, we select 20 compounds along with 68 anticancer compounds for GLIDE extra precision molecular docking. It was discovered in this study that compound epigallocatechin gallate (EGCG) yielded magnificent Gscore with IGF1R (PDB; 1K3A) and VEGFIIR (PDB; 2OH4), and protein-ligand interactions are chart out. Effect of EGCG on biological activity such as mRNA expression of selected protein, cell proliferation, oxidative stress, and cell migration was reported after the 48 h treatments in cancer cell lines. The RT-PCR densitometric bands analysis showed that compound EGCG reduced the mRNA expression of IGF1R, VEGFIIR, and mTOR at 80 μM concentration. Moreover, EGCG significantly reduced cell proliferation and ROS generation after 48 h treatments. Our result also indicated a reduction in the potential for cell migration that might show in vivo anti-metastasis activity of EGCG.

Discovery of biphenyl-aryl ureas as novel VEGFR-2 inhibitors. Part 4: exploration of diverse hinge-binding fragments

VEGFR-2 plays an essential role in angiogenesis and is an important target for cancer therapy. A series of biphenyl-aryl ureas were synthesized and evaluated as novel VEGFR-2 inhibitors. The pyridine, methylamine carbonyl pyridine and pivaloyl amide pyridine were introduced as novel hinge binding fragment. The majority of title compounds displayed potent VEGFR-2 inhibition. In particular, L1, L9, W14 and W15 exhibited significant enzymatic inhibitory activity with IC50 values of 0.36nM, 0.22nM, 0.15nM and 0.14nM. Compounds L1, L9 and W15 displayed potent antiproliferative activity against A549 and SMMC-7721 cells. SAR study suggested that incorporation of 3-trifluoromethyl and methylamine carbonyl on terminal pyridine could improve VEGFR-2 inhibitory activity. Molecular docking illustrated that urea moiety formed two critical hydrogen bonds with the DFG residues of VEGFR-2. The results indicated that these biphenyl-aryl ureas could serve as promising lead compounds for further optimization.

Discovery of potential drugs for human-infecting H7N9 virus containing R294K mutation

Background

After the first epidemic wave from February through May 2013, the influenza A (H7N9) virus emerged and has followed a second epidemic wave since June 2013. As of June 27, 2014, the outbreak of H7N9 had caused 450 confirmed cases of human infection, with 165 deaths included. The case-fatality rate of all confirmed cases is about 36%, making the H7N9 virus a significant threat to people’s health. At present, neuraminidase inhibitors are the only licensed antiviral medications available to treat H7N9 infections in humans. Oseltamivir is the most commonly used inhibitor, and it is also a front-line drug for the threatening H7N9. Unfortunately, it has been reported that patients treated with oseltamivir can induce R294K (Arg294Lys) substitution in the H7N9 virus, which is a rare mutation and can reduce the antiviral efficacy of inhibitors. Even worse, deaths caused by such mutation after oseltamivir treatment have already been reported, indicating that the need to find substitutive neuraminidase inhibitors for currently available drugs to treat drug-resistant H7N9 is really pressing.

Materials and methods

First, the structure of H7N9 containing the R294K substitution was downloaded from the Protein Data Bank, and structural information of approved drugs was downloaded from the ZINC (ZINC Is Not Commercial) database. Taking oseltamivir carboxylate as a reference drug, we then filtered these molecules through virtual screening to find out potential inhibitors targeting the mutated H7N9 virus. For further evaluation, we carried out a 14 ns molecular dynamic simulation for each H7N9–drug complex and calculated the binding energy for each candidate drug.

Results

We found five inhibitors that could be candidate drugs for treating the mutated H7N9 virus. Docking poses showed these drugs could bind to the virus effectively, with the contribution of hydrogen bonds and hydrophobic interactions. With regard to the molecular dynamic simulations, receptor–ligand complexes formed by these candidate drugs were more stable than the one formed by oseltamivir carboxylate. The binding energy of oseltamivir carboxylate was −122.4 kJ/mol, while those for these potential inhibitors were −417.5, −404.7, −372.2, −304.3, and −289.9 kJ/mol, much better than the reference drug.

Conclusion

Given the current and future threat of the mutated H7N9 virus, it is urgent that potent drugs and effective antiviral therapeutics be found. Our study therefore is able to complement currently available drugs for influenza A infectors and helps to prevent the ongoing threat of H7N9 virus.