help button home button Biophys. J.
HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS

Originally published as Biophys J. BioFAST on December 16, 2005.
doi:10.1529/biophysj.105.077982
This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Supplemental
Right arrow All Versions of this Article:
biophysj.105.077982v1
90/4/L33    most recent
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Schlierf, M.
Right arrow Articles by Rief, M.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Schlierf, M.
Right arrow Articles by Rief, M.
Biophysical Journal 90:L33-L35 (2006)
© 2006 The Biophysical Society

Single-Molecule Unfolding Force Distributions Reveal a Funnel-Shaped Energy Landscape

Michael Schlierf and Matthias Rief

Physik-Department E22, Technische Universität München, D-85748 Garching b. Munich, Germany

Correspondence: Address reprint requests and inquiries to Matthias Rief, E-mail: mrief{at}ph.tum.de.


    ABSTRACT
 TOP
 ABSTRACT
 SUPPLEMENTARY MATERIAL
 ACKNOWLEDGEMENTS
 REFERENCES
 
The protein folding process is described as diffusion on a high-dimensional energy landscape. Experimental data showing details of the underlying energy surface are essential to understanding folding. So far in single-molecule mechanical unfolding experiments a simplified model assuming a force-independent transition state has been used to extract such information. Here we show that this so-called Bell model, although fitting well to force velocity data, fails to reproduce full unfolding force distributions. We show that by applying Kramers' diffusion model, we were able to reconstruct a detailed funnel-like curvature of the underlying energy landscape and establish full agreement with the data. We demonstrate that obtaining spatially resolved details of the unfolding energy landscape from mechanical single-molecule protein unfolding experiments requires models that go beyond the Bell model.

In the "new view" of protein folding a funnel-shaped energy landscape steers the folding polypeptide chain toward its native folded conformation and thus facilitates the folding problem. This view has emerged from a host of bulk studies and theoretical work (1Go). Recently single-molecule mechanical techniques have offered the possibility to study the stabilizing forces in protein folding directly (2Go). Up to now all information about the energy landscape extracted from single-molecule mechanical unfolding experiments has been limited to a force-independent barrier position and an unloaded unfolding rate. Although in the field of receptor-ligand mechanics multiple barriers have been reported (3Go,4Go), for forced protein unfolding such effects have never been observed. Experimental limitations, like calibration errors and ill-defined loading rate conditions, have precluded a more detailed insight into the barrier shape. In this letter we demonstrate that the commonly used Bell model to analyze velocity-dependent unfolding force data fails to describe force distributions and yields a misleading view of the unfolding energy landscape.

We investigated unfolding force distributions at various pulling velocities of the Immunoglobulin-like domain 4 (ddFLN4) from the Dictyostelium discoideum F-actin cross-linker filamin (ddFLN) (see also inset in Fig. 1 A). We chose this domain because it offers the possibility to perform experiments under well-defined loading rate conditions in repetitive unfolding-refolding cycles with a truly single domain. The details of these double-jump experiments and the applied methods for obtaining force distributions with ddFLN4 have been described in detail elsewhere (5Go,6Go). It is important to note that for a detailed analysis of force distributions, well-defined loading rate conditions through a constant spacer length are essential. Moreover, it is crucial to ensure that no events get missed in the thermal noise floor.


Figure 1
View larger version (32K):
[in this window]
[in a new window]
 
FIGURE 1  A comparison between the commonly used Bell view and Kramers' view for analysis of protein unfolding data. (A) The pulling velocity data of the native-state unfolding of ddFLN4 show a logarithmic behavior within the experimental error. The inset shows a schematic illustration of the experimental setup. (B) Schematic reconstruction of the energy landscape in the Bell view along the NC-terminal vector. (C) Normalized unfolding probability force distributions (histogram with statistical error) at four different pulling velocities in comparison with the theoretical distributions using the Bell model (lines). (D) Pulling velocity data modeled using Kramers' theory. (E) The reconstructed energy landscape shows detailed curvature along the unfolding/folding pathway. (F) The characteristic behavior of the experimental unfolding force distributions with increasing pulling velocity is reproduced well using the Kramers' model.

 
There are two independent ways of assessing the underlying energy landscape in single-molecule mechanical unfolding experiments. Commonly, the most probable unfolding force is plotted versus the pulling velocity. Typically a logarithmic dependence is observed. Force versus velocity for ddFLN4 unfolding is plotted in Fig. 1 A. The data appear almost perfectly logarithmic. A strict logarithmic dependence has usually served as justification to apply the so-called Bell model, a simple two-state model with a triangular energy profile (Fig. 1 B) exhibiting a force-independent transition state position (7Go). According to this model, force modulates the observed unfolding rate following

Formula
where koff(0) denotes the unloaded unfolding rate, {Delta}x the transition state position, and kBT the thermal energy (7Go). From this force-dependent unfolding rate the probability distributions of rupture forces as well as the most probable unfolding force as a function of the pulling velocity can now be calculated analytically even if the force is applied through a nonlinear polymer spacer chain ((8Go); also Supplementary Material). Our data for the most probable unfolding forces in Fig. 1 A agree perfectly with a fit to the Bell model (dotted line in Fig. 1 A) with koff(0) = 0.35 s–1 and {Delta}x = 5 Å.

However, when trying to fit the full probability distributions for the unfolding force at the various pulling velocities with the Bell model using the above parameters, a strong disagreement between data and Bell model becomes evident (Fig. 1 C). At low pulling velocities the distributions according to the Bell model are much wider, exhibiting a higher fraction of low force events than our data (between 0 and 20 pN, ~25% of all events at 200 nm/s). At this point it is extremely important to be certain about the number of missed low-force events in the experiment. Because our experiment involves repeated unfolding-refolding cycles of a single ddFLN4 domain with well-studied refolding kinetics we can estimate the number of missed events to be <5% (for details see Supplementary Material). In contrast, the probability distribution at high pulling velocities (4000 nm/s) as calculated by the Bell model appears too narrow compared to the data. The changing width of the measured force distributions is a strong indication for a transition state position that moves with force.

How can we resolve the puzzling disagreement between an apparently logarithmic force-velocity behavior and yet a moving transition-state position? Obviously our data call for a more realistic model to describe unfolding force data. Previous theoretical studies have shown that the adaptation of Kramers' diffusion model to forced unbinding can give important information about the underlying energy landscape (9Go,10Go). Following Shillcock and Seifert (11Go) we calculated the mean first passage time for a particle in a given potential at increasing external forces according to:

Formula
where D represents the diffusion constant of the particle and UF(x) represents the energy profile along the pulling direction at force F (for details see Supplementary Material) (11Go). We optimized the energy profile UF(x) shown in Fig. 1 E so that the full distributions were reproduced. It is interesting to note that although the transition state of the potential UF(x) is strongly force dependent the predicted force versus pulling velocity curve is still in perfect agreement with our data (Fig. 1 D). Although a slight curvature of the force-velocity curve can be observed over a large range of four decades, experimental data would still appear almost logarithmic. Moreover, without any additional fit parameter Kramers' theory also reproduces the full experimental probability distributions (Fig. 1 F).

The optimized energy profile has a funnel shape. Broad transition states have been observed in bulk studies, albeit along a very different reaction coordinate (12Go). Such a broad transition state has been suggested to actively steer the protein to its minimum from distant conformations.

Our results demonstrate that merely analyzing unfolding force versus velocity curves may yield an oversimplified and even sometimes misleading picture of the underlying energy landscape. It is also important that not only the transition state position as obtained by the Bell model is drastically smaller (5 vs. 11 Å) but also koff(0) is considerably too high (0.35 s–1 vs. 0.01 s–1). This has important consequences when discussing koff(0) values extracted from force-velocity data analyzed using the Bell model.

To conclude, we could demonstrate that high-resolution recordings of full probability distributions are essential for characterizing in detail the energy landscape governing forced unfolding of proteins.


    SUPPLEMENTARY MATERIAL
 TOP
 ABSTRACT
 SUPPLEMENTARY MATERIAL
 ACKNOWLEDGEMENTS
 REFERENCES
 
An online supplement to this article can be found by visiting BJ Online at http://www.biophysj.org.


    ACKNOWLEDGEMENTS
 TOP
 ABSTRACT
 SUPPLEMENTARY MATERIAL
 ACKNOWLEDGEMENTS
 REFERENCES
 
We thank M. Schleicher and A. A. Noegel for providing the ddFLN1-5 construct.

This project was funded by a SFB413 grant of the Deutsche Forschungsgemeinschaft. M.S. was supported by International Graduate School-Navo Bio Technology (IDK-NBT).

Submitted on November 17, 2005; accepted for publication December 14, 2005.


    REFERENCES
 TOP
 ABSTRACT
 SUPPLEMENTARY MATERIAL
 ACKNOWLEDGEMENTS
 REFERENCES
 
1. Onuchic, J. N., and P. G. Wolynes. 2004. Theory of protein folding. Curr. Opin. Struct. Biol. 14:70–75.[CrossRef][Medline]

2. Carrion-Vazquez, M., A. F. Oberhauser, T. E. Fisher, P. E. Marszalek, H. Li, and J. M. Fernandez. 2000. Mechanical design of proteins studied by single-molecule force spectroscopy and protein engineering. Prog. Biophys. Mol. Biol. 74:63–91.[CrossRef][Medline]

3. Merkel, R., P. Nassoy, A. Leung, K. Ritchie, and E. Evans. 1999. Energy landscapes of receptor-ligand bonds explored with dynamic force spectroscopy. Nature. 397:50–53.[CrossRef][Medline]

4. Pincet, F., and J. Husson. 2005. The solution to the streptavidin-biotin paradox: the influence of history on the strength of single molecular bonds. Biophys. J. 89:4374–81.[Abstract/Free Full Text]

5. Schwaiger, I., M. Schleicher, A. A. Noegel, and M. Rief. 2005. The folding pathway of a fast-folding immunoglobulin domain revealed by single-molecule mechanical experiments. EMBO Rep. 6:46–51.[CrossRef][Medline]

6. Schlierf, M., and M. Rief. 2005. Temperature softening of a protein in single-molecule experiments. J. Mol. Biol. 354:497–503.[CrossRef][Medline]

7. Bell, G. I. 1978. Models for the specific adhesion of cells to cells. Science. 200:618–627.[Abstract/Free Full Text]

8. Evans, E., and K. Ritchie. 1999. Strength of a weak bond connecting flexible polymer chains. Biophys. J. 76:2439–2447.[Abstract/Free Full Text]

9. Evans, E., and K. Ritchie. 1997. Dynamic strength of molecular adhesion bonds. Biophys. J. 72:1541–1555.[Abstract/Free Full Text]

10. Heymann, B., and H. Grubmuller. 2000. Dynamic force spectroscopy of molecular adhesion bonds. Phys. Rev. Lett. 84:6126–6129.[CrossRef][Medline]

11. Shillcock, J., and U. Seifert. 1998. Escape from a metastable well under a time-ramped force. Phys. Rev. E. 57:7301–7304.[CrossRef]

12. Silow, M., and M. Oliveberg. 1997. High-energy channeling in protein folding. Biochemistry. 36:7633–7637.[CrossRef][Medline]




This article has been cited by other articles:


Home page
Biophys. JHome page
S. Guo, C. Ray, A. Kirkpatrick, N. Lad, and B. B. Akhremitchev
Effects of Multiple-Bond Ruptures on Kinetic Parameters Extracted from Force Spectroscopy Measurements: Revisiting Biotin-Streptavidin Interactions
Biophys. J., October 15, 2008; 95(8): 3964 - 3976.
[Abstract] [Full Text] [PDF]


Home page
Biophys. JHome page
E. Eyal and I. Bahar
Toward a Molecular Understanding of the Anisotropic Response of Proteins to External Forces: Insights from Elastic Network Models
Biophys. J., May 1, 2008; 94(9): 3424 - 3435.
[Abstract] [Full Text] [PDF]


Home page
Biophys. JHome page
J. I. Sulkowska and M. Cieplak
Stretching to Understand Proteins A Survey of the Protein Data Bank
Biophys. J., January 1, 2008; 94(1): 6 - 13.
[Abstract] [Full Text] [PDF]


Home page
Proc. Natl. Acad. Sci. USAHome page
R. A. Nome, J. M. Zhao, W. D. Hoff, and N. F. Scherer
Axis-dependent anisotropy in protein unfolding from integrated nonequilibrium single-molecule experiments, analysis, and simulation
PNAS, December 26, 2007; 104(52): 20799 - 20804.
[Abstract] [Full Text] [PDF]


Home page
Biophys. JHome page
M. Schlierf, F. Berkemeier, and M. Rief
Direct Observation of Active Protein Folding Using Lock-in Force Spectroscopy
Biophys. J., December 1, 2007; 93(11): 3989 - 3998.
[Abstract] [Full Text] [PDF]


Home page
Biophys. JHome page
J. Morfill, K. Blank, C. Zahnd, B. Luginbuhl, F. Kuhner, K.-E. Gottschalk, A. Pluckthun, and H. E. Gaub
Affinity-Matured Recombinant Antibody Fragments Analyzed by Single-Molecule Force Spectroscopy
Biophys. J., November 15, 2007; 93(10): 3583 - 3590.
[Abstract] [Full Text] [PDF]


Home page
Biophys. JHome page
M. S. Li
Secondary Structure, Mechanical Stability, and Location of Transition State of Proteins
Biophys. J., October 15, 2007; 93(8): 2644 - 2654.
[Abstract] [Full Text] [PDF]


Home page
Biophys. JHome page
M. Odorico, J.-M. Teulon, T. Bessou, C. Vidaud, L. Bellanger, S.-w. W. Chen, E. Quemeneur, P. Parot, and J.-L. Pellequer
Energy Landscape of Chelated Uranyl: Antibody Interactions by Dynamic Force Spectroscopy
Biophys. J., July 15, 2007; 93(2): 645 - 654.
[Abstract] [Full Text] [PDF]


Home page
Biophys. JHome page
D. E. Makarov
Unraveling Individual Molecules by Mechanical Forces: Theory Meets Experiment
Biophys. J., June 15, 2007; 92(12): 4135 - 4136.
[Full Text] [PDF]


This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Supplemental
Right arrow All Versions of this Article:
biophysj.105.077982v1
90/4/L33    most recent
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Schlierf, M.
Right arrow Articles by Rief, M.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Schlierf, M.
Right arrow Articles by Rief, M.


HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
Copyright © 2006 by the Biophysical Society.