Discovery of new GSK-3β inhibitors through structure-based virtual screening

Glycogen synthase kinase-3β (GSK-3β) is an attractive therapeutic target for human diseases, such as diabetes, cancer, neurodegenerative diseases, and inflammation. Thus, structure-based virtual screening was performed to identify novel scaffolds of GSK-3β inhibitors, and we observed that conserved water molecules of GSK-3β were suitable for virtual screening. We found 14 hits and D1 (IC₅₀ of 0.71?μM) were identified. Furthermore, the neuroprotection activity of D1–D3 was validated on a cellular level. 2D similarity searches were used to find derivatives of high inhibitory compounds and an enriched structure–activity relationship suggested that these skeletons were worthy of study as potent GSK-3β inhibitors.     

Discovery of novel inhibitors of vascular endothelial growth factor-A–Neuropilin-1 interaction by structure-based virtual screening

Neuropilin-1 (NRP-1), one of the most important co-receptors of vascular endothelial growth factor-A (VEGF-A), increases its angiogenic action in several chronic diseases including cancer by increasing the activity of associated tyrosine kinase receptors, VEGFR1 and VEGFR2. Binding of VEGF-A to NRP-1 plays a critical role in pathological angiogenesis and tumor progression. Today, targeting this interaction is a validated approach to fight against angiogenesis-dependent diseases. Only anti-NRP-1 antibodies, peptide and peptidomimetic drug-candidates or hits have been developed thus far. In order to identify potent orally active small organic molecules various experimental and in silico approaches can be used. Here we report, novel promising small drug-like molecules disrupting the binding of VEGF-A₁₆₅ to NRP-1. We carried out structure-based virtual screening experiments using the ChemBridge compound collection on the VEGF-A₁₆₅ binding pocket of NRP-1. After docking and two rounds of similarity search computations, we identified 4 compounds that inhibit the biotinylated VEGF-A₁₆₅ binding to recombinant NRP-1 with K_i of about 10 μM. These compounds contain a common chlorobenzyloxy alkyloxy halogenobenzyl amine scaffold that can serve as a base for further development of new NRP-1 inhibitors.     

Silibinin: a novel inhibitor of Aβ aggregation

Alzheimer's disease (AD) is characterized by the abnormal aggregation of amyloid β peptide (Aβ) into extracellular fibrillar deposits known as amyloid plaque. Inhibition of Aβ aggregation is therefore viewed as a potential method to halt or slow the progression of AD. It is reported that silibinin (silybin), a flavonoid derived from the herb milk thistle (Silybum marianum), attenuates cognitive deficits induced by Aβ25-35 peptide and methamphetamine. However, it remains unclear whether silibinin interacts with Aβ peptide directly and decreases Aβ peptide-induced neurotoxicity. In the present study, we identified, through employing a ThT assay and electron microscopic imaging that silibinin also appears to act as a novel inhibitor of Aβ aggregation and this effect showed dose-dependency. We also show that silibinin prevented SH-SY5Y cells from injuries caused by Aβ(1-42)-induced oxidative stress by decreasing H(2)O(2) production in Aβ(1-42)-stressed neurons. Taken together, these results indicate that silibinin may be a novel therapeutic agent for the treatment of AD.     

Structure–activity relationship of novel DAPK inhibitors identified by structure-based virtual screening

Death-associated protein kinase (DAPK) is a serine/threonine protein kinase implicated in diverse programmed cell death pathways. DAPK is a promising target protein for the treatment of ischemic diseases. We identified novel potent and selective DAPK inhibitors efficiently by structure-based virtual screening, then further developed the hit compounds. In this paper, we describe the development of the hit compounds and the structure–activity relationship studies of the DAPK inhibitors in detail, including calculation of the solvated interaction energy (SIE), and verification of selectivity using a kinase panel.     

Discovery of atorvastatin as a tetramer stabilizer of nuclear receptor RXRα through structure-based virtual screening

Retinoid X receptor alpha (RXRα), a central member of the nuclear receptor superfamily and a key regulator of many signal transduction pathways, has been an attractive drug target. We previously discovered that an N-terminally truncated form of RXRα can be induced by specific ligands to form homotetramers, which, as a result of conformational selection, forms the basis for inhibiting the nongenomic activation of RXRα. Here, we report the identification and characterization of atorvastatin as a new RXRα tetramer stabilizer by using structure-based virtual screening and demonstrate that virtual library screening can be used to aid in identifying RXRα ligands that can induce its tetramerization. In this study, docking was applied to screen the FDA-approved small molecule drugs in the DrugBank 4.0 collection. Two compounds were selected and purchased for testing. We showed that the selected atorvastatin could bind to RXRα to promote RXRα-LBD tetramerization. We also showed that atorvastatin possessed RXRα-dependent apoptotic effects. In addition, we used a chemical approach to aid in the studies of the binding mode of atorvastatin.     

Discovery of a novel hybrid from finasteride and epristeride as 5α-reductase inhibitor

Finasteride and epristeride both inhibit 5α-reductase with high potency via competitive and non-competitive mechanism, respectively. A new hybrid of finasteride and epristeride was designed as a new 5α-reductase inhibitor based on combination principles in medicinal chemistry. Human 5β-reductase was chosen as a plausible surrogate of 5α-reductase type II and the results indicate that although the hybrid compound possesses the main bulk of epristeride, its inhibitory mechanism is same as of finasteride. The hybrid turned out to be a potent 5α-reductase inhibitor in low IC₅₀ ranges.     

Virtual ligand screening of α-glucosidase: Identification of a novel potent noncarbohydrate mimetic inhibitor

5-Thiazoleacetamide derivatives of AR122 and AR125 were screened as α-glucosidase inhibitors by in silico high-throughput screening from commercial drug-like small compound libraries. Inhibition of α-glucosidase with AR122 and AR125 is time dependent: with no preincubation, AR122 and AR125 are relatively moderate inhibitors, but interestingly, after a 120 min incubation, they were 50-fold more potent (AR122: IC(50)=2.47 μM and AR125: IC(50)=27.1 μM). Plots of ln [residual α-glucosidase activity %] versus preincubation time show a pseudo-first order kinetics for both inhibitors. Through dialysis of enzyme-inhibitor complexes, no activity recovery was shown. These results suggest that AR122 and AR125 constitute a new class of noncarbohydrate mimetic inhibitor with an irreversible mechanism.     

Discovery of a novel covalent non-β-lactam inhibitor of the metallo-β-lactamase NDM-1

The inhibition of metallo-β-lactamases (MBL) can prevent the hydrolysis of β-lactam antibiotics and hence is a promising strategy for the treatment of antibiotic resistant infections. In this study, we present a novel reversible covalent inhibitor of the clinically relevant MBL New Delhi metallo-β-lactamase 1 (NDM-1). Electrospray ionization-mass spectrometry (ESI-MS) and single site directed mutagenesis were used to show that the inhibitor forms a covalent bond with Lys224 in the active site of NDM-1. The inhibitor was further characterized using an enzyme inhibition assay, a surface plasmon resonance (SPR) based biosensor assay and covalent docking. The determined inhibition constant (K_I¹) was 580 nM and the inhibition constant for the initial complex (K_I) was 76 μM. To our knowledge, this inhibitor is the first example for a reversible covalent non-β-lactam inhibitor targeting NDM-1 and a promising starting point for the design of potent covalent inhibitors.     

Screening and identification of a novel class of TGF-β type 1 receptor kinase inhibitor

Transforming growth factor β (TGF-β) type I receptor (activin receptor-like kinase 5, ALK5) has been identified as a promising target for fibrotic diseases. To find a novel inhibitor of ALK5, the authors performed a high-throughput screen of a library of 420,000 compounds using dephosphorylated ALK5. From primary hits of 1521 compounds, 555 compounds were confirmed. In total, 124 compounds were then selected for follow-up based on their unique structures and other properties. Repeated concentration-response testing and final interference assays of the above compounds resulted in the discovery of a structurally novel ALK5 inhibitor (compound 8) (N-(thiophen 2-ylmethyl)-3-(3,4,5 trimethoxyphenyl)imidazo[1,2β]pyridazin 6-amine) with a low IC(50) value of 0.7 μM. Compound 8 also inhibited the TGF-β-induced nuclear translocation of SMAD with an EC(50) value of 0.8 μM. Kinetic analysis revealed that compound 8 inhibited ALK5 via mixed-type inhibition, suggesting that it may bind to ALK5 differently than other published adenosine triphosphate site inhibitors.     

Discovery of highly potent, nonsteroidal 17β-hydroxysteroid dehydrogenase type 1 inhibitors by virtual high-throughput screening

17β-Hydroxysteroid dehydrogenase type 1 (17β-HSD1) catalyzes the formation of the potent proliferation-stimulating hormone estradiol, and it is thus involved in the development of hormone-dependent breast cancer. Due to its high substrate specificity and the known relationships between its overexpression and disease incidence, 17β-HSD1 is considered an attractive target for drug development. Here, we have used structure-based virtual high-throughput screening to successfully identify potent nonsteroidal 17β-HSD1 inhibitors. Computational screening of a drug-like database containing 13 million compounds identified hits with a 2-benzylidenebenzofuran-3(2H)-one scaffold that we show to be highly potent 17β-HSD1 inhibitors. The most potent in the series, compound 1, showed an IC(50) of 45nM in our 17β-HSD1 inhibition assay, and also showed good selectivity for 17β-HSD1 over 17β-HSD2.     

Identification of small molecules that inhibit GSK-3β through virtual screening

Glycogen synthase kinase-3β (GSK-3β) is involved in glycogen metabolism, neuronal cell development, osteoblast differentiation. Small molecule inhibitors of GSK-3β have various therapeutic potential for the treatment of diabetes type II, bipolar disorders, stroke and chronic inflammatory disease.To identify GSK-3β inhibitors with novel scaffold from chemical library, we primarily screened out putative inhibitors through computer modeling and subsequently evaluated the inhibitory activity of selected compounds against GSK-3β by in vitro Z’-LYTE? assay. A series of compound KRMs strongly inhibited phosphorylation of its substrate with IC₅₀ value of approximately 0.5 μM. Also, we demonstrated that KRM-189 and KRM-191 competed with ATP for GSK-3β, leading to decreased V_max and constant K_m with increasing concentrations of ATP as determined from Lineweaver–Berk equation. Moreover, they showed the selectivity for GSK-3β over other kinases with IC₅₀ values of 2 to 10 μM or more Incubation of cells with KRM-191 with highly selective and potent inhibitory activity caused accumulation of β-catenin, downstream of GSK-3β signaling pathway, indicating that small molecule can prevent degradation of β-catenin via GSK-3β inhibition. Our results suggest that modeling in combination with in vitro assays can be used for the identification of selective and potent inhibitors.     

The identification of novel PLC-γ inhibitors using virtual high throughput screening

Phospholipase C-γ (PLC-γ) has been identified as a possible biological target for anticancer drug therapy but suitable inhibitors are lacking. Therefore, in order to identify active compounds (hits) virtual high throughput screening was performed. The crystal structure of the PLC-δ isoform was used as a model docking scaffold since no crystallographic data are available on its γ counterpart. A pilot screen was performed using ～9.2 × 10⁴ compounds, where the robustness of the methodology was tested. This was followed by the main screening effort where ～4.4 × 10⁵ compounds were used. In both cases, plausible compounds were identified (virtual hits) and a selection of these was experimentally tested. The most potent compounds were in the single digit micro-molar range as determined from the biochemical (Flashplate) assay. This translated into ～15 μM in a functional assay in cells. About 30% of the virtual hits showed activity against PLC-γ (IC₅₀ < 50 μM).     

E-pharmacophore-based virtual screening to identify GSK-3β inhibitors

Glycogen synthase kinase-3β (GSK-3β) is a serine/threonine kinase which has attracted significant attention during recent years in drug design studies. The deregulation of GSK-3β increased the loss of hippocampal neurons by triggering apoptosis-mediating production of neurofibrillary tangles and alleviates memory deficits in Alzheimer’s disease (AD). Given its role in the formation of neurofibrillary tangles leading to AD, it has been a major therapeutic target for intervention in AD, hence was targeted in the present study. Twenty crystal structures were refined to generate pharmacophore models based on energy involvement in binding co-crystal ligands. Four common e-pharmacophore models were optimized from the 20 pharmacophore models. Shape-based screening of four e-pharmacophore models against nine established small molecule databases using Phase v3.9 had resulted in 1800 compounds having similar pharmacophore features. Rigid receptor docking (RRD) was performed for 1800 compounds and 20 co-crystal ligands with GSK-3β to generate dock complexes. Interactions of the best scoring lead obtained through RRD were further studied with quantum polarized ligand docking (QPLD), induced fit docking (IFD) and molecular mechanics/generalized Born surface area. Comparing the obtained leads to 20 co-crystal ligands resulted in 18 leads among them, lead1 had the lowest docking score, lower binding free energy and better binding orientation toward GSK-3β. The 50?ns MD simulations run confirmed the stable nature of GSK-3β-lead1 docking complex. The results from RRD, QPLD, IFD and MD simulations confirmed that lead1 might be used as a potent antagonist for GSK-3β.     

Identification of chemically diverse, novel inhibitors of 17β-hydroxysteroid dehydrogenase type 3 and 5 by pharmacophore-based virtual screening

17β-Hydroxysteroid dehydrogenase type 3 and 5 (17β-HSD3 and 17β-HSD5) catalyze testosterone biosynthesis and thereby constitute therapeutic targets for androgen-related diseases or endocrine-disrupting chemicals. As a fast and efficient tool to identify potential ligands for 17βHSD3/5, ligand- and structure-based pharmacophore models for both enzymes were developed. The models were evaluated first by in silico screening of commercial compound databases and further experimentally validated by enzymatic efficacy tests of selected virtual hits. Among the 35 tested compounds, 11 novel inhibitors with distinct chemical scaffolds, e.g. sulfonamides and triazoles, and with different selectivity properties were discovered. Thereby, we provide several potential starting points for further 17β-HSD3 and 17β-HSD5 inhibitor development. Article from the Special issue on Targeted Inhibitors.     

Discovery of aminoheterocycles as a novel β-secretase inhibitor class: pH dependence on binding activity part 1

We have developed a novel series of heteroaromatic BACE-1 inhibitors. These inhibitors interact with the enzyme in a unique fashion that allows for potent binding in a non-traditional paradigm. In addition to the elucidation of their binding profile, we have discovered a pH dependent effect on the binding affinity as a result of the intrinsic pK_a of these inhibitors and the pH of the BACE-1 enzyme binding assay.     

Identification of novel PPARα/γ dual agonists by virtual screening,ADMET prediction and molecular dynamics simulations

PPARα and PPARγ have been the most widely studied Peroxisome proliferator-activated receptor (PPAR) subtypes due to their important roles in regulating glucose, lipids, and cholesterol metabolism. By combining the lowering serum triglyceride levels benefit of PPARα agonists (such as fibrates) with the glycemic advantages of the PPARγ agonists (such as TZD), the dual PPAR agonists approach can both improve the metabolic effects and minimize the side effects caused by either agent alone, and hence, has become a promising strategy for designing effective drugs against type-2 diabetes. In this study, by means of virtual screening, ADMET prediction and molecular dynamics (MD) simulations techniques, one compound-ASN15761007 with high binding score, low toxicity were gained. It was observed by MD simulations that ASN15761007 not only possessed the same function as AZ242 did in activating PPARα and BRL did in activating PPARγ, but also had more favorable conformation for binding to the two receptors. Our results provided an approach to rapidly produce novel PPARα/γ dual agonists which might be a potential lead compound to develop against insulin resistance and hyperlipidemia.     

Identification of new antimalarial leads by use of virtual screening against cytochrome bc₁

Cytochrome bc(1) is a validated drug target in malaria parasites. The spread of Plasmodium falciparum strains resistant to multiple antimalarials emphasizes the urgent need for new drugs. We screened in silico the ZINC and MOE databases, using ligand- and structure-based approaches, to identify new leads for development. The most active compound presented an IC(50) value against cultured P. falciparum of 2 μM and a docking pose consistent with its activity.     

Therapeutic effect of a peptide inhibitor of TGF-β on pulmonary fibrosis

Pulmonary fibrosis encompasses several respiratory diseases characterized by epithelial cell injury, inflammation and fibrosis. Transforming growth factor (TGF)-β1 is one of the main profibrogenic cytokines involved in the pathogenesis of lung fibrosis. It induces fibroblast differentiation into myofibroblasts, which produce high levels of collagen and concomitantly loss of lung elasticity and reduction of the respiratory function. In the present study, we have investigated the effects of P17 (a TGF-β inhibitor peptide) on IMR-90 lung fibroblast differentiation in vitro, as well as on the inhibition of the development of bleomycin-induced pulmonary fibrosis in mice. It was found that in IMR-90 cells, P17 inhibited TGF-β1-induced expression of connective tissue growth factor and α-smooth muscle actin. In vivo, treatment of mice with P17 2days after bleomycin administration decreased lung fibrosis, areas of myofibroblast-like cells and lymphocyte infiltrate. P17 also reduced mRNA expression of collagen type I, fibronectin and the fibronectin splice isoform EDA in the lung, and increased the expression of IFN-γ mRNA. Finally, therapeutic treatment with P17 in mice with already established fibrosis was able to significantly attenuate the progression of lung fibrosis. These results suggest that P17 may be useful in the treatment of pulmonary fibrosis.