Exploiting conformational dynamics in drug discovery: design of C-terminal inhibitors of Hsp90 with improved activities

Computing allostery: from the understanding of biomolecular regulation and the discovery of cryptic sites to molecular design

The concept of allostery has become a central tenet in the study of biological systems. In parallel, the discovery of allosteric drugs is generating new opportunities to selectively modulate difficult targets involved in pathologic mechanisms. Molecular simulations can provide atomistically detailed insight into the processes involved in allosteric regulation and signaling, and at the same time, they have the potential to unveil regulatory hotspots or cryptic sites that are not immediately evident from the analysis of static structures. In this context, computational approaches should be able to connect the study of allosteric regulation at different scales to the possibility of leveraging this knowledge to expand the chemical space of new, active drugs. Here, we will discuss recent advances in the study of allosteric regulation via computational methods and connect the mechanistic knowledge generated to the possibility of designing new small molecules that can tweak the functions of their receptors.

Small-molecule dual inhibitors targeting heat shock protein 90 for cancer targeted therapy

Heat shock protein 90, also known as Hsp90, is an extensively preserved molecular chaperone that performs a critical function in organizing various biological pathways and cellular operations. As a potential drug target, Hsp90 is closely linked to cancer. Hsp90 inhibitors are a class of drugs that have been extensively studied in preclinical models and have shown promise in a variety of diseases, especially cancer. However, Hsp90 inhibitors have encountered several challenges in clinical development, such as low efficacy, toxicity, or drug resistance, few Hsp90 small molecule inhibitors have been approved worldwide. Nonetheless, combining Hsp90 inhibitors with other tumor inhibitors, such as HDAC inhibitors, tubulin inhibitors, and Topo II inhibitors, has been shown to have synergistic antitumor effects. Consequently, the development of Hsp90 dual-target inhibitors is an effective strategy in cancer treatment, as it enhances potency while reducing drug resistance. This article provides an overview of Hsp90's domain structure and biological functions, as well as a discussion of the design, discovery, and structure-activity relationships of Hsp90 dual inhibitors, aiming to provide insights into clinical drug research from a medicinal chemistry perspective and discover novel Hsp90 dual inhibitors.

Insight into the Nucleotide Based Modulation of the Grp94 Molecular Chaperone Using Multiscale Dynamics

Grp94, an ER-localized molecular chaperone, is required for the folding and activation of many membrane and secretory proteins. Client activation by Grp94 is mediated by nucleotide and conformational changes. In this work, we aim to understand how microscopic changes from nucleotide hydrolysis can potentiate large-scale conformational changes of Grp94. We performed all-atom molecular dynamics simulations on the ATP-hydrolysis competent state of the Grp94 dimer in four different nucleotide bound states. We found that Grp94 was the most rigid when ATP was bound. ATP hydrolysis or nucleotide removal enhanced mobility of the N-terminal domain and ATP lid, resulting in suppression of interdomain communication. In an asymmetric conformation with one hydrolyzed nucleotide, we identified a more compact state, similar to experimental observations. We also identified a potential regulatory role of the flexible linker, as it formed electrostatic interactions with the Grp94 M-domain helix near the region where BiP is known to bind. These studies were complemented with normal-mode analysis of an elastic network model to investigate Grp94's large-scale conformational changes. SPM analysis identified residues that are important in signaling conformational change, many of which have known functional relevance in ATP coordination and catalysis, client binding, and BiP binding. Our findings suggest that ATP hydrolysis in Grp94 alters allosteric wiring and facilitates conformational changes.

Synthesis and evaluation of 3'- and 4'-substituted cyclohexyl noviomimetics that modulate mitochondrial respiration

KU-32 (2) and KU-596 (3), are first and second generation cytoprotective novologues that are derivatives of novobiocin (1), a heat shock protein 90 (Hsp90) C-terminal inhibitor. Although 2 and 3 improve mitochondrial bioenergetics and have demonstrated considerable cytoprotective activity, they contain a synthetically demanding noviose sugar. This issue was initially addressed by creating noviomimetics, such as KU-1202 (4), which replaced the noviose sugar with ether-linked cyclohexyl derivatives that retained some cytoprotective potential due to their ability to increase mitochondrial bioenergetics. Based on structure-activity relationship (SAR) studies of KU-1202 (4), the current study investigated 3'- and 4'-substituted cyclohexyl scaffolds as noviomimetics and determined their efficacy at increasing mitochondrial bioenergetic as a marker for cytoprotective potential.

In Silico Discovery and Optimisation of a Novel Structural Class of Hsp90 C-Terminal Domain Inhibitors

Hsp90 is a promising target for the development of novel agents for cancer treatment. The N-terminal Hsp90 inhibitors have several therapeutic limitations, the most important of which is the induction of heat shock response, which can be circumvented by targeting the allosteric binding site on the C-terminal domain (CTD) of Hsp90. In the absence of an Hsp90—CTD inhibitor co-crystal structure, the use of structure-based design approaches for the Hsp90 CTD is difficult and the structural diversity of Hsp90 CTD inhibitors is limited. In this study, we describe the discovery of a novel structural class of Hsp90 CTD inhibitors. A structure-based virtual screening was performed by docking a library of diverse compounds to the Hsp90β CTD binding site. Three selected virtual hits were tested in the MCF-7 breast cancer cell line, with compound TVS-23 showing antiproliferative activity with an IC₅₀ value of 26.4 ± 1.1 µM. We report here the optimisation, synthesis and biological evaluation of TVS-23 analogues. Several analogues showed significantly enhanced antiproliferative activities in MCF-7 breast cancer and SK-N-MC Ewing sarcoma cell lines, with 7l being the most potent (IC₅₀ = 1.4 ± 0.4 µM MCF-7; IC₅₀ = 2.8 ± 0.4 µM SK-N-MC). The results of this study highlight the use of virtual screening to expand the structural diversity of Hsp90 CTD inhibitors and provide new starting points for further development.

Expanding the Accessible Chemical Space of SIRT2 Inhibitors through Exploration of Binding Pocket Dynamics

Considerations of binding pocket dynamics are one of the crucial aspects of the rational design of binders. Identification of alternative conformational states or cryptic subpockets could lead to the discovery of completely novel groups of the ligands. However, experimental characterization of pocket dynamics, besides being expensive, may not be able to elucidate all of the conformational states relevant for drug discovery projects. In this study, we propose the protocol for computational simulations of sirtuin 2 (SIRT2) binding pocket dynamics and its integration into the structure-based virtual screening (SBVS) pipeline. Initially, unbiased molecular dynamics simulations of SIRT2:inhibitor complexes were performed using optimized force field parameters of SIRT2 inhibitors. Time-lagged independent component analysis (tICA) was used to design pocket-related collective variables (CVs) for enhanced sampling of SIRT2 pocket dynamics. Metadynamics simulations in the tICA eigenvector space revealed alternative conformational states of the SIRT2 binding pocket and the existence of a cryptic subpocket. Newly identified SIRT2 conformational states outperformed experimentally resolved states in retrospective SBVS validation. After performing prospective SBVS, compounds from the under-represented portions of the SIRT2 inhibitor chemical space were selected for in vitro evaluation. Two compounds, NDJ18 and NDJ85, were identified as potent and selective SIRT2 inhibitors, which validated the in silico protocol and opened up the possibility for generalization and broadening of its application. The anticancer effects of the most potent compound NDJ18 were examined on the triple-negative breast cancer cell line. Results indicated that NDJ18 represents a promising structure suitable for further evaluation.

Comparative molecular dynamics simulation analysis of G20 and C92 mutations in c-di-GMP I riboswitch and the wild type with docked c-di-GMP ligand

Riboswitch, a bacterial regulatory RNA consists of an aptamer (specific ligand binding unit) and an expression platform (gene expression modulation unit), which act as a potential drug target as it regulates critical genes. Therefore, it is of interest to glean information on the binding of c-di-GMP ligand to mutated conserved G20 and C92 residues of cyclic diguanosine monophosphate I (c-di-GMP I) riboswitch using molecular dynamics simulation. The result shows that the binding energy of wild/native type riboswitch-ligand complex (3IRW) is lower than the mutant complexes suggesting that the binding affinity for c-di-GMP ligand decreases in case of mutant riboswitches. The hydrogen bonding interactions analysis also showed a high number of hydrogen bonds formation in the wild type riboswitch-ligand complex as compared to the mutant complexes illustrating stronger interaction of ligand to wild type riboswitch than the mutants. The simulation result shows that the mutations affected riboswitch-ligand interactions. The residues G14, G21, C46, A47, and U92 were identified as the key residues which contributed effectively to the binding of c-di-GMP I riboswitch with the natural ligand.

Structure-Activity Relationships of Benzothiazole-Based Hsp90 C-Terminal-Domain Inhibitors

Heat shock protein 90 (Hsp90) is a chaperone responsible for the maturation of many cancer-related proteins, and is therefore an important target for the design of new anticancer agents. Several Hsp90 N-terminal domain inhibitors have been evaluated in clinical trials, but none have been approved as cancer therapies. This is partly due to induction of the heat shock response, which can be avoided using Hsp90 C-terminal-domain (CTD) inhibition. Several structural features have been shown to be useful in the design of Hsp90 CTD inhibitors, including an aromatic ring, a cationic center and the benzothiazole moiety. This study established a previously unknown link between these structural motifs. Using ligand-based design methodologies and structure-based pharmacophore models, a library of 29 benzothiazole-based Hsp90 CTD inhibitors was prepared, and their antiproliferative activities were evaluated in MCF-7 breast cancer cells. Several showed low-micromolar IC₅₀, with the most potent being compounds 5g and 9i (IC₅₀, 2.8 ± 0.1, 3.9 ± 0.1 μM, respectively). Based on these results, a ligand-based structure-activity relationship model was built, and molecular dynamics simulation was performed to elaborate the binding mode of compound 9i. Moreover, compound 9i showed degradation of Hsp90 client proteins and no induction of the heat shock response.

Machine Learning Prediction of Allosteric Drug Activity from Molecular Dynamics

Allosteric drugs have been attracting increasing interest over the past few years. In this context, it is common practice to use high-throughput screening for the discovery of non-natural allosteric drugs. While the discovery stage is supported by a growing amount of biological information and increasing computing power, major challenges still remain in selecting allosteric ligands and predicting their effect on the target protein's function. Indeed, allosteric compounds can act both as inhibitors and activators of biological responses. Computational approaches to the problem have focused on variations on the theme of molecular docking coupled to molecular dynamics with the aim of recovering information on the (long-range) modulation typical of allosteric proteins.

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.

Targeting the HSP90-CDC37-kinase chaperone cycle: A promising therapeutic strategy for cancer

Heat shock protein 90 (HSP90) is an indispensable molecular chaperone that facilitates the maturation of numerous oncoproteins in cancer cells, including protein kinases, ribonucleoproteins, steroid hormone receptors, and transcription factors. Although over 30 HSP90 inhibitors have steadily entered clinical trials, further clinical advancement has been restricted by their limited efficacy, inevitable heat shock response, and multiple side-effects, likely induced via an ATP inhibition mechanism. Since both ATP and various co-chaperones play essential roles in the HSP90 chaperone cycle to achieve integrated function, optimal therapeutics require an understanding of the dynamic interactions among HSP90, ATP, and cochaperones. To date, continuous research has promoted the exploration of the cochaperone cell division cycle 37 (CDC37) as a kinase-specific recognizer and has shown that the HSP90-CDC37-kinase complex is particularly relevant in cancers. Indeed, disrupting the HSP90-CDC37-kinase complex, rather than totally blocking the ATP function of HSP90, is emerging as an alternative way to avoid the limitations of current inhibitors. In this review, we first briefly introduce the HSP90-CDC37-kinase cycle and present the currently available approaches for inhibitor development targeting this cycle and provide insights into selective regulation of the kinase clients of HSP90 by more directional ways.

Tubocapsenolide A targets C-terminal cysteine residues of HSP90 to exert the anti-tumor effect

Heat shock protein 90 (HSP90) is a chaperone protein that has been shown to regulate cancer progression. As a result, HSP90 has emerged as an attractive target for cancer therapy. Tubocapsenolide A (TA) is an anti-tumor component isolated from Tubocapsicum anomalum. Although the anti-tumor activity of TA was considered to be related to HSP90, the binding site and deep anti-tumor mechanisms still need to be elucidated. In this study, we found that TA is a covalent inhibitor of HSP90, which inhibits HSP90 ATPase activity without blocking ATP binding. Further studies indicated that TA targets the C-terminal Cys521 site, which led to HSP90 partial oligomerization and hindered its anti-aggregation and refolding activity. The damage of the chaperone activity disrupted the interaction between HSP90 and its cochaperone CDC37 as well as its client proteins, thereby inducing cell cycle arrest and apoptosis. Moreover, TA was found to have therapeutic effects on the xenograft tumor model by inducing the degradation of HSP90 client proteins. Together, our results identified HSP90 as the direct target of TA for mediating the anti-tumor activity. TA could serve as a lead compound for developing novel HSP90 C-terminal covalent inhibitors with binding site different from the ATP-binding domain.

Dissecting Molecular Principles of the Hsp90 Chaperone Regulation by Allosteric Modulators Using a Hierarchical Simulation Approach and Network Modeling of Allosteric Interactions: Conformational Selection Dictates the Diversity of Protein Responses and Ligand-Specific Functional Mechanisms

Conformational plasticity of the Hsp90 molecular chaperones underlies the diversity of functional mechanisms that these versatile molecular machines employ to coordinate their vast protein clientele in the cellular environment. Despite a steady progress in studies of the Hsp90 machinery, a great deal remains unknown about molecular principles and ligand-specific functional mechanisms of the Hsp90 regulation by allosteric modulators that attracted significant attention because of their therapeutic potential. Due to structural complexity and dynamic nature of the Hsp90 responses to allosteric modulators, the atomistic details about the mode of action of these small molecules continue to be fairly scarce and controversial. In this work, we employ an integrative strategy that encompassed atomistic simulations of the Hsp90 proteins and hierarchical modeling of Hsp90-ligand binding with network analysis to explore functional mechanisms of the Hsp90 regulation by a panel of allosteric modulators (novobiocin, KU-135, KU-174, and KU-32) with different models of action. The results show that functional mechanisms of allosteric modulation in the Hsp90 proteins may be driven by conformational selection principles in which ligands elicit pre-existing states of the unbound chaperone to drive ligand-specific protein responses and distinct scenarios of Hsp90 regulation. We found that novobiocin can selectively sequester an ensemble of open chaperone conformations and inhibit the progression of the functional cycle through a cascade of cumulative dynamic changes. In contrast, KU-32 displayed unique preferences toward partially closed dynamic states, inducing robust allosteric signaling and stimulation of the ATPase cycle. The proposed model of the Hsp90 regulation by allosteric modulators reconciled diverse experimental data and showed that allosteric modulators may operate via targeted exploitation of dynamic landscapes eliciting vastly different protein responses and diverse mechanisms of action.

Comparative structural and mechanistic studies of 4-hydroxy-tetrahydrodipicolinate reductases from Mycobacterium tuberculosis and Vibrio vulnificus

BACKGROUND

The products of the lysine biosynthesis pathway, meso-diaminopimelate and lysine, are essential for bacterial survival. This paper focuses on the structural and mechanistic characterization of 4-hydroxy-tetrahydrodipicolinate reductase (DapB), which is one of the enzymes from the lysine biosynthesis pathway. DapB catalyzes the conversion of (2S, 4S)-4-hydroxy-2,3,4,5-tetrahydrodipicolinate (HTPA) to 2,3,4,5-tetrahydrodipicolinate in an NADH/NADPH dependent reaction. Genes coding for DapBs were identified as essential for many pathogenic bacteria, and therefore DapB is an interesting new target for the development of antibiotics.

METHODS

We have combined experimental and computational approaches to provide novel insights into mechanism of the DapB catalyzed reaction.

RESULTS

Structures of DapBs originating from Mycobacterium tuberculosis and Vibrio vulnificus in complexes with NAD⁺, NADP⁺, as well as with inhibitors, were determined and described. The structures determined by us, as well as currently available structures of DapBs from other bacterial species, were compared and used to elucidate a mechanism of reaction catalyzed by this group of enzymes. Several different computational methods were used to provide a detailed description of a plausible reaction mechanism.

CONCLUSIONS

This is the first report presenting the detailed mechanism of reaction catalyzed by DapB.

GENERAL SIGNIFICANCE

Structural data in combination with information on the reaction mechanism provide a background for development of DapB inhibitors, including transition-state analogues.

Discovery of Novel Hsp90 C-Terminal Inhibitors Using 3D-Pharmacophores Derived from Molecular Dynamics Simulations

Hsp90 C-terminal domain (CTD) inhibitors are promising novel agents for cancer treatment, as they do not induce the heat shock response associated with Hsp90 N-terminal inhibitors. One challenge associated with CTD inhibitors is the lack of a co-crystallized complex, requiring the use of predicted allosteric apo pocket, limiting structure-based (SB) design approaches. To address this, a unique approach that enables the derivation and analysis of interactions between ligands and proteins from molecular dynamics (MD) trajectories was used to derive pharmacophore models for virtual screening (VS) and identify suitable binding sites for SB design. Furthermore, ligand-based (LB) pharmacophores were developed using a set of CTD inhibitors to compare VS performance with the MD derived models. Virtual hits identified by VS with both SB and LB models were tested for antiproliferative activity. Compounds 9 and 11 displayed antiproliferative activities in MCF-7 and Hep G2 cancer cell lines. Compound 11 inhibited Hsp90-dependent refolding of denatured luciferase and induced the degradation of Hsp90 clients without the concomitant induction of Hsp70 levels. Furthermore, compound 11 offers a unique scaffold that is promising for the further synthetic optimization and development of molecules needed for the evaluation of the Hsp90 CTD as a target for the development of anticancer drugs.

Comparative structure, dynamics and evolution of acyl-carrier proteins from Borrelia burgdorferi, Brucella melitensis and Rickettsia prowazekii

Acyl carrier proteins (ACPs) are small helical proteins found in all kingdoms of life, primarily involved in fatty acid and polyketide biosynthesis. In eukaryotes, ACPs are part of the fatty acid synthase (FAS) complex, where they act as flexible tethers for the growing lipid chain, enabling access to the distinct active sites in FAS. In the type II synthesis systems found in bacteria and plastids, these proteins exist as monomers and perform various processes, from being a donor for synthesis of various products such as endotoxins, to supplying acyl chains for lipid A and lipoic acid FAS (quorum sensing), but also as signaling molecules, in bioluminescence and activation of toxins. The essential and diverse nature of their functions makes ACP an attractive target for antimicrobial drug discovery. Here, we report the structure, dynamics and evolution of ACPs from three human pathogens: Borrelia burgdorferi, Brucella melitensis and Rickettsia prowazekii, which could facilitate the discovery of new inhibitors of ACP function in pathogenic bacteria.

Gaining deeper insights into the surface binding of bedaquiline analogues with the ATP synthase subunit C of Mycobacterium tuberculosis using molecular docking, molecular dynamics simulation and 3D-QSAR techniques

This computational study exclusively illustrates the key molecular features of bedaquiline and its analogues required for binding to mycobacterial ATP synthase.

Ligand Binding, Unbinding, and Allosteric Effects: Deciphering Small-Molecule Modulation of HSP90

The molecular chaperone HSP90 oversees the functional activation of a large number of client proteins. Because of its role in multiple pathways linked to cancer and neurodegeneration, drug discovery targeting HSP90 has been actively pursued. Yet, a number of inhibitors failed to meet expectations due to induced toxicity problems. In this context, allosteric perturbation has emerged as an alternative strategy for the pharmacological modulation of HSP90 functions. Specifically, novel allosteric stimulators showed the interesting capability of accelerating HSP90 closure dynamics and ATPase activities while inducing tumor cell death. Here, we gain atomistic insight into the mechanisms of allosteric ligand recognition and their consequences on the functional dynamics of HSP90, starting from the fully unbound state. We integrate advanced computational sampling methods based on FunnelMetadynamics, with the analysis of internal dynamics of the structural ensembles visited during the simulations. We observe several binding/unbinding events, and from these, we derive an accurate estimation of the absolute binding free energy. Importantly, we show that different binding poses induce different dynamics states. Our work for the first time explicitly correlates HSP90 responses to binding/unbinding of an allosteric ligand to the modulation of functionally oriented protein motions.

Modulation of Human Hsp90α Conformational Dynamics by Allosteric Ligand Interaction at the C-Terminal Domain

Recent years have seen heat shock protein 90 kDa (Hsp90) attract significant interest as a viable drug target, particularly for cancer. To date, designed inhibitors that target the ATPase domain demonstrate potent anti-proliferative effects, but have failed clinical trials due to high levels of associated toxicity. To circumvent this, the focus has shifted away from the ATPase domain. One option involves modulation of the protein through allosteric activation/inhibition. Here, we propose a novel approach: we use previously obtained information via residue perturbation scanning coupled with dynamic residue network analysis to identify allosteric drug targeting sites for inhibitor docking. We probe the open conformation of human Hsp90α for druggable sites that overlap with these allosteric control elements, and identify three putative natural compound allosteric modulators: Cephalostatin 17, 20(29)-Lupene-3β-isoferulate and 3'-Bromorubrolide F. We assess the allosteric potential of these ligands by examining their effect on the conformational dynamics of the protein. We find evidence for the selective allosteric activation and inhibition of Hsp90's conformational transition toward the closed state in response to ligand binding and shed valuable insight to further the understanding of allosteric drug design and Hsp90's complex allosteric mechanism of action.

Synthesis and evaluation of a ring-constrained Hsp90 C-terminal inhibitor that exhibits neuroprotective activity

KU-596 is a second-generation C-terminal heat shock protein 90 KDa (Hsp90) modulator based on the natural product, novobiocin. KU-596 has been shown to induce Hsp70 levels and manifest neuroprotective activity through induction of the heat shock response. A ring-constrained analog of KU-596 was designed and synthesized to probe its binding orientation and ability to induce Hsp70 levels. Compound 2 was found to exhibit comparable or increased activity compared to KU-596, which is under clinical investigation for the treatment of neuropathy.

In silico identification and analysis of the binding site for aminocoumarin type inhibitors in the C-terminal domain of Hsp90

Hsp90 contains two Nucleotide Binding Sites (NBS): one each in its N-terminal domain (NTD) and C-terminal domain (CTD), respectively. Previously we used computational techniques to locate a nucleotide-binding site in the CTD. Nucleotide binding at this site stabilized the structurally labile region within this domain, thus providing a rationale for increased resistance to thermal denaturation and proteolysis. A scan for ligand-binding sites in CTD revealed four potential sites with the requisite volume to accommodate aminocoumarins and -derived inhibitors. Only one of these reproducibly formed docked complexes with inhibitors and showed excellent interactions with residues lining the site. Fortuitously, it was identical to the aforementioned nucleotide-binding site thus providing an explanation for the reported direct competition between inhibitors and nucleotides. Further studies with carefully chosen inhibitors and some inactive analogues provided an explanation for the known Structure-Activity Relationships (SAR) of aminocoumarin and -derived inhibitors. We also performed similar studies of the NTD to discern the reason(s) for its inability to bind aminocoumarins, given the family resemblance to prokaryotic Top-IV and Gyr-B. Our studies permitted the identification of the putative inhibitor binding site in the CTD, an explanation for increased resistance to thermal denaturation and proteolysis upon inhibitor binding as well as direct competition with ATP.

Discovery of new molecular entities able to strongly interfere with Hsp90 C-terminal domain

Heat shock protein 90 (Hsp90) is an ATP dependent molecular chaperone deeply involved in the complex network of cellular signaling governing some key functions, such as cell proliferation and survival, invasion and angiogenesis. Over the past years the N-terminal protein domain has been fully investigated as attractive strategy against cancer, but despite the many efforts lavished in the field, none of the N-terminal binders (termed "classical inhibitors"), currently in clinical trials, have yet successfully reached the market, because of the detrimental heat shock response (HSR) that showed to induce; thus, recently, the selective inhibition of Hsp90 C-terminal domain has powerfully emerged as a more promising alternative strategy for anti-cancer therapy, not eliciting this cell rescue cascade. However, the structural complexity of the target protein and, mostly, the lack of a co-crystal structure of C-terminal domain-ligand, essential to drive the identification of new hits, represent the largest hurdles in the development of new selective C-terminal inhibitors. Continuing our investigations on the identification of new anticancer drug candidates, by using an orthogonal screening approach, here we describe two new potent C-terminal inhibitors able to induce cancer cell death and a considerable down-regulation of Hsp90 client oncoproteins, without triggering the undesired heat shock response.

Therapeutic Strategies for Restoring Tau Homeostasis

Normal tau homeostasis is achieved when the synthesis, processing, and degradation of the protein is balanced. Together, the pathways that regulate tau homeostasis ensure that the protein is at the proper levels and that its posttranslational modifications and subcellular localization are appropriately controlled. These pathways include the enzymes responsible for posttranslational modifications, those systems that regulate mRNA splicing, and the molecular chaperones that control tau turnover and its binding to microtubules. In tauopathies, this delicate balance is disturbed. Tau becomes abnormally modified by posttranslational modification, it loses affinity for microtubules, and it accumulates in proteotoxic aggregates. How and why does this imbalance occur? In this review, we discuss how molecular chaperones and other components of the protein homeostasis (e.g., proteostasis) network normally govern tau quality control. We also discuss how aging might reduce the capacity of these systems and how tau mutations might further affect this balance. Finally, we discuss how small-molecule inhibitors are being used to probe and perturb the tau quality-control systems, playing a particularly prominent role in revealing the logic of tau homeostasis. As such, there is now interest in developing these chemical probes into therapeutics, with the goal of restoring normal tau homeostasis to treat disease.

Computational Modeling of the Hsp90 Interactions with Cochaperones and Small-Molecule Inhibitors

Allosteric interactions of the molecular chaperone Hsp90 with a diverse array of cochaperones and client proteins, such as protein kinases and transcription factors, allow for efficient molecular communication in signal transduction networks. Deregulation of pathways involving these proteins is commonly associated with cancer pathologies and allosteric inhibition of oncogenic clients by targeting Hsp90 provides a powerful therapeutic strategy in cancer research. We review several validated computational approaches and tools used in the studies of the Hsp90 interactions with proteins and small molecules. These methods include experimentally guided docking to predict Hs90-protein interactions, molecular and binding free energy simulations to analyze Hsp90 binding with small molecules, and structure-based network modeling to evaluate allosteric interactions and communications in the Hsp90 regulatory complexes. Through the lens of allosteric-centric view on Hsp90 function and regulation, we discuss newly emerging computational tools that link protein structure modeling with biophysical simulations and network-based systems biology approaches.

Synthesis and Biological Evaluation of Stilbene Analogues as Hsp90 C-Terminal Inhibitors

The design, synthesis, and biological evaluation of stilbene-based novobiocin analogues is reported. Replacement of the biaryl amide side chain with a triazole side chain produced compounds that exhibited good antiproliferative activities. Heat shock protein 90 (Hsp90) inhibition was observed when N-methylpiperidine was replaced with acyclic tertiary amines on the stilbene analogues that also contain a triazole-derived side chain. These studies revealed that ≈24 Å is the optimal length for compounds that exhibit good antiproliferative activity as a result of Hsp90 inhibition.

Design, Synthesis, and Biological Activities of Vibsanin B Derivatives: A New Class of HSP90 C-Terminal Inhibitors

Previously, vibsanin B (ViB) was found to preferentially target HSP90β compared to HSP90α. In this study, multiple experiments, including pull-down assays of biotin-ViB with recombinant HSP90β-NTD, MD, CTD, and full-length HSP90β, molecular docking of ViB and its derivatives to the HSP90 CTD, and a inhibition assay of interaction of the HSP90β CTD with GST-tagged cyclophilin 40 (Cyp40) by ViB derivatives, suggest that ViB can directly bind to the HSP90 C-terminus. On the basis of the docking predictions and primary structure-activity relationships (SARs), a series of ViB analogues devised with focus on the C18 position, along with compounds derivatized at the C4, C7, and C8 positions, were designed and chemically synthesized. Compound 12f (IC₅₀ = 1.12 μM against SK-BR-3) exhibits great potency with drug-like properties. Overall, our findings demonstrate that compounds with the vibsanin B scaffold are a new class of HSP90 C-terminal inhibitors with considerable potential as anticancer agents.

Targeting the Hsp90 C-terminal domain to induce allosteric inhibition and selective client downregulation

BACKGROUND

Inhibition of Hsp90 is desirable due to potential downregulation of oncogenic clients. Early generation inhibitors bind to the N-terminal domain (NTD) but C-terminal domain (CTD) inhibitors are a promising class because they do not induce a heat shock response. Here we present a new structural class of CTD binding molecules with a unique allosteric inhibition mechanism.

METHODS

A hit molecule, NSC145366, and structurally similar probes were assessed for inhibition of Hsp90 activities. A ligand-binding model was proposed indicating a novel Hsp90 CTD binding site. Client protein downregulation was also determined.

RESULTS

NSC145366 interacts with the Hsp90 CTD and has anti-proliferative activity in tumor cell lines (GI₅₀=0.2-1.9μM). NSC145366 increases Hsp90 oligomerization resulting in allosteric inhibition of NTD ATPase activity (IC₅₀=119μM) but does not compete with NTD or CTD-ATP binding. Treatment of LNCaP prostate tumor cells resulted in selective client protein downregulation including AR and BRCA1 but without a heat shock response. Analogs had similar potencies in ATPase and chaperone activity assays and variable effects on oligomerization. In silico modeling predicted a binding site at the CTD dimer interface distinct from the nucleotide-binding site.

CONCLUSIONS

A set of symmetrical scaffold molecules with bisphenol A cores induced allosteric inhibition of Hsp90. Experimental evidence and molecular modeling suggest that the binding site is independent of the CTD-ATP site and consistent with unique induction of allosteric effects.

GENERAL SIGNIFICANCE

Allosteric inhibition of Hsp90 via a mechanism used by the NSC145366-based probes is a promising avenue for selective oncogenic client downregulation.

Design of Allosteric Stimulators of the Hsp90 ATPase as New Anticancer Leads

Allosteric compounds that stimulate Hsp90 adenosine triphosphatase (ATPase) activity were rationally designed, showing anticancer potencies in the low micromolar to nanomolar range. In parallel, the mode of action of these compounds was clarified and a quantitative model that links the dynamic ligand-protein cross-talk to observed cellular and in vitro activities was developed. The results support the potential of using dynamics-based approaches to develop original mechanism-based cancer therapeutics.

In silico identification and computational analysis of the nucleotide binding site in the C-terminal domain of Hsp90

Hsp90 contains two distinct Nucleotide Binding Sites (NBS), in its N-terminal domain (NTD) and C-terminal domain (CTD), respectively. The NTD site belongs to the GHKL super-family of ATPases and has been the subject of extensive characterization. However, a structure of the nucleotide-bound form of CTD is still unavailable. In this study molecular modeling was employed to incorporate experimental data using partial constructs of the CTD, from work published by many research groups, onto existing structural models of its apo- form. Our attempts to locate potential nucleotide ligand-binding sites or cavities yielded one major candidate-a structurally unconventional site-exhibiting the requisite shape and volume for accommodation of tri-phosphate nucleotides. Its structure was refined by molecular dynamics (MD)-based techniques. We reproducibly docked the Mg²⁺ complexed form of ATP, GTP, CTP, TTP and UTP to this putative NBS. These docking simulations and calculated ligand-binding scores are in general agreement with published data about experimentally measured binding to the CTD. The overall pattern of interactions between residues lining the site and docked nucleotides is conserved and broadly similar to that of other nucleotide-binding sites. Our docking simulations suggest that nucleotide binding stabilizes the only structurally labile region, thereby providing a rationale for the increased resistance to thermal denaturation and proteolysis. The docked nucleotides do not intrude onto the surface of residues involved in dimerization or chaperoning. Our molecular modeling permitted recognition of larger structural changes in the nucleotide-bound CTD dimer, including stabilization of helix-2 in both chains and intra- and inter- chain interactions between three residues (I613, Q617, R620).

Synthesis and Biological Evaluation of Novobiocin Core Analogues as Hsp90 Inhibitors

Development of heat shock protein 90 (Hsp90) C-terminal inhibitors has emerged as an exciting strategy for the treatment of cancer. Previous efforts have focused on modifications to the natural products novobiocin and coumermycin. Moreover, variations in both the sugar and amide moieties have been extensively studied, whereas replacements for the coumarin core have received less attention. Herein, 24 cores were synthesized with varying distances and angles between the sugar and amide moieties. Compounds that exhibited good anti-proliferative activity against multiple cancer cell lines and Hsp90 inhibitory activity, were those that placed the sugar and amide moieties between 7.7 and 12.1 Å apart along with angles of 180°.

New dihydropyrimidin-2(1H)-one based Hsp90 C-terminal inhibitors

The inhibition of the C-terminal domain of heat shock protein 90 (Hsp90) is emerging as a novel strategy for cancer therapy, therefore the identification of a new class of C-terminal inhibitors is strongly required.

Targeting the Hsp90 C-terminal domain by the chemically accessible dihydropyrimidinone scaffold

Hsp90 C-terminal ligands are potential new anti-cancer drugs alternative to the more studied N-terminal inhibitors. Here we report the identification of a new dihydropyrimidinone binding the C-terminus, which is not structurally related to other well-known natural and nature-inspired inhibitors of this second druggable Hsp90 site.

Inhibition of succinate dehydrogenase by the mitochondrial chaperone TRAP1 has anti-oxidant and anti-apoptotic effects on tumor cells

TRAP1 is a mitochondrial chaperone highly expressed in many tumor types; it inhibits respiratory complex II, down-modulating its succinate dehydrogenase (SDH) enzymatic activity. SDH inhibition in turn leads to a pseudohypoxic state caused by succinate-dependent HIF1α stabilization and promotes neoplastic growth. Here we report that TRAP1 inhibition of SDH also shields cells from oxidative insults and from the ensuing lethal opening of the mitochondrial permeability transition pore. This anti-oxidant activity of TRAP1 protects tumor cells from death in conditions of nutrient paucity that mimic those encountered in the neoplasm during the process of malignant accrual, and it is required for in vitro tumorigenic growth. Our findings demonstrate that SDH inhibition by TRAP1 is oncogenic not only by inducing pseudohypoxia, but also by protecting tumor cells from oxidative stress.

Structural Stability and Flexibility Direct the Selection of Activating Mutations in Epidermal Growth Factor Receptor Kinase

Herein we investigate the potential of novel methods of molecular dynamics analysis to provide information on the key factors that underlie the preferential localization and the effects of mutations modulating protein activities. Epidermal growth factor receptor (EGFR) kinases are selected as a test case. The combined analysis of protein energetics and internal dynamics indicates a clear polarization in the native protein, whereby a highly stable and ordered scaffold in one domain, namely the C-lobe, is combined to a flexible and loosely stabilized domain, the N-lobe. The subdivision in two portions with different properties directs the presence of point mutations mainly to the N-lobe. This allows modulating protein flexibility so that the protein can more efficiently sample the conformations necessary for substrate recognition, while leaving the stability of the protein unperturbed. In this context, comparative simulations of EGFR in the wild type sequence and in the presence of the activating oncogenic mutation G719S reveal flexibility changes in several key regions, involving in particular the part of the kinase devoted to the regulation of substrate recognition (regulatory core) and an increase in the number of stabilizing interactions in the N-lobe for the activated mutant. Our approaches represent a promising and simple strategy toward rationalizing the effects of mutations in modulating enzymatic activities.

Design, synthesis and biological evaluation of biphenylamide derivatives as Hsp90 C-terminal inhibitors

Modulation of Hsp90 C-terminal function represents a promising therapeutic approach for the treatment of cancer and neurodegenerative diseases. Current drug discovery efforts toward Hsp90 C-terminal inhibition focus on novobiocin, an antibiotic that was transformed into an Hsp90 inhibitor. Based on structural information obtained during the development of novobiocin derivatives and molecular docking studies, scaffolds containing a biphenyl moiety in lieu of the coumarin ring present in novobiocin were identified as new Hsp90 C-terminal inhibitors. Structure-activity relationship studies produced new derivatives that inhibit the proliferation of breast cancer cell lines at nanomolar concentrations, which corresponded directly with Hsp90 inhibition.

Investigating the dynamic aspects of drug-protein recognition through a combination of MD and NMR analyses: implications for the development of protein-protein interaction inhibitors

In this paper, we investigate the dynamic aspects of the molecular recognition between a small molecule ligand and a flat, exposed protein surface, representing a typical target in the development of protein-protein interaction inhibitors. Specifically, we analyze the complex between the protein Fibroblast Growth Factor 2 (FGF2) and a recently discovered small molecule inhibitor, labeled sm27 for which the binding site and the residues mainly involved in small molecule recognition have been previously characterized. We have approached this problem using microsecond MD simulations and NMR-based characterizations of the dynamics of the apo and holo states of the system. Using direct combination and cross-validation of the results of the two techniques, we select the set of conformational states that best recapitulate the principal dynamic and structural properties of the complex. We then use this information to generate a multi-structure representation of the sm27-FGF2 interaction. We propose this kind of representation and approach as a useful tool in particular for the characterization of systems where the mutual dynamic influence between the interacting partners is expected to play an important role. The results presented can also be used to generate new rules for the rational expansion of the chemical diversity space of FGF2 inhibitors.