Efficient conversion of chemical energy into mechanical work by Hsp70 chaperones

  1. Salvatore Assenza
  2. Alberto Stefano Sassi
  3. Ruth Kellner
  4. Benjamin Schuler
  5. Paolo De Los Rios
  6. Alessandro Barducci  Is a corresponding author
  1. ETH Zürich, Switzerland
  2. Universidad Autónoma de Madrid, Spain
  3. École Polytechnique Fédérale de Lausanne (EPFL), Switzerland
  4. Yorktown Heights, United States
  5. University of Zurich, Switzerland
  6. Ecole Polytechnique Fédérale de Lausanne (EPFL), Switzerland
  7. INSERM, CNRS, Université de Montpellier, France

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 (Rg) 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.

Progressive expansion of substrate by multiple Hsp70 binding.

Probability density maps of substrate potential energy and radius of gyration for representative Hsp70/rhodanese complexes with one (left), three (center) and six (right) bound chaperones. The different Hsp70 chaperones have been represented with different colors to ease their discernibility.

Figure 1—source data 1

Text files containing the data used to generate the probability density maps in Figure 1.

https://cdn.elifesciences.org/articles/48491/elife-48491-fig1-data1-v2.zip

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 (Rg>260 Å), 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 ΔG of each distinct chaperone/substrate complex (see Figure 2—figure supplement 1 and Materials and methods).

In Figure 2 (main), we report ΔG 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 ΔΔG 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 ΔG as a function of Rg 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.

Figure 2 with 1 supplement see all
Hsp70-induced swelling increases the conformational free energy of the system.

Conformational free-energy differences ΔG of the Hsp70/rhodanese complexes with respect to the unbound substrate (n = 0) plotted as a function of the corresponding radius of gyration Rg. Each point represents one of the 64 possible binding configurations with color code indicating the number of bound chaperones. The black curve was obtained using the model in Sanchez (1979) (see Appendix 2). (inset) Distribution of ΔΔG corresponding to the free-energy cost for binding an additional Hsp70 to a chaperone/substrate complex.

Figure 2—source data 1

Data from simulation results and Sanchez theory used to generate the plot in Figure 2 and the histogram in the inset of Figure 2.

https://cdn.elifesciences.org/articles/48491/elife-48491-fig2-data1-v2.zip

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 koff, while we expressed the binding rate constant as

(1) kon,ij=kon0exp[-βΔΔGij],

where β=1/kBTkB is the Boltzmann constant, T is the absolute temperature, kon0 is the binding rate measured for a peptide substrate, and ΔΔGij is the conformational free-energy cost of Hsp70 binding, which depends on the specific initial and final binding patterns i and j 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 relevant biochemical reactions of the Hsp70/rhodanese system can be described by a rate model.

Each chaperone binding site on rhodanese (black dots) can be either free or occupied by an Hsp70 (yellow), which in turn can be either ADP- or ATP-bound. We depict here for the sake of clarity only a representative portion of the full model, which takes into account six binding sites. The reaction cycle is governed by the rates for chaperone binding/unbinding to the substrate (konATP,konADP,koffATP,koffADP) and for hydrolysis (khs), synthesis (kss) and exchange (kex,DTeff,kex,TDeff) of nucleotides bound to the chaperones (see Materials and methods for further details). Importantly, the binding rate constants, konATP and konADP, take into account the conformational free energies, according to Equation (1).

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 [ATP]/[ADP] ratio is greatly tilted in favor of ADP ([ATP]eq/[ADP]eq10-9-10-8; Alberty, 2005), whereas in the cell ATP is maintained in excess over ADP by energy-consuming chemostats ([ATP]/[ADP]>1; Milo and Phillips, 2015). The [ATP]/[ADP] 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 n as a function of this nucleotide ratio. In conditions close to equilibrium (very low values of [ATP]/[ADP]), 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 [ATP]/[ADP] between 10−2 and 10−1 and gives way to larger complexes with multiple chaperones for higher values of the nucleotide ratio. For [ATP]/[ADP]>1, most substrates are bound to at least 4 chaperones, with an average stoichiometry n4.9 (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 ([ATP]/[ADP]>10). 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 (khh,ksh=0) without changing any of the other model parameters, efficient chaperone binding is abolished (n1, 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 (koffATP/konATP>koffADP/konADP, see De Los Rios and Barducci, 2014 for further discussion.)

Hsp70 binding is a non-equilibrium process that depends on ATP hydrolysis.

(Top) compound probabilities for Hsp70/substrate complexes with given number of bound chaperones n as a function of [ATP]/[ADP]. (Bottom) Mean value n as a function of [ATP]/[ADP] with (solid line) and without (dashed line) ATP hydrolysis.

Figure 4—source data 1

Data from rate model employed to generate the plots in Figure 4.

https://cdn.elifesciences.org/articles/48491/elife-48491-fig4-data1-v2.zip

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 ([ATP][ADP]) or in non-equilibrium conditions with ATP in large excess over ADP ([ATP]/[ADP]>10). 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 ΔGi.

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 (75<Rg<95 Å) 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.

Simulation results match sm FRET experimental data.

Histograms of the radius of gyration for equilibrium (blue) and non equilibrium (red) values of [ATP]. The black dashed line indicates the average radius of unbound rhodanese. (inset) FRET transfer efficiencies as a function of the sequence separation between the fluorescent dyes. The black circles correspond to the experimental values (Kellner et al., 2014). Calculated efficiencies taking into account uncertainties are reported as blue (equilibrium conditions) and red circles (ATP excess).

Figure 5—source data 1

Molecular simulation data used to generate the plots in Figure 5.

https://cdn.elifesciences.org/articles/48491/elife-48491-fig5-data1-v2.zip

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 [ATP]/[ADP]1, 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, ΔGSwell (Figure 6, top). This quantity measures the excess probability of each complex with respect to equilibrium conditions weighted by its corresponding conformational free energy ΔGi.

(2) ΔGSwell=i[pi([ATP][ADP])-pieq]ΔGi,

where pi([ATP][ADP]) is the probability of complex i for a given value of [ATP]/[ADP] and pieq 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 ΔGSwell and the free energy of hydrolysis of ATP ΔGh,

(3) ΔGh=kBT[ln([ATP][ADP])ln([ATP]eq[ADP]eq)].

The ratio ΔGSwell/ΔGh reports on the effectiveness of the transduction process. We plot in Figure 6 (top) this quantity as a function of the [ATP]/[ADP] 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 ΔGSwell, attain plateaus for [ATP][ADP] (see 4, top panel), whereas ΔGh increases monotonically with the nucleotide ratio. The regime where transduction is maximally efficient intriguingly corresponds to values of [ATP]/[ADP] that are typical of cellular conditions (grey area).

The conversion of chemical energy into mechanical work by Hsp70s is optimal in cellular conditions thanks to the chaperone intrinsic rates.

(Top) Ratio between the conformational free energy and the free energy of ATP hydrolysis, as a function of [ATP]/[ADP]. Dark green curve results from data from molecular simulations and light green curves takes into account normally-distributed uncertainties on calculated ΔGi. (Bottom) Effective dissociation constant in the case of a single binding site normalized with respect to the corresponding value in equilibrium, as a function of [ATP]/[ADP] (solid black line). Ratio between the non-equilibrium excess of binding free energy ΔGb and the free energy of ATP hydrolysis ΔGh, as a function of [ATP]/[ADP] (dashed red line). The gray region indicates the interval corresponding to physiological conditions.

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, Kdneq, normalized with respect to its equilibrium value Kdeq, 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

(4) ΔGb=-kBTln[KdneqKdeq]

that we can compare to the free energy of ATP hydrolysis, ΔGh, as previously done in the case of ΔGSwell. 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 protocol

In all the simulations, rhodanese and Hsp70 were coarse grained at the single-residue level as collections of beads centered on the Cα 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 Vbond=klb(rbl)2/2 and Vbend=12kθα(θαθ0)2, respectively. In the previous formulas, rb denotes bond lengths; θα the bend angles; l=3.9 Å; (kl/kBT)1/2=0.046 Å; θ0=2.12 rad; (kθ/kBT)-12=0.26; and kBT is the thermal energy. Four-body bonded interactions were implemented as Fourier terms, Vdihed=kBTds=14[Ascos(sϕd)+Bssin(sϕd)], where ϕd is the torsion angle and A1=0.705, A2=-0.313, A3=-0.079, A4=0.041, B1=-0.175, B2=-0.093, B3=0.030, B4=0.030. The steric repulsion was implemented through a Weeks-Chandler-Andersen potential, VWCA=ijVr, where

(5) Vr={4kBT[(σrij)12(σrij)6]+kBTif rij216σ0otherwise.

In the previous formula, rij is the distance between beads i and j, while σ=4.8 Å. The hydrophobic part of the potential is specific to the interacting residues and is modeled as the attractive part of the Lennard-Jones potential, Vhydro=ϵhijVh, where

(6) Vh={4ϵij[(σrij)12(σrij)6]if rij216σϵijotherwise.

In the previous formula, ϵh=0.7722kBT sets the overall strength of the hydrophobic interactions, while ϵij depends on the residues i and j involved in the interaction, and is defined as the geometric mean of their hydrophobicities, ϵijϵiϵj. 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 Rg=(23.3±0.1) Å, which is in good agreement with the experimental value Rg=(20.1±0.8) Å (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 protocol

All 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 T=293K 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 r separating the corresponding amino acids as

(7) E=11+(rr0)6,

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 Rg 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 Rg=Rgfin. 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 W performed by the bias potential during the steering process was measured. The free-energy difference δG between the equilibrium starting point and the reference state (corresponding to Rgfin) was then computed via the Jarzynski (1997) equality:

(8) e-δGkBT=e-WkBT,

where denotes statistical average. The error on δG was estimated according to the bootstrap method. The quantity ΔG considered in the main text was finally computed as ΔG=δG0-δG, where δG0 corresponds to the case of rhodanese alone (Figure 2—figure supplement 1, bottom panel). The uncertainty on ΔG was estimated by propagating the error bars on δG 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 Rgfin within the range 260 Å ≤ Rgfin ≤ 290 Å.

Rate model

Request a detailed protocol

For 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 ci of each state evolves in time according to a system of rate equations

(9) dcidt=jkjicj-jkijci

where kij is the transition rate from state i to state j. 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 i, while the second term in the r.h.s. of Equation (9) accounts for the flux of molecules from state i 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

(10) dcidt=0

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, T for sites occupied by an ATP-bound chaperone and D for sites occupied by an ADP-bound chaperone (e.g. (0,T,0,D,0,0), 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

(0,T,0,0,0,0)koffadp[Hsp70ADP]konadpeβΔΔG(0,T,0,D,0,0)(0,T,0,0,0,0)koffatp[Hsp70ADP]konatpeβΔΔG(0,T,0,T,0,0)
  • hydrolysis/synthesis

(0,T,0,T,0,0)ksskhs(0,T,0,D,0,0)
  • nucleotide exchange

(0,T,0,D,0,0)kex,TDeffkex,DTeff(0,T,0,T,0,0).

We further provide, as an example, the equation for a precise configuration, say (0,T,0,D,0,0) (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

(11) ddt(0,T,0,D,0,0)=(0,T,0,D,0,0)(kex,DTeff+kss+kex,TDeff+khs+koffatp+koffadp)++(0,0,0,D,0,0)[Hsp70ATP]konatpeβΔΔG++(0,T,0,0,0,0)[Hsp70ADP]konadpeβΔΔG++(T,T,0,D,0,0)koffatp++(D,T,0,D,0,0)koffadp++(0,T,T,D,0,0)koffatp++(0,T,D,D,0,0)koffadp++(0,T,0,D,T,0)koffatp++(0,T,0,D,D,0)koffadp++(0,T,0,D,0,T)koffatp++(0,T,0,D,0,D)koffadp.

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):

(12) kex,DTeff=kDk+T[ATP][ADP]k+D+k+T[ATP][ADP]
(13) kex,TDeff=kTk+Dk+D+k+T[ATP][ADP],

where k+D, k+T, k-D and k-T 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, kh, is much smaller than the same rate in the presence of the substrate, khs (khkhs). We assumed that the ratio between the rate of hydrolysis kh and the rate of synthesis ks is not altered by the substrate:

(14) khks=khskss.

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 konatp, konadp, koffatp and koffadp are the rate of substrate binding and unbinding from a chaperone in the ATP and ADP states, we must have

(15) konatpkhskoffadpkskonadpksskoffatpkh=konatpkoffadpkonadpkoffatp=1.

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:

(16) konatpkoffadpkonadpkoffatp=r
(17) konatp,koffadpkonatp/r1/4,koffadp/r1/4
(18) konadp,koffatpkonadpr1/4,koffatpr1/4.

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

(19) Hsp70+ADPHsp70ADPHsp70ATPHsp70+ATP.

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)
koffatp  2.31s-1koffadp  2*10-3s-1
k-T  1.33*10-4s-1k-D  0.022s-1
konatp  1.28*106M-1s-1konadp  103M-1s-1
k+T  1.3*105M-1s-1k+D  2.67*105M-1s-1
kh  6*10-4s-1khs  1.8s-1

To test the robustness of the model for the radius of gyration, the average FRET efficiency and the free-energy ΔGswell, 100 realizations were implemented, taking each time the values ΔGi from a Gaussian distribution with σ=0.3 kcal/mol and mean equal to the value obtained with the MD simulations.

Molecular graphics

Request a detailed protocol

Molecular 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
Normalized FRET efficiency histograms of the variant Δ39 (K135C/K174C).

Experiments were performed in the presence of 10 μM DnaK, 500 nM DnaJ, 1 mM ATP and either 5 mM (grey curves) or 100 mM (red curves) of KCl.

Appendix 1—figure 1—source data 1

Experimental data employed to generate the FRET histograms reported in the Appendix 1—figure 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 n monomers characterized by attractive interactions of average magnitude ϵ. Let αRg/Rg,0 be the expansion parameter, where Rg is the radius of gyration of the polymer, while Rg,0 is the value of Rg in the unperturbed case (ϵ=0 and no excluded volume present). Assuming a Flory-Fisk distribution for the unperturbed case, the probability distribution of the expansion parameter P(α) was shown to be (Sanchez, 1979)

(20) P(α)=1Zα6e-72α2+nq(ϵ,α),

where Z is the partition function and

q(ϵ,α)=12ϵϕ0α3-(α3ϕ0-1)ln(1-ϕ0α3),

with ϕ0=19/(27n). The free energy can then be straightforwardly computed as ΔG(α)=-kBTlnP(α), that is making the dependence on Rg explicit,

(21) ΔG(Rg)=-kBTlnP(Rg/Rg,0).

In the present case, rhodanese is constituted by 293 amino acids, that is there are Nb=292 bonds. The bond length is bl=3.8 Å, estimated as the typical distance between Cα 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 bK can be estimated as bK=2lp=8 Å, where lp=4 Å is the persistence length of an unfolded protein under native conditions (Hofmann et al., 2012). The number of monomers n in the equivalent FJC can then be computed by imposing the total contour length of rhodanese: n=Nbbl/bK. From there, both ϕ0 and the unperturbed radius of gyration Rg,0=nbK2/6 can be calculated, leaving ϵ as the only unknown quantity in Equation (20). The latter can be fixed by imposing that ΔG has a minimum at the experimental value Rg=20.1 Å obtained in the absence of chaperones (Kellner et al., 2014; Hofmann et al., 2014), yielding ϵ1.07 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

  1. Book
    1. Milo R
    2. Phillips R
    (2015)
    Cell Biology by the Numbers
    Garland Science.

Article and author information

Author details

  1. Salvatore Assenza

    1. Laboratory of Food and Soft Materials, ETH Zürich, Zürich, Switzerland
    2. Departmento de Física Teórica de la Materia Condensada, Universidad Autónoma de Madrid, Madrid, Spain
    Contribution
    Validation, Investigation, Methodology
    Contributed equally with
    Alberto Stefano Sassi
    Competing interests
    No competing interests declared
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0001-9983-8927
  2. Alberto Stefano Sassi

    1. Institute of Physics, School of Basic Sciences, École Polytechnique Fédérale de Lausanne (EPFL), Lausanne, Switzerland
    2. IBM TJ Watson Research Center, Yorktown Heights, New York, United States
    Contribution
    Validation, Investigation, Methodology
    Contributed equally with
    Salvatore Assenza
    Competing interests
    No competing interests declared
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0002-1269-4746
  3. Ruth Kellner

    Department of Biochemistry, University of Zurich, Zurich, Switzerland
    Contribution
    Validation, Investigation
    Competing interests
    No competing interests declared
  4. Benjamin Schuler

    1. Department of Biochemistry, University of Zurich, Zurich, Switzerland
    2. Department of Physics, University of Zurich, Zurich, Switzerland
    Contribution
    Resources, Methodology
    Competing interests
    No competing interests declared
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0002-5970-4251
  5. Paolo De Los Rios

    1. Institute of Physics, School of Basic Sciences, École Polytechnique Fédérale de Lausanne (EPFL), Lausanne, Switzerland
    2. Institute of Bioengineering, School of Life Sciences, Ecole Polytechnique Fédérale de Lausanne (EPFL), Lausanne, Switzerland
    Contribution
    Resources, Supervision, Methodology
    Competing interests
    No competing interests declared
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0002-5394-5062
  6. Alessandro Barducci

    Centre de Biochimie Structurale (CBS), INSERM, CNRS, Université de Montpellier, Montpellier, France
    Contribution
    Conceptualization, Resources, Supervision, Methodology, Project administration
    For correspondence
    alessandro.barducci@cbs.cnrs.fr
    Competing interests
    No competing interests declared
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0002-1911-8039

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

A two-part list of links to download the article, or parts of the article, in various formats.

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)

  1. Salvatore Assenza
  2. Alberto Stefano Sassi
  3. Ruth Kellner
  4. Benjamin Schuler
  5. Paolo De Los Rios
  6. Alessandro Barducci
(2019)
Efficient conversion of chemical energy into mechanical work by Hsp70 chaperones
eLife 8:e48491.
https://doi.org/10.7554/eLife.48491

Share this article

https://doi.org/10.7554/eLife.48491

Further reading

    1. Computational and Systems Biology
    2. Physics of Living Systems
    Divyoj Singh, Sriram Ramaswamy ... Mohd Suhail Rizvi
    Research Article Updated

    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.

    1. Physics of Living Systems
    Tommaso Amico, Samuel Toluwanimi Dada ... Amos Maritan
    Research Article

    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.