Identification of potential glycogen kinase-3 inhibitors by structure based virtual screening

Glycogen synthase kinase-3 (GSK3) is a serine/threonine kinase that has attracted much drug discovery attention in recent years. Structural crystallography of the kinase has produced several high resolution inhibitor-GSK3 complexes and this is offering valuable information about the important pharmacophoric features present in the inhibitor, the protein target and the bioactive conformation. The availability of several GSK3-inhibitor co-crystals was successfully exploited to derive a pharmacophore query which retains the all important inhibitor-GSK3 interaction chemistry. A hypothesis containing three features: two hydrogen bond donors and one hydrogen acceptor was found to explain much of the inhibitor-GSK3 interaction. Subsequently, the query has been submitted to three databases for electronic screening. The hits obtained were docked into glycogen synthase kinase-3beta active site. A total of 21 novel potential leads were proposed after thorough examination by a combination of methods: (i) visual examination of how well they dock into the glycogen synthase kinase-3beta binding site, (ii) detailed analysis of their FlexX, G_Score, PMF_Score, ChemScore and D_Score values, (iii) comparative investigation of the docking scores of the hits with that of the thus far reported inhibitors (iv) determination of the binding mode and examination of how the hits retain interactions with the important amino acid residues of the kinase binding site. The hydrophobic heterocycles identified in this investigation are expected to be important additions to the armamentarium of GSK3 hyperactivity antagonism. Further more, the present work may further our current knowledge of the molecular basis of activation, inhibition and regulation of this pharmaceutically important kinase.     

A high-throughput turbidometric assay for screening inhibitors of protein disulfide isomerase activity

Protein disulfide isomerase (PDI) plays a key role in protein folding by catalyzing rearrangements of disulfide bonds in substrate proteins following their synthesis in eukaryotic cells. Besides its major role in the processing and maturation of secretory proteins in the endoplasmic reticulum, this enzyme and its homologs have been implicated in multiple important cellular processes; however, they have not served as targets for the development of therapeutic agents. The authors developed a high-throughput screening assay for PDI and its homologous enzymes in 384-well microplates. The method is based on the enzyme-catalyzed reduction of insulin in the presence of dithiothreitol and measures the aggregation of reduced insulin chains at 650 nm. This kinetic assay was converted to an end-point assay by using hydrogen peroxide as a stop reagent. The feasibility of this high-throughput assay for screening chemical libraries was demonstrated in a pilot screen. The authors show that this homogenous turbidometric assay is robust and cost-effective and can be applied to identify PDI inhibitors from chemical libraries, opening this class of enzymes for therapeutic exploration.     

Identification of small molecule inhibitors of beta-amyloid cytotoxicity through a cell-based high-throughput screening platform

Calpain activation is hypothesized to be an early occurrence in the sequence of events resulting in neurodegeneration, as well as in the signaling pathways linking extracellular accumulation of beta-amyloid (Abeta) peptides and intracellular formation of neurofibrillary tangles. In an effort to identify small molecules that prevent neurodegeneration in Alzheimer's disease by early intervention in the cell death cascade, a cell-based assay in differentiated Sh-SY5Y cells was developed using calpain activity as a read-out for the early stages of death in cells exposed to extracellular Abeta. This assay was optimized for high-throughput screening, and a library of approximately 120,000 compounds was tested. It was expected that the compounds identified as calpain inhibitors would include those that act directly on the enzyme and those that prevented calpain activation by blocking an upstream step in the pathway. In fact, of the compounds that inhibited calpain activation by Abeta with IC(50) values of <10 microM and showed little or no toxicity at concentrations up to 30 microM, none inhibit the calpain enzyme directly.     

A high-throughput turbidometric assay for screening inhibitors of Leishmania major protein disulfide isomerase

The use of a high-throughput technique to perform a pilot screen for Leishmania major protein disulfide isomerase (LmPDI) inhibitors identification is reported. In eukaryotic cells, protein disulfide isomerase (PDI) plays a crucial role in protein folding by catalyzing the rearrangement of disulfide bonds in substrate proteins following their synthesis. LmPDI displays similar domain structure organization and functional properties to other PDI family members and is involved in Leishmania virulence. The authors used a method based on the enzyme-catalyzed reduction of insulin in the presence of dithiothreitol. The screen of a small library of 1920 compounds was performed in a 384-well format and led to the identification of 27 compounds with inhibitory activity against LmPDI. The authors further tested the cytotoxicity of these compounds using Jurkat cells as well as their effect on Leishmania donovani amastigotes using high-content analysis. Results show hexachlorophene and a mixture of theaflavin monogallates inhibit Leishmania multiplication in infected macrophages derived from THP-1 cells, although the inhibitory effect on LmPDI enzymatic activity does not necessarily correlate with the antileishmanial activity.     

Identification of diverse natural products as falcipain-2 inhibitors through structure-based virtual screening

Ten natural compounds are successfully identified as falcipain-2 (FP-2) inhibitors from our in-house natural products database using structure-based virtual screening, which show moderate inhibitory activities against FP-2 with IC₅₀ values ranging from 3.18 to 68.19 μM. While one of the inhibitors (compound 5) also exhibits in vitro antiplasmodial activity against chloroquine sensitive strain (3D7) and chloroquine resistant strain (Dd2) of Plasmodium falciparum in the micromolar range (IC₅₀s = 5.54 μM and 4.05 μM against 3D7 cells and Dd2 cells, respectively). Furthermore, the predicted binding poses are analyzed to explain the structure–activity relationships, which will be helpful for further structural modifications.     

Identification of novel acetylcholinesterase inhibitors through e-pharmacophore-based virtual screening and molecular dynamics simulations

Alzheimer’s disease (AD), a progressive neurodegenerative disorder is the most common cause of dementia among elderly people. To date, the successful therapeutic strategy to treat AD is maintaining the levels of acetylcholine via inhibiting acetylcholinesterase (AChE). The present study involves identification of newer AChE inhibitors by dual approach of e-pharmacophore and structure-based virtual screening of Asinex library. Robustness of docking protocol was validated by enrichment calculation with ROC value .71 and BEDROC value .028. Among 11 selected hits, ZINC72338524 with best MM-GBSA dG binding shows optimal range of CNS properties and ligand–AChE complex stability. Further, molecular dynamics study revealed its molecular interactions with Trp86, Phe338, and Tyr341 amino acid residues of catalytic anionic site and Tyr124, Ser125, and Trp286 amino acid residues of peripheral anionic site. Physicochemical properties and ADMET risk prediction indicates their potential in druggability and safety.     

Discovery of novel non-peptide inhibitors of BACE-1 using virtual high-throughput screening

A novel series of isatin-based inhibitors of β-secretase (BACE-1) have been identified using a virtual high-throughput screening approach. Structure–activity relationship studies revealed structural features important for inhibition. Docking studies suggest these inhibitors may bind within the BACE-1 active site through H-bonding interactions involving the catalytic aspartate residues.     

Structure-based prediction of Mycobacterium tuberculosis shikimate kinase inhibitors by high-throughput virtual screening

A structure-based virtual screening protocol was used to predict Mycobacterium tuberculosis shikimate kinase (MtSk) inhibitors. Docking simulations were performed using eHiTS software and 644 drug-like compounds were identified as potential inhibitors. Forty-two percent of such inhibitors had a structural relationship to a triazole or a tetrazole heteroaromatic system which may provide a candidate lead for the discovery of MtSk inhibitors.     

Identification of novel ASK1 inhibitors using virtual screening

Apoptosis signal-regulating kinase 1 (ASK1, also called MAP3K5) is a mitogen-activated protein kinase kinase kinase (MAP3K) that plays important roles in stress-induced cell death and inflammation, and is expected as a new therapeutic target for cancer, cardiovascular diseases, and neurodegenerative diseases. We identified novel ASK1 inhibitors by virtual screening from the public chemical library collected by Chemical Biology Research Initiative (CBRI) at the University of Tokyo.     

Identification of potential leukocyte antigen-related protein (PTP-LAR) inhibitors through 3D QSAR pharmacophore-based virtual screening and molecular dynamics simulation

Abstract

Owing to its negative regulatory role in insulin signaling, protein tyrosine phosphatase of leukocyte antigen-related protein (PTP-LAR) was widely thought as a potential drug target for diabetes. Now, it was urgent to search for potential LAR inhibitors targeting diabetes. Initially, the pharmacophore models of LAR inhibitors were established with the application of the HypoGen module. The cost analysis, test set validation, as well as Fischer’s test was used to verify the efficiency of pharmacophore model. Then, the best pharmacophore model (Hypo-1-LAR) was applied for the virtual screening of the ZINC database. And 30 compounds met the Lipinski’s rule of five. Among them, 10 compounds with better binding affinity than the known LAR inhibitor (BDBM50296375) were discovered by docking studies. Finally, molecular dynamics simulations and post-analysis experiments (RMSD, RMSF, PCA, DCCM and RIN) were conducted to explore the effect of ligands (ZINC97018474 and Compound 1) on LAR and preliminary understand why ZINC97018474 had better inhibitory activity than Compound 1 (BDBM50296375).

Communicated by Ramaswamy H. Sarma     

Identification of new Hsp90 inhibitors by structure-based virtual screening

Structure-based virtual screening identified pyrimidine-2,4,6-trione and 4H-1,2,4-triazole-3-thiol as novel scaffolds of Hsp90 ATPase inhibitors. Their binding modes in the ATP-binding pocket of Hsp90 were analyzed using AutoDoc program combined with molecular dynamics (MD) simulations.     

Identification of novel BRAF kinase inhibitors with structure-based virtual screening

VRAF murine sarcoma viral oncogene homologue B1 (BRAF) kinase has proved to be a promising target for the development of therapeutics for the treatment of a variety of human cancers. Here, we report the first example of a successful application of the structure-based virtual screening to identify novel BRAF inhibitors. These inhibitors have desirable physicochemical properties as a drug candidate, and compound 1 revealed a submicromolar binding affinity (0.7 μM). Therefore, they may serve as promising lead compounds for further development by structure-activity relationship (SAR) studies to optimize the inhibitory activities. Structural features relevant to the stabilization of the newly identified inhibitors in the ATP-binding site of BRAF are discussed in detail.     

Identification of novel PARP-1 inhibitors by structure-based virtual screening

Poly(ADP-ribose)polymerase-1 (PARP-1) is an abundant and ubiquitous chromatin-bound nuclear protein. PARP-1, a DNA repair enzyme, has been in the limelight as a chemotherapeutic target. In this study, we demonstrated the successful use of structure-based virtual screening to identify inhibitors of PARP-1 from Otava databases comprised of nearly 260,000 compounds. Five novel inhibitors belonging to thienopyrimidinone, isoquinolinoquinazolinone, pyrroloquinazolinone, and cyclopentenothienopyrimidinone scaffolds revealed inhibitory potencies with IC₅₀ values ranged from 9.57 μM to 0.72 μM. Structural features relevant to the activity of these novel compounds within the active site of PARP-1 are discussed in detail and will guide future SAR investigation on these scaffolds.     

Discovery of VEGFR inhibitors through virtual screening and energy assessment

Vascular endothelial growth factor receptor-2 (VEGFR-2) is crucial in promoting tumor angiogenesis and cancer metastasis. Thus, inhibition of VEGFR-2 has appeared as a good tactic for cancer treatment. To find out novel VEGFR-2 inhibitors, first, the PDB structure of VEGFR-2, 6GQO, was selected based on atomic nonlocal environment assessment (ANOLEA) and PROCHECK assessment. 6GQO was then further used for structure-based virtual screening (SBVS) of different molecular databases, including US-FDA approved drugs, US-FDA withdrawn drugs, may bridge, MDPI, and Specs databases using Glide. Based on SBVS, receptor fit, drug-like filters, and absorption, distribution, metabolism, excretion, and toxicity (ADMET) analysis of 427877 compounds, the best 22 hits were selected. From the 22 hits, hit 5 complex with 6GQO was put through molecular mechanics/generalized born surface area (MM/GBSA) study and hERG binding. The MM/GBSA study revealed that hit 5 possesses lesser binding free energy with more inferior stability in the receptor pocket than the reference compound. The VEGFR-2 inhibition assay of hit 5 disclosed an IC₅₀ of 165.23 nM against VEGFR-2, which can be possibly enhanced through structural modifications.     

Identification of inhibitors of N5-carboxyaminoimidazole ribonucleotide synthetase by high-throughput screening

The increasing risk of drug-resistant bacterial infections indicates that there is a growing need for new and effective antimicrobial agents. One promising, but unexplored area in antimicrobial drug design is de novo purine biosynthesis. Recent research has shown that de novo purine biosynthesis in microbes is different from that in humans. The differences in the pathways are centered around the synthesis of 4-carboxyaminoimidazole ribonucleotide (CAIR) which requires the enzyme N⁵-carboxyaminoimidazole ribonucleotide (N⁵-CAIR) synthetase. Humans do not require and have no homologs of this enzyme. Unfortunately, no studies aimed at identifying small-molecule inhibitors of N⁵-CAIR synthetase have been published. To remedy this problem, we have conducted high-throughput screening (HTS) against Escherichia coli N⁵-CAIR synthetase using a highly reproducible phosphate assay. HTS of 48,000 compounds identified 14 compounds that inhibited the enzyme. The hits identified could be classified into three classes based on chemical structure. Class I contains compounds with an indenedione core. Class II contains an indolinedione group, and Class III contains compounds that are structurally unrelated to other inhibitors in the group. We determined the Michaelis–Menten kinetics for five compounds representing each of the classes. Examination of compounds belonging to Class I indicates that these compounds do not follow normal Michaelis–Menten kinetics. Instead, these compounds inhibit N⁵-CAIR synthetase by reacting with the substrate AIR. Kinetic analysis indicates that the Class II family of compounds are non-competitive with both AIR and ATP. One compound in Class III is competitive with AIR but uncompetitive with ATP, whereas the other is non-competitive with both substrates. Finally, these compounds display no inhibition of human AIR carboxylase:SAICAR synthetase indicating that these agents are selective inhibitors of N⁵-CAIR synthetase.     

Novel human mPGES-1 inhibitors identified through structure-based virtual screening

Microsomal prostaglandin E synthase-1 (mPGES-1) is an inducible prostaglandin E synthase after exposure to pro-inflammatory stimuli and, therefore, represents a novel target for therapeutic treatment of acute and chronic inflammatory disorders. It is essential to identify mPGES-1 inhibitors with novel scaffolds as new leads or hits for the purpose of drug design and discovery that aim to develop the next-generation anti-inflammatory drugs. Herein we report novel mPGES-1 inhibitors identified through a combination of large-scale structure-based virtual screening, flexible docking, molecular dynamics simulations, binding free energy calculations, and in vitro assays on the actual inhibitory activity of the computationally selected compounds. The computational studies are based on our recently developed three-dimensional (3D) structural model of mPGES-1 in its open state. The combined computational and experimental studies have led to identification of new mPGES-1 inhibitors with new scaffolds. In particular, (Z)-5-benzylidene-2-iminothiazolidin-4-one is a promising novel scaffold for the further rational design and discovery of new mPGES-1 inhibitors. To our best knowledge, this is the first time a 3D structural model of the open state mPGES-1 is used in structure-based virtual screening of a large library of available compounds for the mPGES-1 inhibitor identification. The positive experimental results suggest that our recently modeled trimeric structure of mPGES-1 in its open state is ready for the structure-based drug design and discovery.     

Identification of novel inhibitors of BCR-ABL tyrosine kinase via virtual screening

Inhibition of BCR-ABL tyrosine kinase activity has shown to be essential for the treatment of chronic myelogenous leukemia (CML). However, drug resistance has quickly arisen in recent clinical trials for STI571 (Gleevec), which is the first approved drug of CML by inhibiting ABL tyrosine kinase. It is desirable to develop new types of ABL tyrosine kinase inhibitors that may overcome this drug resistance problem. Here we present the discovery of novel inhibitors targeted at the catalytic domain of ABL tyrosine kinase by using three-dimensional database searching techniques. From a database containing 200,000 commercially available compounds, the top 1000 compounds with the best DOCK energy score were selected and subjected to structural diversity and drug likeness analysis, 15 compounds were submitted for biological assay. Eight out of the 15 showed inhibitory activity against K562 cells with IC(50) value ranging from 10 to 200 microM. Two promising compounds showed inhibition in further ABL tyrosine phosphorylation assay. It is anticipated that those two compounds can serve as lead compounds for further drug design and optimization.     

Discovery of furan carboxylate derivatives as novel inhibitors of ATP-citrate lyase via virtual high-throughput screening

The enzyme ATP citrate lyase (ACL) catalyzes the formation of cytosolic acetyl CoA, the starting material for de novo lipid and cholesterol biosynthesis. The dysfunction and upregulation of ACL in numerous cancers makes it an attractive target for developing anticancer therapies. ACL inhibition by shRNA knockdown limits cancer cell proliferation and reduces cancer stemness. We designed and implemented a dual docking protocol to select virtual ACL inhibitors that were scored among the top 10 percentiles by both the Autodock Vina and the Glamdock algorithms. Via this in silico screens of a focused furoic acid library, we discovered four subtypes of furans and benzofurans as novel ACL inhibitors. The hit rate of our in silico protocol was 45.8% with 11 of 24 virtual hits confirmed as active in an in vitro ACL enzymatic assay. The IC₅₀ of the most potent ACL inhibitor A1 is 4.1 μM. Our results demonstrated remarkable hit rate by the dual docking approach and provided novel chemical scaffolds for the development of ACL inhibitors for the treatment of cancer.