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.     

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β.     

Description and assessment of a model for GSK-3β database virtual screening

Molecular docking was used in order to prioritize organic syntheses or experimental evaluations. Different GSK-3β protein models were generated in silico from a known X-ray structure. A set of 42 known inhibitors were then flexibly docked into each rigid model and re-scored with various functions, which led to different rankings. The biological activities of the chemicals were then compared to each set of results and one of the rigid models emerged in combination with two scoring functions as giving the best correlation. This methodology constitutes an easy and accurate way to generate reliable models for virtual database screening.     

Identification of novel scaffold of benzothiazepinones as non-ATP competitive glycogen synthase kinase-3β inhibitors through virtual screening

Glycogen synthase kinase-3β (GSK-3β) is an important serine/threonine kinase that has been proved as a key target for neurodegenerative diseases and diabetes. Up to date, most of known inhibitors are bound to the ATP-binding pocket of GSK-3β, which might lead widespread effects due to the high homology between kinases. Recently, some of its non-ATP competitive inhibitors had been confirmed having therapeutical effects owing to their high selectivity. This finding opens a new pathway to study hopeful drugs for treatment of these diseases. However, it is still a challenge nowadays on how to efficiently find non-ATP competitors. Here, we successfully discovered a novel scaffold of benzothiazepinones (BTZs) as selective non-ATP competitive GSK-3β inhibitors through virtual screening approach. A 3D receptor model of substrate binding site of GSK-3β was constructed and applied to screen against drug-like Maybridge database through Autodock program. BTZ compounds were top ranked as efficient hits and were then synthesized for further screening. Among them, the representative compound 4j showed activity to GSK-3β (IC₅₀: 25 μM) in non-ATP competitive mechanism, and nearly no inhibitory effect on other 10 related protein kinases. Overall, the results point out that BTZ compounds might be useful in treatment of Alzheimer’s disease and diabetes mellitus as novel GSK-3β inhibitors. It also suggests, on the other hand, that virtual screening would provide a valuable tool in combination with in vitro assays for the identification of novel selective and potent inhibitors.     

Identification of elements in human long 3′ UTRs that inhibit nonsense-mediated decay

The nonsense-mediated mRNA decay (NMD) pathway serves an important role in gene expression by targeting aberrant mRNAs that have acquired premature termination codons (PTCs) as well as a subset of normally processed endogenous mRNAs. One determinant for the targeting of mRNAs by NMD is the occurrence of translation termination distal to the poly(A) tail. Yet, a large subset of naturally occurring mRNAs contain long 3′ UTRs, many of which, according to global studies, are insensitive to NMD. This raises the possibility that such mRNAs have evolved mechanisms for NMD evasion. Here, we analyzed a set of human long 3′ UTR mRNAs and found that many are indeed resistant to NMD. By dissecting the 3′ UTR of one such mRNA, TRAM1 mRNA, we identified a cis element located within the first 200 nt that inhibits NMD when positioned in downstream proximity of the translation termination codon and is sufficient for repressing NMD of a heterologous reporter mRNA. Investigation of other NMD-evading long 3′ UTR mRNAs revealed a subset that, similar to TRAM1 mRNA, contains NMD-inhibiting cis elements in the first 200 nt. A smaller subset of long 3′ UTR mRNAs evades NMD by a different mechanism that appears to be independent of a termination-proximal cis element. Our study suggests that different mechanisms have evolved to ensure NMD evasion of human mRNAs with long 3′ UTRs.     

Identification of small molecule inhibitors of neurite loss induced by Aβ peptide using high content screening

Multiple lines of evidence indicate a strong relationship between Αβ peptide-induced neurite degeneration and the progressive loss of cognitive functions in Alzheimer disease (AD) patients and in AD animal models. This prompted us to develop a high content screening assay (HCS) and Neurite Image Quantitator (NeuriteIQ) software to quantify the loss of neuronal projections induced by Aβ peptide neurons and enable us to identify new classes of neurite-protective small molecules, which may represent new leads for AD drug discovery. We identified thirty-six inhibitors of Aβ-induced neurite loss in the 1,040-compound National Institute of Neurological Disorders and Stroke (NINDS) custom collection of known bioactives and FDA approved drugs. Activity clustering showed that non-steroidal anti-inflammatory drugs (NSAIDs) were significantly enriched among the hits. Notably, NSAIDs have previously attracted significant attention as potential drugs for AD; however their mechanism of action remains controversial. Our data revealed that cyclooxygenase-2 (COX-2) expression was increased following Aβ treatment. Furthermore, multiple distinct classes of COX inhibitors efficiently blocked neurite loss in primary neurons, suggesting that increased COX activity contributes to Aβ peptide-induced neurite loss. Finally, we discovered that the detrimental effect of COX activity on neurite integrity may be mediated through the inhibition of peroxisome proliferator-activated receptor γ (PPARγ) activity. Overall, our work establishes the feasibility of identifying small molecule inhibitors of Aβ-induced neurite loss using the NeuriteIQ pipeline and provides novel insights into the mechanisms of neuroprotection by NSAIDs.     

Identification of 2-(4-pyridyl)thienopyridinones as GSK-3β inhibitors

The discovery of a novel series of 2-(4-pyridyl)thienopyridinone GSK-3β inhibitors is reported. X-ray crystallography reveals its binding mode and enables rationalization of the SAR. The initial optimization of the template for improved cellular activity and predicted CNS penetration is also presented.     

Identification of new γ-hydroxybutenolides that preferentially inhibit the activity of mPGES-1

Microsomal prostaglandin E(2) synthase-1 (mPGES-1) has been recognized as novel, promising drug target for anti-inflammatory and anticancer drugs. mPGES-1 catalyzes the synthesis of the inducible prostaglandin E(2) in response to pro-inflammatory stimuli, rendering this enzyme extremely interesting in drug discovery process owing to the drastic reduction of the severe side effects typical for traditional non-steroidal anti-inflammatory drugs. In the course of our investigations focused on this topic, we identified two interesting molecules bearing the γ-hydroxybutenolide scaffold which potently inhibit the activity of mPGES-1. Notably, the lead compound 2c that inhibited mPGES-1 with IC(50)=0.9μM, did not affect other related enzymes within the arachidonic acid cascade.     

Discovery of a novel IKK-β inhibitor by ligand-based virtual screening techniques

Various inflammatory stimuli that activate the nuclear factor kappa B (NF-κB) signaling pathway converge on a serine/threonine kinase that displays a key role in the activation of NF-κB: the I kappa B kinase β (IKK-β). Therefore, IKK-β is considered an interesting target for combating inflammation and cancer. In our study, we developed a ligand-based pharmacophore model for IKK-β inhibitors. This model was employed to virtually screen commercial databases, giving a focused hit list of candidates. Subsequently, we scored by molecular shape to rank and further prioritized virtual hits by three-dimensional shape-based alignment. One out of ten acquired and biologically tested compounds showed inhibitory activity in the low micromolar range on IKK-β enzymatic activity in vitro and on NF-κB transactivation in intact cells. Compound 8 (2-(1-adamantyl)ethyl 4-[(2,5-dihydroxyphenyl)methylamino]benzoate) represents a novel chemical class of IKK-β inhibitors and shows that the presented model is a valid approach for identification and development of new IKK-β ligands.     

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.     

Discovery of a natural PI3Kδ inhibitor through virtual screening and biological assay study

Phosphoinositide-3-kinase-δ (PI3Kδ) is a key regulator in the process of IgE mediated mast cell degranulation, which directly induces allergic diseases, such as asthma. This study is aimed at discovery of natural PI3Kδ inhibitors from Chinese medicine and evaluating their anti-mast cell degranulation activity. A combined virtual screening based on 3D pharmacophore model and molecular docking was used to screen for bioactive ingredients directly targeting PI3Kδ. Then, an in vitro kinase inhibition assay was conducted to evaluate the PI3Kδ inhibitory activity of the virtual screening hits. Subsequently, a β-hexosaminidase release assay was performed to verify the anti-mast cell degranulation activity of the active compounds. Finally, ginkgoneolic acid was identified as a PI3Kδ inhibitor (IC₅₀?=?2.49?μM) and exhibited anti-mast cell degranulation activity in vitro (IC₅₀?=?2.40?μM). Docking studies showed that Glu826, Val827 and Val828 were key amino acid residues for PI3Kδ inhibitory activity. Ginkgoneolic acid may be a potential lead compound for developing effective and safe PI3Kδ-inhibiting drugs.     

Discovery of potent and bioavailable GSK-3β inhibitors

Here we report on the discovery of a series of maleimides which have high potency and good selectivity for GSK-3β. The incorporation of polar groups afforded compounds with good bioavailability. The most potent compound 34 has an IC₅₀ of 0.6 nM for GSK-3β, over 100-fold selectivity against a panel of other kinases, and shows efficacy in rat osteoporosis models. The X-ray structure of GSK-3β protein with 34 bound revealed the binding mode of the template and provided insights for future optimization opportunities.     

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.     

Beryllium is an inhibitor of cellular GSK-3β that is 1,000-fold more potent than lithium

Glycogen synthase kinase 3β (GSK-3β) is a key regulator in signaling networks that control cell proliferation, metabolism, development, and other processes. Lithium chloride is a GSK-3 family inhibitor that has been a mainstay of in vitro and in vivo studies for many years. Beryllium salt has the potential to act as a lithium-like inhibitor of GSK-3, but it is not known whether this agent is effective under physiologically relevant conditions. Here we show that BeSO₄ inhibits endogenous GSK-3β in cultured human cells. Exposure to 10 µM Be²⁺ produced a decrease in GSK-3β kinase activity that was comparable to that produced by 10 mM Li⁺, indicating that beryllium is about 1,000-fold more potent than the classical inhibitor when treating intact cells. There was a statistically significant dose-dependent reduction in specific activity of GSK-3β immunoprecipitated from cells that had been treated with either agent. Lithium inhibited GSK-3β kinase activity directly, and it also caused GSK-3β in cells to become phosphorylated at serine-9 (Ser-9), a post-translational modification that occurs as part of a well-known positive feedback loop that suppresses the kinase activity. Beryllium also inhibited the kinase directly, but unlike lithium it had little effect on Ser-9 phosphorylation in the cell types tested, suggesting that alternative modes of feedback inhibition may be elicited by this agent. These results indicate that beryllium, like lithium, can induce perturbations in the GSK-3β signaling network of treated cells.     

Role of GSK-3β in the osteogenic differentiation of palatal mesenchyme

Introduction

The function of Glycogen Synthase Kinases 3β (GSK-3β) has previously been shown to be necessary for normal secondary palate development. Using GSK-3ß null mouse embryos, we examine the potential coordinate roles of Wnt and Hedgehog signaling on palatal ossification.

Methods

Palates were harvested from GSK-3β, embryonic days 15.0–18.5 (e15.0–e18.5), and e15.5 Indian Hedgehog (Ihh) null embryos, and their wild-type littermates. The phenotype of GSK-3β null embryos was analyzed with skeletal whole mount and pentachrome stains. Spatiotemporal regulation of osteogenic gene expression, in addition to Wnt and Hedgehog signaling activity, were examined in vivo on GSK-3β and Ihh +/+ and −/− e15.5 embryos using in situ hybridization and immunohistochemistry. To corroborate these results, expression of the same molecular targets were assessed by qRT-PCR of e15.5 palates, or e13.5 palate cultures treated with both Wnt and Hedgehog agonists and anatagonists.

Results

GSK-3β null embryos displayed a 48 percent decrease (*p<0.05) in palatine bone formation compared to wild-type littermates. GSK-3β null embryos also exhibited decreased osteogenic gene expression that was associated with increased Wnt and decreased Hedgehog signaling. e13.5 palate culture studies demonstrated that Wnt signaling negatively regulates both osteogenic gene expression and Hedgehog signaling activity, while inhibition of Wnt signaling augments both osteogenic gene expression and Hedgehog signaling activity. In addition, no differences in Wnt signaling activity were noted in Ihh null embryos, suggesting that canonical Wnt may be upstream of Hedgehog in secondary palate development. Lastly, we found that GSK-3β −/− palate cultures were “rescued” with the Wnt inhibitor, Dkk-1.

Conclusions

Here, we identify a critical role for GSK-3β in palatogenesis through its direct regulation of canonical Wnt signaling. These findings shed light on critical developmental pathways involved in palatogenesis and may lead to novel molecular targets to prevent cleft palate formation.     

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.     

Curcumin ameliorates CKD-induced mitochondrial dysfunction and oxidative stress through inhibiting GSK-3β activity

Curcumin has been reported to attenuate muscle atrophy. However, the underling mechanism remains unclear. The aim of this study was to investigate whether curcumin could improve chronic kidney disease (CKD)-induced muscle atrophy and mitochondrial dysfunction by inhibiting glycogen synthase kinase-3β (GSK-3β) activity. The sham and CKD mice were fed either a control diet or an identical diet containing 0.04% curcumin for 12 weeks. The C2C12 myotubes were treated with H₂O₂ in the presence or absence of curcumin. In addition, wild-type and muscle-specific GSK-3β knockout (KO) CKD model mice were made by 5/6 nephrectomy, and the sham was regarded as control. Curcumin could exert beneficial effects, including weight maintenance and improved muscle function, increased mitochondrial biogenesis, alleviated mitochondrial dysfunction by increasing adenosine triphosphate levels, activities of mitochondrial electron transport chain complexes and basal mitochondrial respiration and suppressing mitochondrial membrane potential. In addition, curcumin modulated redox homeostasis by increasing antioxidant activity and suppressed mitochondrial oxidative stress. Moreover, the protective effects of curcumin had been found to be mediated via inhibiting GSK-3β activity in vitro and in vivo. Importantly, GSK-3β KO contributed to improved mitochondrial function, attenuated mitochondrial oxidative damage and augmented mitochondrial biogenesis in muscle of CKD. Overall, this study suggested that curcumin alleviated CKD-induced mitochondrial oxidative damage and mitochondrial dysfunction via inhibiting GSK-3β activity in skeletal muscle.     

Identification of new benzamide inhibitor against α-subunit of tryptophan synthase from Mycobacterium tuberculosis through structure-based virtual screening,anti-tuberculosis activity and molecular dynamics simulations

Multi-drug-resistant tuberculosis and extensively drug-resistant tuberculosis has emerged as global health threat, causing millions of deaths worldwide. Identification of new drug candidates for tuberculosis (TB) by targeting novel and less explored protein targets will be invaluable for antituberculosis drug discovery. We performed structure-based virtual screening of eMolecules database against a homology model of relatively unexplored protein target: the α-subunit of tryptophan synthase (α-TRPS) from Mycobacterium tuberculosis essential for bacterial survival. Based on physiochemical properties analysis and molecular docking, the seven candidate compounds were selected and evaluated through whole cell-based activity against the H37Rv strain of M. tuberculosis. A new Benzamide inhibitor against α-subunit of tryptophan synthase (α-TRPS) from M. tuberculosis has been identified causing 100% growth inhibition at 25 μg/ml and visible bactericidal activity at 6 μg/ml. This benzamide inhibitor displayed a good predicted binding score (?48.24 kcal/mol) with the α-TRPS binding pocket and has logP value (2.95) comparable to Rifampicin. Further refinement of docking results and evaluation of inhibitor-protein complex stability were investigated through Molecular dynamic (MD) simulations studies. Following MD simulations, Root mean square deviation, Root mean square fluctuation and secondary structure analysis confirmed that protein did not unfold and ligand stayed inside the active pocket of protein during the explored time scale. This identified benzamide inhibitor against the α-subunit of TRPS from M. tuberculosis could be considered as candidate for drug discovery against TB and will be further evaluated for enzyme-based inhibition in future studies.