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

Originally published as Biophys J. BioFAST on June 30, 2006.
doi:10.1529/biophysj.106.089490
This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow All Versions of this Article:
biophysj.106.089490v1
91/5/L51    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 Google Scholar
Google Scholar
Right arrow Articles by West, D. K.
Right arrow Articles by Paci, E.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by West, D. K.
Right arrow Articles by Paci, E.
Biophysical Journal 91:L51-L53 (2006)
© 2006 The Biophysical Society

Prediction of the Translocation Kinetics of a Protein from Its Mechanical Properties

Daniel K. West * {dagger}, David J. Brockwell {dagger} and Emanuele Paci *

* School of Physics and Astronomy, and {dagger} Institute of Molecular and Cellular Biology, University of Leeds, Leeds LS2 9JT, United Kingdom

Correspondence: Address reprint requests and inquiries to E. Paci, E-mail: e.paci{at}leeds.ac.uk.


    ABSTRACT
 TOP
 ABSTRACT
 REFERENCES
 
Proteins are actively unfolded to pass through narrow channels in macromolecular complexes that catalyze protein translocation and degradation. Catalyzed unfolding shares many features that characterize the mechanical unfolding of proteins using the atomic force microscope (AFM). However, simulations of unfolding induced by the AFM and when a protein is translocated through a pore suggest that each process occurs by distinct pathways. The link, if any, between each type of unfolding, therefore, is not known. We show that the mechanical unfolding energy landscape of a protein, obtained using an atomistic molecular model, can be used to predict both the relative mechanical strength of proteins when unfolded using the AFM and when unfolded by translocation into a pore. We thus link the two processes and show that the import rate through a pore not only depends on the location of the initiation tag but also on the mechanical properties of the protein when averaged over all the possible geometries that are relevant for a given translocation initiation site.

When proteins are actively translocated across the mitochondrial membrane or degraded by chambered proteases (such as the proteasome or the Clp proteases), the protein substrate must first be unfolded as the entry pores of these large macromolecular complexes have a diameter that is much smaller than the size of the folded protein (1Go). It is thought that these complexes catalyze unfolding by the application of a mechanical force onto the protein, thereby decreasing the unfolding barrier and exponentially increasing the unfolding rate constant. Proteins targeted for either translocation or degradation by a suitable polypeptide tag are unfolded significantly faster than in solution. The precise mechanism, however, by which these "unfoldases" are able to denature proteins, remains unresolved. In the last 10 years, many proteins have now been mechanically unfolded using the AFM—a process that shares many features of active unfolding catalyzed by the macromolecular complexes described above. For both types of unfolding it has been shown that i), the unfolding pathway is distinct to the intrinsic unfolding pathway that is probed by chemical denaturation experiments (2Go–4Go); ii), the unfolding rates of substrate proteins do not correlate with their global thermodynamic or kinetic stability (3Go,5Go,6Go); and iii), the mechanical strength is determined by the local structure of the polypeptide chain relative to the initiation point of import or degradation (2Go,3Go,7Go,8Go). Although these results suggest a common mechanism, direct comparison of mechanical strength measured by the atomic force microscope (AFM) and import rates (the mitochondrial TOM and TIM complexes (9Go,10Go)) or degradation rates (ClpXP (11Go)) have proved inconclusive. Furthermore, simulation studies that allow the unfolding pathways of proteins to be observed at the molecular level suggest that the unfolding pathways of barnase, when dragged into a very simple geometrical model of a pore, are different to the pathway observed when the two ends of the protein are pulled apart (12Go). An analogous result was also found for ubiquitin (13Go) using an even simpler model of a pore and a coarse-grained model of the protein. Although the mechanisms between pore translocation and AFM pulling were found to be different, both sets of simulations showed that the local stability of the part of the protein that enters the pore first determines its overall resistance, in qualitative agreement with experimental observations.

What then is the relationship between catalyzed mechanical unfolding and unfolding a protein using the AFM? In AFM, unfolding force is applied locally to two distinct locations within the protein (Fig. 1 a), and in this case the mechanical strength of the protein is thought to be dominated by the stability of the contacts between these regions (e.g., the hydrogen bond mechanical clamp in I27 (4Go)). When translocated through a pore, however, force is applied at a precise location, but the areas of the protein resisting the extension are diffuse (those touching the rim of the pore, Fig. 1 b). Therefore, although force is acting as a highly localized peturbant in both processes, the mechanical stability measured by the unfolding force (in AFM experiments) or unfolding rates (translocation experiments) will depend in some manner on the strength of all of the contacts being perturbed, and these differ in each process. If these processes are both accelerated by force, information about both pathways should be contained in the mechanical unfolding landscape of the protein. We have previously used a C{alpha} native-centric or Go-like model to create such a mechanical landscape for a variety of proteins (14Go) in which the mechanical strength of a protein is simulated when extended by pulling every pair of C{alpha} atoms. We showed that these "mechanical unfolding landscapes" capture the anisotropy of the response of a protein to mechanical peturbation: a protein can either resist a mechanical force or yield easily to it (yellow and blue areas in Fig. 2) by altering the pulling geometry.


Figure 1
View larger version (9K):
[in this window]
[in a new window]
 
FIGURE 1  (a) AFM single molecule experiments can pull pairs of atoms apart under a constant loading force F. By repeating the same experiment a number of times, one may determine the unfolding rate for a given residue pair and force. (b) A simple steric model of a pore. Assuming that the protein unravels as a consequence of being dragged through the pore, a special role will be played by the residues that contact the rim of the pore. Their identity depends not only on the site of the translocation initiation site but also on a possible partial unraveling of the protein in the region of the initiation site.

 

Figure 2
View larger version (63K):
[in this window]
[in a new window]
 
FIGURE 2  Mechanical unfolding landscape of I27 generated by estimating unfolding rates pulling apart all pairs of residues at a constant force of 150 pN, T = 300 K. Yellow to blue colors denote high to weak mechanical resistance; black denotes pairs not pulled. Scale is in ps–1.

 
To investigate the relationship between the mechanical unfolding landscape of a protein and its rate of translocation through a pore, it was first necessary to simulate the latter process. A simple model of a pore was designed using a previously proposed model (13Go) of an infinitely long pore interacting repulsively with the protein. The chosen pore diameter was 12 Å, comparable to the diameter of narrowest constriction in the degradation channel of the proteosome (1Go), which ensures that translocation could only occur with the complete unfolding of the domain. We then performed simulations where a selected residue was centered on the pore, the protein oriented so that the vector joining the selected residue and the center of geometry is parallel to the pore axis and a constant force oriented along this axis applied to the selected residue. The resulting translocation rates are shown in Fig. 3 (black line).


Figure 3
View larger version (17K):
[in this window]
[in a new window]
 
FIGURE 3  Unfolding rates measured when pulling I27 by a single residue through a pore (black line). The translocation rate through the pore is related to the mechanical landscape when averaged over suitable geometries. The blue line is the unfolding time averaged over all possible residue pairs; the red from averaging over all pairs with cij determined as described in the text.

 
Interestingly, the rate of translocation changes by more than two orders of magnitude when the force is applied to different residues. Slowest translocation occurs when the initiation site is close to residues 11, 23, 59, 69, and 83, and I27 is translocated relatively rapidly when dragged through the pore from the N-terminus. This is particularly surprising, since it is has been shown by experiment and simulation that I27 is able to withstand large forces when pulled between the N- and C termini using the AFM. This also contrasts with the conclusion of Wilcox et al. (9Go) that the unfolding force when pulling the N- and C termini correlates in general with the translocation rate when the targeting sequence is attached to the N-terminus.

To predict the translocation rates through a pore from a knowledge of the mechanical properties of the protein, a rough approximation consists in assuming that all possible pulling geometries between that residue and other residues contribute equally; i.e., that all mechanical pathways contribute in parallel, which gives the rate

Formula 1(1)
where Formula 1 is the unfolding rate constant for a residue pair i,j. The resulting Formula 1 is shown as a blue line in Fig. 3. Although some features are shared between these data and the pulling through a pore simulation, the correlation between log Formula 1 and log Formula 1 is only 0.48. The primary reason for this moderate correlation comes from the fact that if a single unfolding rate kij is much larger than all the others, it will dominate the sum in Eq. 1. A better approximation consists in considering only those residues that are in contact with the rim of the pore when the force is applied to residue i. Their identity cannot be exactly determined, as the conformation which corresponds to the rate limiting step in the unfolding process (the transition state) is not necessarily the native one (although for mechanically robust proteins we do not expect it to be very topologically different from the native state). If one assumes that the mechanical unfolding transition state is close to the native state, then a reasonable assumption is to consider residues which are close in space to the initiation site and not buried inside the protein. An effective rate is then obtained by summing only over the relevant pairs of residues:

Formula 2(2)
where Formula 2 if residue j is at most 9 Å apart from residue i (slightly larger than the pore radius) and at least 40% of its surface is solvent-accessible in the native state (since buried residues are unlikely to get stuck on the outer rim of the pore) and that are at least 10 residues apart along the polypeptide chain (it makes little sense to define an unfolding rate between residues close in the chain as their extensions might not lead to the unfolding of the protein). The resulting unfolding rate constant (Formula 2) profile calculated in this manner (red line, Fig. 3) agrees more closely with Formula 2 (the correlation increases to 0.68).

One obvious difference between the rate constant profiles between simulated translocation through a pore and that obtained from the mechanical energy landscape is that the latter are about one order of magnitude larger than the former. This is simply due to the fact that the effective constant force F acting on each of the M residues, which are stuck on the rim (see Fig. 1 a), is ~1/M the force acting on each pair when the same constant force F is used to determine Formula 2.

Although the translocation rate profile predicted from the mechanical map and directly estimated from the pore model is similar in the case of I27, this method may not be reliable for all proteins. We performed analogous estimations for E2Lip3 and barnase (both at a constant force of 100 pN) and the best correlation between Formula 2 and Formula 2, which could be found was 0.53 and 0.48, respectively. In particular, for proteins whose mechanical transition state is far from the native state and may have possible unfolding intermediates, this approximation is less accurate. In this case, any partial unraveling of the protein into the pore will make the prediction less accurate. One other factor that might hinder a reliable prediction of the translocation rate from mechanical properties alone is the fact that the geometry of the protein relative to the pore might differ from that probed by the analogous AFM experiments on the various pairs of atoms that might be involved.

Can we predict the translocation rates through a pore from the knowledge of the mechanical properties of the protein? The answer is yes, but it depends on the full mechanical properties of the protein and not only on those probed by extending the protein from two points as measured in AFM experiments. The results presented here could be experimentally tested by measuring the rate at which the ClpXP protease degrades I27 variants in which an ssrA degradation peptide tag is attached to the side chain of unique cysteine residues inserted at specific points (15Go).

Submitted on May 18, 2006; accepted for publication June 16, 2006.


    REFERENCES
 TOP
 ABSTRACT
 REFERENCES
 
1. Prakash, S., and A. Matouschek. 2004. Protein unfolding in the cell. Trends Biochem. Sci. 29:593–600.[CrossRef][Medline]

2. Huang, S. H., K. S. Ratliff, M. P. Schwartz, J. M. Spenner, and A. Matouschek. 1999. Mitochondria unfold precursor proteins by unraveling them from their N-termini. Nat. Struct. Biol. 6:1132–1138.[CrossRef][Medline]

3. Lee, C., M. P. Schwartz, S. Prakash, M. Iwakura, and A. Matouschek. 2001. ATP-dependent proteases degrade their substrates by processively unraveling them from the degradation signal. Mol. Cell. 7:627–637.[CrossRef][Medline]

4. Best, R. B., S. B. Fowler, J. L. T. Herrera, A. Steward, E. Paci, and J. Clarke. 2003. Mechanical unfolding of a titin Ig domain: Structure of transition state revealed by combining atomic force microscopy, protein engineering and molecular dynamics simulations. J. Mol. Biol. 330:867–877.[CrossRef][Medline]

5. Lim, J. H., F. Martin, B. Guiard, N. Pfanner, and W. Voos. 2001. The mitochondrial Hsp70-dependent import system actively unfolds preproteins and shortens the lag phase of translocation. EMBO J. 20:941–950.[CrossRef][Medline]

6. Li, H. B., A. F. Oberhauser, S. B. Fowler, J. Clarke, and J. M. Fernandez. 2000. Atomic force microscopy reveals the mechanical design of a modular protein. Proc. Natl. Acad. Sci. USA. 97:6527–6531.[Abstract/Free Full Text]

7. Brockwell, D. J., E. Paci, R. C. Zinober, G. S. Beddard, P. D. Olmsted, D. A. Smith, R. N. Perham, and S. E. Radford. 2003. Pulling geometry defines the mechanical resistance of a beta-sheet protein. Nat. Struct. Biol. 10:731–737.[CrossRef][Medline]

8. Carrion-Vazquez, M., H. B. Li, H. Lu, P. E. Marszalek, A. F. Oberhauser, and J. M. Fernandez. 2003. The mechanical stability of ubiquitin is linkage dependent. Nat. Struct. Biol. 10:738–743.[CrossRef][Medline]

9. Wilcox, A. J., J. Choy, C. Bustamante, and A. Matouschek. 2005. Effect of protein structure on mitochondrial import. Proc. Natl. Acad. Sci. USA. 102:15435–15440.[Abstract/Free Full Text]

10. Sato, T., M. Esaki, J. M. Fernandez, and T. Endo. 2005. Comparison of the protein-unfolding pathways between mitochondrial protein import and atomic-force microscopy measurements. Proc. Natl. Acad. Sci. USA. 102:17999–18004.[Abstract/Free Full Text]

11. Kenniston, J. A., T. A. Baker, J. M. Fernandez, and R. T. Sauer. 2003. Linkage between ATP consumption and mechanical unfolding during the protein processing reactions of an AAA(+) degradation machine. Cell. 114:511–520.[CrossRef][Medline]

12. Tian, P., and I. Andricioaei. 2005. Repetitive pulling catalyzes co-translocational unfolding of barnase during import through a mitochondrial pore. J. Mol. Biol. 350:1017–1034.[CrossRef][Medline]

13. Huang, L., S. Kirmizialtin, and D. E. Makarov. 2005. Computer simulations of the translocation and unfolding of a protein pulled mechanically through a pore. J. Chem. Phys. 123:126903.

14. West, D. K., D. J. Brockwell, P. D. Olmsted, S. E. Radford, and E. Paci. 2006. Mechanical resistance of proteins explained using simple molecular models. Biophys. J. 90:287–297.[Abstract/Free Full Text]

15. Kenniston, J. A., R. E. Burton, S. M. Siddiqui, T. A. Baker, and R. T. Sauer. 2004. Effects of local protein stability and the geometric position of the substrate degradation tag on the efficiency of ClpXP denaturation and degradation. J. Struct. Biol. 146:130–140.[CrossRef][Medline]





This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow All Versions of this Article:
biophysj.106.089490v1
91/5/L51    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 Google Scholar
Google Scholar
Right arrow Articles by West, D. K.
Right arrow Articles by Paci, E.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by West, D. K.
Right arrow Articles by Paci, E.


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