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.     

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.     

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.     

Combined pharmacophore,virtual screening and molecular dynamics studies to identify Bruton’s tyrosine kinase inhibitors

Bruton’s tyrosine Kinase (BTK) is a cytoplasmic, non-receptor tyrosine kinase expressed in hematopoietic cells. BTK plays a critical role in many cellular signalling pathways making it a potential target to treat autoimmune diseases and cancer. BTK signalling is important for the production of arthritis-associated antibodies, and inhibiting BTK will help the system to block the production of disease-associated antibodies. In this study, we have implemented ligand-based pharmacophore modelling and virtual screening against natural compounds followed by molecular docking, density functional theory and molecular dynamics studies for 50 ns. Four compounds with high affinity towards BTK were identified, and it could be used as a potent lead molecule for designing BTK inhibitor.     

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.     

IKKβ inhibitors identification part I: Homology model assisted structure based virtual screening

Control of NF-κB release through the inhibition of IKKβ has been identified as a potential target for the treatment of inflammatory and autoimmune diseases. We have employed structure based virtual screening scheme to identify lead like molecule from ChemDiv database. Homology models of IKKβ enzyme were developed based on the crystal structures of four kinases. The efficiency of the homology model has been validated at different levels. Docking of known inhibitors library revealed the possible binding mode of inhibitors. Besides, the docking sequence analyses results indicate the responsibility of Glu172 in selectivity. Structure based virtual screening of ChemDiv database has yielded 277 hits. Top scoring 75 compounds were selected and purchased for the IKKβ enzyme inhibition test. From the combined approach of virtual screening followed by biological screening, we have identified six novel compounds that can work against IKKβ, in which 1 compound had highest inhibition rate 82.09% at 10 μM and IC₅₀ 1.76 μM and 5 compounds had 25.35–48.80% inhibition.     

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.     

The inhibitors of cyclin-dependent kinases and GSK-3β enhance osteoclastogenesis

Osteoclasts are multinucleated cells with bone resorption activity that is crucial for bone remodeling. RANK‐RANKL (receptor activator of nuclear factor κB ligand) signaling has been shown as a main signal pathway for osteoclast differentiation. However, the molecular mechanism and the factors regulating osteoclastogenesis remain to be fully understood. In this study, we performed a chemical genetic screen, and identified a Cdks/GSK-3β (cyclin-dependent kinases/glycogen synthase kinase 3β) inhibitor, kenpaullone, and two Cdks inhibitors, olomoucine and roscovitine, all of which significantly enhance osteoclastogenesis of RAW264.7 cells by upregulating NFATc1 (nuclear factor of activated T cells, cytoplasmic 1) levels. We also determined that the all three compounds increase the number of osteoclast differentiated from murine bone marrow cells. Furthermore, the three inhibitors, especially kenpaullone, promoted maturation of cathepsin K, suggesting that the resorption activity of the resultant osteoclasts is also activated. Our findings indicate that inhibition of GSK-3β and/or Cdks enhance osteoclastogenesis by modulating the RANK–RANKL signaling pathway.     

In silico screening and identification of potential GSK3β inhibitors

Abstract

Glycogen synthase kinase-3β (GSK3β) has been reported for its impact on multitude biological processes from cell proliferation to apoptosis. The increase in the ratio of active/inactive GSK3β is the major factor associated in the etiology of several psychiatric diseases, diabetes, muscle hypertrophy, neurodegenerative diseases, and some cancers. These findings made GSK3β a promising therapeutic target, and the interest in the discovery, synthesis of novel drugs to effectively attenuate its function with probably no side effects has been increasing in the chronology of GSK3β drug discovery. In the present study, we applied a combination of computational tools on a chemical library for the virtual discovery of their potency to inhibit GSK3β. The chemical library was screened against a set of filters at different levels. Finally, five compounds in the chemical library were found to potentially inhibit GSK3β with no toxic effects. Furthermore, binding mode analysis revealed that all the compounds bound to the ATP site and most of the hydrogen bonding interactions are conserved as in GSK3β structures deposited in PDB.     

Pharmacophore-based virtual screening of catechol-o-methyltransferase (COMT) inhibitors to combat Alzheimer’s disease

Alzheimer’s disease (AD) is one of the most significant neurodegenerative disorders and its symptoms mostly appear in aged people. Catechol-o-methyltransferase (COMT) is one of the known target enzymes responsible for AD. With the use of 23 known inhibitors of COMT, a query has been generated and validated by screening against the database of 1500 decoys to obtain the GH score and enrichment value. The crucial features of the known inhibitors were evaluated by the online ZINC Pharmer to identify new leads from a ZINC database. Five hundred hits were retrieved from ZINC Pharmer and by ADMET (absorption, distribution, metabolism, excretion, and toxicity) filtering by using FAF-Drug-3 and 36 molecules were considered for molecular docking. From the COMT inhibitors, opicapone, fenoldopam, and quercetin were selected, while ZINC63625100_413 ZINC39411941_412, ZINC63234426_254, ZINC63637968_451, and ZINC64019452_303 were chosen for the molecular dynamics simulation analysis having high binding affinity and structural recognition. This study identified the potential COMT inhibitors through pharmacophore-based inhibitor screening leading to a more complete understanding of molecular-level interactions.     

Molecular dynamics simulation studies of GSK-3β ATP competitive inhibitors: understanding the factors contributing to selectivity

Glycogen synthase kinase-3 is a constitutively acting, multifunctional serine threonine kinase, the role of which has been implicated in several physiological pathways and has emerged as a promising target for the treatment of type-II diabetes and Alzheimer’s disease. In order to provide a detailed understanding of the origin of selectivity determinants of ATP competitive inhibitors, molecular dynamics simulations in combination with MM-PBSA binding energy calculations were performed using crystal structures of GSK-3β and CDK-2 in complex with 12 ATP competitive inhibitors. Analysis of energy contributions indicate that electrostatic interaction energy dictates the selectivity of ATP competitive inhibitors against CDK-2. Key interactions as well as residues that potentially make a major contribution to the binding free energy were identified at the ATP binding site. This analysis stresses the need for the inhibitors to interact with Lys85, Thr138, and Arg141 in the binding site of GSK-3β to show selectivity. The residue-wise energy decomposition analysis further suggested the additional role of Gln185 in determining the selectivity of maleimides. The results obtained in this study can be utilized to design new selective GSK-3 ATP competitive inhibitors.     

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.     

Novel non-peptide β-secretase inhibitors derived from structure-based virtual screening and bioassay

This Letter describes an efficient approach by integrating virtual screening with bioassay technology for finding small organic inhibitors targeting β-secretase (BACE-1). Fifteen hits with inhibitory potencies ranging from 2.8 to 118 μM (IC₅₀) against β-secretase were successfully identified. Compound 12 with IC₅₀ of 2.8 μM is the most potent hit against BACE-1. Docking simulation from gold 3.0 suggests putative binding mode of 12 in BACE-1 and potential key pharmacophore groups for further designing of non-peptide compounds as more powerful inhibitors against BACE-1.     

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

Synthesis of benzimidazole-linked-1,3,4-oxadiazole carboxamides as GSK-3β inhibitors with in vivo antidepressant activity

Recent findings of potential implications of glycogen synthase kinase-3β (GSK-3β) dysfunction in psychiatric disorders like depression, have increased focus for development of GSK-3β inhibitors with possible anti-depressant activity. Keeping this in view, we synthesized a series of benzimidazole-linked-1,3,4-oxadiazole carboxamides and evaluated them for in vitro GSK-3β inhibition. Active compounds were investigated for in vivo antidepressant activity in Wistar rats. Docking studies of active compounds have also been performed. Among nineteen compounds synthesized, compounds 7a, 7r, 7j, and 7d exhibited significant potency against GSK-3β in sub-micromolar range with IC₅₀ values of 0.13 μM, 0.14 μM, 0.20 μM, 0.22 μM respectively and significantly reduced immobility time (antidepressant-like activity) in rats compared to control group. Docking study showed key interactions of these compounds with GSK-3β. These compounds may thus serve as valuable candidates for subsequent development of effective drugs against depression and related disorders.     

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.     

Docking,molecular dynamics,binding energy-MM-PBSA studies of naphthofuran derivatives to identify potential dual inhibitors against BACE-1 and GSK-3β

BACE-1 and GSK-3β both are potential therapeutic drug targets for Alzheimer’s disease. Recently, both these targets received attention for designing dual inhibitors. Till now only two scaffolds (triazinone and curcumin) derivatives have been reported as BACE-1 and GSK-3β dual inhibitors. In our previous work, we have reported first in class dual inhibitor for BACE-1 and GSK-3β. In this study, we have explored other naphthofuran derivatives for their potential to inhibit BACE-1 and GSK-3β through docking, molecular dynamics, binding energy (MM-PBSA). These computational methods were performed to estimate the binding affinity of naphthofuran derivatives towards the BACE-1 and GSK-3β. In the docking results, two derivatives (NS7 and NS9) showed better binding affinity as compared to previously reported inhibitors. Hydrogen bond occupancy of NS7 and NS9 generated from MD trajectories showed good interaction with the flap residues Gln73, Thr72 of BACE-1 and Arg141, Thr138 residues of GSK-3β. MM-PBSA and energy decomposition per residue revealed different components of binding energy and relative importance of amino acid involved in binding. The results showed that the binding of inhibitors was majorly governed by the hydrophobic interactions and suggesting that hydrophobic interactions might be the key to design dual inhibitors for BACE1-1 and GSK-3β. Distance between important pair of amino acid residues indicated that BACE-1 and GSK-3β adopt closed conformation and become inactive after ligand binding. The results suggested that naphthofuran derivatives might act as dual inhibitor against BACE-1 and GSK-3β.     

GSK-3β inhibitors suppressed neuroinflammation in rat cortex by activating autophagy in ischemic brain injury

Previous studies have shown that GSK-3β inhibitor could reduce infarct volume after ischemia brain injury. However, the underlying mechanisms of GSK-3β inhibitor involving neuroprotection remain poorly understood. In the present study, we demonstrated that GSK-3β inhibitor suppressed insult-induced neuroinflammation in rat cortex by increasing autophagy activation in ischemic injury. Male rats were subjected to pMCAO (permanent middle cerebral artery occlusion) followed by treating with SB216763, a GSK-3β inhibitor. We found that insult-induced inflammatory response was significantly decreased by intraperitoneal infusion of SB216763 in rat cortex. A higher level of autophagy was also detected after SB216763 treatment. In the cultured primary microglia, SB216763 activated autophagy and suppressed inflammatory response. Importantly, inhibition of autophagy by Beclin1-siRNA increased inflammatory response in the SB216763-treated microglia. These data suggest that GSK-3β inhibitor suppressed neuroinflammation by activating autophagy after ischemic brain injury, thus offering a new target for prevention of ischemic brain injury.