Visualizing the transiently populated closed-state of human HSP90 ATP binding domain.

HSP90 are abundant molecular chaperones, assisting the folding of several hundred client proteins, including substrates involved in tumor growth or neurodegenerative diseases. A complex set of large ATP-driven structural changes occurs during HSP90 functional cycle. However, the existence of such structural rearrangements in apo HSP90 has remained unclear. Here, we identify a metastable excited state in the isolated human HSP90α ATP binding domain. We use solution NMR and mutagenesis to characterize structures of both ground and excited states. We demonstrate that in solution the HSP90α ATP binding domain transiently samples a functionally relevant ATP-lid closed state, distant by more than 30 Å from the ground state. NMR relaxation enables to derive information on the kinetics and thermodynamics of this interconversion, while molecular dynamics simulations establish that the ATP-lid in closed conformation is a metastable exited state. The precise description of the dynamics and structures sampled by human HSP90α ATP binding domain provides information for the future design of new therapeutic ligands.

NMR assignment of human HSP90 N-terminal domain bound to a long residence time resorcinol ligand.

HSP90 is a major molecular chaperone that helps both folding and stabilization of various client proteins often implicated in growth control and cell survival such as kinases and transcription factors. However, among HSP90 clients are also found numerous oncoproteins and, through its assistance to them, HSP90 has consequently been reported as a promising anticancer target. Several ligand chemotypes, including resorcinol type ligands, were found to inhibit HSP90, most of them in an ATP competitive manner. Binding of some of these ligands modify significantly the NMR spectrum of the HSP90 ATP binding domain compared to the apo protein spectrum, hampering assignment transfer from the previously assigned human HSP90 apo state. Here we report the assignment of the 1HN, 15N, 13C', 13Cα, 13Cß, 1Hmethyl, and 13Cmethyl chemical shifts of the 29 kDa HSP90 N-terminal domain bound to a long residence time resorcinol type inhibitor: 5-[4-(2-Fluoro-phenyl)-5-oxo-4,5-dihydro-1H-[1,2,4]triazol-3-yl]-N-furan-2-ylmethyl-2,4-dihydroxy-N-methyl-benzamide. 92% of the backbone resonances and 100% of the [1H, 13C]-resonances of Aß, Mε, Tγ, Lδ2, Vγ2 and Iδ1 methyl groups were successfully assigned, including for the first time the assignment of the segment covering the nucleotide/drug binding site. Secondary structure predictions based on the NMR assignment reveal a structural rearrangement of HSP90 N-terminal domain upon ligand binding. The long residence time ligand induces the formation of a continuous helix covering the ligand binding site of HSP90 N-terminal domain accounting for the large differences observed in the NMR spectra between the apo and bound proteins.

Optimized precursor to simplify assignment transfer between backbone resonances and stereospecifically labelled valine and leucine methyl groups: application to human Hsp90 N-terminal domain.

Methyl moieties are highly valuable probes for quantitative NMR studies of large proteins. Hence, their assignment is of the utmost interest to obtain information on both interactions and dynamics of proteins in solution. Here, we present the synthesis of a new precursor that allows connection of leucine and valine pro-S methyl moieties to backbone atoms by linear 13C-chains. This optimized 2H/13C-labelled acetolactate precursor can be combined with existing 13C/2H-alanine and isoleucine precursors in order to directly transfer backbone assignment to the corresponding methyl groups. Using this simple approach leucine and valine pro-S methyl groups can be assigned using a single sample without requiring correction of 1H/2H isotopic shifts on 13C resonances. The approach was demonstrated on the N-terminal domain of human HSP90, for which complete assignment of Ala-ß, Ile-δ1, Leu-δ2, Met-ε, Thr-γ and Val-γ2 methyl groups was obtained.

Structure-kinetic relationship reveals the mechanism of selectivity of FAK inhibitors over PYK2.

There is increasing evidence of a significant correlation between prolonged drug-target residence time and increased drug efficacy. Here, we report a structural rationale for kinetic selectivity between two closely related kinases: focal adhesion kinase (FAK) and proline-rich tyrosine kinase 2 (PYK2). We found that slowly dissociating FAK inhibitors induce helical structure at the DFG motif of FAK but not PYK2. Binding kinetic data, high-resolution structures and mutagenesis data support the role of hydrophobic interactions of inhibitors with the DFG-helical region, providing a structural rationale for slow dissociation rates from FAK and kinetic selectivity over PYK2. Our experimental data correlate well with computed relative residence times from molecular simulations, supporting a feasible strategy for rationally optimizing ligand residence times. We suggest that the interplay between the protein structural mobility and ligand-induced effects is a key regulator of the kinetic selectivity of inhibitors of FAK versus PYK2.

Estimation of Protein-Ligand Unbinding Kinetics Using Non-Equilibrium Targeted Molecular Dynamics Simulations.

We here report on nonequilibrium targeted molecular dynamics simulations as a tool for the estimation of protein-ligand unbinding kinetics. Correlating simulations with experimental data from SPR kinetics measurements and X-ray crystallography on two small molecule compound libraries bound to the N-terminal domain of the chaperone Hsp90, we show that the mean nonequilibrium work computed in an ensemble of trajectories of enforced ligand unbinding is a promising predictor for ligand unbinding rates. We furthermore investigate the molecular basis determining unbinding rates within the compound libraries. We propose ligand conformational changes and protein-ligand nonbonded interactions to impact on unbinding rates. Ligands may remain longer at the protein if they exhibit strong electrostatic and/or van der Waals interactions with the target. In the case of ligands with a rigid chemical scaffold that exhibit longer residence times, transient electrostatic interactions with the protein appear to facilitate unbinding. Our results imply that understanding the unbinding pathway and the protein-ligand interactions along this path is crucial for the prediction of small molecule ligands with defined unbinding kinetics.

Ligand-Induced Conformational Changes in HSP90 Monitored Time Resolved and Label Free-Towards a Conformational Activity Screening for Drug Discovery.

Investigation of protein-ligand interactions is crucial during early drug-discovery processes. ATR-FTIR spectroscopy can detect label-free protein-ligand interactions with high spatiotemporal resolution. Here we immobilized, as an example, the heat shock protein HSP90 on an ATR crystal. This protein is an important molecular target for drugs against several diseases including cancer. With our novel approach we investigated a ligand-induced secondary structural change. Two specific binding modes of 19 drug-like compounds were analyzed. Different binding modes can lead to different efficacy and specificity of different drugs. In addition, the kobs values of ligand dissociation were obtained. The results were validated by X-ray crystallography for the structural change and by SPR experiments for the dissociation kinetics, but our method yields all data in a single and simple experiment.

Estimation of Drug-Target Residence Times by τ-Random Acceleration Molecular Dynamics Simulations.

Drug-target residence time (τ), one of the main determinants of drug efficacy, remains highly challenging to predict computationally and, therefore, is usually not considered in the early stages of drug design. Here, we present an efficient computational method, τ-random acceleration molecular dynamics (τRAMD), for the ranking of drug candidates by their residence time and obtaining insights into ligand-target dissociation mechanisms. We assessed τRAMD on a data set of 70 diverse drug-like ligands of the N-terminal domain of HSP90α, a pharmaceutically important target with a highly flexible binding site, obtaining computed relative residence times with an accuracy of about 2.3τ for 78% of the compounds and less than 2.0τ within congeneric series. Analysis of dissociation trajectories reveals features that affect ligand unbinding rates, including transient polar interactions and steric hindrance. These results suggest that τRAMD will be widely applicable as a computationally efficient aid to improving drug residence times during lead optimization.

Ligand Desolvation Steers On-Rate and Impacts Drug Residence Time of Heat Shock Protein 90 (Hsp90) Inhibitors.

Residence time and more recently the association rate constant kon are increasingly acknowledged as important parameters for in vivo efficacy and safety of drugs. However, their broader consideration in drug development is limited by a lack of knowledge of how to optimize these parameters. In this study on a set of 176 heat shock protein 90 inhibitors, structure-kinetic relationships, X-ray crystallography, and molecular dynamics simulations were combined to retrieve a concrete scheme of how to rationally slow down on-rates. We discovered that an increased ligand desolvation barrier by introducing polar substituents resulted in a significant kon decrease. The slowdown was accomplished by introducing polar moieties to those parts of the ligand that point toward a hydrophobic cavity. We validated this scheme by increasing polarity of three Hsp90 inhibitors and observed a 9-, 13-, and 45-fold slowdown of on-rates and a 9-fold prolongation in residence time. This prolongation was driven by transition state destabilization rather than ground state stabilization.

Kinetics for Drug Discovery: an industry-driven effort to target drug residence time.

A considerable number of approved drugs show non-equilibrium binding characteristics, emphasizing the potential role of drug residence times for in vivo efficacy. Therefore, a detailed understanding of the kinetics of association and dissociation of a target-ligand complex might provide crucial insight into the molecular mechanism-of-action of a compound. This deeper understanding will help to improve decision making in drug discovery, thus leading to a better selection of interesting compounds to be profiled further. In this review, we highlight the contributions of the Kinetics for Drug Discovery (K4DD) Consortium, which targets major open questions related to binding kinetics in an industry-driven public-private partnership.

Structural and biophysical characterization of the Syk activation switch.

Syk is an essential non-receptor tyrosine kinase in intracellular immunological signaling, and the control of Syk kinase function is considered as a valuable target for pharmacological intervention in autoimmune or inflammation diseases. Upon immune receptor stimulation, the kinase activity of Syk is regulated by binding of phosphorylated immune receptor tyrosine-based activating motifs (pITAMs) to the N-terminal tandem Src homology 2 (tSH2) domain and by autophosphorylation with consequences for the molecular structure of the Syk protein. Here, we present the first crystal structures of full-length Syk (fl-Syk) as wild type and as Y348F,Y352F mutant forms in complex with AMP-PNP revealing an autoinhibited conformation. The comparison with the crystal structure of the truncated Syk kinase domain in complex with AMP-PNP taken together with ligand binding studies by surface plasmon resonance (SPR) suggests conformational differences in the ATP sites of autoinhibited and activated Syk forms. This hypothesis was corroborated by studying the thermodynamic and kinetic interaction of three published Syk inhibitors with isothermal titration calorimetry and SPR, respectively. We further demonstrate the modulation of inhibitor binding affinities in the presence of pITAM and discuss the observed differences of thermodynamic and kinetic signatures. The functional relevance of pITAM binding to fl-Syk was confirmed by a strong stimulation of in vitro autophosphorylation. A structural feedback mechanism on the kinase domain upon pITAM binding to the tSH2 domain is discussed in analogy of the related family kinase ZAP-70 (Zeta-chain-associated protein kinase 70). Surprisingly, we observed distinct conformations of the tSH2 domain and the activation switch including Tyr348 and Tyr352 in the interdomain linker of Syk in comparison to ZAP-70.

Structure of the leech protein saratin and characterization of its binding to collagen.

The leech protein Saratin from Hirudo medicinalis prevents thrombocyte aggregation by interfering with the first binding step of the thrombocytes to collagen by binding to collagen. We solved the three-dimensional structure of the leech protein Saratin in solution and identified its collagen binding site by NMR titration experiments. The NMR structure of Saratin consists of one alpha-helix and a five-stranded beta-sheet arranged in the topology betabetaalphabetabetabeta. The C-terminal region, of about 20 amino acids in length, adopts no regular structure. NMR titration experiments with collagen peptides show that the collagen interaction of Saratin takes place in a kind of notch that is formed by the end of the alpha-helix and the beta-sheet. NMR data-driven docking experiments to collagen model peptides were used to elucidate the putative binding mode of Saratin and collagen. Mainly, parts of the first and the end of the fifth beta-strand, the loop connecting the alpha-helix and the third beta-strand, and a short part of the loop connecting the fourth and fifth beta-strand participate in binding.

Ubiquitin protein ligase Nedd4 binds to connexin43 by a phosphorylation-modulated process.

Connexin43 is degraded by the proteasomal as well as the lysosomal pathway with ubiquitin playing a role in both degradation pathways. So far, no ubiquitin protein ligase has been identified for any of the connexins. By using pull-down assays, here we show binding of a ubiquitin protein ligase, Nedd4, to the C-terminus of connexin43. This observation was confirmed in vivo by coimmunoprecipitation and immunofluorescence, showing colocalization of Nedd4 and connexin43. Binding of Nedd4 to its interaction partners is generally carried out by its WW domains. Our results indicate that the interaction with connexin43 occurs through all three WW domains of Nedd4. Furthermore, whereas WW1 and WW2 domains mainly interact with the unphosphorylated form of connexin43, WW3 binds phosphorylated and unphosphorylated forms equally. In addition, using the surface plasmon resonance approach we show that only the WW2 domain binds to the PY motif located at the C-terminus of connexin43. Suppression of Nedd4 expression with siRNA resulted in an accumulation of gap junction plaques at the plasma membrane, suggesting an involvement of the ubiquitin protein ligase Nedd4 in gap junction internalization.

Features of the Env leader protein and the N-terminal Gag domain of feline foamy virus important for virus morphogenesis.

Previous studies have shown that foamy virus (FV) particle budding, especially the involvement of the viral env glycoprotein is different from that of other (ortho) retroviruses: the N-terminal Env leader protein Elp is a constituent of released FV particles. A defined sequence in Elp required for particle budding binds to the MA domain of Gag. To extend these findings, we show that feline FV Elp is a membrane-anchored protein with the N-terminus located inside the particle. Thus, the internal/cytoplasmic domain of Elp has the correct topology for interacting with Gag during budding. In addition to Elp, an Elp-related protein of about 9 kDa was shown to be virion associated and is probably generated by cellular signal peptidases. Besides the function of Elp binding, the N-terminal domain of Gag was shown to be required for proper localization of feline FV Gag to the cytoplasm and the perinuclear/nuclear region.

Crystal structure of the extracellular segment of integrin alpha Vbeta3 in complex with an Arg-Gly-Asp ligand.

The structural basis for the divalent cation-dependent binding of heterodimeric alphabeta integrins to their ligands, which contain the prototypical Arg-Gly-Asp sequence, is unknown. Interaction with ligands triggers tertiary and quaternary structural rearrangements in integrins that are needed for cell signaling. Here we report the crystal structure of the extracellular segment of integrin alphaVbeta3 in complex with a cyclic peptide presenting the Arg-Gly-Asp sequence. The ligand binds at the major interface between the alphaV and beta3 subunits and makes extensive contacts with both. Both tertiary and quaternary changes are observed in the presence of ligand. The tertiary rearrangements take place in betaA, the ligand-binding domain of beta3; in the complex, betaA acquires two cations, one of which contacts the ligand Asp directly and the other stabilizes the ligand-binding surface. Ligand binding induces small changes in the orientation of alphaV relative to beta3.