Pharmacophore modelling, validation, 3D virtual screening, docking, design and in silico ADMET simulation study of histone deacetylase class-1 inhibitors

Design, synthesis, and pharmacological evaluation of indazole carboxamides of N-substituted pyrrole derivatives as soybean lipoxygenase inhibitors

In this paper, we attempted to develop a novel class of compounds against lipoxygenase, a key enzyme in the biosynthesis of leukotrienes implicated in a series of inflammatory diseases. Given the absence of appropriate human 5-lipoxygenase crystallographic data, solved soybean lipoxygenase-1 and -3 structures were used as a template to generate an accurate pharmacophore model which was further used for virtual screening purposes. Eight compounds (1-8) have been derived from the in-house library consisting of N-substituted pyrroles conjugated with 5- or 6-indazole moieties through a carboxamide linker. This study led to the discovery of hit molecule 8 bearing a naphthyl group with the IC50 value of 22 μM according to soybean lipoxygenase in vitro assay. Isosteric replacement of naphthyl ring with quinoline moieties and reduction of carbonyl carboxamide group resulted in compounds 9-12 and 13, respectively. Compound 12 demonstrated the most promising enzyme inhibition. In addition, compounds 8 and 12 were found to reduce the carrageenan-induced paw edema in vivo by 52.6 and 49.8%, respectively. In view of the encouraging outcomes concerning their notable in vitro and in vivo anti-inflammatory activities, compounds 8 and 12 could be further optimized for the discovery of novel 5-lipoxygenase inhibitors in future. A structure-based 3D pharmacophore model was used in the virtual screening of in-house library to discover novel potential 5-lipoxygenase inhibitors.

In Silico Exploration of Novel Tubulin Inhibitors: A Combination of Docking and Molecular Dynamics Simulations, Pharmacophore Modeling, and Virtual Screening

Microtubules play a critical role in mitosis and cell division and are regarded as an excellent target for anticancer therapy. Although microtubule-targeting agents have been widely used in the clinical treatment of different human cancers, their clinical application in cancer therapy is limited by both intrinsic and acquired drug resistance and adverse toxicities. In a previous work, we synthesized compound 9IV-c, ((E)-2-(3,4-dimethoxystyryl)-6,7,8-trimethoxy-N-(3,4,5-trimethoxyphenyl)quinoline-4-amine) that showed potent activity against multiple human tumor cell lines, by targeting spindle formation and/or the microtubule network. Accordingly, in this study, to identify potent tubulin inhibitors, at first, molecular docking and molecular dynamics studies of compound 9IV-c were performed into the colchicine binding site of tubulin; then, a pharmacophore model of the 9IV-c-tubulin complex was generated. The pharmacophore model was then validated by Güner-Henry (GH) scoring methods and receiver operating characteristic (ROC) analysis. The IBScreen database was searched by using this pharmacophore model as a screening query. Finally, five retrieved compounds were selected for molecular docking studies. These efforts identified two compounds (b and c) as potent tubulin inhibitors. Investigation of pharmacokinetic properties of these compounds (b and c) and compound 9IV-c displayed that ligand b has better drug characteristics compared to the other two ligands.

Discovery of novel isoform-selective histone deacetylases 5 and 9 inhibitors through combined ligand-based pharmacophore modeling, molecular mocking, and molecular dynamics simulations for cancer treatment

Class IIa histone deacetylases (HDACs) 5 and 9 play crucial roles in several human disorders such as cancer, making them important targets for drug design. Continuous research is pursed to overcome the cytotoxicity side effect that comes with the currently available broad-spectrum HDACs inhibitors. Herein, common features of active HDACs inhibitors in clinical trials and use have been calculated to generate the best pharmacophore hypothesis. Guner-Henry scoring system was used to validate the generated hypotheses. Hypo1 of HDAC5 and Hypo2 of HDAC9 exhibited the most statistically significance hypotheses. Compounds with fit value of 3 and more were examined by QuickVina 2 docking tool to calculate their binding affinity toward all class IIa HDACs. A total of 6 potential selective compounds were subjected to 100 molecular dynamics (MD) simulation to examine their binding modes. The free binding energy calculations were computed according to the MM-PBSA method. Proposed selective compounds displayed good stability with their targets and thus they may offer potent leads for the designing of HDAC5 and HDAC9 isoform selective inhibitors.

Virtual Screening of Natural Compounds as Potential PI3K-AKT1 Signaling Pathway Inhibitors and Experimental Validation

A computational screening for natural compounds suitable to bind the AKT protein has been performed after the generation of a pharmacophore model based on the experimental structure of AKT1 complexed with IQO, a well-known inhibitor. The compounds resulted as being most suitable from the screening have been further investigated by molecular docking, ADMET (Absorption, Distribution, Metabolism, Excretion, and Toxicity) analysis and toxicity profiles. Two compounds selected at the end of the computational analysis, i.e., ZINC2429155 (also named STL1) and ZINC1447881 (also named AC1), have been tested in an experimental assay, together with IQO as a positive control and quercetin as a negative control. Only STL1 clearly inhibited AKT activation negatively modulating the PI₃K/AKT pathway.

Molecular dynamics guided insight, binding free energy calculations and pharmacophore-based virtual screening for the identification of potential VEGFR2 inhibitors

Vascular endothelial growth factor-A (VEGF-A) is a crucial member of the Vascular endothelial growth factor (VEGF) family which mediates the metastasis of tumor by 'angiogenic switch'. Therefore, targeting a VEGF-A mediated VEGFR2 signaling pathway is the most promising approach to repress the angiogenesis of tumor cells. VEGFR2 inhibitors are two types: Type I and Type II. Type II inhibitors have more chemical space to exploit and have better selectivity because of allosteric binding pocket over type I inhibitors. Hence, The present study encompasses identification of potential type II VEGFR2 inhibitors employing pharmacophore based virtual screnning. In this study, ten five featured pharmacophore model were generated from a dataset of 39 biaryl urea analogs.Out of all, ADDRR_1 pharmacophore model were used to screen the library of 5.2 million compounds retrieved from NCI, Maybridge, Asinex and Zinc databases. 7000 hits were filtered out from the pharmacophore-based virtual screening based on the phase fitness score. Among all best ten hits were identified employing extra precision mode of GLIDE module. ZINC00759038 and 211246 were chosen as top hits based on docking score, free binding energy, and ADME profile. They were subjected to molecular-dynamic studies to assess the hits-VEGFR2 binding stability. It suggests that ZINC00759038-VEGFR2 and 211246-VEGFR2 complexes are quite stable for the 20 ns simulation period. The strength of hit-protein complexes were further assessed by thermodynamic analysis of MD simulation studies by MMGBSA. Interestingly, these hits retains 90% similarity with standard VEGFR2 inhibitor (Sorafenib). Hence, these identified hits may led to new lead compounda as VEGFR2 inhibitors.