Introduction
Hsp70 proteins are preeminent among molecular chaperones in that their actions also feed into Hsp60, Hsp90 and Hsp100 chaperone systems and into protein degradation systems (Hartl et al., Reference Hartl, Bracher and Hayer-Hartl2011). They participate in diverse cellular processes; going beyond namesake stress responses, they play crucial roles in normal cells for protein folding, disassembly, degradation and membrane translocation, and they are protective against neurodegenerative diseases (Ciechanover and Kwon, Reference Ciechanover and Kwon2017) and complicit in cancers (Murphy, Reference Murphy2013). Hsp70s are found in all forms of life, excepting certain archaea, and in all ATP-containing cellular compartments of eukaryotes. Their sequences are highly conserved (>40% pairwise amino-acid identity), and especially so within three subfamilies corresponding to the eukaryotic cytosol, to the endoplasmic reticulum (ER), and to mitochondria, chloroplasts and prokaryotes. Although many Hsp70s are stress-induced, others are expressed constitutively from essential genes (Daugaard et al., Reference Daugaard, Rohde and Jäättelä2007).
Hsp70 proteins act in ATP-dependent cycles of binding and release of client substrates, typically exposed hydrophobic polypeptide segments. ATP binding and hydrolysis occurs in the nucleotide-binding domain (NBD), which controls the binding and release of polypeptides in the substrate-binding domain (SBD) (Zuiderweg et al., Reference Zuiderweg, Bertelsen, Rousaki, Mayer, Gestwicki and Ahmad2013; Mayer and Kityk, Reference Mayer and Kityk2015). The binding functions of NBD and SBD are separable, but Hsp70 chaperone activity requires direct, albeit transient, allosteric interactions between these sites as linked together (Hartl et al., Reference Hartl, Bracher and Hayer-Hartl2011; Zuiderweg et al., Reference Zuiderweg, Bertelsen, Rousaki, Mayer, Gestwicki and Ahmad2013; Mayer and Kityk, Reference Mayer and Kityk2015). ATP binding to NBD dramatically decreases SBD affinity for client substrates. Reciprocally, substrate binding stimulates ATP hydrolysis, whereupon substrates are retained tightly bound. Both on and off rates for substrate binding are accelerated in the presence of ATP relative to that with ADP-bound or nucleotide-free states. Hsp40s further stimulate ATP hydrolysis by Hsp70s and help to target them to substrates, and Hsp110s and other nucleotide exchange factors (NEFs) facilitate the release of ADP and rebinding of ATP; nevertheless, the Hsp70 chaperone cycle can proceed in vitro without these cofactors. The current picture of Hsp70 function is consistent with initial suggestions of Pelham (Reference Pelham1986) and Rothman (Reference Rothman1989) that molecular chaperones bind to aggregation-prone surfaces induced by stress and employ the energy of ATP hydrolysis for staged release and folding.
Crystal structures of individual NBD and SBD domains of Hsp70s provide a framework for biochemical understanding of Hsp70 chaperone activity. The prototype NBD structure is that from bovine Hsc70 (bHsc70) (Flaherty et al., Reference Flaherty, DeLuca-Flaherty and McKay1990). It comprises four subdomains (IA, IB, IIA and IIB) built up from two structurally similar lobes (I and II). Adenosine nucleotides bind at the interface between the lobes, making contacts with all four subdomains. The prototype SBD structure is that from Hsp70 DnaK of Escherichia coli (Zhu et al., Reference Zhu, Zhao, Burkholder, Gragerov, Ogata, Gottesman and Hendrickson1996). A substrate peptide is bound in an extended conformation through a channel defined by loops from SBDβ and covered by the SBDα subdomain. Findings from these prototypical structures were extended in numerous biochemical and biophysical studies as reviewed (Hartl et al., Reference Hartl, Bracher and Hayer-Hartl2011; Zuiderweg et al., Reference Zuiderweg, Bertelsen, Rousaki, Mayer, Gestwicki and Ahmad2013; Mayer and Kityk, Reference Mayer and Kityk2015).
Interactions between NBD and SBD are clearly essential for allosteric communication between the nucleotide and peptide-binding sites in an Hsp70 chaperone; however, the contacts are labile and their capture has proved elusive. For example, early efforts to crystallize full-length Hsp70s with ATP led instead to the structure of an NBD–ADP complex after ATP hydrolysis and incidental proteolysis (Sriram et al., Reference Sriram, Osipiuk, Freeman, Morimoto and Joachimiak1997). NBD and SBD are flexibly linked when in ADP or nucleotide-free states (Bertelsen et al., Reference Bertelsen, Chang, Gestwicki and Zuiderweg2009), and crystal structures that are obtained often have the DLLLLD-like NBD–SBD linker segment engaged adventitiously with a lattice mate (Chang et al., Reference Chang, Sun, Wang and Hsiao2008; Adell et al., Reference Adell, Calisto, Fita and Martinelli2018). ATP hydrolysis is typically too facile for the capture of stable Hsp70-ATP complexes, although NBD–linker–SBD interactions are evident when ATP is bound to an Hsp70 (Zhuravleva and Gierasch, Reference Zhuravleva and Gierasch2011), or even simply to an NBD–linker construct (Vogel et al., Reference Vogel, Mayer and Bukau2006; Swain et al., Reference Swain, Dinler, Sivendran, Montgomery, Stotz and Gierasch2007). Fortunately, yeast Hsp110 Sse1 provided us with a plausible Hsp70-ATP model (Liu and Hendrickson, Reference Liu and Hendrickson2007). Hsp110s bind but do not hydrolyze ATP, yet their sequences include unmistakable NBD domains and remote resemblances in SBD domains. The structure of Sse1-ATP showed interfaces between redisposed domains that are extensive, intimately complementary, and engaging of residues strikingly conserved in classic Hsp70 sequences. These features suggested that the Sse1-ATP interfaces might be evolutionary vestiges of functional Hsp70 interfaces. In confirmation, in vivo mutational tests of the inferred interfaces in yeast Hsp70 Ssa1 and E. coli Hsp70 DnaK produced severe phenotypes in each at 8 of 9 tested interfacial contact sites including 10 of 13 specific mutations (Liu and Hendrickson, Reference Liu and Hendrickson2007).
The Hsp110 structure inspired the design of DnaK constructs that could mimic the conformation found in the Sse1 structure and, using the hydrolysis-impaired T199A mutant (Barthel et al., Reference Barthel, Zhang and Walker2001), confirmatory structures were obtained for ATP complexes (Kityk et al., Reference Kityk, Kopp, Sinning and Mayer2012; Qi et al., Reference Qi, Sarbeng, Liu, Le, Xu, Xu, Yang, Wong, Vorvis, Hendrickson, Zhou and Liu2013). Subsequently, this conformational state of Hsp70-ATP was corroborated in structures of human Hsp70-8 (BiP) (Yang et al., Reference Yang, Nune, Zong, Zhou and Liu2015) and yeast SsaB (Gumiero et al., Reference Gumiero, Conz, Gesé, Zhang, Weyer, Lapouge, Kappes, Plehwe, Schermann, Fitzke, Wölfle, Fischer and Sinning2016). Further biochemical tests on selected interface DnaK mutants (Wang et al., Reference Wang, Liu, Liu and Hendrickson2020) led us to speculate that this Hsp110-like conformation must be restrained against ATP hydrolysis and essentially devoid of ATP binding, which in turn begged the questions of how the rebinding of substrate peptides and ATP hydrolysis might occur and how restraints against hydrolysis are effected. Additional structural analyses have shed light on these questions (Wang et al., Reference Wang, Liu, Liu and Hendrickson2020; Wang and Hendrickson, Reference Wang and Hendrickson2020a; Wang and Hendrickson, Reference Wang and Hendrickson2020b).
There is a rich literature from previous theoretical investigations of allosteric interactions in proteins (Monod et al., Reference Monod, Wyman and Changeux1965; Koshland et al., Reference Koshland, Némethy and Filmer1966; Cui and Karplus, Reference Cui and Karplus2008; Motlagh et al., Reference Motlagh, Wrabl, Li and Hilser2014; Cuendet et al., Reference Cuendet, Weinstein and LeVine2016; Thirumalai et al., Reference Thirumalai, Hyeon, Zhuravlev and Lorimer2019). Notably, the influential allosteric theories of Monod, Wyman and Changeux (MWC) (Monod et al., Reference Monod, Wyman and Changeux1965) and of Koshland, Némethy and Filmer (KMF) (Koshland et al., Reference Koshland, Némethy and Filmer1966) have treated cooperativity between similar binding sites in symmetric oligomers such as hemoglobin. The particular treatments of MWC and KMF are not directly applicable to Hsp70s, however; as these are monomeric proteins, predominantly, with distinct binding domains for altogether different ligands. What does apply more broadly is the concept that conformational equilibria between alternative states can govern allosteric regulation, and this has been pursued productively (Cui and Karplus, Reference Cui and Karplus2008; Motlagh et al., Reference Motlagh, Wrabl, Li and Hilser2014; Cuendet et al., Reference Cuendet, Weinstein and LeVine2016; Thirumalai et al., Reference Thirumalai, Hyeon, Zhuravlev and Lorimer2019). Nevertheless, we are not aware of quantitative treatments in the MWC mode for allostery in Hsp70 systems.
In an attempt to understand the distinctive Hsp70 mechanisms for allosteric control, we have developed a theoretical model for equilibria among conformational states in Hsp70 chaperones. This theory explains observations on ATP hydrolysis and polypeptide binding from wild-type (WT) and mutant variant DnaKs by postulating that ATP-bound Hsp70s equilibrate between states with distinct characteristics for the binding of substrate peptides and for the hydrolysis of ATP. Our biochemical results are fitted quantitatively by this allosteric theory, and the postulated but previously uncharacterized stimulating state has now been confirmed by crystal structures (Wang et al., Reference Wang, Liu, Liu and Hendrickson2020).
Theory
We assume that an Hsp70 chaperone protein exists in an equilibrium of states. Its NBD may bind ATP, ADP or be nucleotide-free (which we denote as Apo), and its SBD may bind segments of substrate polypeptides. The population of Hsp70 molecules that are complexed with ATP are in equilibrium between a restraining state, which binds substrate peptides poorly at best and only hydrolyzes ATP at a low basal rate, and a stimulating state, which binds substrate peptide well and hydrolyzes ATP at a substantially elevated rate. After ATP hydrolysis to the ADP state, peptide substrates are retained with high affinity in a state without allosteric coupling between the nucleotide and peptide binding sites. ADP may dissociate to yield the Apo state, remaining allosterically uncoupled and retaining the substrate peptide if present. Hsp70-Apo may rebind ATP to reinitiate the chaperone cycle. Additional intermediate states may also exist.
Overview of equilibrating states and analytic approach
The allosteric interactions between ATP in the NBD domain and a client peptide in the SBD domain can be followed biochemically by measuring peptide binding in the presence ATP and by measuring ATP hydrolysis in the presence of a client peptide. Such biochemical measurements on Sse1-inspired mutants provoked us to contemplate the theoretical basis for Hsp70 allostery (Wang et al., Reference Wang, Liu, Liu and Hendrickson2020), attempting to explain the observations as consequences of an equilibrium between restraining and stimulating states. We identify the restraining state as Hsp70R-ATP, abbreviated as R, and the stimulating state as Hsp70S-ATP, S for short. Both of these are ATP-bound states having NBD and SBD engaged for inter-site communication, as observed for R(Kityk et al., Reference Kityk, Kopp, Sinning and Mayer2012; Qi et al., Reference Qi, Sarbeng, Liu, Le, Xu, Xu, Yang, Wong, Vorvis, Hendrickson, Zhou and Liu2013 ) and proposed for S, whereas the binding domains are flexibly linked and thereby uncoupled in ADP and Apo states (Bertelsen et al., Reference Bertelsen, Chang, Gestwicki and Zuiderweg2009).
We first studied a model featuring an R state that cannot bind client peptides and hydrolyzes ATP at a basal rate being in equilibrium with an S state that can bind peptides to form the SP state, with both S and SP hydrolyzing ATP at a more elevated rate. We found that hydrolysis data were fitted well by this model, but that resulting parameters underestimated the apparent peptide affinity in ATP. We then tested a model that allowed for client peptide binding to R as well as S, but found such binding to be incompatible with the hydrolysis data. Finally, we elaborated an alternative model to include a quasi-intermediate conformation Q that can bind peptides as in S but hydrolyzes ATP at the basal rate as in R, but now with an R that cannot bind peptides.
The equilibria and hydrolytic reactions relating the various states in these alternative models are illustrated schematically in Fig. 1, and mathematical details are described in following sections. The system is complex, even as here stripped of cofactors and partner chaperones. These models are neither fully comprehensive nor singularly unique, and the number of parameters may challenge experimental evaluation. Nevertheless, the theory is providing useful insights as found from experimental validations described below and as used in designing constructs that captured the postulated S state (Wang et al., Reference Wang, Liu, Liu and Hendrickson2020).
Network of Hsp70 equilibria and hydrolytic reactions. Symbols S, R, Q, UD and UZ correspond to the stimulating, restraining, quasi-intermediate, uncoupled ADP and uncoupled Apo states, respectively, and SP, RP, QP UDP and UZP are the corresponding peptide complexes. States are shown connected by equilibria, reactions or single-line designations of identity (––––). Equilibria constants KeqS, KD0S, KD0R, KDU, KDADP, KDATP, KeqQ and KD0Q, are defined by Eqs. (1)–(6) and (34), respectively, and catalytic rate constants k′ and k0 are defined by Eqs. (13) and (14), respectively. Core exchanges of the allosteric system involve only R, S, SP, U and UP states, and these are indicated by red symbols and are contained in the red box, where substrate peptide PS is expelled upon ATP binding to UzP. Exchanges that also include RP are contained in the green box, where ATP binding to UzP yields RP. Exchanges that also include Q and QP are contained in the blue box, where ATP binding to UzP yields QP. The red core is the subset of the blue set with no peptide binding to R (KD0R = ∞).
Our analytic approach is to evaluate the rate of ATP hydrolysis (kcat) as measured in single-turnover reactions at steady state (Wang et al., Reference Wang, Liu, Liu and Hendrickson2020; Davis et al., Reference Davis, Voisine and Craig1999) and the apparent dissociation constant for peptide binding (KDApp) in the presence of ATP, both as functions of hydrolysis rates and intrinsic equilibrium constants (Fig. 1). Because of the hydrolytic reactions, the system cannot be analyzed usefully at equilibrium; however, the steady-state analysis is feasible since Hsp70-catalyzed ATP hydrolysis is relatively slow (k0 = 0.0075 min−1 and k′ = 0.276 min−1; DnaK at 20°C, Wang et al., Reference Wang, Liu, Liu and Hendrickson2020) compared to the kinetics of relevant conformational changes. This allows the distribution of species to equilibrate at any instant. We do not know the dynamics of R-to-S exchange; however, measurements have been reported for kinetics of the uncoupled (U)-to-R conformational changes (498 min−1 for SBDα, 1,656 min−1 for SBDβ, and 6,420 min−1 for the NBD–SBD linker; DnaK at 30°C; Kityk et al., Reference Kityk, Kopp, Sinning and Mayer2012) and for peptide binding and release (kon = 198 min−1 and koff = 1.44 × 106 min−1 M−1, i.e. 2.9 min−1 at 2 μM DnaK; DnaK plus Cro peptide at 25°C; Slepenkov and Witt, Reference Slepenkov and Witt S1998).
Allosteric coupling of substrate-peptide binding to ATP hydrolysis
The equilibria between ATP states are described by Eqs. (1)–(3), where KeqS is the equilibrium constant between R and S states, KD0S is the dissociation constant that relates the complex SP of Hsp70S-ATP with a substrate peptide P to its dissociated constituents, S and P, and KD0R is the dissociation constant that relates the complex RP between Hsp70R-ATP and a substrate peptide P to its dissociated products, R and P. The corresponding kinetic rates for peptide association and dissociation are specified as konS, koffS, konR and koffR.
$$ {\mathrm{K}}_{\mathrm{eqS}}=\frac{\left[\mathrm{R}\right]}{\left[\mathrm{S}\right]}, $$
$$ {{\mathrm{K}}_{\mathrm{D}}}^{0\mathrm{S}}=\frac{\left[\mathrm{S}\right]\ \left[\mathrm{P}\right]}{\left[\mathrm{S}\mathrm{P}\right]}=\frac{{{\mathrm{k}}_{\mathrm{off}}}^{\mathrm{S}}}{{{\mathrm{k}}_{\mathrm{on}}}^{\mathrm{S}}}, $$
$$ {{\mathrm{K}}_{\mathrm{D}}}^{0\mathrm{R}}=\frac{\left[\mathrm{R}\right]\ \left[\mathrm{P}\right]}{\left[\mathrm{R}\mathrm{P}\right]}=\frac{{{\mathrm{k}}_{\mathrm{off}}}^{\mathrm{R}}}{{{\mathrm{k}}_{\mathrm{on}}}^{\mathrm{R}}}. $$Substrate-peptide binding to allosterically uncoupled ADP and Apo Hsp70s
We postulate that Hsp70S-ATP hydrolyzes ATP to Hsp70U-ADP (U) and SP hydrolyzes ATP to yield UP, both at the rate of k′ whether peptide is bound or not, and Hsp70R-ATP states, R and RP, both hydrolyze ATP at the basal rate of k0 to produce the respective ADP states U and UP. U and UP are related by the dissociation equilibrium Eq. (4), which is defined by a dissociation constant KDU that comprises the kinetic constants konU and koffU. We expect the constants of Eq. (4) to be the same for any state of SBD that is uncoupled (hence the superscript U) from NBD, including isolated SBD or nucleotide-free Hsp70 as well as Hsp70U-ADP:
$$ {{\mathrm{K}}_{\mathrm{D}}}^{\mathrm{U}}=\frac{\left[\mathrm{U}\right]\ \left[\mathrm{P}\right]}{\left[\mathrm{U}\mathrm{P}\right]}=\frac{{{\mathrm{k}}_{\mathrm{off}}}^{\mathrm{U}}}{{{\mathrm{k}}_{\mathrm{on}}}^{\mathrm{U}}}. $$This reaction does not affect directly the reactions associated with the ATP-bound states, and it is not needed further for this analysis of allosteric interactions.
ADP dissociation from Hsp70 and ATP rebinding to this Apo state
According with observation, we postulate that ADP can dissociate from Hsp70U-ADP (D) to yield Hsp70U-Apo as related by Eq. (5):
$$ {{\mathrm{K}}_{\mathrm{D}}}^{\mathrm{ADP}}=\frac{\left[\mathrm{Hsp}7{0}_{\mathrm{U}}\harrowextender \mathrm{Apo}\right]\left[\mathrm{ADP}\right]}{\left[\mathrm{Hsp}7{0}_{\mathrm{U}}\harrowextender \mathrm{ADP}\right]} $$and that re-association of ATP with Hsp70U-Apo is governed by Eq. (6):
$$ {{\mathrm{K}}_{\mathrm{D}}}^{\mathrm{ATP}}=\frac{\left[\mathrm{Hsp}7{0}_{\mathrm{U}}\harrowextender \mathrm{Apo}\right]\left[\mathrm{ATP}\right]}{\left[\mathrm{Hsp}7{0}_{\mathrm{U}}\harrowextender \mathrm{ATP}\right]}. $$The inorganic phosphate (Pi) product of ATP hydrolysis remains Hsp70-bound with ADP after ATP hydrolysis (Flaherty et al., Reference Flaherty, DeLuca-Flaherty and McKay1990; Sriram et al., Reference Sriram, Osipiuk, Freeman, Morimoto and Joachimiak1997; Wang and Hendrickson, Reference Wang and Hendrickson2020a); thus, Eq. (5) connotes both ADP and Pi release. Since NBD and SBD are presumed uncoupled in the ADP and Apo states, Eqs. (5) and (6) apply equally to peptide-bound and peptide-free states of Hsp70U-ADP and Hsp70U-Apo. These associations apply when Hsp70 alone can interact freely with the reactants, which typically exist with ATP in excess of ADP; however, the consequent nucleotide exchange reactions can be accelerated substantially by nucleotide exchange factors, notably GrpE for DnaK or Hsp110 for eukaryotic Hsp70s.
Distribution of Hsp70-ATP between restraining and stimulating states
The total concentration of Hsp70-ATP, cT, is given by Eq. (7):
$$ {\mathrm{c}}_{\mathrm{T}}=\left[\mathrm{R}\right]+\left[\mathrm{R}\mathrm{P}\right]+\left[\mathrm{S}\right]+\left[\mathrm{S}\mathrm{P}\right]. $$cT(t) changes as time proceeds because of ATP hydrolysis and possible ATP binding to the apo state or exchange of ATP for ADP; however, we assume that the exchanges governed by equilibria (1)–(3) are sufficiently rapid that cT(t) at any instant is governed by Eq. (8):
$$ {\mathrm{c}}_{\mathrm{T}}=\quad0.5em {{\mathrm{K}}_{\mathrm{D}}}^{0\mathrm{S}}\left\{\left({\mathrm{K}}_{\mathrm{eqS}}\left(1+\left[\mathrm{P}\right]/{{\mathrm{K}}_{\mathrm{D}}}^{0\mathrm{R}}\right)+1+\left[\mathrm{P}\right]/{{\mathrm{K}}_{\mathrm{D}}}^{0\mathrm{S}}\right)\right\}\ \left[\mathrm{SP}\right]/\left[\mathrm{P}\right] $$or
$$ \left[\mathrm{SP}\right]=\frac{{\mathrm{c}}_{\mathrm{T}}\left[\mathrm{P}\right]}{{{\mathrm{K}}_{\mathrm{D}}}^{0\mathrm{S}}\left({\mathrm{K}}_{\mathrm{eqS}}+1\right)+\left[\mathrm{P}\right]\ \{1+\left({{\mathrm{K}}_{\mathrm{D}}}^{0\mathrm{S}}/{{\mathrm{K}}_{\mathrm{D}}}^{0\mathrm{R}}\right)\ {\mathrm{K}}_{\mathrm{eqS}}\}}. $$From Eq. (7), the fraction QS of the total Hsp70-ATP protein that is in the stimulating state is given by Eq. (10):
$$ {\mathrm{Q}}_{\mathrm{S}}=\frac{\left[\mathrm{S}\right]+\left[\mathrm{S}\mathrm{P}\right]}{\left[\mathrm{R}\right]+\left[\mathrm{R}\mathrm{P}\right]+\left[\mathrm{S}\right]+\left[\mathrm{S}\mathrm{P}\right]} $$and in light of Eqs. (1) and (2), this yields
$$ {\mathrm{Q}}_{\mathrm{S}}\left(\left[\mathrm{P}\right]\right)=\frac{\left[\mathrm{P}\right]+{{\mathrm{K}}_{\mathrm{D}}}^{0\mathrm{S}}}{\left[\mathrm{P}\right]\ \left(1+\left({{\mathrm{K}}_{\mathrm{D}}}^{0\mathrm{S}}/{{\mathrm{K}}_{\mathrm{D}}}^{0\mathrm{R}}\right)\ {\mathrm{K}}_{\mathrm{eqS}}\right)+{{\mathrm{K}}_{\mathrm{D}}}^{0\mathrm{S}}\left({\mathrm{K}}_{\mathrm{eqS}}+1\right)}. $$And for the case of KD0R = ∞, the limit of no peptide binding to the restraining state,
$$ {\mathrm{Q}}_{\mathrm{S}}\left(\left[\mathrm{P}\right]\right)=\frac{\left[\mathrm{P}\right]+{{\mathrm{K}}_{\mathrm{D}}}^{0\mathrm{S}}}{\left[\mathrm{P}\right]+{{\mathrm{K}}_{\mathrm{D}}}^{0\mathrm{S}}\left({\mathrm{K}}_{\mathrm{eqS}}+1\right)}. $$ATP hydrolysis in the presence of a peptide substrate
We first wish to understand the allosteric control of ATP hydrolysis by peptide binding. We postulate that Hsp70S-ATP hydrolyzes ATP to ADP and Pi at the same rate, k′, whether complexed with peptide or not, in state SP or S. We similarly postulate that Hsp70R-ATP hydrolyzes ATP to ADP and Pi at its own rate, k0, again whether complexed with peptide or not, that is as RP or R. These reactions are designated in Eqs. (13) and (14):
$$ \mathrm{S}\ \left(\mathrm{ATP}\right)\overset{\mathrm{k}^{\prime }}{\to}\mathrm{U}\left(\mathrm{ADP}\right),\quad3em \mathrm{SP}\ \left(\mathrm{ATP}\right)\overset{\mathrm{k}^{\prime }}{\to}\mathrm{U}\mathrm{P}\ \left(\mathrm{ADP}\right), $$
$$ \mathrm{R}\ \left(\mathrm{ATP}\right)\overset{{\mathrm{k}}^0}{\to}\mathrm{U}\left(\mathrm{ADP}\right),\quad3em \mathrm{RP}\ \left(\mathrm{ATP}\right)\overset{{\mathrm{k}}^0}{\to}\mathrm{U}\mathrm{P}\ \left(\mathrm{ADP}\right). $$For ATP hydrolysis, each state contributes to the observed catalytic rate in proportion to its relative abundance and the associated rate of hydrolysis as described by Eq. (15):
$$ {\mathrm{k}}_{\mathrm{c}\mathrm{at}}=\left\{\ \left[\mathrm{R}\right]\ {\mathrm{k}}^0+\left[\mathrm{R}\mathrm{P}\right]\ {\mathrm{k}}^0+\left[\mathrm{S}\right]\ \mathrm{k}^{\prime }+\left[\mathrm{S}\mathrm{P}\right]\ \mathrm{k}^{\prime}\right\}/{\mathrm{c}}_{\mathrm{T}}. $$Taking Eq. (7) into account, Eq. (15) yields
$$ {\mathrm{k}}_{\mathrm{c}\mathrm{at}}=\left\{\ \left({\mathrm{c}}_{\mathrm{T}}\hbox{--} \left[\mathrm{S}\right]\hbox{--} \left[\mathrm{S}\mathrm{P}\right]\right)\ {\mathrm{k}}^0+\left(\left[\mathrm{S}\right]+\left[\mathrm{S}\mathrm{P}\right]\right)\ \mathrm{k}^{\prime}\right\}/{\mathrm{c}}_{\mathrm{T}} $$
$$ {\mathrm{k}}_{\mathrm{c}\mathrm{at}}={\mathrm{k}}^0+\left(\left[\mathrm{S}\right]+\left[\mathrm{S}\mathrm{P}\right]\right)\ \left(\mathrm{k}^{\prime }-{\mathrm{k}}^0\right)/{\mathrm{c}}_{\mathrm{T}}. $$With reference to Eq. (2), we obtain Eq. (16):
$$ {\mathrm{k}}_{\mathrm{c}\mathrm{at}}={\mathrm{k}}^0+\left(\left[\mathrm{P}\right]+{{\mathrm{K}}_{\mathrm{D}}}^{0\mathrm{S}}\right)\ \left(\mathrm{k}^{\prime }-{\mathrm{k}}^0\right)\ \left[\mathrm{SP}\right]/\left({\mathrm{c}}_{\mathrm{T}}\left[\mathrm{P}\right]\right). $$Then, upon substitution of Eq. (9) into Eq. (16),
$$ {\mathrm{k}}_{\mathrm{cat}}={\mathrm{k}}^0+\frac{\left[\mathrm{P}\right]+{{\mathrm{K}}_{\mathrm{D}}}^{0\mathrm{S}}}{\left[\mathrm{P}\right]\ \left(1+\left({{\mathrm{K}}_{\mathrm{D}}}^{0\mathrm{S}}/{{\mathrm{K}}_{\mathrm{D}}}^{0\mathrm{R}}\right)\ {\mathrm{K}}_{\mathrm{eqS}}\right)+{{\mathrm{K}}_{\mathrm{D}}}^{0\mathrm{S}}\left({\mathrm{K}}_{\mathrm{eqS}}+1\right)}\times \left(\mathrm{k}^{\prime}\hbox{--} {\mathrm{k}}^0\right). $$Equation 17 can be rearranged to yield the mathematically explicit form
$$ {\mathrm{k}}_{\mathrm{c}\mathrm{at}}=\frac{\left[\mathrm{P}\right]\ \left(\mathrm{k}^{\prime }+\left({{\mathrm{K}}_{\mathrm{D}}}^{0\mathrm{S}}/{{\mathrm{K}}_{\mathrm{D}}}^{0\mathrm{R}}\right)\ {\mathrm{K}}_{\mathrm{eqS}}{\mathrm{k}}^0\right)+{{\mathrm{K}}_{\mathrm{D}}}^{0\mathrm{S}}\left({\mathrm{K}}_{\mathrm{eqS}}{\mathrm{k}}^0+{\mathrm{k}}^{\prime}\right)}{\left[\mathrm{P}\right]\ \left(1+\left({{\mathrm{K}}_{\mathrm{D}}}^{0\mathrm{S}}/{{\mathrm{K}}_{\mathrm{D}}}^{0\mathrm{R}}\right)\ {\mathrm{K}}_{\mathrm{eqS}}\right)+{{\mathrm{K}}_{\mathrm{D}}}^{0\mathrm{S}}\left({\mathrm{K}}_{\mathrm{eqS}}+1\right)}=\frac{\mathrm{a}\ \left[\mathrm{P}\right]+\mathrm{b}}{\mathrm{c}\ \left[\mathrm{P}\right]+\mathrm{d}}. $$A particular case, which was in fact motivating to our analysis, arises when the restraining state has no affinity for peptide substrates; that is when KD0R = ∞. Then, Eqs. (17) and (18) reduce to
$$ {\mathrm{k}}_{\mathrm{cat}}={\mathrm{k}}^0+\frac{\left[\mathrm{P}\right]+{{\mathrm{K}}_{\mathrm{D}}}^{0\mathrm{S}}}{\left[\mathrm{P}\right]+{{\mathrm{K}}_{\mathrm{D}}}^{0\mathrm{S}}\left({\mathrm{K}}_{\mathrm{eqS}}+1\right)}\ \mathrm{x}\ \left(\mathrm{k}^{\prime}\hbox{--} {\mathrm{k}}^0\right), $$and
$$ {\mathrm{k}}_{\mathrm{cat}}=\frac{\mathrm{k}^{\prime}\left[\mathrm{P}\right]+{{\mathrm{K}}_{\mathrm{D}}}^{0\mathrm{S}}\left({\mathrm{K}}_{\mathrm{eqS}}{\mathrm{k}}^0+{\mathrm{k}}^{\prime}\right)}{\left[\mathrm{P}\right]+{{\mathrm{K}}_{\mathrm{D}}}^{0\mathrm{S}}\left({\mathrm{K}}_{\mathrm{eqS}}+1\right)}=\frac{\mathrm{a}\ \left[\mathrm{P}\right]+\mathrm{b}}{\left[\mathrm{P}\right]+\mathrm{d}}. $$Degeneracy in allosteric parameters
The mathematical form of Eq. (20) has only three independent parameters: a = k′, b = KD0S (KeqS k0 + k′), and d = KD0S (KeqS + 1) even though the theory is formulated in terms of four physically meaningful parameters, even as simplified by ignoring peptide binding in the restraining state as for Eqs. (17) and (18). Thus, a degeneracy in solutions must arise from the fitting of measurements of hydrolytic rate kcat at varied peptide concentrations [P]; k′ is determined uniquely, but only the b and d combinations of other physical parameters are determined uniquely. In principle, after either KD0S or KeqS is specified, or in certain special cases (e.g. k0 = 0), then the other parameters can be separated. In practice, we break the degeneracy by fixing one parameter from separate measurements. For example, in studies on Hsp70 DnaK, we are able to generate conditions that fix the protein in the stimulating state whereby peptide binding in the presence of ATP serves to define KD0S.
Mathematical fitting
Parameters a, b, c and d of Eq. (18) can only be determined relative to a common factor, which is most conveniently taken as c, since c = 1 corresponds to KD0R = ∞ for no peptide affinity in the R state. Thereby,
$$ {\mathrm{k}}_{\mathrm{cat}}=\frac{\mathrm{a}^{\prime}\left[\mathrm{P}\right]+\mathrm{b}^{\prime }}{\left[\mathrm{P}\right]+\mathrm{d}^{\prime }}, $$where a′ = a/c, b′ = b/c, c′ = 1 and d′ = d/c. The formalism of Eq. (21) permits unique fitting to a set of kcat versus [P] data; however, biochemical interpretations can then be made for any arbitrary value of c using Eq. (18) provided that the mathematical degeneracy is broken by specifying one parameter in formalism Eq. (20), which we take here to be that for KD0S. Then from the coefficients of Eq. (18), we can evaluate the allosteric parameters at arbitrary scalings c:
$$ \mathrm{From}\quad1em \mathrm{d}^{\prime }={{\mathrm{K}}_{\mathrm{D}}}^{0\mathrm{S}}\left({\mathrm{K}}_{\mathrm{eqS}}+1\right),{\mathrm{K}}_{\mathrm{eqS}}=\left(\mathrm{d}/\mathrm{c}{{\mathrm{K}}_{\mathrm{D}}}^{0\mathrm{S}}\right)\hbox{--} 1. $$
$$ \mathrm{From}\quad1em \mathrm{c}=1+\left({{\mathrm{K}}_{\mathrm{D}}}^{0\mathrm{S}}/{{\mathrm{K}}_{\mathrm{D}}}^{0\mathrm{R}}\right)\ {\mathrm{K}}_{\mathrm{eqS}},\quad1.75em {{\mathrm{K}}_{\mathrm{D}}}^{0\mathrm{R}}={\mathrm{K}}_{\mathrm{eqS}}{{\mathrm{K}}_{\mathrm{D}}}^{0\mathrm{S}}/\left(\mathrm{c}\hbox{--} 1\right). $$
$$ {{\mathrm{K}}_{\mathrm{D}}}^{0\mathrm{R}}=\infty \quad1em \mathrm{when}\quad1em \mathrm{c}=1 $$
$$ {{\mathrm{K}}_{\mathrm{D}}}^{0\mathrm{S}}\mathrm{k}^{\prime }+{{\mathrm{K}}_{\mathrm{D}}}^{0\mathrm{S}}{\mathrm{K}}_{\mathrm{eqS}}{\mathrm{k}}^0=\mathrm{b}^{\prime } $$
$$ \mathrm{k}^{\prime }+{\mathrm{K}}_{\mathrm{eqS}}\left({{\mathrm{K}}_{\mathrm{D}}}^{0\mathrm{S}}/{{\mathrm{K}}_{\mathrm{D}}}^{0\mathrm{R}}\right)\ {\mathrm{k}}^0=\mathrm{a}^{\prime } $$
$$ {{\mathrm{K}}_{\mathrm{D}}}^{0\mathrm{S}}\mathrm{k}^{\prime }+{{\mathrm{K}}_{\mathrm{D}}}^{0\mathrm{S}}{\mathrm{K}}_{\mathrm{eqS}}\left({{\mathrm{K}}_{\mathrm{D}}}^{0\mathrm{S}}/{{\mathrm{K}}_{\mathrm{D}}}^{0\mathrm{R}}\right){\mathrm{k}}^0=\mathrm{a}^{\prime }{{\mathrm{K}}_{\mathrm{D}}}^{0\mathrm{S}} $$
$$ {\mathrm{k}}^0=\left(\mathrm{b}-{{\mathrm{K}}_{\mathrm{D}}}^{0\mathrm{S}}\mathrm{a}\right)/\mathrm{c}\ {{\mathrm{K}}_{\mathrm{D}}}^{0\mathrm{S}}{\mathrm{K}}_{\mathrm{eqS}}\left(1-{{\mathrm{K}}_{\mathrm{D}}}^{0\mathrm{S}}/{{\mathrm{K}}_{\mathrm{D}}}^{0\mathrm{R}}\right). $$
$$ {\mathrm{k}}^0=\left(\mathrm{b}-{{\mathrm{K}}_{\mathrm{D}}}^{0\mathrm{S}}\mathrm{a}\right)/{{\mathrm{K}}_{\mathrm{D}}}^{0\mathrm{S}}{\mathrm{K}}_{\mathrm{eqS}}\quad1.25em \mathrm{when}\ \mathrm{c}=1 $$From (25)
$$ \mathrm{k}^{\prime }=\mathrm{a}/\mathrm{c}-{\mathrm{K}}_{\mathrm{eqS}}\left({{\mathrm{K}}_{\mathrm{D}}}^{0\mathrm{S}}/{{\mathrm{K}}_{\mathrm{D}}}^{0\mathrm{R}}\right)\ {\mathrm{k}}^0 $$
$$ \mathrm{k}^{\prime }=\mathrm{a}\quad3em \mathrm{when}\quad.5em {{\mathrm{K}}_{\mathrm{D}}}^{0\mathrm{R}}=\infty\ \left(\mathrm{c}=1\right). $$To summarize, given the fitting with parameters a, b and d, and an arbitrarily chosen value of c and a specified value of KD0S, the other biochemical parameters are determined by Eqs. (22)–(28) for KeqS, KD0R, k0 and k′, respectively.
Alternative allosteric models for hydrolysis
In Eqs. (13) and (14), we postulate a hydrolysis model having rates of k′ and k0 for S and R states, respectively, whether with substrate peptide or not. We can also contemplate an alternative model wherein the S state is in a restraining conformation such that its hydrolysis rate is k0 until peptide binding generates SP with hydrolytic rate k′. More generally, S might have an arbitrary hydrolytic rate kS, not necessarily either k′ or k0. In this case, Eq. 15 will be replaced by Eq. 29:
$$ {\mathrm{k}}_{\mathrm{c}\mathrm{at}}=\left\{\left[\mathrm{R}\right]\ {\mathrm{k}}^0+\left[\mathrm{R}\mathrm{P}\right]\ {\mathrm{k}}^0+\left[\mathrm{S}\right]\ {\mathrm{k}}^{\mathrm{S}}+\left[\mathrm{S}\mathrm{P}\right]\ \mathrm{k}^{\prime}\right\}/{\mathrm{c}}_{\mathrm{T}}. $$Then proceeding as from Eqs. (15) to (17)
$$ {\mathrm{k}}_{\mathrm{cat}}={\mathrm{k}}^0+\frac{\left[\mathrm{P}\right]\ \left(\mathrm{k}^{\prime}\hbox{--} {\mathrm{k}}^0\right)+{{\mathrm{K}}_{\mathrm{D}}}^{0\mathrm{S}}\left({\mathrm{k}}^{\mathrm{S}}\hbox{--} {\mathrm{k}}^0\right)}{\left[\mathrm{P}\right]\ \left(1+\left({{\mathrm{K}}_{\mathrm{D}}}^{0\mathrm{S}}/{{\mathrm{K}}_{\mathrm{D}}}^{0\mathrm{R}}\right)\ {\mathrm{K}}_{\mathrm{eqS}}\right)+{{\mathrm{K}}_{\mathrm{D}}}^{0\mathrm{S}}\left({\mathrm{K}}_{\mathrm{eqS}}+1\right)}, $$which can be recast, as for Eq. (18) from Eq. (17), into
$$ {\mathrm{k}}_{\mathrm{c}\mathrm{at}}=\frac{\left[\mathrm{P}\right]\ \left(\mathrm{k}^{\prime }+\left({{\mathrm{K}}_{\mathrm{D}}}^{0\mathrm{S}}/{{\mathrm{K}}_{\mathrm{D}}}^{0\mathrm{R}}\right)\ {\mathrm{K}}_{\mathrm{eqS}}{\mathrm{k}}^0\right)+{{\mathrm{K}}_{\mathrm{D}}}^{0\mathrm{S}}\left({\mathrm{K}}_{\mathrm{eqS}}{\mathrm{k}}^0+{\mathrm{k}}^{\mathrm{S}}\right)}{\left[\mathrm{P}\right]\ \left(1+\left({{\mathrm{K}}_{\mathrm{D}}}^{0\mathrm{S}}/{{\mathrm{K}}_{\mathrm{D}}}^{0\mathrm{R}}\right)\ {\mathrm{K}}_{\mathrm{eqS}}\right)+{{\mathrm{K}}_{\mathrm{D}}}^{0\mathrm{S}}\left({\mathrm{K}}_{\mathrm{eqS}}+1\right)}=\frac{\mathrm{a}\ \left[\mathrm{P}\right]+{\mathrm{b}}_{\mathrm{S}}}{\mathrm{c}\ \left[\mathrm{P}\right]+\mathrm{d}}. $$Equations (30) and (31) reduce respectively to Eqs. (17) and (18) for kS = k′, and Eq. (31) gives bS = b0 = KD0S (KeqS + 1) k0 for kS = k0.
Structural and biochemical evidence lead us to contemplate another alternative model wherein the R-state itself does not bind substrate peptides, but which is instead in equilibrium with a quasi-intermediate state Q. Q has an S-like SBD conformation, which is peptide associative, and it has an R-like NBD-SBD interface such that its hydrolysis rate is k0. Thus,
$$ {\mathrm{c}}_{\mathrm{T}}=\left[\mathrm{R}\right]+\left[\mathrm{Q}\right]+\left[\mathrm{Q}\mathrm{P}\right]+\left[\mathrm{S}\right]+\left[\mathrm{S}\mathrm{P}\right]. $$and
$$ {\displaystyle \begin{array}{c}{\mathrm{k}}_{\mathrm{c}\mathrm{at}}=\left\{\ \left[\mathrm{R}\right]\ {\mathrm{k}}^0+\left[\mathrm{Q}\right]\ {\mathrm{k}}^0+\left[\mathrm{Q}\mathrm{P}\right]\ {\mathrm{k}}^0+\left[\mathrm{S}\right]\ \mathrm{k}^{\prime }+\left[\mathrm{S}\mathrm{P}\right]\ \mathrm{k}^{\prime}\right\}/{\mathrm{c}}_{\mathrm{T}}.\\ {}\quad-9.25em = {\mathrm{k}}^0+\left(\left[\mathrm{S}\right]+\left[\mathrm{S}\mathrm{P}\right]\right)\ \left(\mathrm{k}^{\prime }-{\mathrm{k}}^0\right)/{\mathrm{c}}_{\mathrm{T}}.\end{array}} $$As for obtaining Eq. (16) from Eq. (15), we refer to Eqs. (1) and (2) and here also add the Q-state conformational and binding equilibria:
$$ {\mathrm{K}}_{\mathrm{eqQ}}=\frac{\left[\mathrm{R}\right]}{\left[\mathrm{Q}\right]};\quad2.5em {{\mathrm{K}}_{\mathrm{D}}}^{0\mathrm{Q}} = \frac{\left[\mathrm{Q}\right]\ \left[\mathrm{P}\right]}{\left[\mathrm{Q}\mathrm{P}\right]}. $$Then,
$$ {\mathrm{c}}_{\mathrm{T}} = {{\mathrm{K}}_{\mathrm{D}}}^{0\mathrm{S}}\left(\left[\mathrm{SP}\right]/\left[\mathrm{P}\right]\right)\ \left\{{\mathrm{K}}_{\mathrm{eqS}}+{\mathrm{K}}_{\mathrm{eqS}}/{\mathrm{K}}_{\mathrm{eqQ}}\left(1+\left[\mathrm{P}\right]/{{\mathrm{K}}_{\mathrm{D}}}^{0\mathrm{Q}}\right)+\left(1+\left[\mathrm{P}\right]/{{\mathrm{K}}_{\mathrm{D}}}^{0\mathrm{S}}\right)\right\} $$and
$$ {\mathrm{k}}_{\mathrm{cat}}={\mathrm{k}}^0+\frac{\left[\mathrm{P}\right]+{{\mathrm{K}}_{\mathrm{D}}}^{0\mathrm{S}}}{\left[\mathrm{P}\right]\ \Big(1+\left({{\mathrm{K}}_{\mathrm{D}}}^{0\mathrm{S}}{\mathrm{K}}_{\mathrm{eqS}}/{{\mathrm{K}}_{\mathrm{D}}}^{0\mathrm{Q}}{\mathrm{K}}_{\mathrm{eqQ}}\right)+{{\mathrm{K}}_{\mathrm{D}}}^{0\mathrm{S}}\left({\mathrm{K}}_{\mathrm{eqS}}+1\right)}\mathrm{x}\ \left(\mathrm{k}^{\prime}\hbox{--} {\mathrm{k}}^0\right). $$
$$ =\frac{\left[\mathrm{P}\right]\ \left(\mathrm{k}^{\prime }+\left({{\mathrm{K}}_{\mathrm{D}}}^{0\mathrm{S}}{\mathrm{K}}_{\mathrm{eqS}}/{{\mathrm{K}}_{\mathrm{D}}}^{0\mathrm{Q}}{\mathrm{K}}_{\mathrm{eqQ}}\right)\ {\mathrm{k}}^0\right)+{{\mathrm{K}}_{\mathrm{D}}}^{0\mathrm{S}}\left({\mathrm{K}}_{\mathrm{eqS}}+{\mathrm{K}}_{\mathrm{eqS}}/{\mathrm{K}}_{\mathrm{eqQ}}\right)\ {\mathrm{k}}^0+{\mathrm{k}}^{\prime}\Big)}{\left[\mathrm{P}\right]\ \left(1+\left({{\mathrm{K}}_{\mathrm{D}}}^{0\mathrm{S}}{\mathrm{K}}_{\mathrm{eqS}}/{{\mathrm{K}}_{\mathrm{D}}}^{0\mathrm{Q}}{\mathrm{K}}_{\mathrm{eqQ}}\right)\right)+{{\mathrm{K}}_{\mathrm{D}}}^{0\mathrm{S}}\left({\mathrm{K}}_{\mathrm{eqS}}+{\mathrm{K}}_{\mathrm{eqS}}/{\mathrm{K}}_{\mathrm{eqQ}}+1\right)}=\frac{\mathrm{a}\ \left[\mathrm{P}\right]+\mathrm{b}}{\mathrm{c}\ \left[\mathrm{P}\right]+\mathrm{d}}. $$Substrate peptide binding in the presence of ATP
Substrate model of Equations (1) to (3)
We also wish to understand the allosteric effect of ATP on substrate peptide binding. The R and S states are not differentiated in typical peptide binding experiments, whereby the apparent dissociation constant that can be observed is (38):
$$ {{\mathrm{K}}_{\mathrm{D}}}^{\mathrm{App}}\left(\mathrm{ATP}\right)=\left[\mathrm{P}\right]\ \left(\ \left[\mathrm{R}\right]+\left[\mathrm{S}\right]\ \right)/\left(\left[\mathrm{R}\mathrm{P}\right]+\left[\mathrm{S}\mathrm{P}\right]\ \right), $$which by Eq. (1) gives
$$ {{\mathrm{K}}_{\mathrm{D}}}^{\mathrm{App}}\left(\mathrm{ATP}\right)=\left[\mathrm{P}\right]\ \left[\mathrm{S}\right]\ \left({\mathrm{K}}_{\mathrm{eqS}}+1\right)/\left(\left[\mathrm{RP}\right]+\left[\mathrm{S}\mathrm{P}\right]\ \right), $$and with reference to Eqs. (2) and (3), the apparent and intrinsic dissociation constants in ATP are then related by (39):
$$ {{\mathrm{K}}_{\mathrm{D}}}^{\mathrm{App}}\left(\mathrm{ATP}\right)=\frac{\left[\mathrm{P}\right]\ \left[\mathrm{S}\right]\ \left({\mathrm{K}}_{\mathrm{eqS}}+1\right)}{\left[\mathrm{P}\right]\ \left[\mathrm{S}\right]\ \left({\mathrm{K}}_{\mathrm{eqS}}/{{\mathrm{K}}_{\mathrm{D}}}^{0\mathrm{R}}+1/{{\mathrm{K}}_{\mathrm{D}}}^{0\mathrm{S}}\right)}=\frac{{{\mathrm{K}}_{\mathrm{D}}}^{0\mathrm{R}}{{\mathrm{K}}_{\mathrm{D}}}^{0\mathrm{S}}\left({\mathrm{K}}_{\mathrm{eqS}}+1\right)}{\left({{\mathrm{K}}_{\mathrm{D}}}^{0\mathrm{S}}{\mathrm{K}}_{\mathrm{eqS}}+{{\mathrm{K}}_{\mathrm{D}}}^{0\mathrm{R}}\right)}{{\mathrm{K}}_{\mathrm{D}}}^{\mathrm{App}}\left(\mathrm{ATP}\right)=\frac{{{\mathrm{K}}_{\mathrm{D}}}^{0\mathrm{S}}\left({\mathrm{K}}_{\mathrm{eqS}}+1\right)}{1+{{\mathrm{K}}_{\mathrm{D}}}^{0\mathrm{S}}{\mathrm{K}}_{\mathrm{eqS}}/{{\mathrm{K}}_{\mathrm{D}}}^{0\mathrm{R}}}. $$In the limit of no peptide binding to the restraining state, KD0R = ∞,
$$ {{\mathrm{K}}_{\mathrm{D}}}^{\mathrm{App}}\left(\mathrm{ATP}\right)={{\mathrm{K}}_{\mathrm{D}}}^{0\mathrm{S}}\left({\mathrm{K}}_{\mathrm{eqS}}+1\right). $$Given a measured value for KDApp, the corresponding value for KD0R can be obtained from Eq. (39):
$$ {{\mathrm{K}}_{\mathrm{D}}}^{0\mathrm{R}}=\frac{{{\mathrm{K}}_{\mathrm{D}}}^{\mathrm{App}}\left(\mathrm{ATP}\right)\ {{\mathrm{K}}_{\mathrm{D}}}^{0\mathrm{S}}{\mathrm{K}}_{\mathrm{eqS}}}{{{\mathrm{K}}_{\mathrm{D}}}^{0\mathrm{S}}\left({\mathrm{K}}_{\mathrm{eqS}}+1\right)-{{\mathrm{K}}_{\mathrm{D}}}^{\mathrm{App}}\left(\mathrm{ATP}\right)}. $$Also from Eq. (39), when the equilibrium is entirely toward the restraining state, KeqS = ∞, KDApp(ATP) = KD0R; and when it is entirely toward the stimulating state, KeqS = 0, KDApp(ATP) = KD0S.
Alternative peptide-binding model
As described in the analysis of ATP hydrolysis, we contemplate an alternative allosteric model in which the quasi-intermediate Q state is in equilibrium with the R and S states and in which Q binds substrate peptides but R does not. In this case, analogous to Eq. (38),
$$ {{\mathrm{K}}_{\mathrm{D}}}^{\mathrm{App}}\left(\mathrm{ATP}\right)=\left[\mathrm{P}\right]\ \left(\ \left[\mathrm{R}\right]+\left[\mathrm{Q}\right]+\left[\mathrm{S}\right]\ \right)/\left(\left[\mathrm{Q}\mathrm{P}\right]+\left[\mathrm{S}\mathrm{P}\right]\ \right), $$and then, after substitutions from Eq. (34) and from Eq. (1),
$$ {{\mathrm{K}}_{\mathrm{D}}}^{\mathrm{App}}\left(\mathrm{ATP}\right) = \left[\mathrm{R}\right]\ \left[\mathrm{P}\right]\ \left(1+1/{\mathrm{K}}_{\mathrm{eqQ}}+1/{\mathrm{K}}_{\mathrm{eqS}}\right)/\left[\mathrm{R}\right]\ \left[\mathrm{P}\right]\ \left(1/{{\mathrm{K}}_{\mathrm{D}}}^{0\mathrm{Q}}{\mathrm{K}}_{\mathrm{eqQ}}+1/{{\mathrm{K}}_{\mathrm{D}}}^{0\mathrm{S}}{\mathrm{K}}_{\mathrm{eqS}}\right) $$or
$$ {{\mathrm{K}}_{\mathrm{D}}}^{\mathrm{App}}\left(\mathrm{ATP}\right) = \frac{\left({\mathrm{K}}_{\mathrm{eqQ}}{\mathrm{K}}_{\mathrm{eqS}}+{\mathrm{K}}_{\mathrm{eqS}}+{\mathrm{K}}_{\mathrm{eqQ}}\right)}{{\mathrm{K}}_{\mathrm{eqQ}}{\mathrm{K}}_{\mathrm{eqS}}}\times \frac{{{\mathrm{K}}_{\mathrm{D}}}^{0\mathrm{S}}{\mathrm{K}}_{\mathrm{eqS}}{{\mathrm{K}}_{\mathrm{D}}}^{0\mathrm{Q}}{\mathrm{K}}_{\mathrm{eqQ}}}{\left({{\mathrm{K}}_{\mathrm{D}}}^{0\mathrm{S}}{\mathrm{K}}_{\mathrm{eqS}}+{{\mathrm{K}}_{\mathrm{D}}}^{0\mathrm{Q}}{\mathrm{K}}_{\mathrm{eqQ}}\right)} $$or
$$ {{\mathrm{K}}_{\mathrm{D}}}^{\mathrm{App}}\left(\mathrm{ATP}\right)\quad1em =\frac{{{\mathrm{K}}_{\mathrm{D}}}^{0\mathrm{S}}\left({\mathrm{K}}_{\mathrm{eqS}}+{\mathrm{K}}_{\mathrm{eqS}}/{\mathrm{K}}_{\mathrm{eqQ}}+1\right)}{\left(1+{{\mathrm{K}}_{\mathrm{D}}}^{0\mathrm{S}}{\mathrm{K}}_{\mathrm{eqS}}/{{\mathrm{K}}_{\mathrm{D}}}^{0\mathrm{Q}}{\mathrm{K}}_{\mathrm{eqQ}}\right)}, $$which reduces to (40) as KeqQ = ∞, that is [Q] = 0.
As for its hydrolysis counterpart Eq. (36), the variables in Eq. (44) are too numerous for independent evaluation; however, we can fix KD0S at a measured value and we can estimate KD0Q in relation to that value. A third specification can come from in-parallel fitting to ATP hydrolysis data by Eq. (37). Thus, using d′ from the fitting of hydrolysis data by Eq. (21) while specifying c as the scaling factor for the evaluation of d from Eq. (37), we obtain f = (cd′/KD0S) = KeqS + KeqS/KeqQ + 1; while letting g = KDApp/KD0S and g′ = KDApp/KD0Q, we obtain g + g′ KeqS/KeqQ = KeqS + KeqS/KeqQ + 1 from Eq. (44). These two observational equations then determine values for the relevant unknowns:
$$ {\mathrm{K}}_{\mathrm{eqS}}/{\mathrm{K}}_{\mathrm{eqQ}}=\left(\mathrm{f}\hbox{--} \mathrm{g}\right)/\mathrm{g}^{\prime }, $$
$$ {\mathrm{K}}_{\mathrm{eqS}}=\left(\mathrm{f}-1\right)-\left(\mathrm{f}\hbox{--} \mathrm{g}\right)/\mathrm{g}^{\prime }. $$Since both KeqS/KeqQ and KeqS must be non-negative for the solution to be physical, it follows from (45) that d ≥ KDApp and from (46) that
$$ {{\mathrm{K}}_{\mathrm{D}}}^{0\mathrm{Q}}\le {{\mathrm{K}}_{\mathrm{D}}}^{\mathrm{App}}\left(\mathrm{f}-1\right)/\left(\mathrm{f}\hbox{--} \mathrm{g}\right)=\left(\mathrm{d}-{{\mathrm{K}}_{\mathrm{D}}}^{0\mathrm{S}}\right)/\left(\mathrm{d}-{{\mathrm{K}}_{\mathrm{D}}}^{\mathrm{App}}\right). $$Distribution of states in the Q-alternative model
The states R, S, Q, SP and QP in the Q-alternative model are mutually exclusive, contributing to the total Hsp70 concentration, CT, as given in Eq. (32). After normalization to [SP] through equilibria defined by Eqs. (1), (2) and 34,
$$ {\mathrm{C}}_{\mathrm{T}}={{\mathrm{K}}_{\mathrm{D}}}^{0\mathrm{S}}\left({\mathrm{K}}_{\mathrm{eqS}}+{\mathrm{K}}_{\mathrm{eqS}}/{\mathrm{K}}_{\mathrm{eqQ}}+{\mathrm{K}}_{\mathrm{eqS}}/{\mathrm{K}}_{\mathrm{eqQ}}\left(\left[\mathrm{P}\right]/{{\mathrm{K}}_{\mathrm{D}}}^{0\mathrm{Q}}\right)\quad-11.9em +1+\left[\mathrm{P}\right]/{{\mathrm{K}}_{\mathrm{D}}}^{0\mathrm{S}}\right)\ \left[\mathrm{SP}\right]/\left[\mathrm{P}\right]. $$Then, following as from Eqs. (7) to (11) for the R/S model, we obtain the fractions QS, QR and QQ in the S, R and Q states, respectively, as a function of peptide concentration [P]:
$$ {\mathrm{Q}}_{\mathrm{S}}\left(\left[\mathrm{P}\right]\right)=\frac{\left[\mathrm{P}\right]+{{\mathrm{K}}_{\mathrm{D}}}^{0\mathrm{S}}}{\left[\mathrm{P}\right]\ \left(1+\left({{\mathrm{K}}_{\mathrm{D}}}^{0\mathrm{S}}{\mathrm{K}}_{\mathrm{eqS}}/{{\mathrm{K}}_{\mathrm{D}}}^{0\mathrm{Q}}{\mathrm{K}}_{\mathrm{eqQ}}\right)\right)+{{\mathrm{K}}_{\mathrm{D}}}^{0\mathrm{S}}\left({\mathrm{K}}_{\mathrm{eqS}}+{\mathrm{K}}_{\mathrm{eqS}}/{\mathrm{K}}_{\mathrm{eqQ}}+1\right)}. $$
$$ {\mathrm{Q}}_{\mathrm{R}}\left(\left[\mathrm{P}\right]\right)=\frac{{{\mathrm{K}}_{\mathrm{D}}}^{0\mathrm{S}}{\mathrm{K}}_{\mathrm{eqS}}}{\left[\mathrm{P}\right]\ \left(1+\left({{\mathrm{K}}_{\mathrm{D}}}^{0\mathrm{S}}{\mathrm{K}}_{\mathrm{eqS}}/{{\mathrm{K}}_{\mathrm{D}}}^{0\mathrm{Q}}{\mathrm{K}}_{\mathrm{eqQ}}\right)\right)+{{\mathrm{K}}_{\mathrm{D}}}^{0\mathrm{S}}\left({\mathrm{K}}_{\mathrm{eqS}}+{\mathrm{K}}_{\mathrm{eqS}}/{\mathrm{K}}_{\mathrm{eqQ}}+1\right)}. $$
$$ {\mathrm{Q}}_{\mathrm{Q}}\left(\left[\mathrm{P}\right]\right)=\frac{\left[\mathrm{P}\right]\ \left({{\mathrm{K}}_{\mathrm{D}}}^{0\mathrm{S}}{\mathrm{K}}_{\mathrm{eqS}}/{{\mathrm{K}}_{\mathrm{D}}}^{0\mathrm{Q}}{\mathrm{K}}_{\mathrm{eqQ}}\right)+{{\mathrm{K}}_{\mathrm{D}}}^{0\mathrm{S}}{\mathrm{K}}_{\mathrm{eqS}}/{\mathrm{K}}_{\mathrm{eqQ}}}{\left[\mathrm{P}\right]\ \left(1+\left({{\mathrm{K}}_{\mathrm{D}}}^{0\mathrm{S}}{\mathrm{K}}_{\mathrm{eqS}}/{{\mathrm{K}}_{\mathrm{D}}}^{0\mathrm{Q}}{\mathrm{K}}_{\mathrm{eqQ}}\right)\right)+{{\mathrm{K}}_{\mathrm{D}}}^{0\mathrm{S}}\left({\mathrm{K}}_{\mathrm{eqS}}+{\mathrm{K}}_{\mathrm{eqS}}/{\mathrm{K}}_{\mathrm{eqQ}}+1\right)}. $$In the limit of [P] = 0,
$$ {\mathrm{Q}}_{\mathrm{S}}(0)=1/\left({\mathrm{K}}_{\mathrm{eqS}}+{\mathrm{K}}_{\mathrm{eqS}}/{\mathrm{K}}_{\mathrm{eqQ}}+1\right). $$
$$ {\mathrm{Q}}_{\mathrm{R}}(0)={\mathrm{K}}_{\mathrm{eqS}}/\left({\mathrm{K}}_{\mathrm{eqS}}+{\mathrm{K}}_{\mathrm{eqS}}/{\mathrm{K}}_{\mathrm{eqQ}}+1\right). $$
$$ {\mathrm{Q}}_{\mathrm{Q}}(0)={\mathrm{K}}_{\mathrm{eqS}}/{\mathrm{K}}_{\mathrm{eqQ}}/\left({\mathrm{K}}_{\mathrm{eqS}}+{\mathrm{K}}_{\mathrm{eqS}}/{\mathrm{K}}_{\mathrm{eqQ}}+1\right). $$Peptide binding profiles
Binding characteristics (KD values) can be evaluated from profiles of the saturation of peptide binding to the protein, and calculated saturation curves can be useful for demonstration purposes. Such analyses can be performed either as a function of peptide concentration at fixed protein concentration or of protein concentration at fixed peptide concentration. The various models of conformational and binding equilibria will have different profiles, and we consider two of these here.
Allosteric model with only S-state binding
For allosteric hydrolysis model (19), R and S are in equilibrium by KeqS = [R]/[S] and peptide P binds to S with intrinsic affinity governed by KD0S, Eq. (2), but with no affinity for R (KD0R = ∞). In this case, the saturation curve y for varied peptide concentrations at a fixed total protein concentration derives from
$$ \mathrm{y}=\frac{\left[\mathrm{S}\mathrm{P}\right]}{\left[\mathrm{R}\right]+\left[\mathrm{S}\right]+\left[\mathrm{S}\mathrm{P}\right]}=\frac{1}{\left(\left[\mathrm{R}\right]+\left[\mathrm{S}\right]\right)/\left[\mathrm{S}\mathrm{P}\right]+1\Big)}. $$In light of the relevant equilibria, the total protein concentration is given by (56):
$$ {\mathrm{c}}_{\mathrm{T}}=\left[\mathrm{R}\right]+\left[\mathrm{S}\right]+\left[\mathrm{S}\mathrm{P}\right]=\left[\mathrm{S}\right]\ \left({\mathrm{K}}_{\mathrm{eqS}}+1+\left[\mathrm{P}\right]/{{\mathrm{K}}_{\mathrm{D}}}^{0\mathrm{S}}\right)=\left\{{{\mathrm{K}}_{\mathrm{D}}}^{0\mathrm{S}}\left({\mathrm{K}}_{\mathrm{eqS}}+1\right)+\left[\mathrm{P}\right]\right\}\ \left[\mathrm{S}\mathrm{P}\right]/\left[\mathrm{P}\right], $$whereby
$$ \left[\mathrm{SP}\right]=\frac{{\mathrm{c}}_{\mathrm{T}}\left[\mathrm{P}\right]}{{{\mathrm{K}}_{\mathrm{D}}}^{0\mathrm{S}}\left({\mathrm{K}}_{\mathrm{eqS}}+1\right)+\left[\mathrm{P}\right]} $$and
$$ \left[\mathrm{R}\right]+\left[\mathrm{S}\right]={\mathrm{c}}_{\mathrm{T}}-\left[\mathrm{S}\mathrm{P}\right]=\frac{{\mathrm{c}}_{\mathrm{T}}\left\{{{\mathrm{K}}_{\mathrm{D}}}^{0\mathrm{S}}\left({\mathrm{K}}_{\mathrm{eqS}}+1\right)+\left[\mathrm{P}\right]\right\}-{\mathrm{c}}_{\mathrm{T}}\left[\mathrm{P}\right]}{{{\mathrm{K}}_{\mathrm{D}}}^{0\mathrm{S}}\left({\mathrm{K}}_{\mathrm{eqS}}+1\right)+\left[\mathrm{P}\right]}. $$From the ratio of Eqs. (58) and (57),
$$ \left(\left[\mathrm{R}\right]+\left[\mathrm{S}\right]\right)/\left[\mathrm{S}\mathrm{P}\right]=\frac{{{\mathrm{K}}_{\mathrm{D}}}^{0\mathrm{S}}\left({\mathrm{K}}_{\mathrm{eqS}}+1\right)}{\left[\mathrm{P}\right]}, $$after substitution of Eq. (59) into Eq. (55) and rearrangement, the desired saturation curve results:
$$ \mathrm{y}=\frac{\left[\mathrm{P}\right]}{\left[\mathrm{P}\right]+{{\mathrm{K}}_{\mathrm{D}}}^{0\mathrm{S}}\left({\mathrm{K}}_{\mathrm{eqS}}+1\right)}. $$Alternatively, as in our peptide-binding experiments, one can obtain saturation curve y from varied protein concentrations at a fixed peptide concentration. The appropriate formulation for this situation derives from
$$ \mathrm{y}=\frac{\left[\mathrm{SP}\right]}{\left[\mathrm{P}\right]+\left[\mathrm{SP}\right]}=\frac{1}{\left(\left[\mathrm{P}\right]/\left[\mathrm{SP}\right]+1\right)}. $$From Eq. (56),
$$ \left[\mathrm{P}\right]/\left[\mathrm{SP}\right]=\left\{{{\mathrm{K}}_{\mathrm{D}}}^{0\mathrm{S}}\left({\mathrm{K}}_{\mathrm{eqS}}+1\right)+\left[\mathrm{P}\right]\right\}/{\mathrm{c}}_{\mathrm{T}}, $$and, on substitution of Eq. (62) into Eq. (61), followed by rearrangement,
$$ \mathrm{y}=\frac{{\mathrm{c}}_{\mathrm{T}}}{{\mathrm{c}}_{\mathrm{T}}+\left({{\mathrm{K}}_{\mathrm{D}}}^{0\mathrm{S}}\left({\mathrm{K}}_{\mathrm{eqS}}+1\right)+\left[\mathrm{P}\right]\right)}. $$When there is no R state, that is KeqS = 0, saturation Eqs. (60) and (63) reduce to the respective single-component counterparts Eqs. (64) and (65):
$$ \mathrm{y}=\frac{\left[\mathrm{P}\right]}{\left[\mathrm{P}\right]+{{\mathrm{K}}_{\mathrm{D}}}^{0\mathrm{S}}}, $$
$$ \mathrm{y}=\frac{{\mathrm{c}}_{\mathrm{T}}}{{\mathrm{c}}_{\mathrm{T}}+\left({{\mathrm{K}}_{\mathrm{D}}}^{0\mathrm{S}}+\left[\mathrm{P}\right]\right)}. $$Q-alternative model
The states R, S, Q, SP and QP in the Q-alternative binding model (38) are mutually exclusive and interrelated by equilibria (1), (2) and (34). The saturation model analogous with Eq. (60) can be developed readily, but it suffices here to develop the model analogous with Eq. (63) for measurements made with fixed peptide concentration and varied total protein concentration, which derives from
$$ \mathrm{y}=\frac{\left[\mathrm{QP}\right]+\left[\mathrm{SP}\right]}{\left[\mathrm{P}\right]+\left[\mathrm{QP}\right]+\left[\mathrm{SP}\right]}. $$Proceeding as from Eqs. (61) to (65), successively evaluating ([QP] + [SP]) and [P]/[SP] and then defining p1 = (KD0S KeqS)/(KD0Q KeqQ), Eq. (66) yields the Q-alternative relationship (67):
$$ \mathrm{y}=\frac{{\mathrm{c}}_{\mathrm{T}}}{{\mathrm{c}}_{\mathrm{T}}+\left\{\Big({{\mathrm{K}}_{\mathrm{D}}}^{0\mathrm{S}}\left({\mathrm{K}}_{\mathrm{eqS}}+{\mathrm{K}}_{\mathrm{eqS}}/{\mathrm{K}}_{\mathrm{eqQ}}+1\right)/\left(1+\mathrm{p}1\right)\right\}+\left[\mathrm{P}\right]}. $$Kinetics of substrate peptide association
The kinetics of substrate peptide binding and release are complicated by having Hsp70-ATP in its two states, whether each binds peptide or not. Of course, the dissociation of prebound substrates must be contemplated even after transition to a binding-deficient state, notably SP to RP when konR = 0. Here, we consider the situation where the peptide concentration [P] is much in excess of [Hsp70-ATP], such that [P] can be considered constant and absorbed into pseudo-first order kon rate constants. Thus,
$$ \mathrm{S}+\mathrm{P}\overset{{{\mathrm{k}}_{\mathrm{on}}}^{\mathrm{S}}}{\to}\mathrm{S}\mathrm{P};\quad3em \mathrm{S}\mathrm{P}\overset{{{\mathrm{k}}_{\mathrm{off}}}^{\mathrm{S}}}{\to}\mathrm{S}+\mathrm{P} $$
$$ \frac{\mathrm{d}\left[\mathrm{S}\right]}{\mathrm{d}\mathrm{t}}=\frac{\mathrm{d}\left[\mathrm{P}\right]}{\mathrm{d}\mathrm{t}}=-{{\mathrm{k}}_{\mathrm{on}}}^{\mathrm{S}}\left[\mathrm{S}\right]+{{\mathrm{k}}_{\mathrm{off}}}^{\mathrm{S}}\left[\mathrm{S}\mathrm{P}\right];\frac{\mathrm{d}\left[\mathrm{S}\mathrm{P}\right]}{\mathrm{d}\mathrm{t}}=-{{\mathrm{k}}_{\mathrm{off}}}^{\mathrm{S}}\left[\mathrm{S}\mathrm{P}\right] $$and
$$ \mathrm{R}+\mathrm{P}\overset{{{\mathrm{k}}_{\mathrm{on}}}^{\mathrm{R}}}{\to}\mathrm{R}\mathrm{P};\quad3em \mathrm{R}\mathrm{P}\overset{{{\mathrm{k}}_{\mathrm{off}}}^{\mathrm{R}}}{\to}\mathrm{R}+\mathrm{P} $$
$$ \frac{\mathrm{d}\left[\mathrm{R}\right]}{\mathrm{d}\mathrm{t}}=\frac{\mathrm{d}\left[\mathrm{P}\right]}{\mathrm{d}\mathrm{t}}=-{{\mathrm{k}}_{\mathrm{on}}}^{\mathrm{R}}\left[\mathrm{R}\right]+{{\mathrm{k}}_{\mathrm{off}}}^{\mathrm{R}}\left[\mathrm{R}\mathrm{P}\right];\frac{\mathrm{d}\left[\mathrm{R}\mathrm{P}\right]}{\mathrm{d}\mathrm{t}}=-{{\mathrm{k}}_{\mathrm{off}}}^{\mathrm{R}}\left[\mathrm{R}\mathrm{P}\right]. $$The rate equations of (68) and (69) cannot be solved analytically, even for either alone. In principle, however, at equilibrium one obtains Eqs. (2) and (3), respectively, for the individual associations; and the two are linked by Eq. (1). Moreover, in the contemplated event of having no peptide binding in the R state (konR = 0), [RP] = 0 at equilibrium. Nevertheless, previously bound peptide P dissociates from RP according to Eq. (69), whereupon re-association must occur via Eq. (68) after equilibration to the S state. Thus, for cases such as WT DnaK at low [P] where KeqS favors the R-state, KDApp(ATP) may be dominated by koffR and konS even when konR → 0 and KD0R →∞.
Considerations on the effect of conformational equilibria apply to reaction kinetics as they do to binding equilibria (38)–(40); however, this is only so for the bimolecular association process and not for the pseudo-first order dissociation process. Moreover, the situation becomes extra complicated if peptide binding occurs to R-state as well as to S-state Hsp70. For the case of binding to a single state, we have
$$ {{\mathrm{K}}_{\mathrm{D}}}^{\mathrm{App}}={\mathrm{k}}_{\mathrm{off}}/{{\mathrm{k}}_{\mathrm{on}}}^{\mathrm{App}}. $$Then with reference to Eq. (40), the case for negligible peptide binding to the R state, and to Eq. (41) for corresponding S-state association, we obtain
$$ {{\mathrm{k}}_{\mathrm{on}}}^{\mathrm{App}}\left(\mathrm{ATP}\right)={{\mathrm{k}}_{\mathrm{on}}}^{\mathrm{S}}\left(\mathrm{ATP}\right)/\left({\mathrm{K}}_{\mathrm{eqS}}+1\right). $$The situation is not so simple in the general case where peptides bind both to R and S states since most experiments will not discriminate.
Kinetics of substrate peptide release
Having to consider substrate peptide binding to both the R and S states also complicates general considerations on peptide disassociation; however, analyses can be made in certain circumstances.
One particular experiment of interest concerns the measurement of koff for peptide release from DnaK in the presence of ATP. Typically, one incubates DnaK with a labelled peptide, P*, in the absence of nucleotides or in ADP, and then measures the release of labelled peptide after mixing with a solution containing ATP and excess unlabelled peptide at concentration [P]. In general, the chaperone will be in equilibrium between S and R states, from which release will occur according to Eqs. (72) and (73):
$$ {\mathrm{S}\mathrm{P}}^{\ast}\left(\mathrm{t}\right)\overset{{{\mathrm{k}}_{\mathrm{off}}}^{\mathrm{S}}}{\to}\mathrm{S}+{\mathrm{P}}^{\ast}\left[{\mathrm{S}\mathrm{P}}^{\ast}\left(\mathrm{t}\right)\right]={\left[{\mathrm{S}\mathrm{P}}^{\ast}\right]}_0\exp \left(-{{\mathrm{k}}_{\mathrm{off}}}^{\mathrm{S}}\mathrm{t}\right) $$
$$ {\mathrm{R}\mathrm{P}}^{\ast}\left(\mathrm{t}\right)\overset{{{\mathrm{k}}_{\mathrm{off}}}^{\mathrm{R}}}{\to}\mathrm{R}+{\mathrm{P}}^{\ast}\left[{\mathrm{R}\mathrm{P}}^{\ast}\left(\mathrm{t}\right)\right]={\left[{\mathrm{R}\mathrm{P}}^{\ast}\right]}_0\exp \left(-{{\mathrm{k}}_{\mathrm{off}}}^{\mathrm{R}}\mathrm{t}\right) $$Thereby, released labelled peptide accumulates as
$$ \left[{\mathrm{P}}^{\ast}\left(\mathrm{t}\right)\right]=\left(\left[{\mathrm{SP}}^{\ast}\left(\mathrm{t}\right)\right]-{\left[{\mathrm{SP}}^{\ast}\right]}_0\right)+\left(\left[{\mathrm{RP}}^{\ast}\left(\mathrm{t}\right)\right]-{\left[{\mathrm{RP}}^{\ast}\right]}_0\right) $$
$$ \left[{\mathrm{P}}^{\ast}\left(\mathrm{t}\right)\right]={\left[{\mathrm{S}\mathrm{P}}^{\ast}\right]}_0\left\{1-\exp \left(-{{\mathrm{k}}_{\mathrm{off}}}^{\mathrm{S}}\mathrm{t}\right)\right\}+{\left[{\mathrm{R}\mathrm{P}}^{\ast}\right]}_0\left\{1-\exp \left(-{{\mathrm{k}}_{\mathrm{off}}}^{\mathrm{R}}\mathrm{t}\right)\right\}. $$The equilibrium between the S and R states is governed by Eqs. (1)–(3), whereby the fraction in the stimulating state, Q, is given by Eq. (10) and the remainder, 1 − Q, is in the restraining state. At t = 0, cT*(0) = [SP*]0 + [RP*]0. As the dissociation proceeds, the products S and R also accumulate along with P*, but perhaps at very different rates. Nevertheless, since we assume that both S and SP equilibrate identically with respect to both R and RP, Q([P], t) is expected by Eq. (11) to be invariant with time:
$$ \left[{\mathrm{P}}^{\ast}\left(\mathrm{t}\right)\right]={{\mathrm{c}}_{\mathrm{T}}}^{\ast }(0)\left\{\mathrm{Q}\left(\left[\mathrm{P}\right]\right)\left(1-\exp \left(-{{\mathrm{k}}_{\mathrm{off}}}^{\mathrm{S}}\mathrm{t}\right)\right)+\left(1\hbox{--} \mathrm{Q}\left(\left[\mathrm{P}\right]\right)\right)\right(1-\exp \left(-{{\mathrm{k}}_{\mathrm{off}}}^{\mathrm{R}}\mathrm{t}\right)\Big\}. $$The sensitivity of Q([P]) to peptide concentration is, of course, governed by the overall peptide concentration, which the excess of unlabelled peptide can be assumed to dominate.
A complication with the experiment for peptide dissociation in the presence of ATP is that ATP hydrolysis by Eqs. (13) and (14) will deplete SP and RP components to yield Hsp70U-ADP, from which the labelled peptide P* will then dissociate by Eq. (4) with koffU. Taking these two steps into account in Eq. (75) yields Eq. (76):
$$ \left[\mathrm{P}\ast \left(\mathrm{t}\right)\right]={\left[\mathrm{SP}\ast \right]}_0\left\{1-\exp \left(-{\mathrm{k}}^{\prime }\ \mathrm{t}\right)\right\}\ \left\{1-\exp \left(-{{\mathrm{k}}_{\mathrm{off}}}^{\mathrm{S}}\mathrm{t}\right)\right\}+{\left[\mathrm{RP}\ast \right]}_0\left\{1-\exp \left(-{\mathrm{k}}^0\mathrm{t}\right)\right\}\ \left\{1-\exp \left(-{{\mathrm{k}}_{\mathrm{off}}}^{\mathrm{R}}\mathrm{t}\right)\right\}+\left\{{\left[\mathrm{SP}\ast \right]}_0\exp \left(-{\mathrm{k}}^{\prime }\ \mathrm{t}\right)+{\left[\mathrm{RP}\ast \right]}_0\exp \left(-{\mathrm{k}}^0\mathrm{t}\right)\ \right\}\ \left\{1-\exp \left(-{{\mathrm{k}}_{\mathrm{off}}}^{\mathrm{U}}\mathrm{t}\right)\right\}, $$and Eq. (64) follows on rearrangement. Notice that while ATP hydrolysis may instantaneously change the relative proportions in S and P states, these will re-equilibrate by (11) to Q([P]) with [P] typically dominated by unlabelled peptide in excess of the labelled P*.
$$ \quad-5em \left[\mathrm{P}\ast \left(\mathrm{t}\right)\right]={\mathrm{c}}_{\mathrm{T}}\ast (0)\ \{\mathrm{Q}\left(\left[\mathrm{P}\right]\right)\ [\left(1-\exp \left(-{\mathrm{k}}^{\prime }\ \mathrm{t}\right)\right)\ \left(1-\exp \left(-{{\mathrm{k}}_{\mathrm{off}}}^{\mathrm{S}}\mathrm{t}\right)\right)\quad-5.25em +\exp \left(-{\mathrm{k}}^{\prime }\ \mathrm{t}\right)\ \left(1-\exp \left(-{{\mathrm{k}}_{\mathrm{off}}}^{\mathrm{U}}\mathrm{t}\right)\right)]+\left(1\hbox{--} \mathrm{Q}\left(\left[\mathrm{P}\right]\right)\right)\ \left[\left(1-\exp \left(-{\mathrm{k}}^0\mathrm{t}\right)\right)\ \left(1-\exp \left(-{{\mathrm{k}}_{\mathrm{off}}}^{\mathrm{R}}\mathrm{t}\right)\right)+\exp \left(-{\mathrm{k}}^0\mathrm{t}\right)\ \left(1-\exp \left(-{{\mathrm{k}}_{\mathrm{off}}}^{\mathrm{U}}\mathrm{t}\right)\right)\right]\}. $$Special cases
One special case of interest is the fraction in the stimulating state at [P] = 0, which from Eq. (12) is
$$ \mathrm{Q}(0)=\frac{1}{\left({\mathrm{K}}_{\mathrm{eqS}}+1\right)}. $$Another special case of interest is the hydrolysis rate at [P] = 0, which from Eq. (17) is
$$ {\mathrm{k}}_{\mathrm{cat}}={\mathrm{k}}^0+\frac{1}{\left({\mathrm{K}}_{\mathrm{eqS}}+1\right)}\times \left(\mathrm{k}^{\prime}\hbox{--} {\mathrm{k}}^0\right). $$Also from Eq. (17), the hydrolysis rate for the fully stimulated case (KeqS = 0) at the limit of KeqS = ∞ is kcat = k′ and that at the fully restrained state (KeqS = ∞) is kcat = k0.
Correction to obtain free peptide concentration
Although the formulations given by Eqs. (14) and (16) depend on [P], the concentration of free peptide, one is only able to control the total concentration of the peptide under study, which at the outset before ATP hydrolysis is given by Eq. (80):
$$ {\mathrm{c}}_{\mathrm{P}}=\left[\mathrm{P}\right]+\left[\mathrm{SP}\right]. $$Then, upon substituting Eq. (9) into Eq. (80) and rearranging the factors, one obtains the quadratic Eq. (81):
$$ {\left[\mathrm{P}\right]}^2+\left\{\left({\mathrm{c}}_{\mathrm{T}}\hbox{--} {\mathrm{c}}_{\mathrm{P}}\right)+{{\mathrm{K}}_{\mathrm{D}}}^{0\mathrm{S}}\left({\mathrm{K}}_{\mathrm{eqS}}+1\right)\right\}\ \left[\mathrm{P}\right]-{{\mathrm{K}}_{\mathrm{D}}}^{0\mathrm{S}}\left({\mathrm{K}}_{\mathrm{eqS}}+1\right)\ {\mathrm{c}}_{\mathrm{P}}=0, $$which can be solved readily for the desired [P] as a function of the experimentally accessible cP. For parameters relevant to WT DnaK, [P] is at 99% of cP even at cP = 5 μM and the fraction increases as cP increases. For the parameters of mutants such as I483D, where KeqS = 0, the [P] fraction is reduced to 82% at cP = 5 μM but reaches 99% by 100 μM. In other words, for practical situations with Hsp70s, measures of cP are reasonably close to the free peptide concentration [P].
Experimental Validation
Our theoretical explication was devised to explain observations that we had made on allosteric phenomena in Hsp70 action, and we use such data to test the formulations. The theory also predicts a previously uncharacterized stimulating state conformation, and we designed constructs that have successfully captured this state.
ATP hydrolysis controlled by peptide binding
A driving motivation for our theoretical development came in explaining the effect of peptide binding on ATP hydrolysis. The rate of hydrolysis is observed to accelerate as a function of substrate peptide concentration and to be affected by certain mutations. In a separate report (Wang et al., Reference Wang, Liu, Liu and Hendrickson2020), we describe single-turnover kinetic measurements of ATP hydrolysis by WT Hsp70 DnaK from E. coli and by selected mutant variants. Here in Fig. 2, we reproduce the resulting data for WT and I483D DnaK together with fittings based on Eq. (21) with alternative interpretations in terms of allosteric parameters as given by Eqs. (18) and (20). The goodness-of-fit to the WT data is excellent (1.02), giving mathematical parameters a′ = 0.276 ± 0.012 min−1, b′ = 1.33 ± 0.14 min−1 μM, and d′ = 115.1 ± 9.0 min−1 μM. Since our single-turnover kcat measurements are highly accurate, these fitting results provide a stringent test of the theory.
Effect of substrate peptide binding on ATP hydrolysis by Hsp70 DnaK. Observed rates of hydrolysis kcat and standard deviations, as reported elsewhere (Wang et al., Reference Wang, Liu, Liu and Hendrickson2020), are plotted as a function of the concentration [P] of NR heptapeptide (sequence NRLLLTG) for WT DnaK (●) and for I483D DnaK (x), which is characterized as fully stimulated constitutively. The smooth curve through WT DnaK data is from the least-squares fitting of measured rates by Eq. (21), which gave a′ = 0.276 min−1, b′ = ATP 1.33 min−1 μM and d′ = 115.1 μM. The straight line through points for I483D is at kcat = k′ = a′, which is the asymptote for the curve fitted to the WT DnaK data. The hydrolysis rates were measured, as reported Wang et al. (Reference Wang, Liu, Liu and Hendrickson2020), in assays of single-turnover kinetics (Schrank et al., Reference Schrank, Bolen and Hilser2009).
As discussed above, necessarily there is a degeneracy in biochemical parameters of the allosteric model since there are four of these in the formulation of Eq. (20), and five in the formulation of Eqs. (18), (31) and (37) as compared to the three intrinsic mathematical parameters of Eq. (21). Either independent experimental information or ad hoc assumptions or approximations are needed to break the degeneracy. We know, for example, that ATP hydrolysis is very slow in absence of peptide; so, if we assume this rate to be negligible, k0 = 0. We also formulated the model based on the premise that peptide affinity in the restraining state is very low; and if we assume it to be negligible, KD0R = ∞. And from our experimental evaluations (Wang et al., Reference Wang, Liu, Liu and Hendrickson2020), we have deduced that certain mutant variants are in defined states wherein particular biochemical parameters for that mutant should also reflect the WT value. For example, I483D and N170D both appear to be fixed in the stimulating state and these mutations are at sites that would not be expected to affect peptide binding; thus, we might assume KD0S (WT) = KDApp (I483D) = KDApp (N170D). Similarly, although N170D does affect the rate of ATP hydrolysis, I483D would not be expected to do so since this residue is exposed on SBD remote from the catalytic center; thus, by the hydrolysis model of Eq. (15) we might then assume k′ (WT) = k′ (I483D).
To break the degeneracy here, we extract allosteric parameters from the fitted WT mathematical parameters by first setting KD0S (WT) = < KDApp (N170D), KDApp (I483D) > = < 1.71 ± 0.26, 1.75 ± 0.16 > = 1.73 ± 0.20 μM. With this specification and assuming c = 1, it follows by (22) that KeqS = 65.5 ± 5.2. With c = 1, KD0R = ∞ by Eq. (23), k0 = 0.0075 ± 0.0020 min−1 by Eq. (27), and k′= a′ = 0.276 ± 0.012 min−1 by (28). This fitted value for k′ is within experimental error of the rate found constitutively for the S-state mutant I483D, k′ = 0.271 ± 0.011 min−1; moreover, the discrepancy is opposite from what could be closed by reducing KD0R in Eq. (28). For example, if c were increased to 1.1, by (23) KD0R would decrease only to 1,249 μM, which still implies negligible peptide affinity, while the k′ discrepancy would increase from 0.005 to 0.033 min−1 (3.0 σ). Indeed, it is fair to conclude that experiments confirm that model (20) is a valid simplification.
Additional hydrolysis experiments allow us to discriminate among the alternative hydrolysis models. DnaK609::NR was constructed with the optimized substrate peptide NRLLLTG fused to DnaK (WT-NR) in a manner disposed for avid binding to the SBD site, thus producing the SP state independent of extrinsic [P] (Wang et al., Reference Wang, Liu, Liu and Hendrickson2020). The rate of ATP hydrolysis measured for DnaK609::NR (kcat = 0.276 min−1) was essentially the same as for the constitutive S-state mutant I483D (kcat = 0.271 min−1), and both are indistinguishable from the WT hydrolytic value, k′ = 0.276 ± 0.012 min−1 (Fig. 2). The observation of k′(S) = k′(SP) is as expected from model (15) and inconsistent with model (29) when kS(S) is distinctly different from k′, such as being k0; however, these results are compatible with the Q-alternative model (33).
Peptide binding controlled by ATP binding and hydrolysis
In accord with many other observations, we find that peptide binding to DnaK is much reduced in the presence of ATP as compared to when with ADP (Wang et al., Reference Wang, Liu, Liu and Hendrickson2020) (Fig. 3). While this is true for WT DnaK, it is not so for mutants that by various biochemical criteria are fixed in the stimulating state, notably N170D and I483D (Wang et al., Reference Wang, Liu, Liu and Hendrickson2020). With these mutations, DnaK binds peptides with similar affinity whether in the presence of ATP or of ADP. Since the N170D mutation is in NBD and I483D is on an SBD surface remote from the peptide-binding site, we assert that these mutants reflect the intrinsic affinity of the site; that is I483D affinity is high because KeqS = 0 whereby KDApp = KD0S. The peptide binding affinity for these and other R-state interface mutants are all nearly the same as for WT when in ADP, KD(ADP) = 1.64 ± 0.08 μM and this is nearly the same as the intrinsic affinity with ATP, KD0S ≡ < KDApp (N170D, ATP), KDApp (I483D, ATP) > = 1.73 ± 0.20 μM.
Effects of ATP and ADP on substrate peptide binding by Hsp70 DnaK. The binding of fluorescein-labelled NR peptide (NRLLLTG, 10 nM) was measured by fluorescence anisotropy as a function of DnaK concentration (cT) as described elsewhere (Wang et al., Reference Wang, Liu, Liu and Hendrickson2020). Measurements are shown for WT DnaK in the presence of ADP (
) and in the presence of ATP (
) and for I483D DnaK in the presence of ATP (
). Peptide binding to I483D in ADP was indistinguishable from that to WT when in ADP.
Our experimental results on peptide binding to WT DnaK are compatible with the allosteric theory as formulated in Eq. (39), but not exactly as given by Eq. (40) for the case of KD0R = ∞. By Eq. (40), we calculate KDApp = 115.0 μM using KeqS = 65.5 from the fitting to hydrolysis data and KD0S = 1.73 μM as defined by binding to S-state mutants; whereas, we actually measure KDApp = 36.7 ± 5.2 μM for WT DnaK-ATP (Wang et al., Reference Wang, Liu, Liu and Hendrickson2020). By Eq. (41), a rearrangement of Eq. (39), we obtain KD0R = 53.0 μM from the measured values for KDApp, KeqS and KD0S; whereas, the fittings to hydrolysis data imply negligible binding to the R-state (KD0R > 1 mM). This dichotomy prompted us to consider the alternative model (42) whereby R is also in equilibrium with a conformation Q that has hydrolysis restrained as for R (33), but which is competent for peptide binding.
The Q-alternative peptide-binding model, Eq. (44), replaces KD0R with parameters KD0Q and KeqQ for the newly postulated conformation. Although Eq. (44) has too many variables for independent evaluation, we can obtain the controlling parameters from KDApp using Eqs. (45) and (46) after assuming values for KD0S and KD0Q. In this instance, with KDApp = 36.7 μM, we take KD0S = 1.73 μM as before and consider two options for KD0Q. At one plausible extreme, KD0Q = KD0S and at another, in light of a Q-like structure with R-like SBD-NBD interfaces and an S-like SBDβ bound to the NR peptide (Wang and Hendrickson, Reference Wang and Hendrickson2020b), we use KD0Q = 5.66 KD0S as observed comparing a lidless construct to WT DnaK (Buczynski et al., Reference Buczynski, Slepenkov, Sehorn and Witt S2001). We obtain KeqS/KeqQ = 2.14 and KeqS = 63.39 for KD0Q = 1.73 μM and KeqS/KeqQ = 12.08 and KeqS = 53.45 for KD0Q = 9.79 μM. The distributions among S, Q and R states follow from Eqs. (49)–(55). The fraction in the S-state stays the same as for the R/S model of Eq. (15), which in absence of substrate peptide is QS(0) = 1.5%; whereas, the predominating remainder is apportioned differently depending on the Q-state affinity: for KD0Q = 1.73 μM, QQ(0) = 3.2% and QR(0) = 95.3% while for KD0Q = 9.79 μM, QQ(0) = 18.2% and QR(0) = 80.3%.
Simulation of allosteric behaviour
The theory for allosteric regulation of Hsp70 activity permits the possibility to simulate the Hsp70 behaviour under varied conditions. It is of particular interest to consider the impact on chaperone properties of variation in the R–S equilibrium, which is governed by KeqS of Eq. (1). The effects of such variation on ATP hydrolysis are shown in Fig. 4a and the effects on peptide binding are shown in Fig. 4b. The respective families of curves for varied fractions in the stimulating state as determined by KeqS can be compared with experimental determinations, including those for WT and I483D DnaK given in Figs. 2 and 3, respectively. WT DnaK is dominantly in restraining state R whereas I483D is an extreme mutant fixed in the stimulating state S; other mutants are intermediate. For the I160D mutant, the fitting to peptide-binding data gave KDApp = 2.9 μM and fitting by Eq. (20) to the hydrolysis data gave KeqS = 4.75 (17% S-state at [P] = 0); however, this implies KDApp = 9.9 μM by Eq. (40). Exact fitting to the Q-alternative model by Eqs. (45) and (46) assuming KD0Q = KD0S is accomplished with a distribution of states QR(0):QQ(0):QS(0) of 40.3%:42.3%:17.4%. By Eq. (47) a physical solution in this case requires that KD0Q ≤ 1.95 KD0S (KD0Q ≤ 3.38 μM), whereby QR(0):QQ(0):QS(0) of 0%:82.6%:17.4%.
Simulations of effects of varied R vs. S distributions on ATP hydrolysis and substrate peptide binding by Hsp70 DnaK. (a) ATP hydrolysis by DnaK as a function of NR substrate peptide concentration. Simulations are by Eq. (20). We assume hydrolytic rate parameters k′ = 0.276 min−1 and k0 = 0.0075 min−1 from the fitting in Fig. 1; we assume the intrinsic dissociation constant KD0S = 1.73 μM, the average for fully stimulating mutants I483D and N170D; and we derive the equilibrium constant from Eq. (12), KeqS = (1 − Q0)/Q0 where Q0 is the specified fraction in the stimulating state at [P] = 0, QS(0). Q0 = 1.5% for the data in Fig. 1a. (b) Peptide saturation as a function of DnaK concentration at fixed NR peptide concentration. Simulations are by Eq. (63), assuming [P] = 10 nM, KD0S = 1.73 μM as for a, and again obtaining KeqS from the specified QS(0) by Eq. (12).
Hsp70 structure
Our initial formulation of a theory on allosteric regulation of Hsp70 molecular chaperones was derived to account for the results of biochemical tests in yeast Ssa1 and E. coli DnaK of interface mutations based on the structure of yeast Hsp110 Sse1 as a prototype for Hsp70s (Liu and Hendrickson, Reference Liu and Hendrickson2007; Wang et al., Reference Wang, Liu, Liu and Hendrickson2020). Structures of hydrolysis-impaired T199A mutants of E. coli DnaK in complexes with ATP corroborated the conjecture that ATP-associated Hsp70s would resemble Hsp110 (Kityk et al., Reference Kityk, Kopp, Sinning and Mayer2012; Qi et al., Reference Qi, Sarbeng, Liu, Le, Xu, Xu, Yang, Wong, Vorvis, Hendrickson, Zhou and Liu2013). Moreover, in keeping with Eq. (19) where KD0R = ∞, the peptide-binding sites in these structures are deformed from those in SBD-peptide complexes (Zhu et al., Reference Zhu, Zhao, Burkholder, Gragerov, Ogata, Gottesman and Hendrickson1996), such as to preclude peptide binding (Kityk et al., Reference Kityk, Kopp, Sinning and Mayer2012; Qi et al., Reference Qi, Sarbeng, Liu, Le, Xu, Xu, Yang, Wong, Vorvis, Hendrickson, Zhou and Liu2013; Wang et al., Reference Wang, Liu, Liu and Hendrickson2020). In addition, in keeping with the low rate of ATP hydrolysis, k0, as fitted to this predominating restraining R state conformation in the absence of peptide substrates, we find from a series of NBD(ATP) structures that the R state has a portion of NDB, which we call the R-to-S switch segment, in a conformation that blocks hydrolysis (Wang and Hendrickson, Reference Wang and Hendrickson2020a) as implied by low k0 in fittings to hydrolysis data.
Perhaps the most important confirmation of the theory is the finding that the alternative S-state conformation, postulated to explain rebinding of peptide substrates and ATP hydrolysis, is found to exist as predicted (Wang et al., Reference Wang, Liu, Liu and Hendrickson2020). These newly discovered S-state structures have molecular features compatible with biochemical activities. NBD in the S-state has its R-to-S switch segment in a conformation permissive of elevated hydrolysis as implied by higher k′ in fittings to hydrolysis data. The NBD-linker construct adopts the S-state conformation when with ATP as expected from the idea that R restrains Hsp70 from hydrolyzing ATP, and when released from SBD interactions it reverts to a potentiated hydrolysis (Wang and Hendrickson, Reference Wang and Hendrickson2020a). SBD in the S-state is receptive to peptide binding and, as seen by the near equivalence of KD0S to peptide KD(ADP), the conformation of SBDβ in the S state is nearly identical to that in the uncoupled U state (Wang et al., Reference Wang, Liu, Liu and Hendrickson2020). SBD in the S-state has the SBDα lid domain flexibly linked in keeping with higher on/off kinetics for WT Hsp70 in ATP as compared to that in ADP or to mutant Hsp70 variants (Wang et al., Reference Wang, Liu, Liu and Hendrickson2020).
Finally, in keeping with the lack of allosteric coupling in the absence of ATP, an NMR analysis of DnaK(ADP) shows NBD and SBD flexibly linked (Bertelsen et al., Reference Bertelsen, Chang, Gestwicki and Zuiderweg2009), and the contacts between NBD and SBD in X-ray structures of Hsp70s in the presence of ADP (Chang et al., Reference Chang, Sun, Wang and Hsiao2008; Adell et al., Reference Adell, Calisto, Fita and Martinelli2018) or without nucleotide (Jiang et al., Reference Jiang, Prasad, Lafer and Sousa2005) appear to be unnatural and nonproductive interactions, for example. SBD-linker lattice contacts or disordered domains.
Discussion
Biochemical properties of structure-inspired mutations of interfaces between domains in Hsp70 DnaK prompted the hypothesis that the state of Hsp70 first recognized by analogy to our structure of Sse1-ATP (Liu and Hendrickson, Reference Liu and Hendrickson2007) and also seen in the structures of DnaKR-ATP (Kityk et al., Reference Kityk, Kopp, Sinning and Mayer2012; Qi et al., Reference Qi, Sarbeng, Liu, Le, Xu, Xu, Yang, Wong, Vorvis, Hendrickson, Zhou and Liu2013) is restrained against its hydrolysis of ATP. The theory that we devised to explain such biochemical behaviour provides a sound basis for understanding allostery in Hsp70s. This theory is reminiscent of the famous MWC allosteric equilibrium model developed to explain oxygen binding by hemoglobin (Monod et al., Reference Monod, Wyman and Changeux1965), and it builds from decades of studies of allostery (Monod et al., Reference Monod, Wyman and Changeux1965; Koshland et al., Reference Koshland, Némethy and Filmer1966; Cui and Karplus, Reference Cui and Karplus2008; Motlagh et al., Reference Motlagh, Wrabl, Li and Hilser2014; Cuendet et al., Reference Cuendet, Weinstein and LeVine2016; Thirumalai et al., Reference Thirumalai, Hyeon, Zhuravlev and Lorimer2019); however, whereas MWC controls one binding activity in an oligomer through alternative quaternary states, here two activities are controlled reciprocally through alternative conformations adopted between domains of a single chain. We postulate an allosteric equilibrium between two ATP states: a restraining state with negligible affinity for polypeptide substrates and very limited ATP hydrolysis, and a stimulating state that hydrolyzes ATP readily and binds substrate peptides with rapid exchange kinetics. In the absence of peptide substrates, the restraining state dominates in the equilibrium, and the apparent peptide affinity is much reduced from its intrinsic value. When substrates are present, the equilibrium is drawn to the stimulating state, enhancing ATP hydrolysis and capturing valid substrates in the ADP state. The DnaKR-ATP structures epitomize the restraining state (Kityk et al., Reference Kityk, Kopp, Sinning and Mayer2012; Qi et al., Reference Qi, Sarbeng, Liu, Le, Xu, Xu, Yang, Wong, Vorvis, Hendrickson, Zhou and Liu2013) and our new structures of DnaKS-ATP model depict the stimulating state (Wang et al., Reference Wang, Liu, Liu and Hendrickson2020).
To be tractable, the theoretical model for peptide-stimulated hydrolysis of ATP needed to be simple. While our model seems to capture the essence of allosteric control in DnaK quite well, the reality may be more complex. For example, whereas we assume that rates of hydrolysis are the same by SP and S, stimulating-state Hsp70 with and without bound peptide, these rates likely differ somewhat. Moreover, the biochemical complexity of the system forced us to consider models, as for ATP hydrolysis experiments, that have more parameters than the resulting data can define unambiguously. In order to break the consequent degeneracy of parameters for such experiments, we have used observations from other experiments to define certain parameters, for example, intrinsic peptide affinity; however, underlying assumptions of equivalence may not hold perfectly.
Our studies have employed cellular and biochemical analyses of mutant variants to test functional hypotheses; however, mutated proteins are imperfect reporters of native function since intended perturbations of activity may extend to unanticipated effects. For example, whereas certain other mutations such as I483D seem to preclude the restraining state and to give a valid picture of stimulating state properties, there may be unanticipated consequences as well. Still, despite shortcomings, such mutation analyses have generated new insights and provided critical tests of Hsp70 function.
The treatment here for allosteric regulation is novel, to the best of my knowledge, both for Hsp70s in particular and for allosteric systems more generally. The result that most closely approaches our Hsp70-ATP equilibrium model (Eqs. 1–3) came from an NMR study (Zhuravleva et al., Reference Zhuravleva, Clerico and Gierasch2012) showing that DnaK-ATP can have two alternative conformations: an “ATP-bound, domain-docked state” modeled on Sse1-ATP (Liu and Hendrickson, Reference Liu and Hendrickson2007), now known to be very similar to restraining-state DnaKR-ATP (Kityk et al., Reference Kityk, Kopp, Sinning and Mayer2012; Qi et al., Reference Qi, Sarbeng, Liu, Le, Xu, Xu, Yang, Wong, Vorvis, Hendrickson, Zhou and Liu2013, Wang et al., Reference Wang, Liu, Liu and Hendrickson2020), and an “allosterically active state” bound to both ATP and substrate peptide. With clarifications from additional experiments (Lai et al., Reference Lai, Clerico, Blackburn, Patel, Robinson, Borbat, Freed and Gierasch2017), the latter likely relates to our stimulating-state DnaKS-ATP; however, this “partially docked state” is insufficiently specified for direct comparison with our crystal structure results (Wang et al., Reference Wang, Liu, Liu and Hendrickson2020). For allosteric systems more generally, the simplified case of Eq. (40) on binding, where KD0R = KD0Q = ∞, has been described before in the context of studies on an adenylate kinase (Schrank et al., Reference Schrank, Bolen and Hilser2009); however, we are not aware of previous formulations comparable to Eqs. (17)–(20), (36) and (37) on hydrolysis or to Eqs. (39) and (44) on binding for any system. An empirical fitting to DnaK hydrolysis data does have a form equivalent to Eq. (20) (Slepenkov and Witt, Reference Slepenkov and Witt S2002).
Methods
Fitting to conformational equilibrium theory
We developed a least-squares program for the fitting of parameters to our theoretical model of peptide-stimulated ATP hydrolysis by Hsp70 proteins. The observed kcat values are weighted by their inverse variances, w = 1/σ2, where σ = (σfit2 + σsys2)½ where σfit is the random error deduced from fitting to the kinetic data and σsys is a systematic error increment to account for added variations in repeated constant measurements, as for N170D, V389D and I483D at varying peptide concentrations. We find σsys = 0.0324 × kcat.
ATP hydrolysis and peptide binding measurements were made and analyzed as described in a companion paper (Wang et al., Reference Wang, Liu, Liu and Hendrickson2020).
Open Peer Review
To view the open peer review materials for this article, please visit http://doi.org/10.1017/qrd.2020.10.
Acknowledgements
I thank Wei Wang, Qun Liu, and Qinglian Liu for discussions about Hsp70 conformations, which were critical to the formulation and testing of the theory, and I thank Qinglian Liu for use of her data on ATP hydrolysis and substrate peptide binding by Hsp70 DnaK in experimental validations of the theory. This work was supported in part by NIH grant GM107462.
Conflict of Interest
The author declares no competing financial interests.
Open access
Comments
Comments to Author: I have no substantial criticism of the theory being advanced. However, the paper would be vastly improved if the authors first submitted the structure to the PDB. The absence of structural information is very strange indeed.