Virtual screening and biophysical studies lead to HSP90 inhibitors

Discovery of New Estrogen-Related Receptor α Agonists via a Combination Strategy Based on Shape Screening and Ensemble Docking

Estrogen-related receptor α (ERRα), a member of nuclear receptors (NRs), plays a role in the regulation of cellular energy metabolism and is reported to be a novel potential target for type 2 diabetes therapy. To date, only a few agonists of ERRα have been identified to improve insulin sensitivity and decrease blood glucose levels. Herein, the discovery of novel potent agonists of ERRα determined using a combined virtual screening approach is described. Molecular dynamics (MD) simulations were used to obtain structural ensembles that can consider receptor flexibility. Then, an efficient virtual screening strategy with a combination of similarity search and ensemble docking was performed against the Enamine, SPECS, and Drugbank databases to identify potent ERRα agonists. Finally, a total of 66 compounds were purchased for experimental testing. Biological investigation of promising candidates identified seven compounds that have activity against ERRα with EC50 values ranging from 1.11 to 21.70 μM, with novel scaffolds different from known ERRα agonists until now. Additionally, the molecule GX66 showed micromolar inverse activity against ERRα with an IC50 of 0.82 μM. The predicted binding modes showed that these compounds were anchored in ERRα-LBP via interactions with several residues of ERRα. Overall, this study not only identified the novel potent ERRα agonists or an inverse agonist that would be the promising starting point for further exploration but also demonstrated a successful molecular dynamics-guided approach applicable in virtual screening for ERRα agonists.

Dual Targeting Strategies On Histone Deacetylase 6 (HDAC6) And Heat Shock Protein 90 (Hsp90)

The design of multi-target drugs acting simultaneously on multiple signaling pathways is a growing field in medicinal chemistry, especially for the treatment of complex diseases such as cancer. Histone deacetylase 6 (HDAC6) is an established anticancer drug target involved in tumor cells transformation. Being an epigenetic enzyme at the interplay of many biological processes, HDAC6 has become an attractive target for polypharmacology studies aimed at improving therapeutic efficacy of anticancer drugs. For example, the molecular chaperone Heat shock protein 90 (Hsp90) is a substrate of HDAC6 deacetylation, and several lines of evidence demonstrate that simultaneous inhibition of HDAC6 and Hsp90 promote synergistic antitumor effects on different cancer cell lines, highlighting the potential benefits of developing a single molecule endowed with multi-target activity. This review will summarize the complex interplay between HDAC6 and Hsp90, providing also useful hints for multi-target drug design and discovery approaches in this field. To this end, crystallographic structures of HDAC6 and Hsp90 complexes will be extensively reviewed in the light of discussing binding pockets features and pharmacophore requirements and providing useful guidelines for the design of dual inhibitors. The few examples of multi-target inhibitors obtained so far, mostly based on chimeric approaches, will be summarized and put into context. Finally, the main features of HDAC6 and Hsp90 inhibitors will be compared, and ligand- and structure-based strategies potentially useful for the development of small molecular weight dual inhibitors will be proposed and discussed.

Designing Specific HSP70 Substrate Binding Domain Inhibitor for Perturbing Protein Folding Pathways to Inhibit Cancer Mechanism

BACKGROUND

HSP70 is a survival factor for tumor cells in transformation and in tumor progression as well as in anti-apoptotic response.

OBJECTIVE

Several inhibitors targeting HSP70 ATPase function displayed off-target effects, but PES, which targets the HSP70 substrate binding domain, prevents tumor cell survival prominently. However, PES may not bind HSP70 in the absence of nucleotide. This research aimed to design a unique inhibitor molecule that works both in the presence and absence of nucleotides to amplify inhibition.

METHODS

A set of chimeric coumarine-pyrazole derivatives were determined by in silico techniques and synthesized to elucidate their inhibitory effects. Cell viability experiments displayed KBR1307 as the most efficient inhibitor. A set of characterization experiments were performed, and the results were compared to that of PES agent. Binding constant, ATP hydrolysis rate, and percent aggregation were determined in the presence and absence of inhibitors.

RESULTS

In silico docking experiments showed that only KBR1307 binds the HSP70 substrate binding domain and interacts with cochaperone interface. Binding experiments indicated that KBR1307 binds HSP70 both in the presence and absence of nucleotides, but PES does not. Both inhibitors significantly lower HSP70 ATPase activity and substrate protein disaggregation activity. However, KBR1307 displays a lower IC₅₀ value at the MCF-7 cell line compared to PES. Both inhibitors do not alter HSP70 secondary structure composition and overall stability.

CONCLUSION

KBR1307 effectively inhibits HSP70 compared to PES and provides a promising template for novel anticancer drug development.

Identification of novel Atg3-Atg8 inhibitors using virtual screening for autophagy modulation

A collection of 9050 natural products, their derivatives, and mimetics, was virtually screened against the human Atg3-Atg8 (Atg - autophagy) binding scaffold. By blocking this interaction, the lipidation of Atg8 does not occur and the formation of autophagosomes is inhibited. Forty-three (43) potential ligands were tested using enhanced Green Fluorescent Protein (eGFP) tagged LC3, the human ortholog of Atg8, in MCF7 breast cancer cells. Three hits showed single digit µM IC₅₀ values with AT110, an isoflavone derivative, being the best at 1.2 ± 0.6 µM. Molecular modelling against Atg8 in conjunction with structural activity relationship (SAR) strongly supports the binding to this target. Testing in a panel of cancer cell lines showed little cytotoxic effect as compared to chloroquine. However, same concentration of AT110 was shown to be toxic to young zebrafish embryos. This can be explained in terms of the autophagy process being very active in the zebrafish embryos rendering them susceptible to AT110 whereas in the cancer cells tested the autophagy is not usually active. Nevertheless, AT110 blocks autophagy flux in the zebrafish confirming that the ligand is modulating autophagy. A small molecule non-cytotoxic autophagy inhibitor would open the door for adjunct therapies to bolster many established anticancer drugs, reducing their efficacious concentration thus limiting undesirable site effects. In addition, since many cancer types rely on the autophagy mechanism to survive a therapeutic regime, recurrence can potentially be reduced. The discovery of AT110 is an important step in establishing such an adjunct therapy.

Assay design and development strategies for finding Hsp90 inhibitors and their role in human diseases

Heat shock protein 90 (Hsp90) is a molecular chaperone that facilitates the maturation of its client proteins including protein kinases, transcription factors, and steroid hormone receptors which are structurally and functionally diverse. These client proteins are involved in various cellular signaling pathways, and Hsp90 is implicated in various human diseases including cancer, inflammation, and diseases associated with protein misfolding; thus making Hsp90 a promising target for drug discovery. Some of its client proteins are well-known cancer targets. Instead of targeting these client proteins individually, however, targeting Hsp90 is more practical for cancer drug development. Efforts have been invested in recognizing potential drugs for clinical use that inhibit Hsp90 activity and result in the prevention of Hsp90 client maturation and dampening of subsequent signaling cascades. Here, we discuss current assays and technologies used to find and characterize Hsp90 inhibitors that include biophysical, biochemical, cell-based assays and computational modeling. This review highlights recent discoveries that N-terminal isoform-selective compounds and inhibitors that target the Hsp90 C-terminus that may offer the potential to overcome some of the detriments observed with pan Hsp90 inhibitors. The tools and assays summarized in this review should be used to develop Hsp90-targeting drugs with high specificity, potency, and drug-like properties that may prove immensely useful in the clinic.

Discovery of novel Hsp90 C-terminal domain inhibitors that disrupt co-chaperone binding

Heat shock protein 90 (Hsp90) is an essential molecular chaperone that performs vital stress-related and housekeeping functions in cells and is a current therapeutic target for diseases such as cancers. Particularly, the development of Hsp90 C-terminal domain (CTD) inhibitors is highly desirable as inhibitors that target the N-terminal nucleotide-binding domain may cause unwanted biological effects. Herein, we report on the discovery of two drug-like novel Hsp90 CTD inhibitors by using virtual screening and intrinsic protein fluorescence quenching binding assays, paving the way for future development of new therapies that employ molecular chaperone inhibitors.

Hsp90 chaperones have an energetic hot-spot for binding inhibitors

Although Hsp90-family chaperones have been extensively targeted with ATP-competitive inhibitors, it is unknown whether high affinity is achieved from a few highly stabilizing contacts or from many weaker contacts within the ATP-binding pocket. A large-scale analysis of Hsp90α:inhibitor structures shows that inhibitor hydrogen-bonding to a conserved aspartate (D93 in Hsp90α) stands out as most universal among Hsp90 inhibitors. Here we show that the D93 region makes a dominant energetic contribution to inhibitor binding for both cytosolic and organelle-specific Hsp90 paralogs. For inhibitors in the resorcinol family, the D93:inhibitor hydrogen-bond is pH-dependent because the associated inhibitor hydroxyl group is titratable, rationalizing a linked-protonation event previously observed by the Matulis group. The inhibitor hydroxyl group pKa associated with the D93 hydrogen-bond is therefore critical for optimizing the affinity of resorcinol derivatives, and we demonstrate that spectrophotometric measurements can determine this pKa value. Quantifying the energetic contribution of the D93 hotspot is best achieved with the mitochondrial Hsp90 paralog, yielding 3-6 kcal/mol of stabilization (35-60% of the total binding energy) for a diverse set of inhibitors. The Hsp90 Asp93➔Asn substitution has long been known to abolish nucleotide binding, yet puzzlingly, native sequences of structurally similar ATPases, such as Topoisomerasese II, have an asparagine at this same crucial site. While aspartate and asparagine sidechains can both act as hydrogen bond acceptors, we show that a steric clash prevents the Hsp90 Asp93➔Asn sidechain from adopting the necessary rotamer, whereas this steric restriction is absent in Topoisomerasese II.

Identification of novel inhibitors for the tyrosyl-DNA-phosphodiesterase 1 (Tdp1) mutant SCAN1 using virtual screening

Spinocerebellar ataxia syndrome with axonal neuropathy (SCAN1) is a debilitating neurological disease that is caused by the mutation the Tyrosyl-DNA phosphodiesterase 1 (TDP1) DNA repair enzyme. The crucial His493 in TDP1's binding site is replaced with an arginine amino acid residue rendering the enzyme dysfunctional. A virtual screen was performed against the homology model of SCAN1 and seventeen compounds were identified and tested in a novel SCAN1 specific biochemical assay. Six compounds showed activity with IC₅₀ values between 3.5 and 25.1 µM. The most active ligand 5 (3.5 µM) is a dicoumarin followed by a close structural analogue 6 at 6.0 µM. A less potent series of β-carbolines (14 and 15) was found with potency in the mid-teens. According to molecular modelling an excellent fit for the active ligands into the binding pocket is predicted. To the best of our knowledge, data on inhibitors of the mutant form of TDP1 has not been reported previously. The virtual hits were also tested for wild type TDP1 activity and all six SCAN1 inhibitors are potent for the former, e.g., ligand 5 has a measured IC₅₀ at 99 nM. In the last decade, TDP1 is considered as a promising target for adjuvant therapy against cancer in combination with Topoisomerase 1 poisons. The active ligands are mostly non-toxic to cancer cell lines A-549, T98G and MCF-7 as well as the immortalized WI-38 human fetal lung cells. Furthermore, ligands 5 and 7, show promising synergy in conjunction with topotecan, a clinically used topoisomerase 1 anticancer drug. The active ligands 5, 7, 14 and 15 have a good balance of the physicochemical properties required for oral bioavailability making the excellent candidates for further development.

Identification of Isoform-Selective Ligands for the Middle Domain of Heat Shock Protein 90 (Hsp90)

The molecular chaperone heat shock protein 90 (Hsp90) is a current inhibition target for the treatment of diseases, including cancer. In humans, there are two major cytosolic isoforms of Hsp90 (Hsp90α and Hsp90β). Hsp90α is inducible and Hsp90β is constitutively expressed. Most Hsp90 inhibitors are pan-inhibitors that target both cytosolic isoforms of Hsp90. The development of isoform-selective inhibitors of Hsp90 may enable better clinical outcomes. Herein, by using virtual screening and binding studies, we report our work in the identification and characterisation of novel isoform-selective ligands for the middle domain of Hsp90β. Our results pave the way for further development of isoform-selective Hsp90 inhibitors.

The Development of Tyrosyl-DNA Phosphodiesterase 1 Inhibitors. Combination of Monoterpene and Adamantine Moieties via Amide or Thioamide Bridges

Eleven amide and thioamide derivatives with monoterpene and adamantine substituents were synthesised. They were tested for their activity against the tyrosyl-DNA phosphodiesterase 1 DNA (Tdp1) repair enzyme with the most potent compound 47a, having an IC50 value of 0.64 M. When tested in the A-549 lung adenocarcinoma cell line, no or very limited cytotoxic effect was observed for the ligands. However, in conjunction with topotecan, a well-established Topoisomerase 1 (Top1) poison in clinical use against cancer, derivative 46a was very cytotoxic at 5 M concentration, displaying strong synergism. This effect was only seen for 46a (IC50—3.3 M) albeit some other ligands had better IC50 values. Molecular modelling into the catalytic site of Tdp1 predicted plausible binding mode of 46a, effectively blocking access to the catalytic site.

A Guide to Run Affinity Screens Using Differential Scanning Fluorimetry and Surface Plasmon Resonance Assays

Over the past 30 years, drug discovery has evolved from a pure phenotypic approach to an integrated target-based strategy. The implementation of high-throughput biochemical and cellular assays has enabled the screening of large compound libraries which has become an important and often times the main source of new chemical matter that serve as starting point for medicinal chemistry efforts. In addition, biophysical methods measuring the physical interaction (affinity) between a low molecular weight ligand and a target protein became an integral part of hit validation/optimization to rule out false positives due to assay artifacts. Recent advances in throughput, robustness, and sensitivity of biophysical affinity screening methods have broadened their application in hit identification and validation such that they can now complement classical functional readouts. As a result, new target classes can be accessed that have not been amenable to functional assays. In this chapter, two affinity screening methods, differential scanning fluorimetry and surface plasmon resonance, which are broadly utilized in both academia and pharmaceutical industry are discussed in respect to their use in hit identification and validation. These methods exemplify how assays which differ in complexity, throughput, and information content can support the hit identification/validation process. This chapter focuses on the practical aspects and caveats of these techniques in order to enable the reader to establish their own affinity-based screens in both formats.

Identification of inhibitors of the E. coli chaperone SurA using in silico and in vitro techniques

SurA is a gram-negative, periplasmic chaperone protein involved in the proper folding of outer membrane porins (OMPs), which protect bacteria against toxins in the extracellular environment by selectively regulating the passage of nutrients into the cell. Previous studies demonstrated that deletion of SurA renders bacteria more sensitive to toxins that compromise the integrity of the outer membrane. Inhibitors of SurA will perturb the folding of OMPs, leading to disruption of the outer membrane barrier and making the cell more vulnerable to toxic insults. The discovery of novel SurA inhibitors is therefore of great importance for developing alternative strategies to overcome antibiotic resistance. Our laboratory has screened over 10,000,000 compoundsin silicoby computationally docking these compounds onto the crystal structure of SurA. Through this screen and a screen of fragment compounds (molecular weight less than 250 g/mol), we found twelve commercially readily available candidate compounds that bind to the putative client binding site of SurA. We confirmed binding to SurA by developing and employing a competitive fluorescence anisotropy-based binding assay. Our results show that one of these compounds, Fmoc-β-(2-quinolyl)-d-alanine, binds the client binding site with high micromolar affinity. Using this compound as a lead, we also discovered that Fmoc-l-tryptophan and Fmoc-l-phenylalanine, but not Fmoc-l-tyrosine, bind SurA with similar micromolar affinity. To our knowledge, this is the first report of a competitive fluorescence anisotropy assay developed for the identification of inhibitors of the chaperone SurA, and the identification of three small molecules that bind SurA at its client binding site.

An update on the discovery and development of selective heat shock protein inhibitors as anti-cancer therapy

INTRODUCTION

Over the years, not a single HSP inhibitor has progressed into the post-market phase of drug development despite the success recorded in various pre-clinical and clinical studies. The inability of existing drugs to specifically target oncogenic HSPs has majorly accounted for these setbacks. Recent combinatorial strategies that incorporated computer-aided drug design (CADD) techniques are geared towards the development of highly specific HSP inhibitors with increased activities and minimal toxicities. Areas covered: In this review, strategic therapeutic approaches that have recently aided the development of selective HSP inhibitors were highlighted. Also, the significant contributions of CADD techniques over the years were discussed in detail. This article further describes promising computational paradigms and their applications towards the discovery of highly specific inhibitors of oncogenic HSPs. Expert opinion: The recent shift towards highly selective and specific HSP inhibition has shown great promise as evidenced by the development of paralog/isoform-selective HSP drugs. It could be further augmented with computer-aided drug design strategies, which incorporate reliable methods that would greatly enhance the design and optimization of novel inhibitors with improved activities and minimal toxicities.

The in silico identification of potent anti-cancer agents by targeting the ATP binding site of the N-domain of HSP90

To identify new HSP90 inhibitors, the ATP binding site of the N-domain of HSP90 was targeted by molecular docking of a library of 23,129,083 compounds (from the ZINC database) to the ATP binding site of the N-domain of HSP90. Structure-based virtual screen (SBVS) was performed using idock software on the istar web platform. Based on idock binding energies, 40 molecules were considered as HSP90 inhibitors. In the next step, the 40 molecules and the compound AT13387 (Onalespib) were docked to the XJX binding site using AutoDock Vina software. By comparing the binding energies of the 40 molecules selected with compound AT13387, 26 molecules were selected. By applying the rule of five, eight molecules were selected as hit compounds. The interactions of these eight compounds with the XJX binding site were obtained and investigated, and two-dimensional interaction maps were provided for the others. Finally, computing the toxicity of these compounds with the ProTox-II webserver shows that three compounds, namely ZINC89453765, ZINC23918431 and ZINC12414793, can be considered as good HSP90 inhibitors. These compounds are inactive for nuclear receptor signalling and stress response pathways including heat shock response, so do not have the limitations of common HSP90 inhibitors. They are also inactive for hepatotoxicity, carcinogenicity, immunotoxicity, mutagenicity and cytotoxicity.

Structural ensemble-based docking simulation and biophysical studies discovered new inhibitors of Hsp90 N-terminal domain

Heat shock protein 90 (Hsp90) is one of the most abundant cellular proteins and plays a substantial role in the folding of client proteins. The inhibition of Hsp90 has been regarded as an attractive therapeutic strategy for treating cancer because many oncogenic kinases are Hsp90 client proteins. In this study, we report new inhibitors that directly bind to N-terminal ATP-binding pocket of Hsp90. Optimized structure-based virtual screening predicted candidate molecules, which was followed by confirmation using biophysical and cell-based assays. Among the reported crystal structures, we chose the two structures that show the most favourable early enrichments of true-positives in the receiver operating characteristic curve. Four molecules showed significant changes in the signals of 2D [¹H, ¹⁵N] correlation NMR spectroscopy. Differential scanning calorimetry analysis supported the results indicating direct binding. Quantified dissociation constant values of the molecules, determined by a series of 2D NMR experiments, lie in the range of 0.1–33 μM. Growth inhibition assay with breast and lung cancer cells confirmed the cellular activities of the molecules. Cheminformatics revealed that the molecules share limited chemical similarities with known inhibitors. Molecular dynamics simulations detailed the putative binding modes of the inhibitors.

Development of NMR and thermal shift assays for the evaluation of Mycobacterium tuberculosis isocitrate lyase inhibitors

The enzymes isocitrate lyase (ICL) isoforms 1 and 2 are essential for Mycobacterium tuberculosis survival within macrophages during latent tuberculosis (TB). As such, ICLs are attractive therapeutic targets for the treatment of tuberculosis. However, there are few biophysical assays that are available for accurate kinetic and inhibition studies of ICL in vitro. Herein we report the development of a combined NMR spectroscopy and thermal shift assay to study ICL inhibitors for both screening and inhibition constant (IC50) measurement. Operating this new assay in tandem with virtual high-throughput screening has led to the discovery of several new ICL1 inhibitors.