Abstract
Theory, simulations and experimental results have suggested an important role of internal friction in the kinetics of protein folding. Recent experiments on spectrin domains provided the first evidence for a pronounced contribution of internal friction in proteins that fold on the millisecond timescale. However, it has remained unclear how this contribution is distributed along the reaction and what influence it has on the folding dynamics. Here we use a combination of single-molecule Förster resonance energy transfer, nanosecond fluorescence correlation spectroscopy, microfluidic mixing and denaturant- and viscosity-dependent protein-folding kinetics to probe internal friction in the unfolded state and at the early and late transition states of slow- and fast-folding spectrin domains. We find that the internal friction affecting the folding rates of spectrin domains is highly localized to the early transition state, suggesting an important role of rather specific interactions in the rate-limiting conformational changes.
Similar content being viewed by others
Introduction
In spite of the complexity of protein folding, which involves the interactions of thousands of atoms, the experimentally observed kinetics of folding can often be described in terms of simple activated reactions that lead to the interconversion between well-defined thermodynamic states. However, experiments have recently started to reveal more of the kinetic complexity underlying the folding process. Examples are the identification of the relaxation dynamics of ‘downhill’ folding1,2, the measurement of transition path times3 or the presence of internal friction4,5,6,7,8. These observations can no longer be accounted for by simple chemical kinetics, but require more detailed physical concepts; in this way, they provide an important link to current theories of protein folding9,10,11,12 and molecular simulations13,14.
A topic that is attracting increasing interest is the role of internal friction in the folding process. The experimental identification of internal friction is commonly based on a deviation from the direct proportionality of the folding time (the inverse folding rate coefficient) to solvent viscosity, resulting in a non-zero value of the folding time when extrapolated to a solvent viscosity of zero5. This observation indicates that the folding (or unfolding) rate is not only influenced by friction caused by collisions with solvent molecules (solvent or ‘external’ viscosity) but also by friction within or between parts of the protein15 (protein or ‘internal’ viscosity). The molecular cause of internal friction is often difficult to identify in detail, but in principle, it can originate from any mechanism that leads to a dissipation of energy into internal degrees of freedom of the protein that do not contribute to progress along the folding reaction. Such mechanisms could range from random collisions with surrounding protein atoms to dihedral angle rotations or the transient formation and breakage of side chain or backbone interactions5,16,17. We use the term ‘internal friction’ to distinguish it from effects on the dynamics that scale with the bulk solvent viscosity and from activation energy terms. Based on recent experimental results, it has been suggested that internal friction effects can be particularly obvious for the kinetics of very fast-folding proteins, where the timescales characterizing the crossing of the transition state region (and even the unfolded state dynamics) can be slow enough to significantly affect the progress of the folding reaction4,6. For slower-folding proteins, however, such effects have largely remained undetectable5,18,19,20,21.
It thus came as a surprise when recent experiments provided evidence for pronounced internal friction in the rate-limiting transition state for folding of two spectrin domains8. Although the folding time of the fast-folding spectrin domain R15 is directly proportional to solvent viscosity, as observed previously for other millisecond folders18,19,20,21, the folding of the structurally very similar domains R16 and R17 is much slower, and the solvent viscosity dependence of their folding times is dramatically reduced. Wensley et al.8 concluded that the slow folding of R16 and R17 is at least in part due to internal friction. However, it has remained unclear whether this unusually large contribution of internal friction in R16 and R17 is present throughout the energy landscape underlying the folding reaction or whether it is confined to the transition state region.
The role of internal friction in protein folding is of great interest also for the theoretical description of protein folding. A model that is frequently used to conceptualize protein-folding dynamics, especially for the analysis of simulations and experiments, is that of diffusion on a free energy surface3,7,9,22,23,24,25,26. In the simplest case, the folding dynamics are represented in terms of diffusion along a single reaction coordinate. Even though the required projection of the large number of degrees of freedom bears the risk of obscuring structural or mechanistic details10,25, suitable reaction coordinates that lead to a faithful description of the folding kinetics and mechanisms can often be identified16,24,25. However, the effective diffusion coefficient can then no longer be assumed to be invariant along the resulting free energy surface16, that is, as the polypeptide chain transits from the highly disordered unfolded state via its transition state ensemble towards the natively folded structure. Conceptually, diffusive models thus allow effects on protein-folding dynamics to be partitioned into differences in barrier heights on the one hand and differences in effective diffusion coefficients on the other16,24,26,27,28 (even though this partitioning will depend on the choice of the reaction coordinate and can thus be somewhat ambiguous). The internal friction observed experimentally is expected to affect the diffusion coefficient, and therefore provides an opportunity to address the question of the position dependence of the diffusion coefficient along the folding reaction.
Here we assess internal friction in the spectrin domains with two complementary approaches. The first is single-molecule Förster resonance energy transfer (FRET) combined with nanosecond fluorescence correlation spectroscopy (nsFCS), which can be used to quantify unfolded state dynamics down to the nanosecond range29,30. The second approach is the analysis of denaturant-dependent protein-folding kinetics, which allows us to differentiate the viscosity dependencies of the folding and unfolding rates over the early and late transition states of the slow folders R16 and R1731,32. By quantitatively comparing the contributions of internal friction to the dynamics in (a) the compact unfolded state, (b) the early transition state and (c) the late transition state, we can establish whether the internal friction is evenly distributed or localized to specific regions along the free energy surface.
Results
Single-molecule FRET experiments
In our single-molecule analysis, we focus on R15 and R17, the domains showing the lowest and the highest contributions of internal friction to their folding kinetics, respectively8. To investigate the spectrin domains with single-molecule FRET, variants of R15 and R17 with two cysteine residues were prepared and labelled with Alexa 488 and 594 as FRET donor and acceptor, respectively (Fig. 1a). Single-molecule FRET by confocal detection of freely diffusing molecules allows the separation of folded and unfolded subpopulations in transfer efficiency histograms (Fig. 1b), and thus the investigation of unfolded state properties even in the presence of a majority of folded molecules. As expected from previous observations on other proteins33,34, unfolded spectrin domains undergo a continuous compaction when the denaturant concentration is reduced and native conditions are approached (see Fig. 1c). The transfer efficiencies and fluorescence lifetimes (Supplementary Fig. S1) of the unfolded subpopulation can be used to quantify the dimensions of the unfolded state based on the intramolecular distance distributions of simple polymer models33,35,36,37,38,39. In a second step, the dynamics of the unfolded polypeptide are described in terms of diffusion on the potential of mean force corresponding to this distance distribution. To determine the relaxation time of the system (that is, the reconfiguration time, τr, of the chain), we use nsFCS29,30,40 (see Supplementary Information for details), which reports directly on the timescale of the distance fluctuations between the donor and acceptor chromophores (Fig. 1d–f).
As in the case of FRET efficiency histograms, the nsFCS experiments require a sufficient fraction of unfolded molecules to be present in the sample. However, the populations of the unfolded state at guanidinium chloride (GdmCl) concentrations far below the unfolding midpoint become too low to be detected reliably. To probe unfolded R17 even in the virtual absence of denaturant, where the conformations are most compact and internal friction is expected to be most pronounced, we employed a microfluidic device that allows millisecond mixing and was optimized for single-molecule fluorescence detection41,42 (Fig. 2). The denaturant is diluted rapidly by mixing GdmCl-unfolded protein with buffer solution, which enables us to populate the unfolded state transiently under native conditions. The resulting compact unfolded state can be probed by confocal single-molecule measurements in the observation channel of the device directly after mixing, where a large fraction of protein molecules are still unfolded (Fig. 2a and b). The unfolded subpopulation will convert to the folded state as the molecules flow down the observation channel (Supplementary Fig. S2). The dynamics of the unfolded state can also be measured in the mixer: Fig. 2d–f shows examples of donor–donor, donor–acceptor and acceptor–acceptor fluorescence intensity correlation curves with a relaxation on the 100 ns timescale, typical of unfolded state dynamics29,30,43,44. The relaxation time is extracted by a global fit to all three correlations (see Supplementary Information for details). Together with the information about the distance distribution from the transfer efficiency histograms, a reconfiguration time or an effective intramolecular diffusion coefficient can be calculated29,30,40. In this way, unfolded state dynamics can be measured over a wide range of denaturant concentrations, and can then be used to quantify the contribution of internal friction.
Internal friction in the unfolded state of spectrin domains
Figure 3 shows the reconfiguration times of unfolded R15 and R17 over the entire accessible denaturant concentration range. With the microfluidic mixing experiments, τr for R17 could be measured down to 0.03 M GdmCl; the folding of R15 in these conditions is too fast for the unfolded state to be populated sufficiently in the mixer. For both R15 and R17, τr ranges between 50 and 200 ns; it shows a significant increase at high GdmCl concentration due to the increasing solvent viscosity, and, more interestingly, an increase below ~2 M GdmCl. The latter observation indicates the presence of internal friction in the compact unfolded state29, but the similarity of the reconfiguration times of both spectrin domains shows that the slow- and fast-folding domains do not exhibit pronounced differences in internal friction in their unfolded states.
To quantify the extent of internal friction, we investigated the solvent viscosity dependence of τr at several GdmCl concentrations (Fig. 3). On extrapolation to zero solvent viscosity, we observed in all cases a value of τr >0, confirming the presence of internal friction5,6,15. Surprisingly, internal friction is present even at 6 M GdmCl, where both domains are very expanded (Fig. 1c), and therefore the interactions within the polypeptide chains should be small. To test the possible influence of the specific viscogen used, both glycerol and glucose were employed. The good agreement of the values for τr extrapolated to zero solvent viscosity (Fig. 3, insets) indicates the absence of specific interactions between protein and viscogen that might affect the results. Note also that transfer efficiency histograms provide a sensitive way of verifying that the distance distributions and thus the energetics of interactions within the unfolded states are independent of the viscogen concentration (see Supplementary Fig. S3), that is, the observed changes in τr are not due to changes in unfolded state dimensions.
Can we employ these results to estimate the magnitude of internal friction in the unfolded state over the entire range of GdmCl concentrations? According to a commonly used relation from polymer dynamics, which can be derived rigorously in the framework of simple models, such as the Rouse model with internal friction45,46,47, the overall reconfiguration time of the unfolded chain, τr, results from the additive contributions of τs, the reconfiguration time in the absence of internal friction, and τi, a relaxation time due to internal friction processes. The common assumption that only the dynamics corresponding to τs depend on solvent viscosity, η, whereas τi does not45,46,47, leads to the relation
where η0 is the viscosity of the solvent in the absence of viscogen44. If we take into account the linear scaling of τs with the mean squared end-to-end distance, r2, of the chain segment, we obtain τs η r2 (in this case, we assume the scaling expected for a Rouse chain, but other common models of polymer dynamics, for example the Zimm model, yield very similar results44). The values of τi measured via the viscosity dependences can be used as constraints for obtaining the GdmCl dependence of both τi and τs (Fig. 3) (see Methods for details). The result shows that internal friction dominates the unfolded state dynamics of the spectrin domains at low GdmCl concentrations, but τi remains in the range of 100–200 ns for both variants, orders of magnitude lower than the folding time. In summary, we find no pronounced difference in τi between fast- and slow-folding spectrin domains, which indicates that the internal friction observed in the folding kinetics of the slow-folding domains is not present in the unfolded state, that is, the part of the reaction preceding the transition state region.
Internal friction in the transition state region
The early transition state, which is rate-limiting at low concentrations of denaturant, has been the subject of previous experiments8. Here we extend our investigation to the late transition state, which becomes rate-limiting at high GdmCl concentrations8,48,49,50, by taking advantage of the nonlinearity in the dependences of the logarithm of the unfolding rate constants on denaturant concentration (the chevron plots)31,32,51 observed for R1648,49 and R178,50 (Fig. 4a). Several models have been proposed to interpret such behaviour31,32,51. For our analysis, we assumed the presence of two sequential transition states separated by a high-energy intermediate49 to quantify to what extent each of these transition states contributes to the kinetics of the overall folding reaction and how sensitive the two relative rates are to solvent viscosity. This approach allows us to assess the presence of internal friction separately for the two transition states (note, however, that our interpretation would also be valid in the framework of a broad energy barrier scenario49,52).
Figure 4b shows the dependence on solvent viscosity of the relative folding and unfolding rates of R15 and those corresponding to the early and late transition states of R16 and R17. Folding and unfolding of R16 and R17 over the early transition state show much less solvent viscosity dependence than folding and unfolding of R158. This behaviour reflects the large internal friction in the early transition state of R16 and R17. Intriguingly, the behaviour resulting from unfolding over the late transition state (Fig. 4a) is very different. In this case, a pronounced solvent viscosity dependence of the folding and unfolding rates is observed for both R16 and R17 (Fig. 4b), with a slope close to unity, as expected for a process without internal friction. The kinetics corresponding to the second transition state of R15 are not accessible by stopped-flow measurements because of the high unfolding rates at high denaturant concentrations. However, considering the fast folding of R15, the linear dependence of its folding and unfolding rates on solvent viscosity, and the mechanistic interpretation suggested on the basis of extensive protein-engineering experiments53, it is very unlikely for internal friction to develop in the context of a hypothetical late transition state. This surprising result clearly indicates that only the early, less native-like, transition states of R16 and R17 show indications for the presence of internal friction; in the second transition states, internal friction is undetectably low. In a next step, can we go beyond such a qualitative statement and compare internal friction at the different stages of the folding reaction more quantitatively?
Internal friction along the folding reaction
A common feature of the viscosity dependencies of both the reconfiguration times in the unfolded state and the folding and unfolding times of the spectrin domains is the existence of a non-zero intercept with the relaxation time axis on extrapolation to η=0. This behaviour is the characteristic signature of friction within the protein that is largely independent of solvent viscosity15,45. For the unfolded state, the overall relaxation time of the chain is well described by equation (1), as expected from polymer theory45,46,47. For the folding dynamics, we use a modified Kramers-type picture originally suggested by Ansari et al.15 in the context of conformational changes in folded proteins. This model assumes that the contribution of friction from collisions, with the solvent on the one hand and other atoms within the protein on the other, results in an additivity of external (that is, solvent-mediated) and internal (that is, protein-mediated) viscosities, η and σ, respectively. The resulting expression for the folding time, τf (the inverse folding rate coefficient), is
where σTS is the value of σ in the transition state region26,54; C contains all contributions to the preexponential factor except η and σTS; Δ is the height of the free energy barrier; kB is the Boltzmann constant and T is temperature. In analogy to equation (1), τf can be considered a sum of a solvent viscosity-dependent relaxation time, τfs, and a solvent viscosity-independent relaxation time, τfi, which would equal zero in the absence of internal friction. An analogous expression is valid for the unfolding time, τu. Assuming that the variation of solvent viscosity does not change the equilibrium constant of the reaction (note that the equilibrium denaturation free energy is kept constant here by the addition of GdmCl), σTS is expected to be the same for folding and unfolding kinetics. Because of microscopic reversibility, the extrapolation of the folding and unfolding times, respectively, to η=0 should thus yield the same value of σTS, that is, the contribution of internal viscosity at the transition state should be independent of the direction of barrier crossing. The viscosity dependences of the folding and unfolding rates observed experimentally (Fig. 5) for R15, R16 and R17 are in good agreement with this prediction (Fig. 5). Note that for Δ→0, that is, in the absence of a free energy barrier, equation (2) converges to equation (1), indicating the consistency of equation (1) for unfolded state dynamics with equation (2) for diffusive barrier crossing. Note also that a non-additive contribution of internal friction (or internal viscosity, equation (2)) to the relaxation dynamics of the system55 would always lead to an overall relaxation time of zero when extrapolated to zero solvent viscosity and is thus not considered here. Likewise, in view of the close structural similarity of the spectrin domains, generic solvent properties56 are unlikely to cause the pronounced differences in experimentally observed internal friction (Fig. 3)8.
Based on equation (2) and Fig. 4b, we can compare the values of τfi corresponding to the early and late transition states of the different spectrin variants (Fig. 6). We obtain values for τfi of (95±11) and (7.3±0.5) ms for the early transition states of R17 and R16, respectively, but only (0.016±0.009) ms for R15, which in part reflects the difference in internal friction between the slow- and fast-folding spectrin domains. For the late transition state, we can only estimate upper bounds for τfi because the underlying internal friction is too low to influence the overall kinetics to a detectable degree; we obtain 28 and 10 μs for R17 and R16, respectively (see Methods for details). For comparison, the values of τi in the unfolded states of R17 and R15 under native conditions are only about 110 and 200 ns, respectively. For the rate-limiting early transition states of R15, R16 and R17, we obtain values for σTS of (0.14±0.08), (3.8±1.0) and (5.6±2.0) mPa s, respectively (Fig. 5, Supplementary Table S1). According to equation (2), we expect the folding and unfolding times to be increased compared with the case without internal friction (but assuming the same values of Δ and C) by (σTS+η)/η, that is, by factors of approximately 5 and 7 for R16 and R17, respectively, (Supplementary Table S1). Remarkably, the relaxation rates of R15 and the recently identified E18F/K25V variant of R16, which folds and unfolds 40-fold faster than wild-type R16 but has an unchanged value of σTS, differ by about this factor of five57. This agreement supports the hypothesis that the 40-fold increase in rate for the variant is caused by a decrease in barrier height, but the remaining difference in rates between R15 and R16 E18F/K25V is caused by the difference in σTS, that is, in internal friction57.
Discussion
We find that the marked effects of internal friction on the folding and unfolding rates of R16 and R17 are present only in the early transition state; in the late transition state, internal friction remains undetectable; in the compact unfolded states of both slow- and fast-folding spectrin domains, even though they show a pronounced influence of internal friction on their reconfiguration dynamics, the internal friction times are in the sub-microsecond range, orders of magnitude too low for having a measurable effect on the millisecond folding times. The effect of internal friction on the free energy surface of the spectrin domains thus seems to be far from evenly distributed along the reaction; instead, it is highly localized to the early transition state (Fig. 6).
This finding is surprising, given the reasonable assumption that internal friction should increase monotonically with the number of interactions within a polypeptide chain, which would imply that it should increase on moving from the early to the late transition state. Such a scenario, which assumes that internal friction results from averaging over rather uniform properties of the polypeptide chain, and describes, for example, the behaviour we observe in the unfolded state (Fig. 3)29,44, does not seem to apply to the barrier-crossing process. Instead, our results are more consistent with a group of rather specific interactions that need to be correctly established during passage of the early transition state, in keeping with a nucleation-type mechanism58. Additional support for this interpretation comes from the results of extensive protein-engineering experiments53, which show that changing only five amino acids in a key region of the A helix is sufficient to increase the folding and unfolding rates of the slow folders R16 and R17 into the same range as the much faster-folding R15 and to greatly reduce the influence of internal friction. The Φ-value analysis of the corresponding variants is indicative of a concomitant change in folding mechanism, from a diffusion–collision-type in R16 and R17 to the nucleation–condensation-type reaction characteristic of R15. This seems to confirm the hypothesis that the origin of the internal friction that leads to the insensitivity to solvent viscosity observed for R16 and R17 may lie in the search for the correct register in the packing of preformed helices8,53. The low solvent viscosity dependence would indicate, however, that only small-amplitude motions against the solvent are involved.
What do these findings implicate for a theoretical description of protein-folding dynamics? In the simplest one-dimensional diffusive model based on Kramers theory54, the rate of barrier crossing, kf, is given by
where ωU and ωTS characterize the curvatures of the potential in the unfolded well and on top of the barrier (that is, in the transition state region), respectively, and DTS is the effective diffusion coefficient in the barrier region. In the case of experiments based on single-molecule FRET or force spectroscopy, for instance, the reaction coordinate would correspond to an intramolecular distance, and the corresponding shape and diffusion parameters are starting to become available from optical and force-based single-molecule experiments29,36,39,44,59. For the spectrin domains, where the results suggest that interactions between a few residues have a key role for internal friction in folding kinetics, it is plausible that a relatively small change in the number of contacts, resulting from conformational rearrangements involved in crossing the transition state barrier, could give rise to a comparatively large change in free energy. In the Kramers-type picture (equation (3)), such behaviour would correspond to a larger curvature of the free energy surface at the transition state compared with that for unfolded state, resulting in a substantial difference between ωTS and ωU23. The results shown in Fig. 5 suggest that a diffusive description of the folding dynamics is applicable, in spite of the importance of highly specific interactions at the transition state. A promising strategy for testing the importance of internal friction on the dynamics in the transition state region and the applicability of a diffusive description for the folding of spectrin domains might be measurements of transition path times, which are starting to become accessible in single-molecule experiments3; for folding reactions involving internal friction, the transition path time should be longer by a factor of approximately (σTS+η)/η compared with the reaction without internal friction3. We hope that the increasing amount of detail available from experimental results like ours will stimulate simulations and the development of quantitative theoretical concepts that go beyond the simple one-dimensional picture used here, and will allow us to better describe the complexity of folding dynamics and understand the role and mechanisms of internal friction.
Methods
Preparation and labelling of proteins
The labelling sites were selected to minimize the effects on the conformational stability and folding mechanism, based on the extensive previous protein-engineering studies48,50. The separations of the dyes in the sequence and the folded structure were chosen such that (a) the transfer efficiencies expected in the unfolded state would be in a range resulting in good sensitivity for the determination of intramolecular distance distributions and dynamics, and that (b) a pronounced separation of folded and unfolded subpopulations in FRET efficiency histograms could be achieved. For labelling of the spectrin domains, cysteine residues were introduced by site-directed mutagenesis at positions 39 and 99, and the proteins were expressed and purified as described previously60. In R17, an endogenous cysteine at position 68 was exchanged to alanine to avoid multiple labelling. For labelling, a 1.3:1 molar excess of reduced protein was incubated with Alexa Fluor 488 maleimide (Invitrogen) at 4 °C for ~10 h. Unreacted dye was removed by gel filtration (G25 desalting; GE Healthcare Biosciences AB, Uppsala, Sweden), and the protein was incubated with Alexa Fluor 594 maleimide at room temperature for ~2 h. Differently labelled variants were separated by ion-exchange chromatography (MonoQ HR 5/5; GE Healthcare Biosciences AB).
Ensemble folding and unfolding kinetics of R15, R16 and R17 were monitored by the change in tryptophan fluorescence on an SX20 stopped-flow spectrometer (Applied Photophysics) and fitted as described previously8. Briefly, chevron plots for R15 were fitted to a two-state model and those for R16 and R17 to a sequential transition state model31 as described in detail by Scott and Clarke49. The calculation of folding and unfolding times (τf and τu) at ΔG=1.5 kcal mol−1 and of the relaxation time τmp at ΔG=0 has previously been described for the first transition state8. Those for the second transition state were calculated in an analogous manner using the fitted second transition state parameters.
Single-molecule fluorescence spectroscopy
Measurements were performed using both a custom-built confocal microscope and a MicroTime 200 confocal microscope equipped with a HydraHarp 400 counting module (PicoQuant, Berlin, Germany). The donor dye was excited with a diode laser at 485 nm (dual mode: continuous wave and pulsed; LDH-D-C-485, PicoQuant) at an average power of 100 μW at the sample. Single-molecule FRET efficiency histograms were acquired in samples with a protein concentration of about 20–50 pM, with the laser in either continuous-wave mode or pulsed mode at a repetition rate of 64 MHz; photon counts were recorded with a resolution of 16 ps at the counting electronics (time resolution was thus limited by the timing jitter of the detectors). The nsFCS measurements were performed in samples with a protein concentration of ~1 nM, over a period of 16 h, at an average power of 70 μW with the laser in continuous-wave mode. For rapid mixing experiments, microfluidic mixers fabricated by replica moudling in polydimethylsiloxane (PDMS) were used as described previously41,42. The experiments were performed with pressures of 13.8 kPa (2.0 psi) applied to the side channels (no. 3 and 4 in Fig. 2c) and 6.9 kPa (1.0 psi) to the middle channel. The flow stability was monitored by analysing the time evolution of the fluorescence intensity cross-correlation. The nsFCS measurements in the mixer were performed over a period of 2–10 h and analysed as detailed below (‘nsFCS measurements’). Controls of the contribution of dye quenching to the acceptor autocorrelation were carried out with a continuous wave HeNe laser (594 nm; CVI Melles Griot, Albuquerque, NM, USA) at a power of 16 μW at the sample. All measurements were performed in 50 mM sodium phosphate buffer, pH 7.0, 140 mM β-mercaptoethanol (Sigma), 20mM cysteamine hydrochloride (Sigma) and 0.001% Tween 20 (Pierce) with varying concentrations of GdmCl (Pierce). Tween 20 was used to prevent surface adhesion of the proteins, and the photoprotective agents β-mercaptoethanol and cysteamine hydrochloride were employed to minimize chromophore damage and enhance brightness. For experiments in the microfluidic device, the Tween 20 concentration was increased to 0.01%.
Guanidine concentrations were measured with an Abbe refractometer (Krüss, Germany), and viscosities of the solutions were measured with a digital viscometer (DV-I+; Brookfield Engineering, Middleboro, MA, USA) with a CP40 spindle at 30 or 60 r.p.m., which allows determination of viscosity with an uncertainty of 0.05–0.1 mPa·s.
nsFCS measurements
Fluorescence intensities of donor and acceptor fluorophores were recorded with a four-channel HydraHarp 400 Picosecond Event Timer (PicoQuant) and correlated with a binning time of 1 ns. To avoid effects of detector dead times and afterpulsing on the correlation functions, the signal was recorded with two detectors each for donor and acceptor and cross-correlated between detectors. Auto- and cross-correlation curves were fitted over a time window of 2.5 μs with
where i and j correspond to donor or acceptor fluorescence emission, N is the effective mean number of molecules in the confocal volume, cab, τab, ccd and τcd are the amplitudes and time constants of the antibunching and chain dynamics, respectively29,30, whereas cT and τT refer to the triplet blinking of the fluorophores. The three correlation curves were fitted globally with the same values of τcd and a fixed value of τT estimated from the fit of the whole acceptor–acceptor autocorrelation curve from nanoseconds to seconds. The amplitude and the lifetime of the antibunching and the triplet amplitude were fit with a free independent decay component for each correlation curve. The errors associated with the global fit were estimated by resampling the experimental dataset via a bootstrap method.
Viscosity dependence of τ r
The reconfiguration times τr measured as a function of viscosity at different GdmCl concentrations (Fig. 3 insets) were fitted to a linear function, yielding the reconfiguration time at 0 M GdmCl, which corresponds to the τi, that is, the internal friction contribution, at every set of conditions. To quantify the contribution of internal friction to the observed reconfiguration time, we subtract τi (6 M), the average of the two values obtained for the internal friction time in 6 M GdmCl with the two viscogens, from the value of the reconfiguration time calculated at the same GdmCl concentration, τr (6 M). The resulting time, τs (6 M), represents the reconfiguration time of an ideal polypeptide chain in the absence of internal friction. τs depends on both the polypeptide chain dimensions and the solvent viscosity, therefore its value can be calculated at any GdmCl concentration, provided that the dimension of the chains are known, according to:
where η is the solvent viscosity, and and are the GdmCl-dependent dye-to-dye distance and the dye-to-dye distance at 6 M GdmCl, respectively, obtained from the transfer efficiency histograms. The value of τi is then calculated at all GdmCl concentrations by subtracting τs from τr. Because the quantification of τi and τs relies on the value of τi at 6 M GdmCl, an estimation of the error associated with these quantities can be obtained by propagating the uncertainty of τi (6 M). The s.d. of τi (6 M) is calculated from the values of τi extrapolated at η=0 from the viscosity-dependent experiment carried out with glycerol and glucose. In Fig. 3, the grey- and blue-shaded areas for τi and τs are calculated with equation (5), assuming a distribution of τi (6 M) within one s.d.
Calculation of friction parameters for the transition states
From the definition of σ at the transition state (equation (2) and Supplementary Fig. S5),
where τf (η=0) is the folding time extrapolated to zero solvent viscosity, and τfs is the solvent-dependent contribution to the folding time, calculated as τf −τfi, where τfi is the solvent-independent internal friction time. Given that the value of σ seems to be invariant with GdmCl concentration (‘Internal friction along the folding reaction’ in main text), we calculate σ for the early transition state from a global fit of refolding, unfolding and midpoint kinetics (Fig. 5). τfi can then be calculated for any denaturant concentration from
Individual and global linear fits of the three data sets corresponding to refolding, unfolding and midpoint kinetics over the late transition state in R16 and R17 (Fig. 4) extrapolate to small negative values on the τobs axis (Supplementary Fig. S4, corresponding to small negative values on the k0/kobs axis of Fig. 4b). Instead of just assuming that internal friction is negligible within experimental error for these processes, we calculated confidence intervals (confidence level of 90%) for the individual fits of each of these three data sets for each domain and used the average of the highest values within these confidence intervals to estimate an upper bound on σ within experimental uncertainty (Supplementary Fig. S4). The maximum σ values were then used to calculate an upper bound for the expected internal friction times for the two domains according to equation (7) (Supplementary Table S1).
Additional information
How to cite this article: Borgia, A. et al. Localizing internal friction along the reaction coordinate of protein folding by combining ensemble and single-molecule fluorescence spectroscopy. Nat. Commun. 3:1195 doi: 10.1038/ncomms2204 (2012).
References
Yang W. Y., Gruebele M. Folding at the speed limit. Nature 423, 193–197 (2003).
Li P., Oliva F. Y., Naganathan A. N., Munoz V. Dynamics of one-state downhill protein folding. Proc. Natl Acad. Sci. USA 106, 103–108 (2009).
Chung H. S., McHale K., Louis J. M., Eaton W. A. Single-molecule fluorescence experiments determine protein folding transition path times. Science 335, 981–984 (2012).
Qiu L. L., Hagen S. J. A limiting speed for protein folding at low solvent viscosity. J. Am. Chem. Soc. 126, 3398–3399 (2004).
Hagen S. J., Qiu L. L., Pabit S. A. Diffusional limits to the speed of protein folding: fact or friction? J. Phys. Cond. Mat. 17, S1503–S1514 (2005).
Cellmer T., Henry E. R., Hofrichter J., Eaton W. A. Measuring internal friction of an ultrafast-folding protein. Proc. Natl Acad. Sci. USA 105, 18320–18325 (2008).
Liu F., Nakaema M., Gruebele M. The transition state transit time of WW domain folding is controlled by energy landscape roughness. J. Chem. Phys. 131, 195101 (2009).
Wensley B. G. et al. Experimental evidence for a frustrated energy landscape in a three-helix-bundle protein family. Nature 463, 685–U122 (2010).
Onuchic J. N., Luthey Schulten Z., Wolynes P. G. Theory of protein folding: the energy landscape perspective. Annu. Rev. Phys. Chem. 48, 545–600 (1997).
Shakhnovich E. Protein folding thermodynamics and dynamics: where physics, chemistry, and biology meet. Chem. Rev. 106, 1559–1588 (2006).
Thirumalai D., O'Brien E. P., Morrison G., Hyeon C. Theoretical perspectives on protein folding. Annu. Rev. Biophys. 39, 159–183 (2010).
Dill K. A., Ozkan S. B., Shell M. S., Weikl T. R. The protein folding problem. Annu. Rev. Biophys. 37, 289–316 (2008).
Snow C. D., Sorin E. J., Rhee Y. M., Pande V. S. How well can simulation predict protein folding kinetics and thermodynamics? Annu. Rev. Biophys. Biomol. Struct. 34, 43–69 (2005).
Shaw D. E. et al. Atomic-level characterization of the structural dynamics of proteins. Science 330, 341–346 (2010).
Ansari A., Jones C. M., Henry E. R., Hofrichter J., Eaton W. A. The role of solvent viscosity in the dynamics of protein conformational changes. Science 256, 1796–1798 (1992).
Best R. B., Hummer G. Coordinate-dependent diffusion in protein folding. Proc. Natl Acad. Sci. USA 107, 1088–1093 (2010).
Schulz J. C., Schmidt L., Best R. B., Dzubiella J., Netz R. R. Peptide chain dynamics in light and heavy water: zooming in on internal friction. J. Am. Chem. Soc. 134, 6273–6279 (2012).
Chrunyk B. A., Matthews C. R. Role of diffusion in the folding of the alpha subunit of tryptophan synthase from Escherichia coli. Biochemistry 29, 2149–2154 (1990).
Plaxco K. W., Baker D. Limited internal friction in the rate-limiting step of a two-state protein folding reaction. Proc. Natl Acad. Sci. USA 95, 13591–13596 (1998).
Jacob M., Geeves M., Holtermann G., Schmid F. X. Diffusional barrier crossing in a two-state protein folding reaction. Nat. Struct. Biol. 6, 923–926 (1999).
Ladurner A. G., Fersht A. R. Upper limit of the time scale for diffusion and chain collapse in chymotrypsin inhibitor 2. Nat. Struct. Biol. 6, 28–31 (1999).
Socci N. D., Onuchic J. N., Wolynes P. G. Diffusive dynamics of the reaction coordinate for protein folding funnels. J. Chem. Phys. 104, 5860–5868 (1996).
Klimov D. K., Thirumalai D. Viscosity dependence of the folding rates of proteins. Phys. Rev. Lett. 79, 317–320 (1997).
Best R. B., Hummer G. Diffusive model of protein folding dynamics with Kramers turnover in rate. Phys. Rev. Lett. 96, 228104 (2006).
Krivov S. V., Karplus M. Diffusive reaction dynamics on invariant free energy profiles. Proc. Natl Acad. Sci. USA 105, 13841–13846 (2008).
Best R. B., Hummer G. Diffusion models of protein folding. Phys. Chem. Chem. Phys. 13, 16902–16911 (2011).
Chahine J., Oliveira R. J., Leite V. B. P., Wang J. Configuration-dependent diffusion can shift the kinetic transition state and barrier height of protein folding. Proc. Natl. Acad. Sci. USA 104, 14646–14651 (2007).
Wang J., Oliveira R. J., Whitford P. C., Chahine J., Leite V. B. P. Coordinate and time-dependent diffusion dynamics in protein folding. Methods 52, 91–98 (2010).
Nettels D., Gopich I. V., Hoffmann A., Schuler B. Ultrafast dynamics of protein collapse from single-molecule photon statistics. Proc. Natl Acad. Sci. USA 104, 2655–2660 (2007).
Nettels D., Hoffmann A., Schuler B. Unfolded protein and peptide dynamics investigated with single-molecule FRET and correlation spectroscopy from picoseconds to seconds. J. Phys. Chem. B 112, 6137–6146 (2008).
Sánchez I. E., Kiefhaber T. Evidence for sequential barriers and obligatory intermediates in apparent two-state protein folding. J. Mol. Biol. 325, 367–376 (2003).
Walkenhorst W. F., Green S. M., Roder H. Kinetic evidence for folding and unfolding intermediates in staphylococcal nuclease. Biochemistry 36, 5795–5805 (1997).
Schuler B., Eaton W. A. Protein folding studied by single-molecule FRET. Curr. Opin. Struct. Biol. 18, 16–26 (2008).
Ziv G., Thirumalai D., Haran G. Collapse transition in proteins. Phys. Chem. Chem. Phys. 11, 83–93 (2009).
Schuler B., Lipman E. A., Eaton W. A. Probing the free-energy surface for protein folding with single-molecule fluorescence spectroscopy. Nature 419, 743–747 (2002).
Hoffmann A. et al. Mapping protein collapse with single-molecule fluorescence and kinetic synchrotron radiation circular dichroism spectroscopy. Proc. Natl Acad. Sci. USA 104, 105–110 (2007).
Sherman E., Haran G. Coil-globule transition in the denatured state of a small protein. Proc. Natl Acad. Sci. USA 103, 11539–11543 (2006).
Merchant K. A., Best R. B., Louis J. M., Gopich I. V., Eaton W. A. Characterizing the unfolded states of proteins using single-molecule FRET spectroscopy and molecular simulations. Proc. Natl Acad. Sci. USA 104, 1528–1533 (2007).
Hofmann H. et al. Polymer scaling laws of unfolded and intrinsically disordered proteins quantified with single-molecule spectroscopy. Proc. Natl Acad. Sci. USA 109, 16155–16160 (2012).
Gopich I. V., Nettels D., Schuler B., Szabo A. Protein dynamics from single-molecule fluorescence intensity correlation functions. J. Chem. Phys. 131, 095102 (2009).
Hofmann H et al. Single-molecule spectroscopy of protein folding in a chaperonin cage. Proc. Natl Acad. Sci. USA 107, 11793–11798 (2010).
Pfeil S. H., Wickersham C. E., Hoffmann A., Lipman E. A. A microfluidic mixing system for single-molecule measurements. Rev. Sci. Instrum. 80, 055105 (2009).
Hillger F. et al. Probing protein-chaperone interactions with single molecule fluorescence spectroscopy. Angew. Chem. Int. Ed. 47, 6184–6188 (2008).
Soranno A. et al. Quantifying internal friction in unfolded and intrinsically disordered proteins with single molecule spectroscopy. Proc. Natl Acad. Sci. USA early edition, www.pnas.org/cgi/doi/10.1073/pnas.1117368109 (2012).
De Gennes P. G. Scaling concepts in polymer physics Cornell University Press: Ithaca and London, (1979).
Portman J. J., Takada S., Wolynes P. G. Microscopic theory of protein folding rates. II. Local reaction coordinates and chain dynamics. J. Chem. Phys. 114, 5082–5096 (2001).
Khatri B. S., McLeish T. C. B. Rouse model with internal friction: A coarse grained framework for single biopolymer dynamics. Macromolecules 40, 6770–6777 (2007).
Scott K. A., Randles L. G., Clarke J. The folding of spectrin domains II: phi-value analysis of R16. J. Mol. Biol. 344, 207–221 (2004).
Scott K. A., Clarke J. Spectrin R16: Broad energy barrier or sequential transition states? Prot. Sci. 14, 1617–1629 (2005).
Scott K. A., Randles L. G., Moran S. J., Daggett V., Clarke J. The folding pathway of spectrin R17 from experiment and simulation: using experimentally validated MD simulations to characterize states hinted at by experiment. J. Mol. Biol. 359, 159–173 (2006).
Otzen D. E., Kristensen O., Proctor M., Oliveberg M. Structural changes in the transition state of protein folding: Alternative interpretations of curved chevron plots. Biochemistry 38, 6499–6511 (1999).
Gianni S., Brunori M., Jemth P., Oliveberg M., Zhang M. Distinguishing between smooth and rough free energy barriers in protein folding. Biochemistry 48, 11825–11830 (2009).
Wensley B. G., Kwa L. G., Shammas S. L., Rogers J. M., Clarke J. Protein folding: adding a nucleus to guide helix docking reduces landscape roughness. J. Mol. Biol. 423, 273–283 (2012).
Hänggi P., Talkner P., Borkovec M. Reaction-rate theory—50 years after Kramers. Rev. Mod. Phys. 62, 251–341 (1990).
Zwanzig R. Diffusion in a rough potential. Proc. Natl Acad. Sci. USA 85, 2029–2030 (1988).
Frauenfelder H., Fenimore P. W., Chen G., McMahon B. H. Protein folding is slaved to solvent motions. Proc. Natl Acad. Sci. USA 103, 15469–15472 (2006).
Wensley B. G. et al. Separating the effects of internal friction and transition state energy to explain the slow, frustrated folding of spectrin domains. Proc. Natl Acad. Sci. USA early edition doi:10.1073/pnas.1201793109 (2012).
Abkevich V. I., Gutin A. M., Shakhnovich E. I. Specific nucleus as the transition state for protein folding: evidence from the lattice model. Biochemistry 33, 10026–10036 (1994).
Yu H. et al. Energy landscape analysis of native folding of the prion protein yields the diffusion constant, transition path time, and rates. Proc. Natl Acad. Sci. USA 109, 14452–14457 (2012).
Scott K. A., Batey S., Hooton K. A., Clarke J. The folding of spectrin domains I: wild-type domains have the same stability but very different kinetic properties. J. Mol. Biol. 344, 195–205 (2004).
Acknowledgements
We thank Robert Best, Hagen Hofmann, Dmitrii Makarov, Eugene Shakhnovich, Devarajan Thirumalai and Peter Wolynes for discussion. This work was supported by the Swiss National Science Foundation (B.S.), the Swiss National Centre of Competence in Research for Structural Biology (B.S.), a Starting Grant of the European Research Council (SingleMolFolding; B.S.), Marie Curie Actions (A.B.), the Human Frontier Science Programme (B.S. and E.A.L.), the Wellcome Trust (J.C. and B.G.W., grant number WT095195MA) and the National Science Foundation of the USA (E.A.L.). J.C. is a Wellcome Trust Senior Research Fellow.
Author information
Authors and Affiliations
Contributions
A.B., J.C. and B.S. designed the research. A.B., B.G.W. and M.B.B. prepared protein samples. D.N. developed instrumentation and assisted with data analysis. A.B. and B.G.W. performed the measurements. A.B., B.G.W., A.S., J.C. and B.S. analysed and interpreted the data. S.H.P. and E.A.L. provided the microfluidic mixing system. A.H. helped with microfluidics experiments. A.B., J.C. and B.S. wrote the paper with the help of the other authors.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing financial interests.
Supplementary information
Supplementary Figure
Supplementary Figures S1-S5, Supplementary Table S1, Supplementary Methods and Supplementary References (PDF 4048 kb)
Rights and permissions
This work is licensed under a Creative Commons Attribution-NonCommercial-No Derivative Works 3.0 Unported License. To view a copy of this license, visit http://creativecommons.org/licenses/by-nc-nd/3.0/
About this article
Cite this article
Borgia, A., Wensley, B., Soranno, A. et al. Localizing internal friction along the reaction coordinate of protein folding by combining ensemble and single-molecule fluorescence spectroscopy. Nat Commun 3, 1195 (2012). https://doi.org/10.1038/ncomms2204
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/ncomms2204
This article is cited by
-
Transition path times of coupled folding and binding reveal the formation of an encounter complex
Nature Communications (2018)
-
Hypothesis: structural heterogeneity of the unfolded proteins originating from the coupling of the local clusters and the long-range distance distribution
Biophysical Reviews (2018)
-
Cotranslational folding of spectrin domains via partially structured states
Nature Structural & Molecular Biology (2017)
-
Molecular origins of internal friction effects on protein-folding rates
Nature Communications (2014)
-
Microfluidic mixer designed for performing single-molecule kinetics with confocal detection on timescales from milliseconds to minutes
Nature Protocols (2013)
Comments
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.