Efficient conversion of chemical energy into mechanical work by Hsp70 chaperones
Abstract
Hsp70 molecular chaperones are abundant ATP-dependent nanomachines that actively reshape non-native, misfolded proteins and assist a wide variety of essential cellular processes. Here, we combine complementary theoretical approaches to elucidate the structural and thermodynamic details of the chaperone-induced expansion of a substrate protein, with a particular emphasis on the critical role played by ATP hydrolysis. We first determine the conformational free-energy cost of the substrate expansion due to the binding of multiple chaperones using coarse-grained molecular simulations. We then exploit this result to implement a non-equilibrium rate model which estimates the degree of expansion as a function of the free energy provided by ATP hydrolysis. Our results are in quantitative agreement with recent single-molecule FRET experiments and highlight the stark non-equilibrium nature of the process, showing that Hsp70s are optimized to effectively convert chemical energy into mechanical work close to physiological conditions.
Introduction
Even though in vitro most proteins can reach their native structure spontaneously (Anfinsen, 1973), this is not always the case in cellular conditions and proteins can populate misfolded states which can form cytotoxic aggregates (Dobson, 2003). In order to counteract misfolding and aggregation, cells employ specialized proteins, called molecular chaperones, which act on non-native protein substrates by processes that stringently depend on ATP hydrolysis for most chaperone families (Hartl, 1996). Among them, the ubiquitous 70 kDa heat-shock proteins (Hsp70s) play a special role because they assist a plethora of fundamental cellular processes beyond prevention of aggregation (Clerico et al., 2019; Rosenzweig et al., 2019).
Decades of biochemical and structural studies have clarified the core elements of the Hsp70 functional cycle at the molecular level (Mayer, 2013). Hsp70s consist of two domains: the substrate binding domain (SBD) interacts with disparate substrate proteins, whereas the nucleotide binding domain (NBD) is responsible for the binding and hydrolysis of ATP. The two domains are allosterically coupled, and the nature of the nucleotide bound to the NBD affects the structure of the SBD and as a consequence the affinity for the substrate and its association/dissociation rates. More precisely, when the chaperone is in the ATP-bound state, the SBD is open and easily accessible to the substrate, whereas the SBD is closed when ADP is bound. These structural differences result in substrate binding and unbinding rates when ATP is bound that are orders of magnitude faster than when ADP is bound (Mayer et al., 2000). Furthermore, the coupling is bidirectional: the substrate, together with a co-localized J-domain protein (JDP) that serves as cochaperone (Kampinga and Craig, 2010; Kampinga et al., 2019), greatly accelerates the hydrolysis of ATP. Substrate binding thus benefits from the fast association rate of the ATP-bound state and the slow dissociation rate of the ADP-bound state, resulting in a non-equilibrium affinity (ultra-affinity) that can be enhanced beyond the maximum limit allowed by thermodynamic equilibrium, namely the affinity of the ADP-bound state (De Los Rios and Barducci, 2014; Barducci and De Los Rios, 2015).
More recently, the consequences of Hsp70 binding on the conformational ensembles of its substrates have also been investigated. Several lines of evidence indicate that the binding of Hsp70s to a polypeptide induces its expansion. Biochemical assays revealed that binding of Hsp70 increases the sensitivity of misfolded Luciferase to proteolysis and decreases its propensity to bind Thioflavin-T, strongly suggesting a loss of compactness (Sharma et al., 2010). Nuclear Magnetic Resonance (NMR) measurements have shown that Hsp70s destabilize the tertiary structure of several different substrates (Lee et al., 2015; Sekhar et al., 2015). Moreover, a single-molecule study based on Förster resonance energy transfer (FRET) spectroscopy quantified the considerable expansion of unfolded rhodanese in native conditions upon binding of multiple Hsp70 chaperones. In particular, this study revealed that the expansion is stringently ATP-dependent, because upon ATP exhaustion the system relaxes to the expansion values observed in the absence of chaperones (Kellner et al., 2014).
Despite these advances in the characterization of Hsp70 functioning, the mechanistic understanding of how the energy of ATP hydrolysis is used to expand a substrate has lagged behind. Our goal here is precisely to fill this gap between the molecular and functional characterization of Hsp70. To this aim, we first explore the structural and energetic features of Hsp70-bound rhodanese using Molecular Dynamics (MD) simulations. We next integrate this molecular information into a rate model that explicitly includes the Hsp70-rhodanese interactions and the chaperone ATPase cycle, thus elucidating how Hsp70s convert the chemical energy of ATP into mechanical work necessary to expand their substrates.
Results
Structural and thermodynamic characterization of chaperone-substrate complexes
To characterize the main features of chaperone-induced expansion, we performed MD simulations of the Hsp70/rhodanese complexes. We relied on a one-bead-per-residue Coarse Grained (CG) force field (Smith et al., 2014), which has been tailored to match experimental FRET data of intrinsically disordered proteins and satisfactorily reproduces the compactness of unfolded rhodanese in native conditions without any further tuning (see Materials and methods). Hsp70 chaperones were modeled with a structure-based potential to account for their excluded volume and they were artificially restrained to binding sites on the substrate. We identified six binding sites on the rhodanese sequence using two distinct bioinformatic algorithms (Rüdiger et al., 1997; Van Durme et al., 2009). Considering that each binding site could be either free or bound to a Hsp70 protein, we thus took into account a total of 26 = 64 distinct chaperone/substrate complexes, which were exhaustively simulated. In Figure 1, we report the distributions of the substrate potential energy and of the radius of gyration () for three representative complexes with one (left), three (center) and six (right) bound chaperones. Consistently with FRET results (Kellner et al., 2014), chaperone binding leads to larger radii of gyration and higher potential energies, implying that the excluded-volume interactions due to the large Hsp70s progressively expand the complex and disrupt the attractive intra-chain interactions in rhodanese.
We then calculated the conformational free energy of all the possible chaperone/rhodanese complexes to obtain a quantitative picture of the energy landscape governing the chaperone-induced expansion. To this aim, we performed extensive sets of non-equilibrium steering MD trajectories for each complex, and measured the work needed to steer it to a completely extended reference structure ( Å), whose conformational free energy is not affected by chaperone binding. Equilibrium free-energy differences with respect to this reference state were then estimated from non-equilibrium work distributions via the Jarzynski equality (Jarzynski, 1997), thus allowing the determination of the conformational free energy of each distinct chaperone/substrate complex (see Figure 2—figure supplement 1 and Materials and methods).
In Figure 2 (main), we report for each complex as a function of its mean radius of gyration using different colors for different stoichiometries. The conformational free energy increased with the swelling of the substrate due to the progressive binding of the chaperones. The increase in substrate potential energy due to the loss of intra-chain interactions upon Hsp70 binding is therefore only marginally compensated by the gain in conformational entropy. Notably, the conformational free energy is not uniquely determined by the stoichiometry, and is significantly affected by the specific binding pattern. The conformational free-energy cost of adding a single chaperone (inset in Figure 2) is positive for all complexes, but it varies from 2 kcal/mol up to 7 kcal/mol depending on the stoichiometry of the complex and on the particular choice of the binding sites. The increase of as a function of is quantitatively captured by Sanchez (1979) theory for the coil-to-globule collapse transition in polymers (see Figure 2 and Materials and methods). Remarkably, the excellent agreement is not the outcome of a fitting procedure since all the parameters were extracted from experiments (see Appendix 2). This result further reinforces the reliability of our simulations as well as the general applicability of the present setup beyond the particular system considered in this work.
ATP hydrolysis promotes multiple chaperone binding
The structural and thermodynamic characterization obtained by molecular simulations can be profitably complemented by a kinetic model encompassing relevant biochemical processes in order to determine the probability of each chaperone/substrate complex as a function of the chemical conditions. Notably, a model of the Hsp70 biochemical cycle based on experimental rates was previously used to illustrate how ATP-hydrolysis may result into non-equilibrium ultra-affinity for peptide substrates (De Los Rios and Barducci, 2014). Here, we extend this result to the more complex case of Hsp70-induced expansion by taking into account multiple chaperone binding events and their consequences on the conformational free energy of the substrate.
In our model, each state corresponds to a single configuration of the chaperone/substrate complex, which is defined by the occupation state of the six Hsp70 binding sites on rhodanese. Each site can be either free or occupied by an ADP- or ATP-bound chaperone for a total of 36 = 729 different states. All the relevant molecular processes corresponding to transitions between these states are explicitly modeled, including chaperone binding/unbinding, nucleotide exchange and ATP hydrolysis (see Figure 3). We took advantage of available biochemical data for determining the rate constants associated to all the relevant reactions (see Materials and methods). Importantly, kinetic rates for Hsp70 binding were modulated by the conformational free energies determined by CG MD simulations. Indeed, the unbinding rates of Hsp70 from large-sized protein substrates were observed to be similar to the ones from small peptides, whereas the binding rates can be up to two orders of magnitude smaller (Schmid et al., 1994; Mayer et al., 2000; Kellner et al., 2014). This evidence was further corroborated by a recent NMR study (Sekhar et al., 2018) suggesting a conformational selection scenario where the energetic cost due to substrate expansion mostly affects the Hsp70/rhodanese binding rate. Accordingly, we thus considered a substrate-independent unbinding rate constant , while we expressed the binding rate constant as
where , is the Boltzmann constant, T is the absolute temperature, is the binding rate measured for a peptide substrate, and is the conformational free-energy cost of Hsp70 binding, which depends on the specific initial and final binding patterns and in the rhodanese/chaperone complex (see Figure 2, inset). The interactions with JDP cochaperones were not explicitly modeled but the cochaperones were assumed to be colocalized with the substrate, so that their effect was implicitly taken into account in the choice of the rate constants for ATP hydrolysis (Kampinga and Craig, 2010; Hu et al., 2006).
The analytical solution of the model provides the steady-state probability of each binding configuration and allows the exploration of their dependence on external conditions. It is particularly instructive to investigate the system behavior as a function of the ratio between the concentrations of ATP and ADP, which is intimately connected to the energy released by ATP hydrolysis. At thermodynamic equilibrium, the ratio is greatly tilted in favor of ADP (; Alberty, 2005), whereas in the cell ATP is maintained in excess over ADP by energy-consuming chemostats (; Milo and Phillips, 2015). The ratio hence determines how far the system is from equilibrium, thus representing a natural control parameter for the non-equilibrium biochemical cycle. We thus report in Figure 4 (top panel) the compound probabilities for complexes with the same stoichiometry as a function of this nucleotide ratio. In conditions close to equilibrium (very low values of ), the vast majority of the substrate proteins are free and only about 10% of them are bound to a single chaperone. The population of equimolar complexes increases for between 10−2 and 10−1 and gives way to larger complexes with multiple chaperones for higher values of the nucleotide ratio. For , most substrates are bound to at least 4 chaperones, with an average stoichiometry (solid line in Figure 4), bottom panel). Further increase of the nucleotide ratio does not significantly change this scenario indicating an almost constant behaviour in large excess of ATP . It is important here to underscore that the binding of the chaperones in these conditions is a non-equilibrium effect, driven by the Hsp70-induced hydrolysis of ATP, and it is not a mere consequence of the excess of ATP over ADP or of the large difference between the substrate association rates to the ATP- and ADP-bound chaperones. Indeed, if we neglect Hsp70 ATPase activity () without changing any of the other model parameters, efficient chaperone binding is abolished (, as shown in the bottom panel of Figure 4, dashed line). As a matter of fact, in such equilibrium scenario an excess of ATP over ADP actually slightly disfavors chaperone binding, because the Hsp70 affinity for the substrate is slightly lower in the ATP-bound state than in the ADP-bound state (, see De Los Rios and Barducci, 2014 for further discussion.)
Combining the steady-state probabilities derived from the rate model with the results of the MD simulations, we can now exhaustively characterize the structural properties of the system. This provides the opportunity to directly compare our model with the results from FRET experiments both in equilibrium and non-equilibrium conditions. To this aim, we first focused on the average radius of gyration of the system at thermodynamic equilibrium () or in non-equilibrium conditions with ATP in large excess over ADP (). In order to probe the robustness of our results with respect to inaccuracies in the molecular model, we also took into account normally distributed errors on the conformational free energies .
The results are reported as histograms in Figure 5 and they suggest that at equilibrium the average radius of gyration is extremely close to what would be measured in the case of free substrate (dashed line). This is in agreement with the experimental observation that the formation of rhodanese–DnaK complexes is strictly dependent on the hydrolysis of ATP and that ADP cannot trigger the expansion of the substrate (Kellner et al., 2014). Conversely, in large excess of ATP we observe a substantial swelling of the substrate ( Å) due to the ultra-affine binding of Hsp70s. This finding is fully compatible with the size of DnaK/DnaJ/rhodanese complexes determined by single-molecule FRET experiments in excess of ATP (Kellner et al., 2014). In this regime, the limited effects of cochaperone binding on substrate conformations, which are not explicitly included in the model, play a minor role in determining the global expansion of the complex.
A more quantitative comparison between the model and the FRET results can be achieved by back-calculating the transfer efficiencies that were experimentally measured for five distinct pairs of fluorescent dyes (Kellner et al., 2014). In equilibrium conditions, namely when , the calculated FRET efficiency is ~0.8 for all considered pairs of fluorescent dyes (inset of Figure 5, blue circles) and it matches the experimental results for the compact unbound rhodanese (∼0.8). A dramatic difference is instead observed in excess of ATP (red circles), where the expansion of the substrate leads to a significant decrease of the calculated efficiency , in excellent agreement with the experimental values measured in similar conditions (black circles). Remarkably, the results correctly captured the non-monotonic behaviour of FRET efficiency as a function of the sequence separation between the dyes, which was not reproduced in previous calculations (Kellner et al., 2014). This agreement corroborates the prediction of the DnaK binding sites on the rhodanese sequence and the overall reliability of our model.
Energy balance and thermodynamic efficiency
Molecular chaperones consume energy via ATP hydrolysis in order to expand rhodanese. It is hence important to determine how effective they are as molecular machines, as well as to assess how favorable the physiological conditions are to perform their biological task.
To this aim, we calculated the global increase in the overall conformational free energy of the substrate with respect to equilibrium conditions, (Figure 6, top). This quantity measures the excess probability of each complex with respect to equilibrium conditions weighted by its corresponding conformational free energy .
where is the probability of complex for a given value of and is the same quantity computed at equilibrium conditions. In order to investigate the conversion of chemical energy into mechanical work, it is instructive to focus on the ratio between and the free energy of hydrolysis of ATP ,
The ratio reports on the effectiveness of the transduction process. We plot in Figure 6 (top) this quantity as a function of the ratio considering the estimated inaccuracies of the model as previously done for the gyration radius. Not surprisingly, all these curves exhibit a maximum because the probabilities of the different states, and thus also , attain plateaus for (see 4, top panel), whereas increases monotonically with the nucleotide ratio. The regime where transduction is maximally efficient intriguingly corresponds to values of that are typical of cellular conditions (grey area).
We highlight that in our model Hsp70 functioning encompasses two distinct yet intertwined processes: the ATP-dependent binding of the chaperones to the substrate, and its consequent expansion. In this two-step mechanism, the amount of energy available for the mechanical expansion is limited by that provided by non-equilibrium Hsp70 binding, which does not explicitly depend on the overall conformational properties of the substrate. To further dissect the energetic determinants of Hsp70 functioning and obtain more general conclusions, we thus analyzed the energy balance of chaperone binding to a model substrate, such as a peptide with a single binding site. To this aim, we focused on a simplified reaction cycle, which essentially corresponds to a single triangle within the overall scheme in Figure 3 and does not imply any conformational free-energy variation upon Hsp70 binding. We report in Figure 6 (bottom panel, black solid line) the non-equilibrium dissociation constant, , normalized with respect to its equilibrium value , as a function of [ATP]/[ADP]. When the ratio between the concentrations of ATP and ADP approaches the physiological regime, the dissociation constant drops significantly until it settles at a value that is two orders of magnitude lower than its equilibrium counterpart, as already discussed in De Los Rios and Barducci (2014). Here, we convert the dissociation constant into a binding free-energy excess with respect to equilibrium
that we can compare to the free energy of ATP hydrolysis, , as previously done in the case of . Interestingly, also in this case the energy ratio is maximal in cellular conditions (red dashed line in Figure 6, bottom panel), suggesting that the optimality of the overall expansion process does not depend on specific features of the substrate but it is a direct consequence of the intrinsic kinetic parameters of Hsp70 chaperones.
Discussion
Integrating molecular simulations, polymer theory, single-molecule experimental data and non-equilibrium rate models, we have developed a comprehensive framework that provides a quantitative picture of Hsp70-induced expansion of substrate proteins and offers a broad insight into the cellular functioning of this versatile chaperone machine.
We relied on molecular simulations for characterizing the structural and thermodynamic features of the complexes formed by the bacterial chaperone DnaK and its unfolded substrate rhodanese. Notably, we investigated a large variety of possible chaperone-substrate complexes for determining their conformational free energy as a function of stoichiometry and chaperone binding patterns. This computational strategy based on an enhanced-sampling protocol confirmed that excluded-volume interactions upon chaperone binding can greatly perturb the conformational ensemble of the unfolded substrate leading to its expansion. Remarkably, simulation results were found to be in excellent agreement with the predictions of Sanchez theory for globule to coil transition, thus providing another example of how polymer theory can be successfully used to decipher the behaviour of disordered proteins (Sherman and Haran, 2006; Hofmann et al., 2012; Schuler et al., 2016). We then combined conformational free energies with available biochemical data to develop an analytical rate model of the chaperone/substrate reaction cycle, which included both chaperone binding/unbinding and nucleotide hydrolysis/exchange processes.
This model fully takes into account non-equilibrium effects due to ATP hydrolysis and represents a natural extension of the ultra-affinity framework originally developed for peptide substrates with a single Hsp70 binding site (De Los Rios and Barducci, 2014). We could thus investigate the population of each complex and the average structural properties of the system as a function of the ATP/ADP nucleotide ratio, which measures how far the system is from thermodynamic equilibrium. The reliability of the model was corroborated by a quantitative comparison with recent single-molecule FRET data, indicating that our non-equilibrium framework accurately captures the salient features of the ATP-dependent expansion. We then used this unprecedented access to the thermodynamics details of this complex molecular process to compare the free-energy cost associated with substrate swelling with the chemical energy released by ATP-hydrolysis. Remarkably, this analysis revealed that energy transduction is maximally efficient for ATP/ADP values in cellular conditions. This result hints at the possibility that Hsp70 chaperones have been tuned by evolution to optimize the conversion of chemical energy into mechanical work for substrate expansion. Further analysis indicated that this optimality is likely inherited from the intrinsic properties of Hsp70 chaperones, which can convert up to 20% of the ATP chemical energy into non-equilibrium, excess binding energy at physiological conditions (Figure 6, bottom panel).
From a broader perspective, the ATP-driven action of Hsp70s induces a non-equilibrium redistribution of their protein-substrates over their structural ensemble. In particular, thanks to the fine-tuning of the process by co-chaperones (J-domain proteins and Nucleotide Exchange Factors), the expansion process highlighted here, followed by substrate release, may result in the enhancement of the native state population beyond the predictions of thermodynamic equilibrium, as recently observed even under otherwise denaturing conditions (Goloubinoff et al., 2018). Consistently, Zhao and coworkers have recently observed that Hsp70 chaperones crucially contribute in vivo to the solubility and functionality of a sizeable fraction of the E. coli proteome that, in their absence, would instead spontaneously misfold and aggregate (Zhao et al., 2019). Remarkably, a similar effect has been observed in vitro for the GroEL chaperonin (Chakrabarti et al., 2017; Goloubinoff et al., 2018), hinting at the possibility that multiple chaperone families might reshape the equilibrium conformational distribution of proteins through energy-consuming processes. These considerations might have important consequences for our ability to translate results from in vitro experiments to the active cellular context (Bershtein et al., 2013). Likewise, they raise fundamental questions about the evolution of protein sequences: indeed, since chaperones are ubiquitous and very much conserved across the different kingdoms of life, their ability to favor native states might have partially relieved the selection pressure for strong equilibrium thermodynamic stability, thus allowing evolution to proceed faster and to be more tolerant for slightly destabilizing mutations, as suggested in Rutherford and Lindquist (1998); Tokuriki and Tawfik (2009).
Besides the unfolding of non-native substrates discussed in this work, Hsp70s are highly versatile machines that play a fundamental role in a variety of diverse cellular functions such as protein translocation, protein translation, and disassembly of protein complexes. All these processes share basic analogies from the mechanistic point of view: in all these cases, Hsp70 binding to flexible substrates in constrained environments requires the energy of ATP hydrolysis (ultra-affinity) and results in the generation of effective forces due to excluded volume effects (entropic pulling), which ultimately drive protein translocation into mitochondria (De Los Rios et al., 2006; Assenza et al., 2015), clathrin cage disassembly (Sousa et al., 2016) and/or prevention of ribosome stalling (Liu et al., 2013). Here, by detailing how energy flows from ATP hydrolysis to mechanical work due to entropic pulling, we have elucidated a general force-generating mechanism of Hsp70 chaperones. This mechanism does not rely on any power-stroke conformational change but it rather depends on the efficient conversion of ATP chemical energy into ultra-affinity.
Materials and methods
Molecular model
Request a detailed protocolIn all the simulations, rhodanese and Hsp70 were coarse grained at the single-residue level as collections of beads centered on the atom of each amino acid. The unfolded state of bovine rhodanese (PDB:2RHS) was modeled according to the force field for disordered proteins from Smith et al. (2014). Two- and three-body bonded interactions along the substrate backbone were included via harmonic potentials, namely and , respectively. In the previous formulas, denotes bond lengths; the bend angles; Å; Å; rad; ; and is the thermal energy. Four-body bonded interactions were implemented as Fourier terms, , where is the torsion angle and , , , , , , , . The steric repulsion was implemented through a Weeks-Chandler-Andersen potential, , where
In the previous formula, is the distance between beads and , while Å. The hydrophobic part of the potential is specific to the interacting residues and is modeled as the attractive part of the Lennard-Jones potential, , where
In the previous formula, sets the overall strength of the hydrophobic interactions, while depends on the residues and involved in the interaction, and is defined as the geometric mean of their hydrophobicities, . The values of the hydrophobicities considered are based on a shifted and normalized Monera hydrophobicity scale (Smith et al., 2014). Electrostatic interactions were neglected based on control FRET experiments (see section 5.3). Without further tuning, this force field gives a radius of gyration of unbound rhodanese equal to Å, which is in good agreement with the experimental value Å (Kellner et al., 2014; Hofmann et al., 2014).
We modeled Hsp70 by means of a simple structure-based potential (Assenza et al., 2015) built on the conformation of ADP-bound Hsp70 (Bertelsen et al., 2009, PDB:2KHO). We described both the NBD (residues 4–680) and the SBD (residues 690–603) as rigid bodies whereas the interdomain linker (residues 681–689) was modeled according to the potential for flexible proteins described above. Importantly, non-bonded interactions of Hsp70 residues were limited to excluded volume effects and described by WCA potential (see Equation (5)). Electrostatic interactions were not explicitly included due to their marginal role in Hsp70/rhodanese complexes evidenced by FRET experiments (Appendix 1).
The binding sites for DnaK on the substrate were identified by applying the algorithms by Rüdiger et al. (1997) and Van Durme et al., 2009 on rhodanese and selecting only fragments for which at least partial consensus between the two predictions was obtained. Following this procedure, we identified six binding sites roughly centered on residues 10, 118, 131, 162, 188, 260 of the rhodanese sequence. The residues of the binding site were aligned to a SBD-bound peptide reported in the literature (PDB:1DKX Zhu et al., 1996) and constrained to move rigidly with the corresponding SBD, thus ensuring that each chaperone was irreversibly bound to the substrate. Following this procedure, 26 = 64 different chaperone/substrate complexes were built depending on the occupancy of each binding site.
Simulation protocols
Request a detailed protocolAll the simulations were performed with a version of LAMMPS (Plimpton, 1995) patched with the open-source, community-developed PLUMED library (Bonomi et al., 2019), version 2.1 (Tribello et al., 2014). The temperature was controlled through a Langevin thermostat with damping parameter 16 ns-1. The time step was set equal to 1 fs, and each residue had a mass equal to 1 Da.
In order to obtain conformational properties, for each of the 64 chaperone/rhodanese complexes we performed at least 10 independent simulations of 2 · 107 timesteps. To ensure that full equilibration was achieved, only the last 107 timesteps of the obtained trajectories were considered for analysis. Statistical errors on the computed quantities were estimated as standard errors of the mean computed across independent realizations and are smaller than the size of symbols reported in the figures. The FRET efficiency E for a given couple of dyes was computed starting from the distance separating the corresponding amino acids as
where r0 = 54 Å, as in Kellner et al. (2014). For each realization, the time average of E was computed. The final values employed to compute the results reported in the inset of Figure 5 in the main text were obtained as the average between independent realizations.
The conformational free energies were computed by means of steered simulations, where for each complex rhodanese was pulled from equilibrium until an elongated conformation was obtained (Figure 2—figure supplement 1, top panel). Due to the large intermolecular distances, the effect of chaperones on the conformational properties of fully-stretched rhodanese is negligible, so that this state can be used as a reference to compute the free-energy differences between different chaperone/rhodanese complexes. The pulling was implemented by adding a harmonic potential acting on the radius of gyration of rhodanese. The equilibrium position of the harmonic trap was increased at a constant pulling speed v = 10-5Å/fs from the equilibrium value up to . For each chaperone/substrate complex, 100 independent pulling simulations were performed, starting from uncorrelated initial snapshots extracted from the equilibrium distribution. For each realization, the work performed by the bias potential during the steering process was measured. The free-energy difference between the equilibrium starting point and the reference state (corresponding to ) was then computed via the Jarzynski (1997) equality:
where denotes statistical average. The error on was estimated according to the bootstrap method. The quantity considered in the main text was finally computed as , where corresponds to the case of rhodanese alone (Figure 2—figure supplement 1, bottom panel). The uncertainty on was estimated by propagating the error bars on and is always smaller than the size of symbols. In order to enhance the robustness of the results, the final values reported in the main text were obtained as a further average over the values of within the range 260 Å ≤ ≤ 290 Å.
Rate model
Request a detailed protocolFor the kinetic model we consider a system in which each of the six binding sites can either be occupied by a chaperone in the ATP or ADP state, or it can be free, so that in total there are 36 = 729 possible configurations. The concentration of each state evolves in time according to a system of rate equations
where is the transition rate from state to state . The first term in the right hand side (r.h.s.) of Equation (9) represents the total flux of molecules from the other states toward state , while the second term in the r.h.s. of Equation (9) accounts for the flux of molecules from state to any other state. We focused on the steady-state, when the concentrations of the various states do not change over time, which is defined by
Here, we provide a list of the relevant reactions that must be taken into account and of their corresponding rates. Each configuration is labelled by means of six symbols: 0 for empty sites, for sites occupied by an ATP-bound chaperone and for sites occupied by an ADP-bound chaperone (e.g. , where the first, third, fifth and sixth Hsp70 binding sites are unoccupied, the second binding site is associated to a chaperone in the ATP-bound state while the fourth binding site is associated with a chaperone in the ADP-bound state). With this notation, the rates corresponding to every reaction are easily determined. Examples of the reactions that need to be considered are
binding/unbinding
hydrolysis/synthesis
nucleotide exchange
We further provide, as an example, the equation for a precise configuration, say (here the label stands for the concentration of the configuration). The two binding sites that are occupied can undergo chaperone unbinding, ATP hydrolysis/synthesis or nucleotide exchange. The remaining unoccupied binding sites can bind either an ATP- or an ADP-bound chaperone. We thus have
Below we further detail the rates of our model.
It is possible to move from an ATP-state to an ADP-state either via hydrolysis/synthesis or via nucleotide exchange. In the case of exchange, effective constants are used, which take into account the unbinding of one nucleotide species and the binding of the different one. The effective exchange rates are thus a function of the ratio [ATP]/[ADP] (see also De Los Rios and Barducci, 2014):
where , , and are the binding and unbinding rates for ADP and ATP respectively.
The rates of binding between the chaperone and single peptides have been previously determined experimentally (Mayer et al., 2000), and they were corrected in order to take into account the conformational change of the full polypeptide substrate upon binding, as we illustrated in the main text.
Substrate binding enhances the chaperone ATPase activity. Furthermore, the stimulation of ATP hydrolysis always takes place in cooperation with JDP co-chaperones. In our model, we did not consider them explicitly but their contribution was implicitly included through the choice of the rate constants.
In particular, the hydrolysis rate in the absence of the substrate, , is much smaller than the same rate in the presence of the substrate, (). We assumed that the ratio between the rate of hydrolysis and the rate of synthesis is not altered by the substrate:
The substrate binding/unbinding rates, the rates of nucleotide exchange and the hydrolysis and synthesis rates are collectively constrained by thermodynamic relations. Indeed, when the ratio between the concentrations of ATP and ADP is equal its equilibrium value (when the spontaneous hydrolysis and synthesis reactions are at steady state and compensate each other), detailed balance must be satisfied (Ge et al., 2012). As a consequence, for every closed cycle in the reaction network the product of the rates in one direction must be equal to the product of the rates in the opposite direction. Therefore, if , , and are the rate of substrate binding and unbinding from a chaperone in the ATP and ADP states, we must have
Remarkably, taking the rates as provided in Mayer et al. (2000); Hu et al. (2006); Kellner et al. (2014), this relation is not satisfied, and we had thus to modify them. We thus calculated the product in the formula above and then corrected the rates in the following way:
The concentration of free chaperones in the ATP and in the ADP states was obtained, at the leading order, by solving a three-state system whose reactions have the form
Since we worked in the assumption of excess of chaperones in the system, once these concentrations were obtained, they remained fixed once for all, without being considered as a variable of the biochemical network.
We report in the following table the rates used in the model.
Parameters of the model (Mayer et al., 2000; Hu et al., 2006; Kellner et al., 2014) | |
---|---|
To test the robustness of the model for the radius of gyration, the average FRET efficiency and the free-energy , 100 realizations were implemented, taking each time the values from a Gaussian distribution with kcal/mol and mean equal to the value obtained with the MD simulations.
Molecular graphics
Request a detailed protocolMolecular graphics in Figures 1 and 6 have been generated with UCSF Chimera, developed by the Resource for Biocomputing, Visualization, and Informatics at the University of California, San Francisco, with support from NIH P41-GM103311 (Pettersen et al., 2004).
Appendix 1
Limited role of electrostatic interactions in chaperone/substrate complexes
The choice of neglecting electrostatic interactions in the CG simulation of DnaK/rhodanese complexes was motivated by the results of two control experiments, which were performed exactly at the same conditions except for the salt concentration, namely 5 mM and 100 mM KCl. A double-cysteine variant of rhodanese (K135C/K174C), which was produced by site-directed mutagenesis and prepared as described before for wild-type rhodanese (Miller et al., 1992), was labeled with Alexa Fluor 488 C5 maleimide and Alexa Fluor 594 C5 maleimide (Invitrogen, Molecular Probes) (Hillger et al., 2007). A twofold molar excess of the dyes was added to the protein and incubated for 1 hr at room temperature. Unreacted dye was removed by gel filtration followed by anion exchange chromatography to reduce the amount of incorrectly labelled protein using a MonoQ 5/50 GL column (GE Healthcare) installed on an ÄKTA purifier FPLC system equilibrated in 50 mM Tris·HCl, pH 7.0, and eluted with a gradient from 0 to 500 mM sodium chloride over 60 mL (12 column volumes). The chaperone proteins DnaK and DnaJ (stock solution concentration 100 μM in 50 mM Tris HCl, pH 7.7, 100 mM NaCl) were gifts from H.-J. Schönfeld (Hoffmann-La Roche Ltd., Basel). Labelled rhodanese was denatured in 4 M guanidinium chloride in buffer (50 mM Tris HCl, 10 mM MgCl2, 200 mM β-mercaptoethanol, and 0.001% Tween 20) with either 5 or 100 mM KCl added. The denatured rhodanese was diluted 100x into buffer containing 10 μM DnaK, 0.5 μM DnaJ, 1 mM ATP and 5 or 100 mM KCl to form chaperone rhodanese complexes at a final concentration of 50 pM rhodanese. Single-molecule Förster Resonance Energy Transfer (FRET) measurements were started immediately after dilution and data recorded at 22°C for 30 min to construct the FRET efficiency histograms. Data were recorded with a MicroTime 200 confocal microscope (PicoQuant) and on a custom-built confocal microscope. All measurements were obtained with pulsed interleaved excitation (Müller et al., 2005). The instrument set up and data reduction were the same as described before (Kellner et al., 2014). The collected FRET histograms are plotted in Appendix 1—figure 1. The average FRET efficiencies were 0.38 ± 0.01 with 5 mM KCl and 0.40 ± 0.01 with 100 mM KCl (errors indicate uncertainties estimated from the standard deviation of the larger data set collected at 5 mM salt). The very similar values obtained at different salt concentrations suggest that electrostatic interactions do not play a significant role in determining the conformational properties of rhodanese/chaperone complexes.
-
Appendix 1—figure 1—source data 1
- https://cdn.elifesciences.org/articles/48491/elife-48491-app1-fig1-data1-v2.zip
Appendix 2
Sanchez theory for coil-to-globule transition
In his theory on polymer coil-to-globule transition Sanchez (1979) considers a Freely-Jointed Chain (FJC) made of monomers characterized by attractive interactions of average magnitude . Let be the expansion parameter, where is the radius of gyration of the polymer, while is the value of in the unperturbed case ( and no excluded volume present). Assuming a Flory-Fisk distribution for the unperturbed case, the probability distribution of the expansion parameter was shown to be (Sanchez, 1979)
where is the partition function and
with . The free energy can then be straightforwardly computed as , that is making the dependence on explicit,
In the present case, rhodanese is constituted by 293 amino acids, that is there are bonds. The bond length is Å, estimated as the typical distance between atoms belonging to consecutive residues (Hofmann et al., 2012). In order to apply the Sanchez theory, one needs to consider the FJC equivalent to rhodanese. To this aim, the Kuhn length can be estimated as Å, where Å is the persistence length of an unfolded protein under native conditions (Hofmann et al., 2012). The number of monomers in the equivalent FJC can then be computed by imposing the total contour length of rhodanese: . From there, both and the unperturbed radius of gyration can be calculated, leaving as the only unknown quantity in Equation (20). The latter can be fixed by imposing that has a minimum at the experimental value Å obtained in the absence of chaperones (Kellner et al., 2014; Hofmann et al., 2014), yielding kcal/mol. Plugging the values of the parameters derived above into Equation (21) gives the black line reported in Figure 2 in the main text, in excellent agreement with the simulation results. We stress that neither the force field nor the MD results were employed in the derivation of the parameters used Equation (20).
Data availability
All the source data used for generating relevant figures (Figures 1, 2, 4, 5, 6, Figure 2—figure supplement 1, Appendix 1—figure 1) have been provided as supporting files. All the information necessary for reproducing the molecular simulations have been deposited in github (https://github.com/saassenza/Hsp70Unfoldase; also provided as Source data 1) and PLUMED NEST (plumID:19.076) repositories.
References
-
Quantifying the role of chaperones in protein translocation by computational modelingFrontiers in Molecular Biosciences 2:8.https://doi.org/10.3389/fmolb.2015.00008
-
Non-equilibrium conformational dynamics in the function of molecular chaperonesCurrent Opinion in Structural Biology 30:161–169.https://doi.org/10.1016/j.sbi.2015.02.008
-
Hsp70 molecular chaperones: multifunctional allosteric holding and unfolding machinesBiochemical Journal 476:1653–1677.https://doi.org/10.1042/BCJ20170380
-
Chaperones convert the energy from ATP into the nonequilibrium stabilization of native proteinsNature Chemical Biology 14:388–395.https://doi.org/10.1038/s41589-018-0013-8
-
Modeling Hsp70-mediated protein foldingBiophysical Journal 91:496–507.https://doi.org/10.1529/biophysj.106.083394
-
Nonequilibrium equality for free energy differencesPhysical Review Letters 78:2690–2693.https://doi.org/10.1103/PhysRevLett.78.2690
-
Function, evolution, and structure of J-domain proteinsCell Stress and Chaperones 24:7–15.https://doi.org/10.1007/s12192-018-0948-4
-
The HSP70 chaperone machinery: j proteins as drivers of functional specificityNature Reviews Molecular Cell Biology 11:579–592.https://doi.org/10.1038/nrm2941
-
Multistep mechanism of substrate binding determines chaperone activity of Hsp70Nature Structural Biology 7:586–593.https://doi.org/10.1038/76819
-
Hsp70 chaperone dynamics and molecular mechanismTrends in Biochemical Sciences 38:507–514.https://doi.org/10.1016/j.tibs.2013.08.001
-
Recombinant bovine rhodanese: purification and comparison with bovine liver rhodaneseBiochimica Et Biophysica Acta (BBA) - Protein Structure and Molecular Enzymology 1121:286–292.https://doi.org/10.1016/0167-4838(92)90158-A
-
Pulsed interleaved excitationBiophysical Journal 89:3508–3522.https://doi.org/10.1529/biophysj.105.064766
-
UCSF Chimera?A visualization system for exploratory research and analysisJournal of Computational Chemistry 25:1605–1612.https://doi.org/10.1002/jcc.20084
-
Fast parallel algorithms for Short-Range molecular dynamicsJournal of Computational Physics 117:1–19.https://doi.org/10.1006/jcph.1995.1039
-
The Hsp70 chaperone networkNature Reviews Molecular Cell Biology 20:665–680.https://doi.org/10.1038/s41580-019-0133-3
-
Phase transition behavior of the isolated polymer chainMacromolecules 12:980–988.https://doi.org/10.1021/ma60071a040
-
The kinetic parameters and energy cost of the Hsp70 chaperone as a polypeptide unfoldaseNature Chemical Biology 6:914–920.https://doi.org/10.1038/nchembio.455
-
Clathrin-coat disassembly illuminates the mechanisms of Hsp70 force generationNature Structural & Molecular Biology 23:821–829.https://doi.org/10.1038/nsmb.3272
-
PLUMED 2: new feathers for an old birdComputer Physics Communications 185:604–613.https://doi.org/10.1016/j.cpc.2013.09.018
-
Accurate prediction of DnaK-peptide binding via homology modelling and experimental dataPLOS Computational Biology 5:e1000475.https://doi.org/10.1371/journal.pcbi.1000475
Article and author information
Author details
Funding
Agence Nationale de la Recherche (ANR-14-ACHN-0016)
- Alessandro Barducci
Schweizerischer Nationalfonds zur Förderung der Wissenschaftlichen Forschung (200020_163042)
- Paolo De Los Rios
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Copyright
© 2019, Assenza et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 2,291
- views
-
- 328
- downloads
-
- 26
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Computational and Systems Biology
- Physics of Living Systems
Planar cell polarity (PCP) – tissue-scale alignment of the direction of asymmetric localization of proteins at the cell-cell interface – is essential for embryonic development and physiological functions. Abnormalities in PCP can result in developmental imperfections, including neural tube closure defects and misaligned hair follicles. Decoding the mechanisms responsible for PCP establishment and maintenance remains a fundamental open question. While the roles of various molecules – broadly classified into ‘global’ and ‘local’ modules – have been well-studied, their necessity and sufficiency in explaining PCP and connecting their perturbations to experimentally observed patterns have not been examined. Here, we develop a minimal model that captures the proposed features of PCP establishment – a global tissue-level gradient and local asymmetric distribution of protein complexes. The proposed model suggests that while polarity can emerge without a gradient, the gradient not only acts as a global cue but also increases the robustness of PCP against stochastic perturbations. We also recapitulated and quantified the experimentally observed features of swirling patterns and domineering non-autonomy, using only three free model parameters - rate of protein binding to membrane, the concentration of PCP proteins, and the gradient steepness. We explain how self-stabilizing asymmetric protein localizations in the presence of tissue-level gradient can lead to robust PCP patterns and reveal minimal design principles for a polarized system.
-
- Physics of Living Systems
Many proteins have been recently shown to undergo a process of phase separation that leads to the formation of biomolecular condensates. Intriguingly, it has been observed that some of these proteins form dense droplets of sizeable dimensions already below the critical concentration, which is the concentration at which phase separation occurs. To understand this phenomenon, which is not readily compatible with classical nucleation theory, we investigated the properties of the droplet size distributions as a function of protein concentration. We found that these distributions can be described by a scale-invariant log-normal function with an average that increases progressively as the concentration approaches the critical concentration from below. The results of this scaling analysis suggest the existence of a universal behaviour independent of the sequences and structures of the proteins undergoing phase separation. While we refrain from proposing a theoretical model here, we suggest that any model of protein phase separation should predict the scaling exponents that we reported here from the fitting of experimental measurements of droplet size distributions. Furthermore, based on these observations, we show that it is possible to use the scale invariance to estimate the critical concentration for protein phase separation.