Chemical validation of a druggable site on Hsp27/HSPB1 using in silico solvent mapping and biophysical methods

https://doi.org/10.1016/j.bmc.2020.115990Get rights and content

Abstract

Destabilizing mutations in small heat shock proteins (sHsps) are linked to multiple diseases; however, sHsps are conformationally dynamic, lack enzymatic function and have no endogenous chemical ligands. These factors render sHsps as classically “undruggable” targets and make it particularly challenging to identify molecules that might bind and stabilize them. To explore potential solutions, we designed a multi-pronged screening workflow involving a combination of computational and biophysical ligand-discovery platforms. Using the core domain of the sHsp family member Hsp27/HSPB1 (Hsp27c) as a target, we applied mixed solvent molecular dynamics (MixMD) to predict three possible binding sites, which we confirmed using NMR-based solvent mapping. Using this knowledge, we then used NMR spectroscopy to carry out a fragment-based drug discovery (FBDD) screen, ultimately identifying two fragments that bind to one of these sites. A medicinal chemistry effort improved the affinity of one fragment by ~50-fold (16 µM), while maintaining good ligand efficiency (~0.32 kcal/mol/non-hydrogen atom). Finally, we found that binding to this site partially restored the stability of disease-associated Hsp27 variants, in a redox-dependent manner. Together, these experiments suggest a new and unexpected binding site on Hsp27, which might be exploited to build chemical probes.

Introduction

Small heat shock proteins (sHsps) are a recondite class of ATP-independent molecular chaperones1, 2 that aid in suppressing protein aggregation under stress conditions.3 It is thought that these chaperones bind to unfolded proteins, keeping them soluble until they can be refolded. Structurally, these chaperones are defined by a conserved α-crystallin core domain (ACD), flanked by disordered N- and C-terminal extensions.4 Each of these motifs are involved in extensive protein-protein interactions (PPIs). Specifically, the ACD is composed of six β-strands that form an immunoglobulin fold and two of these β-strands, associated in an antiparallel fashion, form a stable dimer interface. Additional contacts involving the N- and C-terminal extensions mediate higher-order assembly of these dimers into larger oligomers.5, 6, 7, 8 In theory, this architecture creates potential opportunities for drug discovery, such as targeting the intra- or inter-molecular PPI surfaces. However, the sHSP oligomers are polydisperse (typically composed of 10–20 dimers) and dynamic in response to cellular stress.5, 9 Thus, it is not clear which of the putative, transient pockets might allow chemical ligands to bind with good affinity.

To explore these questions, we focused on heat shock protein 27 (Hsp27/HSBP1). This member of the sHSP family is of interest because mutations in Hsp27’s ACD, including R127W and S135F, are closely linked to diseases, such as distal hereditary motor neuropathy and Charcot-Marie-Tooth diease.10, 11 These mutations are known to destabilize Hsp27 oligomers and aberrantly increase its chaperone activity.12 This observation is consistent with an emerging model in which smaller, less stable Hsp27 oligomers are more active,5, 9 at least in some contexts. This idea is further supported by the fact that an intermolecular disulfide bond formed at Hsp27’s dimer interface is known to limit its conformational flexibility and, consistent with the model, decrease chaperone activity.13, 14, 15 Together, these observations suggest that one therapeutic objective, at least in a subset of Hsp27-related diseases, may be to identify drug-like molecules that bind and partially restore the stability of disease-associated mutations, such as R127W and S135F. However, there are no known binding pocket(s) for small molecules on Hsp27.

A number of techniques have been developed to identify binding sites on other “undruggable” targets, including computational methods for suggesting sites16 and fragment-based screening approaches.17 Here, we combine computational and experimental platforms towards the discovery of proof-of-principle ligands for Hsp27. Using this combination of approaches, we identified a promising and previously unexplored binding site within Hsp27’s ACD and found that a chemical fragment (Kd > 1 mM) could be matured through a medicinal chemistry campaign to bind with good affinity (Kd ~ 16 µM) and ligand efficiency (0.32 kcal/mol/non-hydrogen atom). Importantly, the resulting probe partially stabilized a disease-associated Hsp27 mutant in vitro, suggesting that this site might be tractable for further development. While additional work is required to increase affinity, these studies provide a starting point for Hsp27 drug discovery and, potentially, for exploration of other sHSPs.

Section snippets

Identification of three binding sites on Hsp27c.

A key challenge with Hsp27 is that binding sites for ligands are not known. To explore this question, we focused on the dimeric, Hsp27 core domain (Hsp27c; residues 79–176), because high-resolution structures are available6, 8 and this region includes the disease-associated mutations, R127W and S135F. In this manuscript, we will use the numbering system of the Hsp27c domain to discuss these hereditary mutations, such that R127W in full length Hsp27 is referred to as R49W and S135F is S57F. In

Summary and conclusion

The identification of druggable sites is a challenge in complex, dynamic systems, such as Hsp27. Indeed, many promising drug targets, including protein complexes and proteins with high intrinsic disorder, share this issue.38 Both computational and experimental methods have been developed to identify potentially druggable sites in these types of systems systems. Here, we sought to combine these methods to identify and exploit putative sites on Hsp27c. Of the HSPs, Hsp27 is a particularly

Mixed-solvent molecular dynamics (MixMD)

MixMD simulations were conducted as described previously.6, 20, 21 Briefly, the Hsp27c NMR structure (PDB 2N3J) was surrounded with a 9-Å layer of cosolvent molecules (acetonitrile, isopropanol, or pyrimidine) and placed in a box of TIP3P water.50 The number of water molecules in each cosolvent simulation was adjusted such that a ratio of 5% cosolvent:water was achieved.

The force field parameters for the probes were obtained from our previous work51 and simulations were carried out in AMBER 1252

Declaration of Competing Interest

The authors declare the following financial interests/personal relationships which may be considered as potential competing interests: J. E. Gestwicki and L. N. Makley are co-founders of Viewpoint Therapeutics.

Acknowledgements

The authors thank Rachel Klevit (University of Washington – Seattle) and Erik R. P. Zuiderweg (University of Michigan, Ann Arbor) for helpful advice. This work was supported by grants from the NIH (NS059690) to J.E.G., (GM065372) to H.A.C. and participation on the PSTP training grant (GM007767) to L.N.M.

References (63)

  • S.D. Bembenek et al.

    Ligand efficiency and fragment-based drug discovery

    Drug Discovery Today

    (2009)
  • H. Shao et al.

    Differential scanning fluorimetry (DSF) screen to identify inhibitors of Hsp60 protein–protein interactions

    Org Biomol Chem

    (2020)
  • I.V. Dudich

    Dimer structure as a minimum cooperative subunit of small heat-shock proteins

    Biochim Biophys Acta (BBA) - Protein Struct Mol Enzymol

    (1995)
  • A. Zavialov

    The effect of the intersubunit disulfide bond on the structural and functional properties of the small heat shock protein Hsp25

    Int J Biol Macromol

    (1998)
  • M. Michishita

    Phosphorylation of the Stress Protein Hsp27 is an Early Event in Murine Myelomoncytic Leukemic Cell Differentiation Induced by Leukemia Inhibitory Factor/D-Factor

    Biochem Biophys Res Commun

    (1991)
  • J.-P. Ryckaert et al.

    Numerical integration of the cartesian equations of motion of a system with constraints: molecular dynamics of n-alkanes

    J Comput Phys

    (1977)
  • U. Jakob et al.

    Assiting spontaneity: the role of Hsp90 and small Hsps as molecular chaperones

    Trends Biochem Sci

    (1994)
  • Montfort RV, Slingsby C, Vierling E. Structure and function of the small heat shock protein/a-crystallin family of...
  • G. Kappé

    The human genome encodes 10 alpha-crystallin-related small heat shock proteins: HspB1-10

    Cell Stress Chaperones

    (2003)
  • T. Kriehuber

    Independent evolution of the core domain and its flanking sequences in small heat shock proteins

    FASEB J

    (2010)
  • S. Jehle

    Solid-state NMR and SAXS studies provide a structural basis for the activation of αB-crystallin oligomers

    Nat Struct Mol Biol

    (2010)
  • P. Rajagopal et al.

    Structure of the α-crystallin domain from the redox-sensitive chaperone, HSPB1

    J Biomol NMR

    (2015)
  • H. Ito

    Hsp27 suppresses the formation of inclusion bodies induced by expression of R120GαB-crystallin, a cause of desmin-related myopathy

    CMLS Cell Mol Life Sci

    (2003)
  • L. Almeida-Souza

    Small Heat-Shock Protein HSPB1 Mutants Stabilize Microtubules in Charcot-Marie-Tooth Neuropathy

    J Neurosci

    (2011)
  • T.R. Alderson

    Local unfolding of the HSP27 monomer regulates chaperone activity

    Nat Commun

    (2019)
  • M.A. Marti-Renom

    The AnnoLite and AnnoLyze programs for comparative annotation of protein structures

    BMC Bioinf

    (2007)
  • C.W. Murray et al.

    The rise of fragment-based drug discovery

    Nat Chem

    (2009)
  • A.-P. Arrigo

    Hsp27 Consolidates Intracellular Redox Homeostasis by Upholding Glutathione in Its Reduced Form and by Decreasing Iron Intracellular Levels

    Antioxid Redox Signal

    (2005)
  • P.M.U. Ung et al.

    Identifying binding hot spots on protein surfaces by mixed-solvent molecular dynamics: HIV-1 protease as a test case: Identifying Binding Hot Spots on Protein Surfaces by MixMD

    Biopolymers

    (2016)
  • P. Ghanakota et al.

    Large-Scale Validation of Mixed-Solvent Simulations to Assess Hotspots at Protein-Protein Interaction Interfaces

    J Chem Inf Model

    (2018)
  • P. Ghanakota et al.

    Moving Beyond Active-Site Detection: MixMD Applied to Allosteric Systems

    J Phys Chem B

    (2016)
  • 1

    These authors contributed equally to the work.

    View full text