Expansion of the scaffold diversity for the development of highly selective butyrylcholinesterase (BChE) inhibitors: Discovery of new hits through the pharmacophore model generation,virtual screening and molecular dynamics simulation

Butyrylcholinesterase (BChE) is recently considered as a new target for the treatment of Alzheimer’s disease (AD). There is an increasing interest in the development of BChE inhibitors. In the present study, a set of pharmacophore models for BChE was developed and validated. Based on the models, virtual screening was performed on five compound collections, from which seventeen potential hits were retained for biological investigation. In total, eight of these seventeen potential hits showed selective BChE inhibitory activity. Moreover, four compounds displayed IC₅₀ values in sub-micromolar range on eqBChE and three displayed IC₅₀ values < 2 μM on huBChE. The diverse scaffolds of the active compounds provided good starting point further development of selective BChE inhibitors. As far as we concerned, here we disclose the first selective pharmacophore model targeting BChE. The high rate of the model in the identification of active hits indicates it is a valuable tool for the development of selective BChE inhibitors, which may benefit the treatment of AD.     

Fragment-based discovery of hydroxy-indazole-carboxamides as novel small molecule inhibitors of Hsp90

Inhibitors of the Hsp90 molecular chaperone are showing considerable promise as potential molecular therapeutic agents for the treatment of cancer. Here we describe the identification of novel small molecular weight inhibitors of Hsp90 using a fragment based approach. Fragments were selected by docking, tested in a biochemical assay and the confirmed hits were crystallized. Information gained from X-ray structures of these fragments and other chemotypes was used to drive the fragment evolution process. Optimization of these high μM binders resulted in 3-benzylindazole derivatives with significantly improved affinity and anti-proliferative effects in different human cancer cell lines.     

The hexameric structures of human heat shock protein 90

Background

The human 90-kDa heat shock protein (HSP90) functions as a dimeric molecular chaperone. HSP90 identified on the cell surface has been found to play a crucial role in cancer invasion and metastasis, and has become a validated anti-cancer target for drug development. It has been shown to self-assemble into oligomers upon heat shock or divalent cations treatment, but the functional role of the oligomeric states in the chaperone cycle is not fully understood.

Principal Findings

Here we report the crystal structure of a truncated HSP90 that contains the middle segment and the carboxy-terminal domain, termed MC-HSP90. The structure reveals an architecture with triangular bipyramid geometry, in which the building block of the hexameric assembly is a dimer. In solution, MC-HSP90 exists in three major oligomer states, namely dimer, tetramer and hexamer, which were elucidated by size exclusion chromatography and analytical ultracentrifugation. The newly discovered HSP90 isoform HSP90N that lacks the N-terminal ATPase domain also exhibited similar oligomerization states as did MC-HSP90.

Conclusions

While lacking the ATPase domain, both MC-HSP90 and HSP90N can self-assemble into a hexameric structure, spontaneously. The crystal structure of MC-HSP90 reveals that, in addition to the C-terminal dimerization domain, the residue W320 in the M domain plays a critical role in its oligomerization. This study not only demonstrates how the human MC-HSP90 forms a hexamer, but also justifies the similar formation of HSP90N by using 3D modeling analysis.     

Inhibitors of human tyrosyl-DNA phospodiesterase (hTdp1) developed by virtual screening using ligand-based pharmacophores

Human tyrosyl-DNA phosphodiesterase (hTdp1) inhibitors have become a major area of drug research and structure-based design since they have been shown to work synergistically with camptothecin (CPT) and selectively in cancer cells. The pharmacophore features of 14 hTdp1 inhibitors were used as a filter to screen the ChemNavigator iResearch Library of about 27 million purchasable samples. Docking of the inhibitors and hits obtained from virtual screening was performed into a structural model of hTdp1 based on a high resolution X-ray crystal structure of human Tdp1 in complex with vanadate, DNA and a human topoisomerase I (TopI)-derived peptide (PDB code: 1NOP). A total of 46 compounds matching the three-dimensional arrangement of the pharmacophoric features were assayed. Using a high-throughput screening assay, we have identified an 1H-indol-3-yl-acetic acid derivative as a potent Tdp1 inhibitor with an IC₅₀ value of 7.94 μM. The obtained novel chemotype may provide a new scaffold for developing inhibitors of Tdp1.     

Synthesis of novel dual target inhibitors of PARP and HSP90 and their antitumor activities

Poly (ADP-ribose) polymerase (PARP) inhibitors have achieved great success in clinical application, especially for the prolonged survival of cisplatin-sensitive ovarian cancer patients. However, there are still many patients who do not respond to PARP inhibitors. Novel PARP inhibitors with higher activity are urgently needed. Herein we report a series of compounds by molecular hybridization PARP-1 inhibitor Olaparib (Ola) with HSP90 inhibitor C0817 (one curcumin derivative). All synthesized compounds were evaluated for their antiproliferative activity in vitro, and some were further assessed for their inhibitory activities of the PARP enzyme and HSP90 affinity. Our results indicated that compound 4 could bind to HSP90 and cause static quenching, indicating that compound 4 was able to bind to HSP90, moreover, downstream molecular breast cancer 1 (BRAC-1) was reduced. In conclusion, dual target inhibitors of PARP and HSP90 exhibited stronger selective cytotoxicities against cancer.     

Identification of selective MMP-9 inhibitors through multiple e-pharmacophore,ligand-based pharmacophore,molecular docking,and density functional theory approaches

Matrix metalloproteinase-9 (MMP-9) is a significant target for the development of drugs for the treatment of arthritis, CNS disorders, and cancer metastasis. The structure-based and ligand-based methods were used for the virtual screening (VS) of database compounds to obtain potent and selective MMP-9 inhibitors. Experimentally known MMP-9 inhibitors were used to grow up ligand-based three pharmacophore models utilizing Schrodinger suite. The X-ray crystallographic structures of MMP-9 with different inhibitors were used to develop five energy-optimized structure-based (e-pharmacophore) models. All developed pharmacophores were validated and applied to screen the Zinc database. Pharmacophore matched compounds were subjected to molecular docking to retrieve hits with novel scaffolds. The molecules with diverse structures, high docking scores and low binding energies for various crystal structures of MMP-9, were selected as final hits. The Induced fit docking (IFD) analysis provided significant information about the driving of inhibitor to approve a suitable bioactive conformational position in the active site of protein. Since charge transfer reaction occurs during receptor–ligand interaction, therefore, electronic features of hits (ligands) are interesting parameters to explain the binding interactions. Density functional theory (DFT) at B3LYP/6-31G* level was utilized to explore electronic features of hits. The docking study of hits using AutoDock was helpful to establish the binding interactions. The study illustrates that the combined pharmacophore approach is advantageous to identify diverse hits which have better binding affinity to the active site of the enzyme for all possible bioactive conformations. The approach used in the study is worthy to design drugs for other targets.     

Virtual screening using ligand-based pharmacophores for inhibitors of human tyrosyl-DNA phospodiesterase (hTdp1)

Human tyrosyl-DNA phosphodiesterase (hTdp1) inhibitors have become a major area of drug research and structure-based design since they have been shown to work synergistically with camptothecin (CPT) and selectively in cancer cells. The pharmacophore features of 14 hTdp1 inhibitors were used as a filter to screen the ChemNavigator iResearch Library of about 27 million purchasable samples. Docking of the inhibitors and hits obtained from virtual screening was performed into a structural model of hTdp1 based on a high resolution X-ray crystal structure of human Tdp1 in complex with vanadate, DNA and a human topoisomerase I (Top1)-derived peptide (PDB code: 1NOP). We present and discuss in some detail 46 compounds matching the three-dimensional arrangement of the pharmacophoric features. The presented novel chemotypes may provide new scaffolds for developing inhibitors of Tdp1.     

Quinazoline Based HSP90 Inhibitors: Synthesis,Modeling Study and ADME Calculations Towards Breast Cancer Targeting

A new 2-thioquinazolinones series was designed and synthesized as HSP90 inhibitors based on the structure of hit compound VII obtained by virtual screening approach. Their in vitro anti-proliferative activity was evaluated against three human cancer cell lines rich in HSP90 namely; colorectal carcinoma (HCT-116), and cervical carcinoma (Hela), breast carcinoma (MCF-7). Compounds 5a, 5d, 5e and 9h showed a significant broad spectrum anti-proliferative activity against all tested cell lines. They were characterized by potent effect against breast cancer in particular with IC₅₀ of 11.73, 8.56, 7.35 and 9.48 μM, respectively against Doxorubicin (IC₅₀ 4.17 μM). HSP90 ATPase activity inhibition assay were conducted where compound 5d exhibited the best IC₅₀ with 1.58 μM compared to Tanespimycin (IC₅₀ = 2.17 μM). Compounds 5a and 9h showed higher IC₅₀ values of 3.21 and 3.41 μM, respectively. The effects of 5a, 5d and 9h on Her2 (a client proteins of HSP90) and HSP70 were evaluated in MCF-7 cells. All tested compounds were found to reduce Her2 protein expression levels and induce Hsp70 protein expression levels significantly, emphasizing that antibreast cancer effect is a consequence of HSP90 chaperone inhibition. Cell cycle analysis of MCF-7 cells treated with 5d showed cell cycle arrest at G2/M phase 38.89% and pro-apoptotic activity as indicated by annexin V-FITC staining by 22.42%. Molecular docking studies suggested mode of interaction to HSP90 via hydrogen bonding. ADME properties prediction of the active compounds suggested that they could be used as orally absorbed anticancer drug candidates.     

A Combination of Receptor-Based Pharmacophore Modeling & QM Techniques for Identification of Human Chymase Inhibitors

Inhibition of chymase is likely to divulge therapeutic ways for the treatment of cardiovascular diseases, and fibrotic disorders. To find novel and potent chymase inhibitors and to provide a new idea for drug design, we used both ligand-based and structure-based methods to perform the virtual screening(VS) of commercially available databases. Different pharmacophore models generated from various crystal structures of enzyme may depict diverse inhibitor binding modes. Therefore, multiple pharmacophore-based approach is applied in this study. X-ray crystallographic data of chymase in complex with different inhibitors were used to generate four structure–based pharmacophore models. One ligand–based pharmacophore model was also developed from experimentally known inhibitors. After successful validation, all pharmacophore models were employed in database screening to retrieve hits with novel chemical scaffolds. Drug-like hit compounds were subjected to molecular docking using GOLD and AutoDock. Finally four structurally diverse compounds with high GOLD score and binding affinity for several crystal structures of chymase were selected as final hits. Identification of final hits by three different pharmacophore models necessitates the use of multiple pharmacophore-based approach in VS process. Quantum mechanical calculation is also conducted for analysis of electrostatic characteristics of compounds which illustrates their significant role in driving the inhibitor to adopt a suitable bioactive conformation oriented in the active site of enzyme. In general, this study is used as example to illustrate how multiple pharmacophore approach can be useful in identifying structurally diverse hits which may bind to all possible bioactive conformations available in the active site of enzyme. The strategy used in the current study could be appropriate to design drugs for other enzymes as well.     

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

Computational approach to the identification of novel Aurora-A inhibitors

Aurora kinase A has been emerging as a key therapeutic target for the design of anticancer drugs. For the purpose of finding biologically active and novel compounds and providing new ideas for drug-design, we performed virtual screening using commercially available databases. A three-dimensional common feature pharmacophore model was developed with the HipHop program provided in the Catalyst software package, and this model was used as a query for screening the databases. A recursive partitioning (RP) model was developed as a filtering system, which was able to classify active and inactive compounds. Eventually, a step-wise virtual screening procedure was conducted by applying the common feature pharmacophore and the RP model in succession to discover novel potent Aurora-A inhibitors. A total of 68 compounds were selected for testing of their in vitro anticancer activities against various human cancer cell lines. Based on the activity data, we have identified fifteen compounds that warrant further investigation. Several compounds have a high inhibition rate (above 80% at 10 ??M) and a GI₅₀ lower than 5 ??M for the cell lines DU145 and HT29. Enzyme assay for these compounds identified hits with micro molar activity. Compound C11 has the highest activity (IC₅₀ = 5.09 ??M). The hits obtained from this screening scheme could be potential drug candidates after further optimization.     

Synthesis and antitumor activity of 1,2,4-triazoles having 1,4-benzodioxan fragment as a novel class of potent methionine aminopeptidase type II inhibitors

A series of 1,2,4-triazole derivatives containing 1,4-benzodioxan (5a-5q) have been designed, synthesized, structurally determined, and their biological activities were evaluated as potential MetAP2 inhibitors. All the synthesized compounds were first reported. Among the compounds, compound 5k showed the most potent biological activity against HEPG2 cancer cell line (IC(50)=0.81 μM for HEPG2 and IC(50)=0.93 μM for MetAP2), which was comparable to the positive control. Docking simulation by positioning compound 5k into the MetAP2 structure active site was performed to explore the possible binding model. The results of apoptosis and Western-blot assay demonstrated that compound 5k possessed good antitumor activity against HEPG2 cancer cell line. Therefore, compound 5k with potent inhibitory activity in tumor growth inhibition may be a potential antitumor agent against HEPG2 cancer cell.     

Reducing agents affect inhibitory activities of compounds: results from multiple drug targets

High-throughput screening (HTS) of large compound libraries has become a commonly used method for the identification of drug leads, and nonphysiological reducing agents have been widely used for HTS. However, a comparison of the difference in the HTS results based on the choice of reducing agent used and potency comparisons of selected inhibitors has not been done with the physiological reducing agent reduced glutathione (GSH). Here, we compared the effects of three reducing agents-dithiothreitol (DTT), β-mercaptoethanol (β-MCE), and tris(2-carboxyethyl)phosphine (TCEP)-as well as GSH against three drug target proteins. Approximately 100,000 compounds were computationally screened for each target protein, and experimental testing of high-scoring compounds (~560 compounds) with the four reducing agents surprisingly produced many nonoverlapping hits. More importantly, we found that various reducing agents altered inhibitor potency (IC(50)) from approximately 10 μM with one reducing agent to complete loss (IC(50)>200 μM) of inhibitory activity with another reducing agent. Therefore, the choice of reducing agent in an HTS is critical because this may lead to the pursuit of falsely identified active compounds or failure to identify the true active compounds. We demonstrate the feasibility of using GSH for in vitro HTS assays with these three target enzymes.     

Ligand-based modelling followed by synthetic exploration unveil novel glycogen phosphorylase inhibitory leads

Glycogen phosphorylase (GP) is a valid anti-diabetic target. Accordingly, we applied a drug discovery workflow to unveil novel inhibitory GP leads via combining pharmacophore modeling, QSAR analysis and in silico screening, followed by synthetic exploration of active hits. Virtual screening identified six low micromolar inhibitory leads from the National Cancer Institute (NCI) list of compounds. The most potent hits exhibited anti-GP IC(50) values of 3.2 and 4.1 μM. Synthetic exploration of hit 59 (IC(50)=4.1 μM) yielded 25 lead inhibitors with the best illustrating IC(50) of 3.0 μM. Interestingly, we prepared several novel mixed oxalyl amide anti-GP leads employing new chemical reaction involving succinic acid-based adducts.     

Ligand- and structure-based in silico studies to identify kinesin spindle protein (KSP) inhibitors as potential anticancer agents

Kinesin spindle protein (KSP) belongs to the kinesin superfamily of microtubule-based motor proteins. KSP is responsible for the establishment of the bipolar mitotic spindle which mediates cell division. Inhibition of KSP expedites the blockade of the normal cell cycle during mitosis through the generation of monoastral MT arrays that finally cause apoptotic cell death. As KSP is highly expressed in proliferating/cancer cells, it has gained considerable attention as a potential drug target for cancer chemotherapy. Therefore, this study envisaged to design novel KSP inhibitors by employing computational techniques/tools such as pharmacophore modelling, virtual database screening, molecular docking and molecular dynamics. Initially, the pharmacophore models were generated from the data-set of highly potent KSP inhibitors and the pharmacophore models were validated against in house test set ligands. The validated pharmacophore model was then taken for database screening (Maybridge and ChemBridge) to yield hits, which were further filtered for their drug-likeliness. The potential hits retrieved from virtual database screening were docked using CDOCKER to identify the ligand binding landscape. The top-ranked hits obtained from molecular docking were progressed to molecular dynamics (AMBER) simulations to deduce the ligand binding affinity. This study identified MB-41570 and CB-10358 as potential hits and evaluated these experimentally using in vitro KSP ATPase inhibition assays.     

High-throughput screening assay for inhibitors of heat-shock protein 90 ATPase activity

The molecular chaperone heat-shock protein 90 (HSP90) plays a key role in the cell by stabilizing a number of client proteins, many of which are oncogenic. The intrinsic ATPase activity of HSP90 is essential to this activity. HSP90 is a new cancer drug target as inhibition results in simultaneous disruption of several key signaling pathways, leading to a combinatorial approach to the treatment of malignancy. Inhibitors of HSP90 ATPase activity including the benzoquinone ansamycins, geldanamycin and 17-allylamino-17-demethoxygeldanamycin, and radicicol have been described. A high-throughput screen has been developed to identify small-molecule inhibitors that could be developed as therapeutic agents with improved pharmacological properties. A colorimetric assay for inorganic phosphate, based on the formation of a phosphomolybdate complex and subsequent reaction with malachite green, was used to measure the ATPase activity of yeast HSP90. The Km for ATP determined in the assay was 510+/-70 microM. The known HSP90 inhibitors geldanamycin and radicicol gave IC(50) values of 4.8 and 0.9 microM respectively, which compare with values found using the conventional coupled-enzyme assay. The assay was robust and reproducible (2-8% CV) and used to screen a compound collection of approximately 56,000 compounds in 384-well format with Z' factors between 0.6 and 0.8.     

Synthesis, biological evaluation, and molecular docking studies of 1,3,4-oxadiazole derivatives possessing 1,4-benzodioxan moiety as potential anticancer agents

In present study, a series of new 1,3,4-oxadiazole derivatives containing 1,4-benzodioxan moiety (6a-6s) as potential telomerase inhibitors were synthesized. The bioassay tests demonstrated that compounds 6k, 6l, 6m, 6n and 6s exhibited broad-spectrum antitumor activity with IC(50) concentration range from 7.21 μM to 25.87 μM against the four cancer cell lines, HEPG2, HELA, SW1116 and BGC823. Moreover, all the title compounds were assayed for telomerase inhibition using the TRAP-PCR-ELISA assay. The results showed compound 6k possessed the most potent telomerase activity (IC(50)=1.27 ± 0.05 μM). Docking simulation was performed to position compound 6k into the active site of telomerase (3DU6) to determine the probable binding model.     

Discovery of potent inhibitors for interleukin-2-inducible T-cell kinase: structure-based virtual screening and molecular dynamics simulation approaches

In our study, a structure-based virtual screening study was conducted to identify potent ITK inhibitors, as ITK is considered to play an important role in the treatment of inflammatory diseases. We developed a structure-based pharmacophore model using the crystal structure (PDB ID: 3MJ2) of ITK complexed with BMS-50944. The most predictive model, SB-Hypo1, consisted of six features: three hydrogen-bond acceptors (HBA), one hydrogen-bond donor (HBD), one ring aromatic (RA), and one hydrophobic (HY). The statistical significance of SB-Hypo1 was validated using wide range of test set molecules and a decoy set. The resulting well-validated model could then be confidently used as a 3D query to screen for drug-like molecules in a database, in order to retrieve new chemical scaffolds that may be potent ITK inhibitors. The hits retrieved from this search were filtered based on the maximum fit value, drug-likeness, and ADMET properties, and the hits that were retained were used in a molecular docking study to find the binding mode and molecular interactions with crucial residues at the active site of the protein. These hits were then fed into a molecular dynamics simulation to study the flexibility of the activation loop of ITK upon ligand binding. This combination of methodologies is a valuable tool for identifying structurally diverse molecules with desired biological activities, and for designing new classes of selective ITK inhibitors.

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.     

Thiosemicarbazone fragment embedded within 1,2,4-triazole ring as inhibitors of Entamoeba histolytica

A series of 1,2,4-triazole derivatives containing thiosemicarbazone linkage was synthesized and evaluated for their in vitro antiamoebic activity against HM1:IMSS strain of Entamoeba histolytica. All the compounds were capable of inhibiting the growth of E. histolytica out of which four compounds (IC(50)=0.28-1.38 μM) were found to have better efficacy than the standard drug Metronidazole (IC(50)=1.8 μM). Cytotoxicity of the active compounds was assessed by MTT assay using human breast cancer MCF-7 cell line, which revealed that all the compounds were low cytotoxic in the concentration range of 2.5-250 μM.