Skip to main page content
U.S. flag

An official website of the United States government

Dot gov

The .gov means it’s official.
Federal government websites often end in .gov or .mil. Before sharing sensitive information, make sure you’re on a federal government site.

Https

The site is secure.
The https:// ensures that you are connecting to the official website and that any information you provide is encrypted and transmitted securely.

Access keys NCBI Homepage MyNCBI Homepage Main Content Main Navigation
. 2021 Nov 4;17(11):e1009567.
doi: 10.1371/journal.pcbi.1009567. eCollection 2021 Nov.

Physics-based modeling provides predictive understanding of selectively promiscuous substrate binding by Hsp70 chaperones

Affiliations

Physics-based modeling provides predictive understanding of selectively promiscuous substrate binding by Hsp70 chaperones

Erik B Nordquist et al. PLoS Comput Biol. .

Abstract

To help cells cope with protein misfolding and aggregation, Hsp70 molecular chaperones selectively bind a variety of sequences ("selective promiscuity"). Statistical analyses from substrate-derived peptide arrays reveal that DnaK, the E. coli Hsp70, binds to sequences containing three to five branched hydrophobic residues, although otherwise the specific amino acids can vary considerably. Several high-resolution structures of the substrate -binding domain (SBD) of DnaK bound to peptides reveal a highly conserved configuration of the bound substrate and further suggest that the substrate-binding cleft consists of five largely independent sites for interaction with five consecutive substrate residues. Importantly, both substrate backbone orientations (N- to C- and C- to N-) allow essentially the same backbone hydrogen-bonding and side-chain interactions with the chaperone. In order to rationalize these observations, we performed atomistic molecular dynamics simulations to sample the interactions of all 20 amino acid side chains in each of the five sites of the chaperone in the context of the conserved substrate backbone configurations. The resulting interaction energetics provide the basis set for deriving a predictive model that we call Paladin (Physics-based model of DnaK-Substrate Binding). Trained using available peptide array data, Paladin can distinguish binders and nonbinders of DnaK with accuracy comparable to existing predictors and further predicts the detailed configuration of the bound sequence. Tested using existing DnaK-peptide structures, Paladin correctly predicted the binding register in 10 out of 13 substrate sequences that bind in the N- to C- orientation, and the binding orientation in 16 out of 22 sequences. The physical basis of the Paladin model provides insight into the origins of how Hsp70s bind substrates with a balance of selectivity and promiscuity. The approach described here can be extended to other Hsp70s where extensive peptide array data is not available.

PubMed Disclaimer

Conflict of interest statement

The authors have declared that no competing interests exist.

Figures

Fig 1
Fig 1. Structural overview of DnaK in substrate-bound and -unbound states.
(A) Structure of DnaK in an unassociated (ADP-bound, substrate bound) conformation (PDB ID: 2KHO [9]). The nucleotide-binding domain (NBD) is shown in blue. The substrate-binding domain (SBD) has two subdomains: the α-helical lid is shown in red, and the βSBD is shown in teal. The substrate from the canonical NR peptide (PDB ID: 1DKZ [10]) is overlaid in the binding cleft of the βSBD in magenta. (B) Structure of DnaK in an associated (ATP-bound, no substrate) conformation (PDB ID: 4B9Q [11]). The ATP molecule is shown based on PDB ID: 3AY9 [12].
Fig 2
Fig 2. Conserved βSBD-substrate binding conformation.
(A) Overlay of 18 bound substrate structures (S1 Table), aligned using the backbone atoms of the βSBD. The SBDs are show in gray cartoons, and the substrates in magenta backbone traces. Only the forward orientation substrates are included. (B) Overlay of two forward and reverse orientation peptide substrates (forward: NRLLLTG, PDB: 1DKZ [10]; reverse: NRLILTG, PDB: 4EZY [19]; underlines indicate residues at Site 0). The forward orientation is shown in thick bonds and the reverse in thin bonds. The backbone hydrogen bonds between βSBD and the substrate are shown in magenta dashed lines. Note that backbone hydrogen bonds involving M404, S427 and A429 completely overlap in two orientations. Side chains are shown only for the central L at Site 0. For more viewing angles of panel B, see S1 Movie.
Fig 3
Fig 3. Substrate-interacting sites of DnaK.
The center image shows a NRLLLTG substrate (magenta sticks) bound to βSBD (gray cartoon) (PDB: 1DKX). Each zoom-in shows the surface of one of the five binding sites on DnaK, colored by atom type: O (red), N (blue), C/H (white), S (yellow). The side-chain conformations shown were selected to clearly illustrate the extent of each site. These additional side-chain and substrate backbone conformations are taken from: (Site -2) 4JWD, (Site -1) 1DKX, (Site 0) 1DKX (for more viewing angles, see S2 Movie), (Site +1) 1DKZ, 4F00 and (Site +2) 4JWI, 4EZN, 4EZO. The binding site surface was generated using DnaK residues within 5 Å of the substrate side chain(s) at each site, with some atoms excluded for clarity.
Fig 4
Fig 4. Individual, unscaled physical interaction terms for all possible substrate amino acid side chains occupying the site 0 in the SBD.
Error bars are calculated as standard error of the mean. The residues are ordered by Wimley-White interfacial hydrophobicity scale to facilitate easy reading [37]. Error bars report standard error of the mean from the MD simulations.
Fig 5
Fig 5. Receiver operating characteristic curves of various predictors.
The raw score was calculated using the energy terms with uniform weights. (A) Discrimination of both strong binders and binders from nonbinders, and (B) Discrimination of strong binders from nonbinders. All solid lines were derived from the training set and dotted lines from the validation set. The red dashed diagonal line represents a random prediction.
Fig 6
Fig 6. Registry predictions for forward-binding substrates.
For each substrate, the Paladin scores are shown for all possible binding registries that allow all five sites on DnaK to be occupied. Both backbone orientations are considered (forward: blue; reverse: orange). The registries as observed in the PDB structures are shown with blue stars. Note that the NRLLLTG peptide has been crystallized in two registries (PDB: 1DKZ, 4EZW).
Fig 7
Fig 7. Prediction of substrate binding orientation by Paladin.
The difference between the lowest forward and reverse scores will be negative when the for5tward orientation is favored, and vice versa. Substrates that actually bind in the forward orientation are labeled on the left of the red dividing line, and those that actually bind reverse are on the right.

Similar articles

Cited by

References

    1. Hartl FU, Bracher A, Hayer-Hartl M. Molecular chaperones in protein folding and proteostasis. Nature. 2011. Jul;475(7356):324–32. doi: 10.1038/nature10317 - DOI - PubMed
    1. Kim YE, Hipp MS, Bracher A, Hayer-Hartl M, Ulrich Hartl F. Molecular Chaperone Functions in Protein Folding and Proteostasis. Annu Rev Biochem. 2013. Jun 2;82(1):323–55. doi: 10.1146/annurev-biochem-060208-092442 - DOI - PubMed
    1. Morimoto RI, Cuervo AM. Proteostasis and the Aging Proteome in Health and Disease. J Gerontol A Biol Sci Med Sci. 2014. Jun 1;69(Suppl 1):S33–8. doi: 10.1093/gerona/glu049 - DOI - PMC - PubMed
    1. Clerico EM, Meng W, Pozhidaeva A, Bhasne K, Petridis C, Gierasch LM. Hsp70 molecular chaperones: multifunctional allosteric holding and unfolding machines. Biochem J. 2019. Jun 14;476(11):1653–77. doi: 10.1042/BCJ20170380 - DOI - PMC - PubMed
    1. Mayer MP. Intra-molecular pathways of allosteric control in Hsp70s. Philos Trans R Soc B Biol Sci. 2018. Jun 19;373(1749):20170183. doi: 10.1098/rstb.2017.0183 - DOI - PMC - PubMed

Publication types

MeSH terms

Substances

LinkOut - more resources

-