DnaK as Antibiotic Target: Hot Spot Residues Analysis for Differential Inhibition of the Bacterial Protein in Comparison with the Human HSP70

Federica Chiappori; Marco Fumian; Luciano Milanesi; Ivan Merelli

doi:10.1371/journal.pone.0124563

Abstract

DnaK, the bacterial homolog of human Hsp70, plays an important role in pathogens survival under stress conditions, like antibiotic therapies. This chaperone sequesters protein aggregates accumulated in bacteria during antibiotic treatment reducing the effect of the cure. Although different classes of DnaK inhibitors have been already designed, they present low specificity. DnaK is highly conserved in prokaryotes (identity 50–70%), which encourages the development of a unique inhibitor for many different bacterial strains. We used the DnaK of Acinetobacter baumannii as representative for our analysis, since it is one of the most important opportunistic human pathogens, exhibits a significant drug resistance and it has the ability to survive in hospital environments. The E.coli DnaK was also included in the analysis as reference structure due to its wide diffusion. Unfortunately, bacterial DnaK and human Hsp70 have an elevated sequence similarity. Therefore, we performed a differential analysis of DnaK and Hsp70 residues to identify hot spots in bacterial proteins that are not present in the human homolog, with the aim of characterizing the key pharmacological features necessary to design selective inhibitors for DnaK. Different conformations of DnaK and Hsp70 bound to known inhibitor-peptides for DnaK, and ineffective for Hsp70, have been analysed by molecular dynamics simulations to identify residues displaying stable and selective interactions with these peptides. Results achieved in this work show that there are some residues that can be used to build selective inhibitors for DnaK, which should be ineffective for the human Hsp70.

Citation: Chiappori F, Fumian M, Milanesi L, Merelli I (2015) DnaK as Antibiotic Target: Hot Spot Residues Analysis for Differential Inhibition of the Bacterial Protein in Comparison with the Human HSP70. PLoS ONE 10(4): e0124563. https://doi.org/10.1371/journal.pone.0124563

Academic Editor: Andrew R. Dalby, University of Westminster, UNITED KINGDOM

Received: December 30, 2014; Accepted: March 16, 2015; Published: April 23, 2015

Copyright: © 2015 Chiappori et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited

Data Availability: All relevant data are within the paper and its Supporting Information files.

Funding: This paper has been supported by the Italian Ministry of Education and Research (MIUR) through the Flagship (PB05) InterOmics, HIRMA (RBAP11YS7 K), and the European MIMOMICS projects. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Competing interests: The authors have declared that no competing interests exist.

Introduction

Heat Shock Proteins (Hsp) are essential for the survival cells and their expression levels rely on cellular conditions. In particular, proteins belonging to the Hsp70 family are involved, under stress conditions, in signal transduction, cell cycle regulation, and programmed cell death. Other conditions that involve these proteins are principally native protein folding, refolding and prevention of protein aggregation [1]. Their essential role for pathogenic microorganisms growing in a host is of particular interest for drug discovery. DnaK belongs to the Hsp70 family and is the bacterial homolog of human Hsp70. In particular, DnaK displays up to 70% of sequence identity with respect to the other eukaryotic proteins of this family [2].

DnaK has been characterized in several pathogenic bacteria and seems to have important functions in stress resistance and pathogenicity in multiple-drug-resistant bacteria, such as Acinetobacter baumannii [3–4], which is one of the most important opportunistic human pathogens displaying several antibiotic resistances. The heat-shock response, and in particular the DnaK machinery, is involved in the antibiotic resistance mechanism of A baumannii [5]. In detail, it results necessary for bacteria survival in unfavourable conditions, such as exposure to oxidative stress, nutrient limitation, extreme temperatures, and presence of heavy metals or antibiotics [6–8]. DnaK mutations increase the bacterial sensitivity to fluoroquinolones, oxacillin and methicillin in normally resistant strains [6–8], since this protein sequesters the aggregates that accumulate in cells exposed to these antibiotics [6] and assists the refolding of proteins misfolded after a stress event [9].

Like all other Hsp70 proteins, DnaK is composed of about 650 residues, arranged in two domains: the nucleotide binding domain (NBD) and the substrate binding domain (SBD), these are connected by a highly flexible linker involved in the allosteric communication between the two domains. When the NBD domain hydrolyses an ATP molecule, the SBD domain assumes a closed conformation, which binds a short extended hydrophobic polypeptide sequence [10]. Therefore, DnaK displays two extreme conformations (Fig 1): in the open state, ATP is bound to the NDB cavity and the substrate affinity is low, while in the closed conformation, after the ATP hydrolysis, the affinity for the substrate is high. Furthermore, the binding of peptides to the SBD induces the ATP hydrolysis in the NDB and the ADP presence induces the SBD rearrangement to the closed conformation, which correspond to a ~10 fold affinity increase for the peptides [11]. The nucleotide exchange from ADP to ATP induces the SBD opening and the substrate release, this brings back the protein to the open conformation (Fig 1). The open/closed state rearrangement depends on the nucleotide that is bound to the NBD domain: in particular, the ATP-bound DnaK is characterized by a low affinity and a fast exchange rate for the substrate, while the ADP-bound form displays a high affinity and low exchange rates. To complete the allosteric cycle, two classes of cochaperone help DnaK/Hsp70 proteins, respectively the DnaJ/Hsp40 and the GrpE/Hip. DnaJ has a J-domain that presents the substrate to DnaK and induces its ATPase activity, resulting essential for the DnaK functionality, while GrpE is a nucleotide exchange factor that increases the basal ADP/ATP exchange rate of DnaK [10].

Download:

Fig 1. Schematic representation of the Hsp70 allosteric cycle.

NBD is in blue, βSBD is in green, helices HA and HB are in yellow and helices HC-HE are in red.

https://doi.org/10.1371/journal.pone.0124563.g001

SDB consists of two sub-domains, βSBD and αSBD. The former is a β-sandwich of two antiparallel β–sheets composed of four strands each one and connected by seven loops. Among these loops, L_1,2 and L_3,4 in combination with strands 1 and 2 conform the deep peptide binding cavity. αSBD consists of five helices, one of which, indicated as HB, forms a lid over the polypeptide binding site [12]. The binding cavity is mainly hydrophobic and regions surrounding the cavity are negatively charged [9].

Short hydrophobic peptides (~7 residues) were identified by experimentally screening libraries of compounds able to bind SBD [13]. The hydrophobic core of these peptides is composed of 2–4 residues containing Leu and excluding acid residues, also from the flanking regions [9, 13]. A general core motif ±HyHyHyHyHy± (where ± stands for positively charged residues and Hy refers to hydrophobic residues Leu, Ile, Val, Ala, Gly) was identified [13–15]. Interactions between the peptide-substrate and the SBD binding site involve, respectively, the side-chains of the peptide in van der Waals and hydrophobic contacts, and the peptide backbone in hydrogen bonds with the cavity-forming loops [9, 13]. Once the lid is closed on the substrate, other interactions with the helices (electrostatic and van der Waals) stabilize the conformation [14].

Recently some peptides with antibacterial characteristics have been isolated in insects. The optimized forms of these peptides able to bind DnaK are relatively stable to protease digestion, due to the high proline content, are nontoxic to human cell lines, and are predominantly active against Gram-negative bacteria [2, 16]. In particular, there are few well-known peptides that inhibit DnaK, which are Api88, drosocin, heliocin, oncocin and pyrrhocoricin (Fig 2). These peptides are composed of 18–22 amino acids capable of binding the SBD and keeping the lid in a closed position, preventing the substrate release [8]. Unfortunately, these peptides display low specificity for DnaK: for example, api88 shows a pro-apoptotic activity in DnaK-null cells and it is supposed to bind GroEL [17]. Moreover, peptide drugs generally suffer from their intrinsic metabolic instability [18].

Download:

Fig 2. Structures and sequences of the five peptides employed in this study.

Atoms are coloured by atom type.

https://doi.org/10.1371/journal.pone.0124563.g002

In this study, we identify SBD hot spot residues that can be exploited for the differential inhibition of the bacterial DnaK in comparison to the human Hsp70. The key residues for the interaction with the inhibitor-peptides were recognized, selecting residues conserved in bacteria and not-conserved in the human Hsp70 family. The study is based on sequence analysis, evaluation of the binding site through molecular dynamics (MD) and interaction analysis. MD is a good approximation for studying protein-peptide binding with atomistic detail, and to take into account the dynamic behaviour of the whole system [14, 19]. In detail, three proteins, A.baumannii DnaK, representative for the opportunistic human pathogens, E.coli DnaK, included as reference and control structure, and human Hsp70-2 in closed conformation, were simulated in complex with Api88, drosocin, heliocin, oncocin, and pyrrhocoricin bound to the SBD. Results achieved for the bacterial complexes were compared to data obtained for the human model. Our analysis resulted in a set of key residues of the DnaK binding site that are involved in stable interactions with the ligands. Moreover, we developed a pharmacophore model, which can be employed in a further virtual screening study to identify ligands with high specificity and selectivity.

Materials and Methods

Sequence analysis

Protein sequences were downloaded from Uniprot [20]. Bacterial sequences were obtained searching for the “DnaK” gene name in the “reference proteome” text box and only reviewed sequences within the taxonomy group “Bacteria” were chosen. Sequences were aligned with Clustal Omega [21] using default settings. Moreover, the SIAS web server [22] was employed to calculate pairwise identities and similarities of multiple sequence alignments, taking gaps into account. WebLogo [23] was employed to obtain a graphical representation of residue frequencies in the βSBD and HA- HB-SBD positions.

Structure

The closed conformation of E.coli DnaK in complex with ADP-Mg²⁺ (pdb ID: 2KHO) [24] was obtained from the Protein Data Bank [25]. Structural models of A.baumannii DnaK and human Hsp70-2 were obtained relying on E.coli DnaK in closed conformation, using the Swiss-model server [26].

Coordinates of the inhibitors were extracted, respectively, from the X-ray structure of E.coli SBD bound to Api88 (pdb ID: 4E81) [16], Drosocin (pdb ID: 4EZR) [11], Heliocin (pdb ID: 4EZT) [11], Oncocin (pdb ID: 3QNJ) [5] and pyrrhocoricin (pdb ID: 4EZN) [11]. These peptides were superimposed to the whole structures of E.coli and A.baumannii DnaK and human Hsp70 to obtain the complexes to simulate. Peptides are composed of 18–22 amino acids, but for each inhibitor only the 8 core-residues were co-crystallized and therefore included in our complexes.

Molecular dynamics

MD simulations were performed with GROMACS 4.0 [27], by employing the GROMOS96 (ff43a1) force field [28]. Each complex was included in a rhombic dodecahedric box filled with water (TIP3P model), setting the distance between the complex and the box walls to 1nm. Na⁺ ions were added to complexes for neutralizing the net charge of the system. All complexes were energy relaxed with 1000 steps of “steepest-descendent” energy minimization. MD simulations were performed using the LINCS algorithm [29] to constrain bond lengths, while periodic boundary conditions were applied in all directions. Long-range electrostatic forces were treated using the Fast Particle-Mesh Ewald method (PME) [30]. Van der Waals forces and Coulomb potential were treated using a cut-off distance of 0.9nm and the neighbour list was updated every 5 steps. The simulation time step was set to 2fs and coordinates were saved every 10 ps. An initial velocity was given to all atoms according to a Maxwell distribution at 300K. All simulations were run in an NVT environment employing V-rescale as temperature coupling algorithm, with reference temperature set at 300K.

Molecular dynamics analysis

Analyses of MD simulations were performed using different software from the GROMACS package. System stability was verified analysing the energy components and the RMSD of the trajectory conformations compared to the starting structure. The cluster analysis was performed using the g_cluster module of GROMACS to obtain representative conformations of protein complexes. Clustering of the trajectories was achieved fitting C_α atoms with the gromos method [31], with a cut-off distance of 0.5 nm. Trajectories were visualized with VMD [32].

Contact maps were generated using the g_mdmat tool of GROMACS, which identifies the minimum distance between residues by calculating the smallest distance between any pair of their atoms. A truncation distance of 0.5 nm was employed to obtain the contact map. This analysis was useful for identifying residues in the neighbourhood of inhibitor-peptides during simulations.

The interaction analysis was performed with the g_hbond module of GROMACS, applied to MD simulations, while in-house developed scripts were employed to identify key residues for the interaction with inhibitor-peptides [33]. In particular, both hydrogen bonds and electrostatic contacts were evaluated during the whole trajectories. Starting from the H-bonds existence matrix, the fraction of time that a particular H-bond or contact is present during the trajectory was calculated using the script plot_hbmap.pl [34]. Residues involved in hydrogen bonds for at least 25% of the trajectory (contacts for at least 50%) were considered for the analysis. Among these residues, only amino acids conserved in the bacterial DnaK and not in the human Hsp70 were filtered. For contacts analysis, cut-offs were set to 0.39 and 0.41, respectively, as set in a previous work [35].

Computational alanine-scanning was accomplished by employing the Robetta server [36–38], which uses Rosetta as backend software. In particular, Rosetta alanine-scanning individually substitutes all residues at the protein-protein (or protein-peptide) interface with alanine and estimates the change in binding free energy.

Structure-based pharmacophore models were generated with LigandScout [39]. Interactions between ligands and proteins were automatically analysed considering both electrostatic and geometric (distances, size, and angles) properties of the binding regions. In particular, hydrophobic anchoring points, H-bond acceptors/donors, positively or negatively ionizable/charged groups, and aromatic rings involved in aromatic interactions were evaluated for building the pharmacophore. As a result, the spatial distribution of the pharmacophoric feature points, which represents the protein residues involved in the interaction, was obtained. We achieved a pharmacophore model for each DnaK-peptide complex representative conformation identified in the cluster analysis of the MD trajectories. Pharmacophore models were aligned and only common characteristics were considered in the final model.

A previously developed approach was applied to estimate the enthalpy variation caused by the ligand binding [35] to the protein, which can also be decomposed in the individual contribution of each residue include in the binding sites. In particular, the Lennard-Jones and Coulomb terms were estimated using the Gromacs tool g_energy. The interaction energy of the three proteins and of residue Met404, Val425, Ser427, Ala429 and Ala435 for the ligands, were analyzed.

Results

A preliminary analysis of the available bacterial DnaK sequences was performed to evaluate the reliability of employing the DnaK of E.coli and A.baumannii as representative for all the possible bacterial proteins. Moreover, this analysis was useful for verifying residue conservation in positions identified as hot spots at the end of the study. The Hsp70-2 was chosen as a representative for human Hsp70s since it is stress-induced, it localizes in the cytoplasm, and it is one of the most diffused in cells [40].

A comparative analysis was carried out by simulating complexes of the X-ray structures of E.coli DnaK in closed conformation, with the SBD lid closed over the peptide binding cavity, and the homology model of A.baumannii DnaK compared to the homology model of the human Hsp70 bound to Api88, drosocin, heliocin, oncocin and pyrrhocoricin. For each of the five inhibitors, closed conformations of both DnaK and Hsp70 in complex with ADP were simulated for 100 ns. Overall, we obtained about 1,3 μs of equilibrated simulation. The equilibrated parts of MD trajectories were evaluated by interaction analysis, conformation clustering, and were employed for the identification of residues belonging to the peptide binding cavity. Computational alanine scanning and per residue binding energy decomposition were carried out to evaluate the role of residues in protein binding. Moreover, multiple alignments were employed to assess conservation of residues, with particular attention to the SBD. At last, identified residues allowed the design of a pharmacophore model, which can be employed in virtual screening studies for the identification of selective and specific DnaK inhibitors.

SBD residue conservation

The 64 sequences obtained from Uniprot (querying the database for DnaK genes in bacteria, see the Experimental section for details) and the DnaK of A.baumannii were aligned (S1 Fig) to identify conserved residues. For the same reason, the 13 human Hsp70 sequences were aligned (S2 Fig). DnaK is highly conserved in bacteria, the residue identity is 73% and similarity is about 80% along the whole sequence. Taking into account NBD, conserved residues are distributed mainly outside the secondary structures: for instance, catalytic residues and residues involved in contacts with the nucleotide belong to conserved loops. βSBD displays a high level of residue conservation, compared to αSBD. In detail, NBD displays 20% of fully conserved and about 20% of strongly conserved residues; βSBD displays 22% of fully conserved residues and 26% of residue similarity; on the other hand, only 2% and 5%, respectively, of residues are fully and strongly conserved in αSBD. Similar results were obtained for the human Hsp70 proteins: considering the whole sequence, NBD and βSBD display about 50–55% of residue identity, while only 35% of αSBD residues are strictly conserved. Comparing bacterial DnaK and human Hsp70 proteins, about 45% of residues are strongly conserved along the whole sequence; similar results are displayed for NBD and βSBD (45–40%, respectively); on the other hand, αSBD displays only 18% of conserved residues. Considering only Hsp70-2, the identities increase slightly by, 47, 52 and 58%, respectively, for the whole sequences, NBD and βSBD, while αSBD remains at 20% of residue identity.

Both identities and similarities between sequences analysed in this study have been evaluated. As shown in Table 1, the global similarity between bacterial sequences is about 80%, while between bacterial and human sequences is about 57%. Identity values are about 10% lower than similarity values.

Download:

Table 1. Identity and similarity percentages of DnaK and Hsp70 protein sequences.

https://doi.org/10.1371/journal.pone.0124563.t001

Identification of interaction key residues

In order to identify stable DnaK residues anchoring the inhibitor-peptides, the temporal evolution of the interactions was evaluated. The presence of hydrogen bonds and contacts was analysed along the whole MD trajectories. Residues involved in hydrogen bonds with the inhibitor-peptide for at least 25% of the simulation (50% for contacts) were taken into account. Among these residues, only those conserved in DnaK and not in Hsp70 were further considered for identifying selective amino acids for the bacterial protein (Table 2 and Table 3). Four residues were identified as involved in hydrogen bonds: Met404, Ser427, Ala429 and Gln433. The first three residues are substituted by amino acids with different characteristics in Hsp70, respectively Ala404, Thr427 and Tyr429, while Gln433 is maintained in the human protein, but it is rarely involved in hydrogen bonds. Moreover, Met404 and Ala429 hydrogen bonds are conserved in different DnaK complexes, while only two hydrogen bonds were observed in this position for Hsp70 complexes. In detail, the hydrogen bond involving Ser427 is maintained in Hsp70 with a different amino acid, Thr427. In addition to residues involved in hydrogen bonds, Val436 was identified as involved in contacts. Val435 is conserved, but it is rarely involved in contacts in human Hsp70.

Download:

Table 2. Residues forming hydrogen bonds with the bound peptides.

https://doi.org/10.1371/journal.pone.0124563.t002

Download:

Table 3. Residues involved in contacts with the bound peptides.

https://doi.org/10.1371/journal.pone.0124563.t003

From the peptide point of view, amino acids involved in interactions with the binding site were analysed. In particular, Met 404 frequently interacts with both side-chain and backbone of residues in the position 4 of peptides, while Ala429 and Gln433 interact with the backbone of residues in positions 5 and 6, respectively. Moreover, Ser427 and Val436 establish contacts with the backbone of residues in positions 2–5 and with the core residue side-chains of the inhibitor-peptides (5 and 6). Concerning residues conservation in peptides, only Proline and Arginine enrichment can be reported.

Definition of the binding cavity

The binding cavity of DnaK consists of βSBD strands and loops and is closed by the HB helix of αSBD. Residues belonging to the peptide binding cavity were identified from MD simulations of DnaK and Hsp70 bound to the inhibitor-peptides. To this end, the distance matrix between peptide residues and SBD was obtained and residues within 0.5nm from the peptides were considered as pertaining to the binding cavity. Identified residues belong to β₁—L_1,2 - β₂, β₃—L_3,4 - β₄, L_5,6 - β₆ and HB (Table 4). All the inhibitors display similar results: differences between DnaK and Hsp70 include Met404, Val425, Ser427, Ala429, Ala435, Ile438 and Gly468 substituted by an Ala, Thr, Thr, Tyr, Ser and Val in the human protein, respectively, while Gly468 (469 in human) is conserved, but point towards the outside of the cavity.

Download:

Table 4. Binding cavity residues.

https://doi.org/10.1371/journal.pone.0124563.t004

Computational alanine scanning

Alanine scanning is a computational approach to predict the importance of each residue for a specific protein-protein interface, relying on an energy criterion. Alanine scanning of residues composing the binding cavity was performed on representative conformations achieved through MD simulations of the two DnaK proteins bound to ADP and the inhibitor-peptides (Api88, drosocin, heliocin, oncocin, and pyrrhocoricin). Representative conformations were obtained by cluster analysis, see Experimental section for details. A binding energy threshold (ΔΔGbind) of 1 kcal/mol for residues replaced with alanine was applied to identify residues responsible of protein-protein binding as specified in Kortemme et al. [41]. Residues with a ΔΔGbind value greater than 1 kcal/mol for at least three inhibitor-peptides in both bacterial DnaK proteins were considered. Overall relevant residues are: Thr403, Thr409, Phe426, Gln433, Gln534 and His541 (Fig 3). The first four residues belong to β-SBD and are highly conserved in the bacterial DnaK, while less conserved in the human Hsp70. On the other end, the last two residues belong to the HB helix and are slightly conserved both in bacterial and in human Hsp70 protein family.

Download:

Fig 3. ΔΔGbind results of alanine scanning on binding cavity residues.

Interesting residues are displayed for E.coli, A.baumannii DnaK and human Hsp70 in complex with known inhibitors.

https://doi.org/10.1371/journal.pone.0124563.g003

Per residues binding energy decomposition

We evaluated the energetic contribution made by the sum of Lennard-Jones and Coulomb terms to the binding enthalpy for each residue of the protein binding cavity displaying differences between DnaK and Hsp70 (Table 5). For each residue the relative fraction of the total enthalpy was calculated. Residues that mainly contribute to the enthalpy of binding are Met404, Ser427 and Ala429. In particular, Ser427 is a key residue for the interaction of all the analysed ligands bound to the three proteins, both DnaK proteins and hHsp70. On the other hand, Met404 displays a relevant enthalpic contribution only for api88 bound to all proteins, for Drosocin bound to both DnaK proteins and for Pyrrhocoricin bound to E.coli DnaK. Moreover, Ala429 shows a large contribution for all Api88 and Pyrrhocoricin complexes, for hHsp70-Drosocin, for A.baumannii DnaK and for hHsp70-Oncocin.

Download:

Table 5. Enthalpic contribution to the binding energy.

https://doi.org/10.1371/journal.pone.0124563.t005

Evaluation of residue conservation

At last, the conservation of residues identified with the previous analysis was evaluated. In particular, residue conservation among the 65 DnaK bacterial sequences (Fig 4A) and among the 13 Hsp70 human sequences (Fig 4B) was analysed. Overall, residues belonging to the binding cavities of both proteins are quite conserved. In detail, positions 401–411 (β₁—L_1,2 – β₂) are highly conserved in bacteria, while in human proteins only positions 405, 409 and 411 are well conserved and 402, 403 are quite conserved. Considering residues of β₃—L_3,4 – β_4, positions 426–429, 431–433, 436, 438 and 440 are highly conserved in bacterial DnaK, while considering human Hsp70 only positions 426, 431 and 436 are conserved. Among residues belonging to L_{5,6 –} β₆ of DnaK, only positions 458 and 469 are not conserved, while Arg468 and Ile474 are quite conserved also in Hsp70 sequences. Concerning residues of HB, only Glu531, Asn537 and Asp540 (respectively Asn540 and Glu543 for Hsp70) are highly conserved.

Download:

Fig 4. WebLogo representation obtained from the alignment of the 65 DnaK (a) and of the 13 human Hsp70 (b) sequences.

Residues belonging to the binding site are underlined, residues identified with alanine scanning analysis are in red boxes, residues obtained by interaction analysis are in blue boxes and residues identified in the pharmacophoric analysis are in green boxes. Residues double underlined are identified hot spots.

https://doi.org/10.1371/journal.pone.0124563.g004

Analysis of structure-based pharmacophoric model

The structure-based pharmacophore model comprises the most important pharmacological features for inhibitors binding. Relying on the identified residues, a pharmacophore was obtained for each DnaK-peptide complex (Fig 5). The final pharmacophore, shown in Fig 6, includes the common characteristics of all the models achieved for the different peptides. Several hydrophobic groups (yellow spheres) can be found in the central part of the pharmacophore model. On the other hand, one positive and several negative chargeable groups (red and blue spheres, respectively) localized at both ends of the model. Moreover, along the whole pharmacophore, hydrogen bond donors and acceptors (red and green arrows, respectively) are distributed. This model is in agreement with the inhibitor model identified by Rüdiger et al. [9, 13].

Download:

Fig 5. Structure-based pharmacophore.

Pharmacophore built on residues important for the interaction between DnaK and inhibitor-peptides and not identified in human proteins. The yellow spheres represent the hydrophobic groups, red and green arrows identify, respectively, H-bonds acceptors and donors. Red and blue spheres are, respectively, negative and positive ionizable groups.

https://doi.org/10.1371/journal.pone.0124563.g005

Download:

Fig 6. Key residues identified in DnaK SBD.

Residue side-chains are in ball and stick, C_α atoms of residues previously exploited by known inhibitors are in light blue, C_α of the new identified hot spot residues are in orange. The Api88 peptide is represented in red.

https://doi.org/10.1371/journal.pone.0124563.g006

Concerning protein residues, these principally correspond to amino acids identified in previous analysis. In particular, Thr403 and Phe426 are involved only in hydrophobic contacts; Ser427 and Asp/Thr437 are donors or acceptors for hydrogen bonds, while Gln433 results only as an acceptor of hydrogen bond. Moreover, Met404, Thr409 and Ala429 are occupied both in contacts and hydrogen bonds as acceptors or donors.

Discussion

DnaK is the bacterial homolog of the human Hsp70 and it contributes to bacterial multi-drug resistance, which makes this protein a very interesting drug target. For this reason we propose a detailed study of DnaK binding cavity and its interaction capability in comparison to five recently discovered peptide inhibitors (Api88, drosocin, heliocin, oncocin and pyrrhocoricin) in comparison to the human Hsp70, though structure analysis and molecular dynamics simulations. These methods allow the analysis of complexes in static and dynamic situations. The analysis of known complexes returns residues stably involved in the interaction with inhibitor-peptides, conserved in bacterial DnaK proteins, but not in the human homolog, while more general analyses on the binding cavity residues allow to identify amino acids potentially useful for the design of high affinity and selective ligands.

Overall, residues of βSBD are quite conserved, while residues belonging to helices HA and HB are poorly conserved in bacterial DnaK sequences. Considering only the binding cavity residues, several non-conserved amino acids between bacterial DnaK and human Hsp70 were identified: Met404, Val425, Ser427, Ala429, Ala435, Ile438 and Gly468. These residues are non-conserved in human Hsp70 or are not included in its binding cavity. Our prediction concerning residues, Met404, Ser427, Ala429, Ala435 and Ile438, included in L_3,4 and L_1,2 cavity-forming loops, was experimentally confirmed in contacts with the bound peptide [9]. These residues are also involved in the specificity of peptide sequence recognition though van der Waals interactions [9]. Some of these residues were already discussed in the literature: Ile438 was identified among the residues involved in the hydrophobic pocket binding of DnaK [11], while the single mutant S427P displays an ATPase activity similar to the wild type DnaK, but with a weakened peptide binding affinity [42–44]. None of these residues occupy a key position in comparison to the alanine scanning, although Ser427, Ala429, Ile438 and Gly468 result highly conserved in bacterial DnaK proteins, suggesting a selective pressure on these positions. These data were obtained relying on a static conformation, which can explain the different results achieved with our analysis.

Among the identified residues, Met404, Ser427 and Ala429 (Fig 6) are involved in conserved hydrogen bonds and contacts with inhibitor-peptides and are prevalent in the enthalpic contribution to the binding energy. As confirmed by the pharmacophore model, Met404, Ala429 and Ser427 are anchoring points for hydrogen bonds and Met404, Ala429 for hydrophobic contacts, too. Residue Ser427 localized outward the cavity, while Met404 and Ala429 side-chains are closed up to the binding cavity. Looking at these residues, Ser427 (β₃) and Ala429 (L_3,4) are highly conserved in bacterial DnaK proteins, while in the human Hsp70 they are substituted by non-conserved residues. On the other hand, Met404 is quite conserved in the bacterial DnaK family and no conserved in the human family. Ala429 forms with Met404 the hydrophobic arch over the peptide backbone [43] and is considered particularly important for the peptide binding [44]. Experimental data confirms that mutations in these positions impair the capability of the proteins to bind peptides [43–45]. Moreover, positions 427 and 429 are included in the binding cavity of Hsp70 with a Thr and a Tyr, respectively.

The binding cavity analysis returned four other residues not involved in the binding of known peptides, which could furthermore be exploited for drug design. Two of these residues, Ile438 and Gly468 are highly conserved both in bacterial DnaK proteins and in the human Hsp70. However Val425 and Ala435 are quite conserved in bacterial, but not in the human protein family, although position 425 presents a conserved Thr residue in human proteins. These two differential conserved residues can be useful anchoring points for the design of selective peptides (Fig 6).

The SBD domain is partially involved in co-chaperones binding, mainly with proteins carrying the J-domain [45]. In particular, the C-terminal helical subdomain is involved in this interaction [46]. Residues identified in this study are not involved in interactions with J-domain proteins.

Our analysis allowed to identify three residues, Met404, Ser427 and Ala429 previously exploited by known inhibitors, and two new anchoring points, Val425 and Ala435, belonging to the binding cavity, probably involved in the inhibition of DnaK binding and selective for bacterial proteins compared to the human Hsp70. In particular, the three known residues localized in the binding cavity anchor the central part of the ligand, while the two new identified residues are on the binding cavity border and can be employed to recognize the external part of the ligand. Moreover, from the binding energy analysis, residues Val425 and Ala435 can be better exploited as anchoring points.

Rational drug design of compounds capable of binding these residues may increase the chance of delivering selective inhibitors. On the other hand, a virtual screening strategy, followed by the selection of ligands that bind critical residues, could also be applied for identifying inhibitors capable of selectively impairing the bacterial DnaK with high affinity.

Supporting Information

S1 Fig. Alignment of bacterial DnaK representative sequences.

https://doi.org/10.1371/journal.pone.0124563.s001

(PDF)

S2 Fig. Human Hsp70 isoforms alignment.

https://doi.org/10.1371/journal.pone.0124563.s002

(PDF)

Acknowledgments

We thank Dr. Marco Moscatelli of Institute of Biomedical Technologies (CNR-ITB) for providing technical support for molecular dynamics simulations and John Hatton of Institute of Biomedical Technologies (CNR-ITB) for proofreading the manuscript.

Author Contributions

Conceived and designed the experiments: FC IM LM. Performed the experiments: FC MF. Analyzed the data: FC MF. Contributed reagents/materials/analysis tools: IM LM. Wrote the paper: FC IM. Critically revised the manuscript: LM MF.

References

1. Hesterkamp T, Bukau B (1998) Role of the DnaK and HscA homologs of Hsp70 chaperones in protein folding in E.coli. The EMBO Journal 17: 4818–4828. pmid:9707441
- View Article
- PubMed/NCBI
- Google Scholar
2. Henderson B, Allan E, Coates ARM (2006) Stress Wars: the Direct Role of Host and bacterial Molecular Chaperones in Bacterial Infection. Infect. Immun. 74: 3693–3706. pmid:16790742
- View Article
- PubMed/NCBI
- Google Scholar
3. Tomoyasu T, Tabata A, Imaki H, Tsuruno K, Miyazaki A, Sonomoto K, et al. (2012) Role of Streptococcus intermedius DnaK chaperone system in stress tolerance and pathogenicity. Cell Stress and Chaperones 17: 41–55. pmid:21822788
- View Article
- PubMed/NCBI
- Google Scholar
4. Patury S, Miyata Y, Gestwicki JE (2009) Pharmacological Targeting of the Hsp70 Chaperone. Curr Top Med Chem. 9: 1337–1351. pmid:19860737
- View Article
- PubMed/NCBI
- Google Scholar
5. Knappe D, Zahn M, Sauer U, Schiffer G, Sträter N, Hoffmann R (2011) Rational Design of Oncocin Derivatives with Superior Protease Stabilities and Antibacterial Activities Based on the High-Resolution Structure of the Oncocin-DnaK Complex. ChemBioChem 12: 874–876. pmid:21387510
- View Article
- PubMed/NCBI
- Google Scholar
6. Yamaguchi Y, Tomoyasu T, Takaya A, Morioka M, Yamamoto T (2003) Effects of disruption of heat shock genes on susceptibility of Escherichia coli to fluoroquinolones. BMC Microbiology 3: 16. pmid:12911840
- View Article
- PubMed/NCBI
- Google Scholar
7. Singh VK, Utaida S, Jackson LS, Jayaswal RK, Wilkinson BJ, Chamberlain NR (2007) Role for dnaK locus in tolerance of multiple stresses in Staphylococcus aureus. Microbiology 153: 3162–3173. pmid:17768259
- View Article
- PubMed/NCBI
- Google Scholar
8. Evans CG, Chang L, Gestwicki JE (2010) Heat Shock Protein 70 (Hsp70) as an Emerging Drug Target. J. Med. Chem. 53: 4585–4602. pmid:20334364
- View Article
- PubMed/NCBI
- Google Scholar
9. Rüdiger S, Buchberger A, Bukau B (1997) Interaction of Hsp70 chaperones with substrates. Nature Structural Biology 4: 342–349. pmid:9145101
- View Article
- PubMed/NCBI
- Google Scholar
10. Genevaux P, Georgopoulos C, Kelley WL (2007) The Hsp70 chaperone machines of Escherichia coli: a paradigm for the repartition of chaperone functions. Molecular Microbiology 66: 840–857. pmid:17919282
- View Article
- PubMed/NCBI
- Google Scholar
11. Zahn M, Berthold N, Kieslich B, Knappe D, Hoffmann R, Sträter N (2013) Structural Studies on the Forward and Reverse Binding Modes of Peptides to the Chaperone DnaK. J. Mol. Biol. 425: 2463–2479. pmid:23562829
- View Article
- PubMed/NCBI
- Google Scholar
12. Bertelsen EB, Chang L, Gestwicki JE, Zuiderweg ERP (2009) Solution conformation of wild-type E. coli Hsp70 (DnaK) chaperone complexed with ADP and substrate. Proc Natl Acad Sci U SA 106: 8471–8476. pmid:19439666
- View Article
- PubMed/NCBI
- Google Scholar
13. Rüdiger S, Germeroth L, Schneider-Mergener J, Bukau B (1997) Substrate specificity of the DnaK chaperone determined by screening cellulose-bound peptide libraries. The EMBO Journal 16: 1501–1507. pmid:9130695
- View Article
- PubMed/NCBI
- Google Scholar
14. Azoulay I, Kucherenko N, Nachliel E, Gutman M, Azem A, Tsfadia Y (2013) Tracking the Interplay between Bound Peptide and the Lid Domain of DnaK, Using Molecular Dynamics. Int. J. Mol. Sci. 14: 12675–12695. pmid:23774839
- View Article
- PubMed/NCBI
- Google Scholar
15. Srinivasan SR, Gillies AT, Chang L, Thompsona AD, Gestwicki JE (2012) Molecular chaperones DnaK and DnaJ share predicted binding sites on most proteins in the E. coli proteome. Mol. BioSyst. 8: 2323–2333. pmid:22732719
- View Article
- PubMed/NCBI
- Google Scholar
16. Czihal P, Knappe D, Fritsche S, Zahn M, Berthold N, Piantavigna S, et al. (2012) Api88 Is a Novel Antibacterial Designer Peptide To Treat Systemic Infections with Multidrug-Resistant Gram-Negative Pathogens. ACS Chem. Biol. 7: 1281–1291. pmid:22594381
- View Article
- PubMed/NCBI
- Google Scholar
17. Hansen A, Schäfer I, Knappe D, Seibel P, Hoffmann R (2012) Intracellular Toxicity of Proline-Rich Antimicrobial Peptides Shuttled into Mammalian Cells by the Cell-Penetrating Peptide Penetratin. Antimicrob. Agents Chemother. 56: 5194–5201. pmid:22850523
- View Article
- PubMed/NCBI
- Google Scholar
18. Craik DJ, Fairlie DP, Liras S, Price D (2013) The Future of Peptide-based Drugs. Chemical Biology & Drug Design. 81: 136–147.
- View Article
- Google Scholar
19. Shan Y, Kim ET, Eastwood MP, Dror RO, Seeliger MA, Shaw DE (2011) How does a drug molecule find its target binding site? J. Am. Chem. Soc. 133: 9181–9183. pmid:21545110
- View Article
- PubMed/NCBI
- Google Scholar
20. UniProt website. Available: http://www.uniprot.org/
21. Sievers F, Wilm A, Dineen DG, Gibson TJ, Karplus K, Li W, et al. (2011) Fast, scalable generation of high-quality protein multiple sequence alignments using Clustal Omega. Molecular Systems Biology 7: 539 pmid:21988835
- View Article
- PubMed/NCBI
- Google Scholar
22. SIAS. Available: http://imed.med.ucm.es/Tools/sias.html
23. Crooks GE, Hon G, Chandonia JM, Brenner SE (2004) WebLogo: A sequence logo generator, Genome Research 14: 1188–1190. pmid:15173120
- View Article
- PubMed/NCBI
- Google Scholar
24. Bertelsen EB, Chang L, Gestwicki JE, Zuiderweg ER (2009) Solution conformation of wild-type E. coli Hsp70 (DnaK) chaperone complexed with ADP and substrate. Proc. Natl. Acad. Sci. USA 106: 8471–8476. pmid:19439666
- View Article
- PubMed/NCBI
- Google Scholar
25. Berman HM, Westbrook J, Feng Z, Gilliland G, Bhat TN, Weissig H, et al. (2000) The Protein Data Bank. Nucl. Acids Res. 28: 235–242. pmid:10592235
- View Article
- PubMed/NCBI
- Google Scholar
26. Arnold K, Bordoli L, Kopp J, Schwede T (2006) The SWISS-MODEL Workspace: A web-based environment for protein structure homology modelling. Bioinformatics 22: 195–201. pmid:16301204
- View Article
- PubMed/NCBI
- Google Scholar
27. Hess B, Kutzner C, Van der Spoel D, Lindahl E (2008) GROMACS 4: Algorithms for Highly Efficient, Load-Balanced, and Scalable Molecular Simulations. J. Chem. Theory Comput. 4: 435–447.
- View Article
- Google Scholar
28. van Gunsteren WF, Billeter SR, Eising AA, Hunenberger PH, Krüger P, Mark AE, et al. (1996) Biomolecular Simulation: The GROMOS96 manual and user guide. Zurich, Switzerland: Hochschulverlag AG an der ETH Zurich.
29. Hess B, Bekker H, Berendsen HJC, Fraaije JGEM (1997) LINCS. A linear constraint solver for molecular simulations. J Comp Chem. 18: 1463–1472.
- View Article
- Google Scholar
30. Essmann U, Perera L, Berkowitz ML, Darden T, Lee H, Pedersen LGA (1995) Smooth Particle Mesh Ewald Potential. J. Chem. Phys. 103: 8577–8593.
- View Article
- Google Scholar
31. Daura X, Gademann K, Jaun B, Seebach D, van Gunsteren WF, Mark AE (1999) Peptide folding: When simulation meets experiment, Angew. Chem. Int. Ed. Engl. 38: 236–240.
- View Article
- Google Scholar
32. Humphrey W, Dalke A, Schulten K (1996) VMD—Visual Molecular Dynamics. J. Molec. Graphics 14: 33–38.
- View Article
- Google Scholar
33. Chiappori F, Merelli I, Colombo G, Milanesi L, Morra G (2012) Molecular Mechanism of Allosteric Communication in Hsp70 Revealed by Molecular Dynamics Simulations. PLoS Comput Biol 8: e1002844. pmid:23300424
- View Article
- PubMed/NCBI
- Google Scholar
34. Bevan Lab website. Available: http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin/Scripts/plot_hbmap.pl
35. Chiappori F, Merelli I, Milanesi L, Marabotti A (2013) Static and dynamic interactions between GALK enzyme and known inhibitors: Guidelines to design new drugs for galactosemic patients. European Journal of Medicinal Chemistry 63: 423–434. pmid:23517731
- View Article
- PubMed/NCBI
- Google Scholar
36. Robetta. Available: http://robetta.bakerlab.org/
37. Kortemme T, Kim DE, Baker D (2004) Computational Alanine Scanning of Protein-Protein Interfaces. Sci. STKE 219: pl2.
- View Article
- Google Scholar
38. Kortemme T, Baker D (2002) A simple physical model for binding energy hot spots in protein-protein complexes. Proc Natl Acad Sci U S A. 99: 14116–21. pmid:12381794
- View Article
- PubMed/NCBI
- Google Scholar
39. Wolber G, Langer T (2005) LigandScout: 3-D pharmacophores derived from protein-bound ligands and their use as virtual screening filters, J. Chem. Inf. Model. 45: 160–169. pmid:15667141
- View Article
- PubMed/NCBI
- Google Scholar
40. Milner CM, Campbell RD (1990) Structure and expression of the three MHC-linked HSP70 genes. Immunogenetics 32: 242–51 pmid:1700760
- View Article
- PubMed/NCBI
- Google Scholar
41. Kortemme T, Baker D (2002) A simple physical model for binding energy hot spots in protein–protein complexes. Proc Natl Acad Sci U SA 99: 14116–14121. pmid:12381794
- View Article
- PubMed/NCBI
- Google Scholar
42. Montgomery DL, Morimoto RI, Gierasch LM (1999) Mutations in the Substrate Binding Domain of the Escherichia coli 70 kDa Molecular Chaperone, DnaK, which Alter Substrate Affinity or Interdomain Coupling. J. Mol. Biol. 286: 915–932. pmid:10024459
- View Article
- PubMed/NCBI
- Google Scholar
43. Rüdiger S, Mayer MP, Schneider-Mergener J, Bukau B (2000) Modulation of Substrate Specificity of the DnaK Chaperone by Alteration of a Hydrophobic Arch. J. Mol. Biol. 304: 245–251. pmid:11090270
- View Article
- PubMed/NCBI
- Google Scholar
44. Qi R, Sarbeng EB, Liu Q, Le KQ, Xu X, Xu H, et al. (2013) Allosteric opening of the polypeptide-binding site when an Hsp70 binds ATP. Nat Struct Mol Biol. 20: 900–907. pmid:23708608
- View Article
- PubMed/NCBI
- Google Scholar
45. Suh W, Burkholder W, Lu CZ, Zhao X, Gottesman ME, Gross CA (1998) Interaction of the Hsp70 molecular chaperone, DnaK, with its cochaperone DnaJ. Proc Natl Acad Sci USA. 95: 15223–15228. pmid:9860950
- View Article
- PubMed/NCBI
- Google Scholar
46. Gao X, Zhou C, Zhou Z, Wu M, Cao C, Hu H (2012) The C-terminal Helices of Heat Shock Protein 70 Are Essential for J-domain Binding and ATPase Activation. J Biol Chem. 287: 6044–6052. pmid:22219199
- View Article
- PubMed/NCBI
- Google Scholar

[ref1] 1. Hesterkamp T, Bukau B (1998) Role of the DnaK and HscA homologs of Hsp70 chaperones in protein folding in E.coli. The EMBO Journal 17: 4818–4828. pmid:9707441
View Article
PubMed/NCBI
Google Scholar

[2] View Article

[3] PubMed/NCBI

[4] Google Scholar

[ref2] 2. Henderson B, Allan E, Coates ARM (2006) Stress Wars: the Direct Role of Host and bacterial Molecular Chaperones in Bacterial Infection. Infect. Immun. 74: 3693–3706. pmid:16790742
View Article
PubMed/NCBI
Google Scholar

[6] View Article

[7] PubMed/NCBI

[8] Google Scholar

[ref3] 3. Tomoyasu T, Tabata A, Imaki H, Tsuruno K, Miyazaki A, Sonomoto K, et al. (2012) Role of Streptococcus intermedius DnaK chaperone system in stress tolerance and pathogenicity. Cell Stress and Chaperones 17: 41–55. pmid:21822788
View Article
PubMed/NCBI
Google Scholar

[10] View Article

[11] PubMed/NCBI

[12] Google Scholar

[ref4] 4. Patury S, Miyata Y, Gestwicki JE (2009) Pharmacological Targeting of the Hsp70 Chaperone. Curr Top Med Chem. 9: 1337–1351. pmid:19860737
View Article
PubMed/NCBI
Google Scholar

[14] View Article

[15] PubMed/NCBI

[16] Google Scholar

[ref5] 5. Knappe D, Zahn M, Sauer U, Schiffer G, Sträter N, Hoffmann R (2011) Rational Design of Oncocin Derivatives with Superior Protease Stabilities and Antibacterial Activities Based on the High-Resolution Structure of the Oncocin-DnaK Complex. ChemBioChem 12: 874–876. pmid:21387510
View Article
PubMed/NCBI
Google Scholar

[18] View Article

[19] PubMed/NCBI

[20] Google Scholar

[ref6] 6. Yamaguchi Y, Tomoyasu T, Takaya A, Morioka M, Yamamoto T (2003) Effects of disruption of heat shock genes on susceptibility of Escherichia coli to fluoroquinolones. BMC Microbiology 3: 16. pmid:12911840
View Article
PubMed/NCBI
Google Scholar

[22] View Article

[23] PubMed/NCBI

[24] Google Scholar

[ref7] 7. Singh VK, Utaida S, Jackson LS, Jayaswal RK, Wilkinson BJ, Chamberlain NR (2007) Role for dnaK locus in tolerance of multiple stresses in Staphylococcus aureus. Microbiology 153: 3162–3173. pmid:17768259
View Article
PubMed/NCBI
Google Scholar

[26] View Article

[27] PubMed/NCBI

[28] Google Scholar

[ref8] 8. Evans CG, Chang L, Gestwicki JE (2010) Heat Shock Protein 70 (Hsp70) as an Emerging Drug Target. J. Med. Chem. 53: 4585–4602. pmid:20334364
View Article
PubMed/NCBI
Google Scholar

[30] View Article

[31] PubMed/NCBI

[32] Google Scholar

[ref9] 9. Rüdiger S, Buchberger A, Bukau B (1997) Interaction of Hsp70 chaperones with substrates. Nature Structural Biology 4: 342–349. pmid:9145101
View Article
PubMed/NCBI
Google Scholar

[34] View Article

[35] PubMed/NCBI

[36] Google Scholar

[ref10] 10. Genevaux P, Georgopoulos C, Kelley WL (2007) The Hsp70 chaperone machines of Escherichia coli: a paradigm for the repartition of chaperone functions. Molecular Microbiology 66: 840–857. pmid:17919282
View Article
PubMed/NCBI
Google Scholar

[38] View Article

[39] PubMed/NCBI

[40] Google Scholar

[ref11] 11. Zahn M, Berthold N, Kieslich B, Knappe D, Hoffmann R, Sträter N (2013) Structural Studies on the Forward and Reverse Binding Modes of Peptides to the Chaperone DnaK. J. Mol. Biol. 425: 2463–2479. pmid:23562829
View Article
PubMed/NCBI
Google Scholar

[42] View Article

[43] PubMed/NCBI

[44] Google Scholar

[ref12] 12. Bertelsen EB, Chang L, Gestwicki JE, Zuiderweg ERP (2009) Solution conformation of wild-type E. coli Hsp70 (DnaK) chaperone complexed with ADP and substrate. Proc Natl Acad Sci U SA 106: 8471–8476. pmid:19439666
View Article
PubMed/NCBI
Google Scholar

[46] View Article

[47] PubMed/NCBI

[48] Google Scholar

[ref13] 13. Rüdiger S, Germeroth L, Schneider-Mergener J, Bukau B (1997) Substrate specificity of the DnaK chaperone determined by screening cellulose-bound peptide libraries. The EMBO Journal 16: 1501–1507. pmid:9130695
View Article
PubMed/NCBI
Google Scholar

[50] View Article

[51] PubMed/NCBI

[52] Google Scholar

[ref14] 14. Azoulay I, Kucherenko N, Nachliel E, Gutman M, Azem A, Tsfadia Y (2013) Tracking the Interplay between Bound Peptide and the Lid Domain of DnaK, Using Molecular Dynamics. Int. J. Mol. Sci. 14: 12675–12695. pmid:23774839
View Article
PubMed/NCBI
Google Scholar

[54] View Article

[55] PubMed/NCBI

[56] Google Scholar

[ref15] 15. Srinivasan SR, Gillies AT, Chang L, Thompsona AD, Gestwicki JE (2012) Molecular chaperones DnaK and DnaJ share predicted binding sites on most proteins in the E. coli proteome. Mol. BioSyst. 8: 2323–2333. pmid:22732719
View Article
PubMed/NCBI
Google Scholar

[58] View Article

[59] PubMed/NCBI

[60] Google Scholar

[ref16] 16. Czihal P, Knappe D, Fritsche S, Zahn M, Berthold N, Piantavigna S, et al. (2012) Api88 Is a Novel Antibacterial Designer Peptide To Treat Systemic Infections with Multidrug-Resistant Gram-Negative Pathogens. ACS Chem. Biol. 7: 1281–1291. pmid:22594381
View Article
PubMed/NCBI
Google Scholar

[62] View Article

[63] PubMed/NCBI

[64] Google Scholar

[ref17] 17. Hansen A, Schäfer I, Knappe D, Seibel P, Hoffmann R (2012) Intracellular Toxicity of Proline-Rich Antimicrobial Peptides Shuttled into Mammalian Cells by the Cell-Penetrating Peptide Penetratin. Antimicrob. Agents Chemother. 56: 5194–5201. pmid:22850523
View Article
PubMed/NCBI
Google Scholar

[66] View Article

[67] PubMed/NCBI

[68] Google Scholar

[ref18] 18. Craik DJ, Fairlie DP, Liras S, Price D (2013) The Future of Peptide-based Drugs. Chemical Biology & Drug Design. 81: 136–147.
View Article
Google Scholar

[70] View Article

[71] Google Scholar

[ref19] 19. Shan Y, Kim ET, Eastwood MP, Dror RO, Seeliger MA, Shaw DE (2011) How does a drug molecule find its target binding site? J. Am. Chem. Soc. 133: 9181–9183. pmid:21545110
View Article
PubMed/NCBI
Google Scholar

[73] View Article

[74] PubMed/NCBI

[75] Google Scholar

[ref20] 20. UniProt website. Available: http://www.uniprot.org/

[ref21] 21. Sievers F, Wilm A, Dineen DG, Gibson TJ, Karplus K, Li W, et al. (2011) Fast, scalable generation of high-quality protein multiple sequence alignments using Clustal Omega. Molecular Systems Biology 7: 539 pmid:21988835
View Article
PubMed/NCBI
Google Scholar

[78] View Article

[79] PubMed/NCBI

[80] Google Scholar

[ref22] 22. SIAS. Available: http://imed.med.ucm.es/Tools/sias.html

[ref23] 23. Crooks GE, Hon G, Chandonia JM, Brenner SE (2004) WebLogo: A sequence logo generator, Genome Research 14: 1188–1190. pmid:15173120
View Article
PubMed/NCBI
Google Scholar

[83] View Article

[84] PubMed/NCBI

[85] Google Scholar

[ref24] 24. Bertelsen EB, Chang L, Gestwicki JE, Zuiderweg ER (2009) Solution conformation of wild-type E. coli Hsp70 (DnaK) chaperone complexed with ADP and substrate. Proc. Natl. Acad. Sci. USA 106: 8471–8476. pmid:19439666
View Article
PubMed/NCBI
Google Scholar

[87] View Article

[88] PubMed/NCBI

[89] Google Scholar

[ref25] 25. Berman HM, Westbrook J, Feng Z, Gilliland G, Bhat TN, Weissig H, et al. (2000) The Protein Data Bank. Nucl. Acids Res. 28: 235–242. pmid:10592235
View Article
PubMed/NCBI
Google Scholar

[91] View Article

[92] PubMed/NCBI

[93] Google Scholar

[ref26] 26. Arnold K, Bordoli L, Kopp J, Schwede T (2006) The SWISS-MODEL Workspace: A web-based environment for protein structure homology modelling. Bioinformatics 22: 195–201. pmid:16301204
View Article
PubMed/NCBI
Google Scholar

[95] View Article

[96] PubMed/NCBI

[97] Google Scholar

[ref27] 27. Hess B, Kutzner C, Van der Spoel D, Lindahl E (2008) GROMACS 4: Algorithms for Highly Efficient, Load-Balanced, and Scalable Molecular Simulations. J. Chem. Theory Comput. 4: 435–447.
View Article
Google Scholar

[99] View Article

[100] Google Scholar

[ref28] 28. van Gunsteren WF, Billeter SR, Eising AA, Hunenberger PH, Krüger P, Mark AE, et al. (1996) Biomolecular Simulation: The GROMOS96 manual and user guide. Zurich, Switzerland: Hochschulverlag AG an der ETH Zurich.

[ref29] 29. Hess B, Bekker H, Berendsen HJC, Fraaije JGEM (1997) LINCS. A linear constraint solver for molecular simulations. J Comp Chem. 18: 1463–1472.
View Article
Google Scholar

[103] View Article

[104] Google Scholar

[ref30] 30. Essmann U, Perera L, Berkowitz ML, Darden T, Lee H, Pedersen LGA (1995) Smooth Particle Mesh Ewald Potential. J. Chem. Phys. 103: 8577–8593.
View Article
Google Scholar

[106] View Article

[107] Google Scholar

[ref31] 31. Daura X, Gademann K, Jaun B, Seebach D, van Gunsteren WF, Mark AE (1999) Peptide folding: When simulation meets experiment, Angew. Chem. Int. Ed. Engl. 38: 236–240.
View Article
Google Scholar

[109] View Article

[110] Google Scholar

[ref32] 32. Humphrey W, Dalke A, Schulten K (1996) VMD—Visual Molecular Dynamics. J. Molec. Graphics 14: 33–38.
View Article
Google Scholar

[112] View Article

[113] Google Scholar

[ref33] 33. Chiappori F, Merelli I, Colombo G, Milanesi L, Morra G (2012) Molecular Mechanism of Allosteric Communication in Hsp70 Revealed by Molecular Dynamics Simulations. PLoS Comput Biol 8: e1002844. pmid:23300424
View Article
PubMed/NCBI
Google Scholar

[115] View Article

[116] PubMed/NCBI

[117] Google Scholar

[ref34] 34. Bevan Lab website. Available: http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin/Scripts/plot_hbmap.pl

[ref35] 35. Chiappori F, Merelli I, Milanesi L, Marabotti A (2013) Static and dynamic interactions between GALK enzyme and known inhibitors: Guidelines to design new drugs for galactosemic patients. European Journal of Medicinal Chemistry 63: 423–434. pmid:23517731
View Article
PubMed/NCBI
Google Scholar

[120] View Article

[121] PubMed/NCBI

[122] Google Scholar

[ref36] 36. Robetta. Available: http://robetta.bakerlab.org/

[ref37] 37. Kortemme T, Kim DE, Baker D (2004) Computational Alanine Scanning of Protein-Protein Interfaces. Sci. STKE 219: pl2.
View Article
Google Scholar

[125] View Article

[126] Google Scholar

[ref38] 38. Kortemme T, Baker D (2002) A simple physical model for binding energy hot spots in protein-protein complexes. Proc Natl Acad Sci U S A. 99: 14116–21. pmid:12381794
View Article
PubMed/NCBI
Google Scholar

[128] View Article

[129] PubMed/NCBI

[130] Google Scholar

[ref39] 39. Wolber G, Langer T (2005) LigandScout: 3-D pharmacophores derived from protein-bound ligands and their use as virtual screening filters, J. Chem. Inf. Model. 45: 160–169. pmid:15667141
View Article
PubMed/NCBI
Google Scholar

[132] View Article

[133] PubMed/NCBI

[134] Google Scholar

[ref40] 40. Milner CM, Campbell RD (1990) Structure and expression of the three MHC-linked HSP70 genes. Immunogenetics 32: 242–51 pmid:1700760
View Article
PubMed/NCBI
Google Scholar

[136] View Article

[137] PubMed/NCBI

[138] Google Scholar

[ref41] 41. Kortemme T, Baker D (2002) A simple physical model for binding energy hot spots in protein–protein complexes. Proc Natl Acad Sci U SA 99: 14116–14121. pmid:12381794
View Article
PubMed/NCBI
Google Scholar

[140] View Article

[141] PubMed/NCBI

[142] Google Scholar

[ref42] 42. Montgomery DL, Morimoto RI, Gierasch LM (1999) Mutations in the Substrate Binding Domain of the Escherichia coli 70 kDa Molecular Chaperone, DnaK, which Alter Substrate Affinity or Interdomain Coupling. J. Mol. Biol. 286: 915–932. pmid:10024459
View Article
PubMed/NCBI
Google Scholar

[144] View Article

[145] PubMed/NCBI

[146] Google Scholar

[ref43] 43. Rüdiger S, Mayer MP, Schneider-Mergener J, Bukau B (2000) Modulation of Substrate Specificity of the DnaK Chaperone by Alteration of a Hydrophobic Arch. J. Mol. Biol. 304: 245–251. pmid:11090270
View Article
PubMed/NCBI
Google Scholar

[148] View Article

[149] PubMed/NCBI

[150] Google Scholar

[ref44] 44. Qi R, Sarbeng EB, Liu Q, Le KQ, Xu X, Xu H, et al. (2013) Allosteric opening of the polypeptide-binding site when an Hsp70 binds ATP. Nat Struct Mol Biol. 20: 900–907. pmid:23708608
View Article
PubMed/NCBI
Google Scholar

[152] View Article

[153] PubMed/NCBI

[154] Google Scholar

[ref45] 45. Suh W, Burkholder W, Lu CZ, Zhao X, Gottesman ME, Gross CA (1998) Interaction of the Hsp70 molecular chaperone, DnaK, with its cochaperone DnaJ. Proc Natl Acad Sci USA. 95: 15223–15228. pmid:9860950
View Article
PubMed/NCBI
Google Scholar

[156] View Article

[157] PubMed/NCBI

[158] Google Scholar

[ref46] 46. Gao X, Zhou C, Zhou Z, Wu M, Cao C, Hu H (2012) The C-terminal Helices of Heat Shock Protein 70 Are Essential for J-domain Binding and ATPase Activation. J Biol Chem. 287: 6044–6052. pmid:22219199
View Article
PubMed/NCBI
Google Scholar

[160] View Article

[161] PubMed/NCBI

[162] Google Scholar

Figures

Abstract

Introduction

Materials and Methods

Sequence analysis

Structure

Molecular dynamics

Molecular dynamics analysis

Results

SBD residue conservation

Identification of interaction key residues

Definition of the binding cavity

Computational alanine scanning

Per residues binding energy decomposition

Evaluation of residue conservation

Analysis of structure-based pharmacophoric model

Discussion

Supporting Information

S1 Fig. Alignment of bacterial DnaK representative sequences.

S2 Fig. Human Hsp70 isoforms alignment.

Acknowledgments

Author Contributions

References