Artificial Accelerators of the Molecular Chaperone Hsp90 Facilitate Rate‐Limiting Conformational Transitions

The molecular chaperone Hsp90 undergoes an ATP‐driven cycle of conformational changes in which large structural rearrangements precede ATP hydrolysis. Well‐established small‐molecule inhibitors of Hsp90 compete with ATP‐binding. We wondered whether compounds exist that can accelerate the conformational cycle. In a FRET‐based screen reporting on conformational rearrangements in Hsp90 we identified compounds. We elucidated their mode of action and showed that they can overcome the intrinsic inhibition in Hsp90 which prevents these rearrangements. The mode of action is similar to that of the co‐chaperone Aha1 which accelerates the Hsp90 ATPase. However, while the two identified compounds influence conformational changes, they target different aspects of the structural transitions. Also, the binding site determined by NMR spectroscopy is distinct. This study demonstrates that small molecules are capable of triggering specific rate‐limiting transitions in Hsp90 by mechanisms similar to those in protein cofactors.     

Hierarchical dynamics in allostery following ATP hydrolysis monitored by single molecule FRET measurements and MD simulations

We report on a study that combines advanced fluorescence methods with molecular dynamics (MD) simulations to cover timescales from nanoseconds to milliseconds for a large protein. This allows us to delineate how ATP hydrolysis in a protein causes allosteric changes at a distant protein binding site, using the chaperone Hsp90 as test system. The allosteric process occurs via hierarchical dynamics involving timescales from nano- to milliseconds and length scales from Ångstroms to several nanometers. We find that hydrolysis of one ATP is coupled to a conformational change of Arg380, which in turn passes structural information via the large M-domain α-helix to the whole protein. The resulting structural asymmetry in Hsp90 leads to the collapse of a central folding substrate binding site, causing the formation of a novel collapsed state (closed state B) that we characterise structurally. We presume that similar hierarchical mechanisms are fundamental for information transfer induced by ATP hydrolysis through many other proteins.

We report on a study that combines advanced fluorescence methods with molecular dynamics simulations to cover timescales from nanoseconds to milliseconds for a large protein, the chaperone Hsp90.     

Identification of Limonol Derivatives as Heat Shock Protein 90 (Hsp90) Inhibitors through a Multidisciplinary Approach

The identification of inhibitors of Hsp90 is currently a primary goal in the development of more effective drugs for the treatment of various types of multidrug resistant malignancies. In an attempt to identify new small molecules modulating the activity of Hsp90, we screened a small library of tetranortriterpenes. A high‐affinity interaction with Hsp90 inducible form was uncovered for eight of these compounds, five of which are described here for the first time. By monitoring the ATPase activity and the citrate synthase thermal induced aggregation, compound 1 (cedrelosin A), 3 (7α‐limonylacetate), and 5 (cedrelosin B), containing a limonol moiety, were found to be the most effective in compromising the Hsp90α chaperone activity. Consistent with these findings, the three compounds caused a depletion of c‐Raf and pAkt Hsp90 client proteins in HeLa and MCF/7 cell lines. Induced fit docking protocol and molecular dynamics were used to rationalize the structural basis of the biological activity of the limonol derivatives. Taken together, these results point to limonol‐derivatives as promising scaffolds for the design of novel Hsp90α inhibitors.     

Selective Inhibitors of FKBP51 Employ Conformational Selection of Dynamic Invisible States

The recently discovered SAFit class of inhibitors against the Hsp90 co‐chaperone FKBP51 show greater than 10 000‐fold selectivity over its closely related paralogue FKBP52. However, the mechanism underlying this selectivity remained unknown. By combining NMR spectroscopy, biophysical and computational methods with mutational analysis, we show that the SAFit molecules bind to a transient pocket in FKBP51. This represents a weakly populated conformation resembling the inhibitor‐bound state of FKBP51, suggesting conformational selection rather than induced fit as the major binding mechanism. The inhibitor‐bound conformation of FKBP51 is stabilized by an allosteric network of residues located away from the inhibitor‐binding site. These residues stabilize the Phe67 side chain in a dynamic outward conformation and are distinct in FKBP52, thus rationalizing the basis for the selectivity of SAFit inhibitors. Our results represent a paradigm for the selective inhibition of transient binding pockets.     

Reinventing Hsp90 Inhibitors: Blocking C‐Terminal Binding Events to Hsp90 by Using Dimerized Inhibitors

Heat shock protein 90 (Hsp90) is a molecular chaperone (90 kDa) that functions as a dimer. This protein facilitates the folding, assembly, and stabilization of more than 400 proteins that are responsible for cancer development and progression. Inhibiting Hsp90’s function will shut down multiple cancer‐driven pathways simultaneously because oncogenic clients rely heavily on Hsp90, which makes this chaperone a promising anticancer target. Classical inhibitors that block the binding of adenine triphosphate (ATP) to the N‐terminus of Hsp90 are highly toxic to cells and trigger a resistance mechanism within cells. This resistance mechanism comprises a large increase in prosurvival proteins, namely, heat shock protein 70 (Hsp70), heat shock protein 27 (Hsp27), and heat shock factor 1 (HSF‐1). Molecules that modulate the C‐terminus of Hsp90 are effective at inducing cancer‐cell death without activating the resistance mechanism. Herein, we describe the design, synthesis, and biological binding affinity for a series of dimerized C‐terminal Hsp90 modulators. We show that dimers of these C‐terminal modulators synergistically inhibit Hsp90 relative to monomers.     

Large scale free energy calculations for blind predictions of protein–ligand binding: the D3R Grand Challenge 2015

We describe binding free energy calculations in the D3R Grand Challenge 2015 for blind prediction of the binding affinities of 180 ligands to Hsp90. The present D3R challenge was built around experimental datasets involving Heat shock protein (Hsp) 90, an ATP-dependent molecular chaperone which is an important anticancer drug target. The Hsp90 ATP binding site is known to be a challenging target for accurate calculations of ligand binding affinities because of the ligand-dependent conformational changes in the binding site, the presence of ordered waters and the broad chemical diversity of ligands that can bind at this site. Our primary focus here is to distinguish binders from nonbinders. Large scale absolute binding free energy calculations that cover over 3000 protein–ligand complexes were performed using the BEDAM method starting from docked structures generated by Glide docking. Although the ligand dataset in this study resembles an intermediate to late stage lead optimization project while the BEDAM method is mainly developed for early stage virtual screening of hit molecules, the BEDAM binding free energy scoring has resulted in a moderate enrichment of ligand screening against this challenging drug target. Results show that, using a statistical mechanics based free energy method like BEDAM starting from docked poses offers better enrichment than classical docking scoring functions and rescoring methods like Prime MM-GBSA for the Hsp90 data set in this blind challenge. Importantly, among the three methods tested here, only the mean value of the BEDAM binding free energy scores is able to separate the large group of binders from the small group of nonbinders with a gap of 2.4 kcal/mol. None of the three methods that we have tested provided accurate ranking of the affinities of the 147 active compounds. We discuss the possible sources of errors in the binding free energy calculations. The study suggests that BEDAM can be used strategically to discriminate binders from nonbinders in virtual screening and to more accurately predict the ligand binding modes prior to the more computationally expensive FEP calculations of binding affinity.     

Design, synthesis, and structure--activity relationships for chimeric inhibitors of Hsp90

Inhibition of the 90 kDa heat shock protein (Hsp90) family of molecular chaperones represents a promising new chemotherapeutic approach toward the treatment of several cancers. Previous studies have demonstrated that the natural products, radicicol and geldanamycin, are potent inhibitors of the Hsp90 N-terminal ATP binding site. The cocrystal structures of these molecules bound to Hsp90 have been determined, and through molecular modeling and superimposition of these ligands, hybrids of radicicol and geldanamycin have been designed. A series of macrocylic chimeras of radicicol and geldanamycin and the corresponding seco-agents have been prepared and evaluated for both antiproliferative activity and their ability to induce Hsp90-dependent client protein degradation.     

Design and Synthesis of Hsp90 Inhibitors with B-Raf and PDHK1 Multi-Target Activity

The design of multi-target ligands has become an innovative approach for the identification of effective therapeutic treatments against complex diseases, such as cancer. Recent studies have demonstrated that the combined inhibition of Hsp90 and B-Raf provides synergistic effects against several types of cancers. Moreover, it has been reported that PDHK1, which presents an ATP-binding pocket similar to that of Hsp90, plays an important role in tumor initiation, maintenance and progression, participating also to the senescence process induced by B-Raf oncogenic proteins. Based on these premises, the simultaneous inhibition of these targets may provide several benefits for the treatment of cancer. In this work, we set up a design strategy including the assembly and integration of molecular fragments known to be important for binding to the Hsp90, PDHK1 and B-Raf targets, aided by molecular docking for the selection of a set of compounds potentially able to exert Hsp90-B-Raf-PDHK1 multi-target activities. The designed compounds were synthesized and experimentally validated in vitro. According to the in vitro assays, compounds 4 a , 4 d and 4 e potently inhibited Hsp90 and moderately inhibited the PDHK1 kinase. Finally, molecular dynamics simulations were performed to provide further insights into the structural basis of their multi-target activity.     

Parameterization of Org27569: An allosteric modulator of the cannabinoid CB1 G protein‐coupled receptor

The cannabinoid CB₁ receptor is a class A G protein‐coupled receptor (GPCR) that is the most widely expressed GPCR in the brain. Many GPCRs contain allosteric binding sites for endogenous and/or synthetic ligands, which are topographically distinct from the agonist‐binding site that is known as the orthosteric site. While both endogenous and synthetic ligands that act at the CB₁ orthosteric site have been known for some time, compounds that act at a CB₁ allosteric site have only recently been discovered. The most studied of these is 5‐chloro‐3‐ethyl‐1H‐indole‐2‐carboxylic acid [2‐(4‐piperidin‐1‐ylphenyl)ethyl]amide (Org27569). Because allosteric ligands are thought to act through conformational changes in the receptor that are transmitted from the allosteric to the orthosteric site, computational studies of the structural and dynamic interactions of Org27569 with the CB₁ receptor are crucial to achieve a molecular level understanding of the basis of action of this important new class of compounds. To date, such computational studies have not been possible due to the lack of a complete set of molecular mechanics force field parameters for Org27569. Here, we present the development of missing CHARMM force field parameters for Org27569 using previously published methods and the validation and application of these new parameters using normal mode analysis and molecular dynamics simulations combined with experimental infrared measurements. © 2011 Wiley Periodicals, Inc. J Comput Chem, 2011     

Exploring Weak Ligand–Protein Interactions by Long‐Lived NMR States: Improved Contrast in Fragment‐Based Drug Screening

Ligands that have an affinity for protein targets can be screened very effectively by exploiting favorable properties of long‐lived states (LLS) in NMR spectroscopy. In this work, we describe the use of LLS for competitive binding experiments to measure accurate dissociation constants of fragments that bind weakly to the ATP binding site of the N‐terminal ATPase domain of heat shock protein 90 (Hsp90), a therapeutic target for cancer treatment. The LLS approach allows one to characterize ligands with an exceptionally wide range of affinities, since it can be used for ligand concentrations [L] that are several orders of magnitude smaller than the dissociation constants K_D. This property makes the LLS method particularly attractive for the initial steps of fragment‐based drug screening, where small molecular fragments that bind weakly to a target protein must be identified, which is a difficult task for many other biophysical methods.     

Fuel‐Responsive Allosteric DNA‐Based Aptamers for the Transient Release of ATP and Cocaine

We show herein that allostery offers a key strategy for the design of out‐of‐equilibrium systems by engineering allosteric DNA‐based nanodevices for the transient loading and release of small organic molecules. To demonstrate the generality of our approach, we used two model DNA‐based aptamers that bind ATP and cocaine through a target‐induced conformational change. We re‐engineered these aptamers so that their affinity towards their specific target is controlled by a DNA sequence acting as an allosteric inhibitor. The use of an enzyme that specifically cleaves the inhibitor only when it is bound to the aptamer generates a transient allosteric control that leads to the release of ATP or cocaine from the aptamers. Our approach confirms that the programmability and predictability of nucleic acids make synthetic DNA/RNA the perfect candidate material to re‐engineer synthetic receptors that can undergo chemical fuel‐triggered release of small‐molecule cargoes and to rationally design non‐equilibrium systems.     

A Bioresistant Nitroxide Spin Label for In‐Cell EPR Spectroscopy: In Vitro and In Oocytes Protein Structural Dynamics Studies

Approaching protein structural dynamics and protein–protein interactions in the cellular environment is a fundamental challenge. Owing to its absolute sensitivity and to its selectivity to paramagnetic species, site‐directed spin labeling (SDSL) combined with electron paramagnetic resonance (EPR) has the potential to evolve into an efficient method to follow conformational changes in proteins directly inside cells. Until now, the use of nitroxide‐based spin labels for in‐cell studies has represented a major hurdle because of their short persistence in the cellular context. The design and synthesis of the first maleimido‐proxyl‐based spin label (M‐TETPO) resistant towards reduction and being efficient to probe protein dynamics by continuous wave and pulsed EPR is presented. In particular, the extended lifetime of M‐TETPO enabled the study of structural features of a chaperone in the absence and presence of its binding partner at endogenous concentration directly inside cells.     

Rational design of allosteric ribozymes

Background: Efficient operation of cellular processes relies on the strict control that each cell exerts over its metabolic pathways. Some protein enzymes are subject to allosteric regulation, in which binding sites located apart from the enzyme's active site can specifically recognize effector molecules and alter the catalytic rate of the enzyme via conformational changes. Although RNA also performs chemical reactions, no ribozymes are known to operate as true allosteric enzymes in biological systems. It has recently been established that small-molecule receptors can readily be made of RNA, as demonstrated by the in vitro selection of various RNA aptamers that can specifically bind corresponding ligand molecules. We set out to examine whether the catalytic activity of an existing ribozyme could be brought under the control of an effector molecule by designing conjoined aptamer-ribozyme complexes.Results: By joining an ATP-binding RNA to a self-cleaving ribozyme, we have created the first example of an allosteric ribozyme that has a catalytic rate that can be controlled by ATP. A 180-fold reduction in rate is observed upon addition of either adenosine or ATP, but no inhibition is detected in the presence of dATP or other nucleoside triphosphates. Mutations in the aptamer domain that are expected to eliminate ATP binding or that increase the distance between aptamer and ribozyme domains result in a loss of ATP-specific allosteric control. Using a similar design approach, allosteric hammerhead ribozymes that are activated in the presence of ATP were created and another ribozyme that can be controlled by theophylline was created.Conclusions: The catalytic features of these conjoined aptamer-ribozyme constructs demonstrate that catalytic RNAs can also be subject to allosteric regulation — a key feature of certain protein enzymes. Moreover, by using simple rational design strategies, it is now possible to engineer new catalytic polynucleotides which have rates that can be tightly and specifically controlled by small effector molecules.     

Insights into Conformational Dynamics and Allostery in DNMT1-H3Ub/USP7 Interactions

DNA methyltransferases (DNMTs) including DNMT1 are a conserved family of cytosine methylases that play crucial roles in epigenetic regulation. The versatile functions of DNMT1 rely on allosteric networks between its different interacting partners, emerging as novel therapeutic targets. In this work, based on the modeling structures of DNMT1-ubiquitylated H3 (H3Ub)/ubiquitin specific peptidase 7 (USP7) complexes, we have used a combination of elastic network models, molecular dynamics simulations, structural residue perturbation, network modeling, and pocket pathway analysis to examine their molecular mechanisms of allosteric regulation. The comparative intrinsic and conformational dynamics analysis of three DNMT1 systems has highlighted the pivotal role of the RFTS domain as the dynamics hub in both intra- and inter-molecular interactions. The site perturbation and network modeling approaches have revealed the different and more complex allosteric interaction landscape in both DNMT1 complexes, involving the events caused by mutational hotspots and post-translation modification sites through protein-protein interactions (PPIs). Furthermore, communication pathway analysis and pocket detection have provided new mechanistic insights into molecular mechanisms underlying quaternary structures of DNMT1 complexes, suggesting potential targeting pockets for PPI-based allosteric drug design.     

Exhaustive search of the configurational space of heat-shock protein 90 with its inhibitor by multicanonical molecular dynamics based dynamic docking

Multicanonical molecular dynamics based dynamic docking was used to exhaustively search the configurational space of an inhibitor binding to the N-terminal domain of heat-shock protein 90 (Hsp90). The obtained structures at 300 K cover a wide structural ensemble, with the top two clusters ranked by their free energy coinciding with the native binding site. The representative structure of the most stable cluster reproduced the experimental binding configuration, but an interesting conformational change in Hsp90 could be observed. The combined effects of solvation and ligand binding shift the equilibrium from a preferred loop-in conformation in the unbound state to an α-helical one in the bound state for the flexible lid region of Hsp90. Thus, our dynamic docking method is effective at predicting the native binding site while exhaustively sampling a wide configurational space, modulating the protein structure upon binding.     

Designing Molecular Spanners to Throw in the Protein Networks

Proteins govern most aspects of cellular life and, through specific interfaces, are typically involved in intricate protein–protein interaction (PPI) networks and signaling pathways. Subtle up- or downregulation of key protein functions and PPIs results in disease; still, the preferred option to contrast the role of a protein in disease and healthy conditions alike remains its outright shutdown through orthosteric ligands that block its active site. Here, we explore subtler alternatives to modulate proteins and PPIs. Driven by a view of proteins as dynamic entities, we discuss ways to identify allosteric binding sites, which, when targeted by tailored ligands, can induce significant changes in the active site of a protein, and lead to agonistic or antagonistic effects. We also summarize the selective regulation of specific PPIs—either direct or allosteric—and show that effects can be stabilizing as well as destabilizing, depending on how the conformational equilibrium of a protein is shifted.     

A small molecule designed to bind to the adenine nucleotide pocket of Hsp90 causes Her2 degradation and the growth arrest and differentiation of breast cancer cells

BACKGROUND: The Hsp90s contain a conserved pocket that binds ATP/ADP and plays an important role in the regulation of chaperone function. Occupancy of this pocket by several natural products (geldanamycin (GM) and radicicol) alters Hsp90 function and results in the degradation of a subset of proteins (i.e. steroid receptors, Her2, Raf). We have used the structural features of this pocket to design a small molecule inhibitor of Hsp90. RESULTS: The designed small molecule PU3 competes with GM for Hsp90 binding with a relative affinity of 15-20 microM. PU3 induces degradation of proteins, including Her2, in a manner similar to GM. Furthermore, PU3 inhibits the growth of breast cancer cells causing retinoblastoma protein hypophosphorylation, G1 arrest and differentiation. CONCLUSIONS: PU3 is representative of a novel class of synthetic compounds that binds to Hsp90 and inhibits the proliferation of cancer cells. These reagents could provide a new strategy for the treatment of cancers.     

Anticancer Gold(III) Porphyrins Target Mitochondrial Chaperone Hsp60

Identification of the molecular target(s) of anticancer metal complexes is a formidable challenge since most of them are unstable toward ligand exchange reaction(s) or biological reduction under physiological conditions. Gold(III) meso‐tetraphenylporphyrin ( gold‐1 a ) is notable for its high stability in biological milieux and potent in vitro and in vivo anticancer activities. Herein, extensive chemical biology approaches employing photo‐affinity labeling, click chemistry, chemical proteomics, cellular thermal shift, saturation‐transfer difference NMR, protein fluorescence quenching, and protein chaperone assays were used to provide compelling evidence that heat‐shock protein 60 (Hsp60), a mitochondrial chaperone and potential anticancer target, is a direct target of gold‐1 a in vitro and in cells. Structure–activity studies with a panel of non‐porphyrin gold(III) complexes and other metalloporphyrins revealed that Hsp60 inhibition is specifically dependent on both the gold(III) ion and the porphyrin ligand.     

Inhibition of Hsp90 with Resorcylic Acid Macrolactones: Synthesis and Binding Studies

A series of resorcylic acid macrolactones, analogues of the natural product radicicol has been prepared by chemical synthesis, and evaluated as inhibitors of heat shock protein 90 (Hsp90), an emerging attractive target for novel cancer therapeutic agents. The synthesis involves acylation of an ortho‐toluic acid dianion, esterification, followed by a ring‐closing metathesis to form the macrocycle. Subsequent manipulation of the protected hydroxymethyl side chain allows access to a range of new analogues following deprotection of the two phenolic groups. Co‐crystallization of one of the new macrolactones with the N‐terminal domain of yeast Hsp90 confirms that it binds in a similar way to the natural product radicicol and to our previous synthetic analogues, but that the introduction of the additional hydroxymethyl substituent appears to result in an unexpected change in conformation of the macrocyclic ring. As a result of this conformational change, the compounds bound less favorably to Hsp90.     

Natural compounds as potential Hsp90 inhibitors for breast cancer-Pharmacophore guided molecular modelling studies

Breast cancer is one of the major impediments affecting women globally. The ATP-dependant heat shock protein 90 (Hsp90) forms the central component of molecular chaperone machinery that predominantly governs the folding of newly synthesized peptides and their conformational maturation. It regulates the stability and function of numerous client proteins that are frequently upregulated and/or mutated in cancer cells, therefore, making Hsp90 inhibition a promising therapeutic strategy for the development of new efficacious drugs to treat breast cancer. In the present in silico investigation, a structure-based pharmacophore model was generated with hydrogen bond donor, hydrogen bond acceptor and hydrophobic features complementary to crucial residues Ala55, Lys58, Asp93, Ile96, Met98 and Thr184 directed at inhibiting the ATP-binding activity of Hsp90. Subsequently, the phytochemical dataset of 3210 natural compounds was screened to retrieve the prospective inhibitors after rigorous validation of the model pharmacophore. The retrieved 135 phytocompounds were further filtered by drug-likeness parameters including Lipinski’s rule of five and ADMET properties, then investigated via molecular docking-based scoring. Molecular interactions were assessed using Genetic Optimisation for Ligand Docking program for 95 drug-like natural compounds against Hsp90 along with two clinical drugs as reference compounds – Geldanamycin and Radicicol. Docking studies revealed three phytochemicals are better than the investigated clinical drugs. The reference and hit compounds with dock scores of 48.27 (Geldanamycin), 40.90 (Radicicol), 73.04 (Hit1), 72.92 (Hit2) and 68.12 (Hit3) were further validated for their binding stability through molecular dynamics simulations. We propose that the non-macrocyclic scaffolds of three identified phytochemicals might aid in the development of novel therapeutic candidates against Hsp90-driven cancers.