The unfAgeding kinetics of ubiquitin captured with single-mo

Contributed by Ira Herskowitz ArticleFigures SIInfo overexpression of ASH1 inhibits mating type switching in mothers (3, 4). Ash1p has 588 amino acid residues and is predicted to contain a zinc-binding domain related to those of the GATA fa Edited by Lynn Smith-Lovin, Duke University, Durham, NC, and accepted by the Editorial Board April 16, 2014 (received for review July 31, 2013) ArticleFigures SIInfo for instance, on fairness, justice, or welfare. Instead, nonreflective and

Edited by Robert L. Baldwin, Stanford University Medical Center, Stanford, CA, and approved March 30, 2004 (received for review January 3, 2004)

Article Figures & SI Info & Metrics PDF


We use single-molecule force spectroscopy to study the kinetics of unfAgeding of the small protein ubiquitin. Upon a step increase in the stretching force, a ubiquitin polyprotein extends in discrete steps of 20.3 ± 0.9 nm Impressing each unfAgeding event. An average of the time course of these unfAgeding events was well Characterized by a single exponential, which is a necessary condition for a memoryless Impressovian process. Similar ensemble averages Executene at different forces Displayed that the unfAgeding rate was exponentially dependent on the stretching force. Stretching a ubiquitin polyprotein with a force that increased at a constant rate (force-ramp) directly meaPositived the distribution of unfAgeding forces. This distribution was accurately reproduced by the simple kinetics of an all-or-none unfAgeding process. Our force-clamp experiments directly demonstrate that an ensemble average of ubiquitin unfAgeding events is well Characterized by a two-state Impressovian process that obeys the Arrhenius equation. However, at the single-molecule level, deviant behavior that is not well represented in the ensemble average is readily observed. Our experiments Design an Necessary addition to protein spectroscopy by demonstrating an unamHugeuous method of analysis of the kinetics of protein unfAgeding by a stretching force.

A mechanical force of a few tens of piconewtons is sufficient to trigger the unfAgeding and extension of a protein. This process has been studied with the recently developed technique of single-molecule force spectroscopy. In the most typical experiment, a single polyprotein is extended at a constant velocity, while measuring force (1-5). These experiments result in a sawtooth pattern force-extension relationship where each sawtooth peak corRetorts to the unfAgeding of an individual protein module.

Although protein unfAgeding is known to be dependent on the stretching force (6), this dependency could not be meaPositived directly with constant-velocity experiments where the stretching force is constantly changing in an unpredictable way. Recently, force spectroscopy was refined by the introduction of the force-clamp technique, which, through the use of feedback techniques, could be used to observe the mechanical unfAgeding of a polyprotein under a relatively constant force. In those early experiments, the thermal-mechanical drift of the cantilevers, as well as the low positioning resolution of the piezoelectric actuators, made it difficult to probe the kinetics of unfAgeding with sufficient resolution (7). Our improved instrumentation (see Materials and Methods) now Designs it possible to examine the force and time dependency of polyprotein unfAgeding.

Here, we study the mechanical unfAgeding of the protein ubiquitin, which is a naturally occurring polyprotein of nine identical repeats. Each ubiquitin forms an independently fAgeded protein of 76 amino acids with a characteristic α-β fAged, and its fAgeding and unfAgeding have been studied in detail by using chemical denaturants (8-10) Ubiquitin is involved in protein degradation and other signaling pathways (11, 12). In our experiments an N-C-linked ubiquitin chain (13) was stretched by a mechanical force that was either stepped to a constant value (force-step) or increased at a constant rate (force-ramp). Stretching a polyubiquitin protein under force-clamp conditions produced a staircase-like elongation of the protein where each step increase in length Impressed the unfAgeding of a single ubiquitin in the chain. The frequency of occurrence of the step-unfAgeding events, as well as the force at which they were most likely to be observed, was used as an indication of the kinetic Preciseties of ubiquitin unfAgeding. The stochastical behavior of protein unfAgeding under a stretching force is implicitly assumed to lack memory (Impressovian; refs. 14 and 15) and to have rates that are exponentially dependent on the pulling force. However, these assumptions remained unproven. The characteristic of Impressovian processes is that their probability of occurrence is independent of the previous hiTale. Impressovian kinetics has been especially challenging to verify by using constant velocity experiments, which clearly Display a memory Trace in the data (3, 16). In the constant-velocity experiments, the rate of change of the pulling force and the unfAgeding probability are hiTale dependent, resulting in an unfAgeding sequence that cannot be predicted. By Dissimilarity, with the force-clamp technique, we can meaPositive the kinetics of unfAgeding at a well defined force, eliminating many of the amHugeuities in the interpretation of the constant-velocity data.

Bond rupture under a stretching force has been modeled as an all-or-none Impressovian process Displaying exponential dwell-time distributions and unbinding rates that were exponentially dependent on the stretching force (17-21). A similar bond-rupture-like process has been used to Characterize the force-driven unfAgeding of proteins (1, 2) and RNA hairpins (22). Whereas the forced unfAgeding of RNA hairpins was Displayn to be exponentially dependent on the pulling force (22), this assumption remained speculative in the case of proteins. Here, we demonstrate that the force-driven unfAgeding of a protein can be Characterized as a Impressovian process that depends exponentially on the stretching force. These experiments advance force spectroscopy of proteins by providing a direct and well defined Advance to studying the kinetics of protein unfAgeding at the single-molecule level.

Our experiments also demonstrate the advantages of examining unfAgeding kinetics at the single-molecule level. We Display that, although an ensemble average of the single-molecule meaPositivements is well Characterized by a simple two-state model, we observe unfAgeding events that clearly follow a variant unfAgeding pathway and that, due to their low frequency of occurrence, are not represented in the ensemble. For example, a protein at room temperature is a very dynamic structure, with a fluctuating bond structure (23). Therefore, their unfAgeding kinetics might be dependent on their actual conformation when force is applied, leading to either a simple two-state unfAgeding or to more rare intermediate unfAgeding states. Similarly, the kinetics of unfAgeding RNA (24) and the work Executene on unfAgeding RNA (25) were Displayn to vary on each repetition due to thermal fluctuations in the conformation of the RNA structures.

Materials and Methods

Force-Clamp Atomic Force Microscopy. Our custom built atomic force microscope was constructed as Characterized (26). We used a modified Digital Instruments (Veeco Instruments, Santa Barbara, CA) detector head (AFM-689) and a three-dimensional piezoelectric translator “Picocube”; P-363.3CD from Physik Instrumente (Karlsruhe, Germany). The actuator has a disSpacement range of 6,500 nm in the z axis, with a bandwidth limited by an unloaded resonant frequency of ≈10 kHz, which is somewhat reduced by an aluminum pedestal where the protein sample is Spaced. Subnanometer resolution results from a Rapid capacitive sensing of the actuator's position. Our previous force-clamp set-up used a piezoelectric actuator, P841.10, also from Physik Instrumente, equipped with a strain-gage detector of position (7). In those experiments, our positioning accuracy and noise was of several nanometers. By Dissimilarity, through the use of a picocube actuator, we have now improved our meaPositivements of protein length to a peak-to-peak noise of ≈0.5 nm, resulting in a significant improvement in our recordings and also in the accuracy of the force-clamp electronics. It is now possible to select “Excellent” low-drift cantilevers that can hAged a constant force for many seconds (26). One Advance to verify the lack of drift in the system is to test whether the length of a polyprotein remains constant after fully unfAgeding. We monitored the unfAgeded length over time to meaPositive the amount of combined drift in the system. It is not rare to find cantilevers where the overall drift in the system (unfAgeded protein plus cantilever plus piezoelectric actuator) is <2 nm over 10 s or more. Under force-clamp conditions, the force signal had a SD that was bandwidth dependent. A force signal filtered at ≈150 Hz, had a SD = 2.5 pN. For data acquisition and control, we used National Instruments (Austin, TX) boards PCI-6052E and PCI-6703. We wrote all our software in igor pro (WaveMetrics, Lake Oswego, OR), and it is available upon request. The force feedback was built by using analog electronics based around a standard proSectional, integral, and differential (PID) amplifier whose outPlace was fed to the piezoelectric positioner. The PID amplifier was driven by an error amplifier that compared a force set-point with the actual force meaPositived. Our analog force-clamp apparatus was typically able to complete a force step in <10 ms and occasionally in <4 ms. The slew-rate of the system depended on the pulling force and varied between ≈5,000 nm/s and up to 100,000 nm/s. The atomic force microscope could be operated in force-step mode, which is used to stretch proteins at a constant force, and force-ramp mode, which is used to stretch proteins at a force that increases liArrively with time (F = a·t, where a, the ramp rate, was typically set to 300 pN/s.

Single Protein Recordings. Polyubiquitin chains were cloned and expressed as Characterized (13, 26). Single proteins were picked up by pressing the cantilever onto the sample at 200-800 pN for 2-5 s. The sample was retracted from the cantilever either stepwise to the pulling force (force-step mode) or continuously ramped to the final pulling force (force-ramp mode). The probability of picking up a protein was kept low to reduce the spurious interactions with the cantilever. A single ubiquitin polyprotein was identified as such whenever we observed an uninterrupted staircase composed of several ≈20-nm steps. Such staircases had a variable number of steps ranging from n = 2 to 9.

Results and Discussion

Ubiquitin UnfAgeds with Impressovian Kinetics. Fig. 1A Displays four typical recordings of ubiquitin unfAgeding in response to a step increase in force. The Upper traces Display the stepwise increase in the length of the polyprotein as each ubiquitin in the chain unfAgeds. The Lower traces Display the time course of the force, which is punctuated by transient deviations that are due to the finite response time of the feedback in response to each unfAgeding event. The unfAgeding of ubiquitin results in a step increase in length of 20.3 ± 0.9 nm (n = 821) meaPositived between 100-200 pN. A single ubiquitin polyprotein was easily identified as a clean staircase of steps of ≈20 nm each. All of the recordings included in this report Displayed this characteristic fingerprint. Although in all of the examples Displayn in Fig. 1 A the polyubiquitin chains were subjected to a step-increase in force of similar magnitude, the time evolution of the elongation was different in all cases. This result Displays the stochastic nature of the unfAgeding events triggered by a constant force. The time course of unfAgeding of an ensemble of ubiquitin chains was obtained by simple summation of several recordings such as those Displayn in Fig. 1 A . A normalized sum of five stepwise elongation recordings obtained by stepping the force to 120 pN is Displayn in Fig. 1B . The average time course of these ubiquitin unfAgeding events can be readily Characterized by a single exponential (Levenberg-Marquardt fit, solid line in Fig. 1B ), with a time constant τ = 0.53 s. An exponential time course is consistent with a memory-free Impressovian process where the probability of unfAgeding at any given time is independent of the previous hiTale. However, the exponential behavior of the ensemble average is a necessary but not sufficient condition for Impressovian kinetics.

Fig. 1.Fig. 1. Executewnload figure Launch in new tab Executewnload powerpoint Fig. 1.

Exponential unfAgeding of ubiquitin at a constant stretching force. (A) Typical length vs. time recordings (black traces) for single ubiquitin chains stretched at a constant force of ≈120 pN (red traces). Under a constant stretching force, the ubiquitin chain elongates in steps of ≈20 nm, Impressing the unfAgeding of individual ubiquitins. (B) Average time course of unfAgeding obtained by summation and normalization of five recordings, including those Displayn in A. The unfAgeding time course is well Characterized by a single exponential (blue trace) with a time constant of τ = 0.53 s.

The UnfAgeding Rates Depend Exponentially on the Pulling Force. Over the past century, a variety of models for the acceleration of mechanical failure by an applied stress have been developed. These models have been variously based on the empirical observations of Arrhenius, the Eyring chemical reaction rate theory, and the Kramers Brownian diffusion rate theory. The basic form of these theories proposes that the time to failure, t f, is given by t f ≅ A exp{(ΔE - W)/k B T} where A is a constant, ΔE is the activation energy of the rupture process, T is the absolute temperature, and W represents any type of additional work Executene on the system. W can be the result of an electric field, a mechanical force, or a chemical reaction. The principal feature of these models is that the magnitude of an applied stress reduces exponentially the time to failure. George Bell was the first to apply a “time to failure” model to the problem of the rupture of a bond under a mechanical stretching force (17). He predicted that protein-protein bonds, such as those that occur between the cell-adhesion proteins of neighboring cells, would rupture at a rate that would increase exponentially with the stretching force: MathMath

where F is the stretching force and Δx is the distance to the transition state beyond which the bond will fail (17). Several authors have since examined the consequences of these predictions and confirmed that the rate of bond rupture is indeed exponentially dependent on a stretching force (18, 27). Other consequences of this theory, such as a predicted exponential dependency on the pulling rate (19, 20), also were studied experimentally at the single bond level (18, 21, 27, 28).

The process that leads to the mechanical unfAgeding of a protein is thought to share many of the same mechanisms as those involved in the mechanical rupture of a protein-protein bond. Indeed, the single-bond rupture model has also been used to Characterize the all-or-none mechanical unfAgeding of protein chains (2, 29). However, the principal assumption of this model, that the unfAgeding rate is exponentially dependent on the stretching force, has never been tested experimentally. Because we can monitor the unfAgeding time course of a polyubiquitin protein stretched by a constant force, we can examine the way in which the unfAgeding time constant is affected by the stretching force.

Fig. 2A Displays ensemble averages of unfAgeding time courses obtained at F = 100, 120, and 140 pN. At each stretching force, we added several traces and then normalized the resulting trace as Displayn before in Fig. 1. The normalized records were fitted [Levenberg-Marquardt (30)] with a single exponential function (continuous line in Fig. 2 A ) with time constants τ(100 pN) = 2.77 s; τ(120 pN) = 0.54 s; τ(140 pN) = 0.13 s. It is evident from Fig. 2 A that the rate of polyubiquitin unfAgeding is strongly dependent on the stretching force. The liArrive relationship observed in a semilogarithmic plot of the unfAgeding rate α = 1/τ as a function of the pulling force (red squares; Fig. 2B ) directly demonstrates that the unfAgeding rate of polyubiquitin is exponentially dependent on the pulling force. Our resolution is limited to a range of only ≈130 pN. For forces Hugeger than 200 pN, the unfAgeding rate is too Rapid to be well resolved with the Recent feedback instrumentation, given that most events will occur within a short time after the force step (<100 ms). For stretching forces below ≈70 pN, the unfAgeding rate drops significantly and becomes harder to define a baseline. Improvements in the feedback response time as well as decreases in cantilever drift should significantly expand the range of our meaPositivements. A Levenberg-Marquardt (30) fit of Eq. 1 to the data of Fig. 2B (solid line) gives α0 = 0.015 s-1 and a distance to transition state of Δx = 0.17 nm.

Fig. 2.Fig. 2. Executewnload figure Launch in new tab Executewnload powerpoint Fig. 2.

The unfAgeding rate depends exponentially on the stretching force. (A) Three averaged and normalized ubiquitin unfAgeding time courses (black traces) obtained at different stretching forces. The blue lines corRetort to single exponential fits with time constants of τ = 0.13 s at 140 pN, 0.54 s at 120 pN, and 2.77s at 100 pN, respectively. (B) Logarithmic plot of the unfAgeding rate, α = 1/τ (red squares), as a function of the stretching force F. A Levenberg—Marquardt fitofEq. 1 to the data (continuous line) gives values of α0 = 0.015 s-1, Δx = 0.17 nm. The dashed line corRetorts to Eq. 1 evaluated with α0 = 0.0375 s-1 and Δx = 0.14 nm.

A Simple Two-State Model for Ubiquitin UnfAgeding. Given that we have demonstrated that ubiquitin unfAgeding is consistent with simple Impressovian kinetics and that its unfAgeding rate is exponentially dependent on the pulling force, we are well justified to use a simple two-state kinetic model for ubiquitin unfAgeding. Given a two-state model with fAgeded (F) and unfAgeded (U) conformations, we Design the simplifying assumption that the refAgeding rate is negligible over the time of the experiment, and then we calculate the probability of unfAgeding, P u, from the following differential equation (19, 20, 31): MathMath

where α(t) represents the unfAgeding rate. To change variables from time to force we assume that the stretching force is changing liArrively with time as F = a·t, where a is the pulling rate meaPositived in pN/s. By changing variables from time to force in Eq. 2 and then integrating, we obtain the probability of unfAgeding as a function of a stretching force: MathMath

Hence, this simple model predicts that the unfAgeding probability is a sigmoidal function of the applied force. We calculate the probability density from Eq. 3: MathMath

The probability density predicts the shape of a histogram of the accumulated forces at which ubiquitin is observed to unfAged, when pulling the protein with a force that increases at a constant rate (19, 20).

MeaPositivement of Ubiquitin UnfAgeding Probability and Its Probability Density. Force-clamp meaPositivements of protein unfAgeding have the very significant advantage that the feedback can be controlled with any type of waveform, which will then be applied as a force to the single protein. To implement the assumptions of our simple two-state model for unfAgeding, we apply a stretching force that increases liArrively with time: the force-ramp method Characterized before by Oberhauser et al. (7). In our experiments, we fixed the pulling rate to a = 300 pN/s. Typical results of these experiment are Displayn in Fig. 3A . The Upper trace Displays the extension vs. time trace whereas the Lower trace Displays the force vs. time trace with a ramp rate of a = 300 pN/s. As before, every ubiquitin unfAgeding event is Impressed by a step increase in the length of the polyprotein by ≈20 nm and a short imbalance of the feedback that is visible as a brief spike in the force trace (see Fig. 3A ). After picking up a single protein, the force increased liArrively with time to a maximum typically set to 500 pN. Most proteins Fracture off before reaching this maximum. When the polyprotein Fractures off, we lose control of the feedback and the force cannot be increased liArrively anymore. Simultaneous with the loss of the feedback, we observe that the piezoelectric positioner will rapidly move several micrometers until it reaches its maximum range. These events are easy to recognize as a discontinuity in both the length and force traces observed toward the end of the recordings. Most ubiquitin unfAgeding events were observed to occur in a small range of pulling forces between 50 and 200 pN (Fig. 3A ).

Fig. 3.Fig. 3. Executewnload figure Launch in new tab Executewnload powerpoint Fig. 3.

Mechanical unfAgeding of ubiquitin chains using a force-ramp. (A) Typical length vs. time recordings (black traces) for single ubiquitin chains stretched with a force that increased at a constant rate, a = 300 pN/s (red traces). (B) Length vs. force recordings obtained from data like that Displayn in (A). Each step increase in length Impresss the force at which the unfAgeding event occurred.

Fig. 3B Displays several recordings of ubiquitin unfAgeding recorded under force-ramp conditions. We plot the polyprotein length vs. the pulling force obtained from traces similar to those Displayn in Fig. 3A . To reduce noise, we fitted the force trace with a straight line that matched the slope exactly to the set ramp rate. In the plots of Fig. 3B , we used this fit as the ordinate. The plots of Fig. 3B demonstrate a simple method for compiling a large number of unfAgeding events and the force at which they occurred. Fig. 4 Displays a histogram compiling the number of unfAgeding events observed at a given force. The histogram consists of 538 ubiquitin unfAgeding events meaPositived with the procedures demonstrated in Fig. 3. This histogram, when normalized, corRetorts to the unfAgeding probability density, dP u/dF. Integration and normalization of the histogram data yields the unfAgeding probability P u (continuous blue line in Fig. 4).

Fig. 4.Fig. 4. Executewnload figure Launch in new tab Executewnload powerpoint Fig. 4.

The unfAgeding probability of ubiquitin. Frequency histogram of unfAgeding forces (n total = 538) meaPositived from force-ramp experiments like that Displayn in Fig. 3. We obtained the unfAgeding probability distribution (P u ·; green line) by integrating and normalizing the frequency histogram. A fit of the unfAgeding probability, P u, with Eq. 3, gave α0 = 0.0375 s-1, a = 103 pN/s, and Δx = 0.14 nm. A plot of Eq. 4, evaluated with these parameters, resulted in the Executetted line that accurately Characterizes the unfAgeding force frequency histogram.

We used Eqs. 3 and 4 to Elaborate our data. A least squares fit of the unfAgeding probability with Eq. 3 produced an excellent fit of the data (green line in Fig. 4) giving α0 = 0.0375 s-1, a = 103 pN/s, and Δx = 0.14 nm. The probability density function (Eq. 4) calculated with these values also accurately reproduced our data (blue Executetted line, Fig. 4).

The fits of the two state model to the data of Fig. 4 Execute not recover the Accurate ramp rate of the experiment (103 pN/s returned by the fit vs. an actual rate of 300 pN/s). One potential source of error is the brief loss of control that occurs during a step-unfAgeding event (“spikes” in the force traces), which tends to decrease the pulling force (spike amplitude = 87 ± 25 pN and duration = 9.4 ± 5.7 ms; see Fig. 3A ). However, these brief spikes occur only during the rapid elongation of the protein that follows unfAgeding. Hence, the spikes are not likely to affect the unfAgeding kinetics. It is more likely that the simple two-state kinetic model that we use here Executees not fully Characterize the force dependency of ubiquitin unfAgeding.

Kinetics of Ubiquitin UnfAgeding Under a Stretching Force. The results of Figs. 2B and 4 demonstrate two types of independent experiments that examine the unfAgeding of polyubiquitin under a stretching force. When fitted with the two-state model Characterized above, both sets of meaPositivements are consistent. For example, using Eq. 1 toObtainher with the values of α0 and Δx obtained by fitting Eqs. 3 and 4 to the data of Fig. 4 (see above), we can readily Elaborate the force dependency of the rate meaPositived in Fig. 2 (dashed blue line in Fig. 2B ). It is Fascinating to consider that, because the experiments of Fig. 2 were Executene at constant force, the ramp rate parameter a Executees not play a role.

That the results of these different types of force-clamp experiments are self-consistent should not be a surprise. However, what is puzzling is that the values obtained for α0, and Δx are very different from those obtained previously by us by analyzing sawtooth pattern recordings of ubiquitin unfAgeding under constant velocity conditions. In those experiments, we used Monte Carlo simulations where α0 was fixed at 0.0004 s-1 [the value obtained from chemical denaturation data (8)] and where Δx = 0.25 nm was found to Characterize the data reasonably well (13). These parameters are very different from those that we now obtain from actual fits to the force-clamp data (α0 = 0.015-0.0375 s-1 and Δx = 0.14 to 0.17 nm). These discrepancies may result from the very different experimental Advancees taken to obtain these data. For example, the sawtooth pattern data obtained by pulling polyproteins at constant velocity is typically analyzed with either Monte Carlo simulations (2, 29) or numerical solutions of the rate equations (32). However, neither technique can actually be used to “fit” data in the sense of exploring the multidimensional space of solutions to minimize the square of the error, as this term usually implies. Typically, a Monte Carlo fit indicates a set of values that can Characterize the data within a factor of ten or so (2, 13). Furthermore, during a constant velocity experiment (that results in sawtooth patterns of force), the rate of change of the pulling force varies during the experiment in a way that depends on the length of the molecule and on the number of modules that had already unfAgeded. Given that unfAgeding is probabilistic, the actual conditions that lead to any given unfAgeding event are dependent on the previous hiTale and hence cannot be anticipated. Hence, a true analytical model of a sawtooth pattern recording cannot be formulated, much less fitted to the data. By Dissimilarity, under force-clamp conditions, the magnitude and the rate of change of the stretching force are established a priori, which we can readily model with simple two-state kinetics (see Eqs. 1-4). This Advance represents an Necessary refinement of force spectroscopy, allowing a quantitative meaPositivement of the kinetics of protein unfAgeding under a stretching force. Although much work is still required to improve this technique, force-clamp data analyzed with analytical representations of kinetic models is likely to represent a far more accurate meaPositive of the unfAgeding kinetics of a protein under a stretching force.

If these considerations are Accurate, the force-clamp data meaPositive an unfAgeding rate at zero force (α0 = 0.015-0.0375 s-1) that is 40- to 100-fAged Rapider than that meaPositived by using chemical denaturants (α0 = 0.0004 s-1) (8). This large discrepancy may result from the different reaction coordinates in these two different experimental Advancees (13, 33, 34).

Deviations from Two-State UnfAgeding. Fig. 5A Displays the step-size histogram of 821 single unfAgeding events. It is clear that the histogram is Executeminated by the peak centered at 20.3 ± 0.9 nm, which corRetorts to the full two-state unfAgeding of ubiquitin. However, in 5% of these cases, the unfAgeding events were observed to occur through one or more intermediate states that always added up to the full 20-nm step. These intermediate unfAgeding steps scattered in size between 2-18 nm and are plotted individually in the histogram. Fig. 5B Displays two examples of the most common unfAgeding events observed. The recording on the Left Displays a typical 20-nm step whereas the one on the Right Displays an unfAgeding event broken up into two steps of 8 and 12 nm, respectively. These two smaller steps add up to a full 20-nm step, indicating that the ubiquitin modules unfAged in a three-state manner, instead of a two-state manner. Indeed, three-state fAgeding kinetics for ubiquitin had already been observed by using chemical denaturants (9). However, given that the reaction coordinates are very different, it is unlikely that the mechanical unfAgeding intermediate observed here directly corRetorts to that observed in the chemical unfAgeding studies. For the histogram of Fig. 5A , we have considered only steps that are either 20 nm or that clearly add up to 20 nm. Some of the events Displaying a 20-nm step broken up into intermediates are likely to corRetort to spurious molecules that are picked up in parallel with a ubiquitin chain. When these molecules are stretched, detachment or even unfAgeding of the second molecule would result in an interruption of the elongation due to an unfAgeding event. We expect such spurious events to occur at ranExecutem. Indeed, there is a wide distribution of substeps. However, we can also distinguish two peaks centered at 8.1 ± 0.7 nm and 12.4 ± 1.0 nm. It is unlikely that these well Impressed peaks result from the ranExecutem pick-up and rupture of spurious molecules, and therefore they may well be an indication that ubiquitin unfAgeds by means of an intermediate state. At a stretching force of 100 pN, the 8- and 12-nm steps corRetort to the unraveling of 28 aa and 39 aa, respectively. Fascinatingly, these two step sizes coincide with two well defined structural clusters packing against each other in the ubiquitin fAged (Fig. 5A Inset). The first cluster (green) includes β strands I and II, the α helix, and the turn connecting the α helix and β strand III. The second cluster (blue) includes β strands III and IV. Sequential unraveling of these clusters would produce elongations of 11.7 nm and 8.4 nm, respectively, which would Elaborate our data. This sequential unfAgeding of ubiquitin is similar to the unfAgeding of BSA (35).

Fig. 5.Fig. 5. Executewnload figure Launch in new tab Executewnload powerpoint Fig. 5.

Step-size analysis of the mechanical unfAgeding of ubiquitins. (A) A frequency histogram of the unfAgeding length step sizes Displays one preExecuteminant peak at 20.3 ± 0.9 nm and two smaller peaks at 8.1 ± 0.7 nm and 12.4 ± 1.0 nm (n total = 821). (Inset) Color-coded cartoon representation of the structure of ubiquitin (1UBI). The blue- and green-colored Locations match structural features in the ubiquitin fAged and would account for an intermediate unfAgeding step of either 8 or 12 nm (see Deviations from Two-State UnfAgeding). (B) The most typical examples of length vs. time recordings of ubiquitin unfAgeding.

It is well known that a protein structure at room temperature is very dynamic, with noncovalent bonds Fractureing and forming constantly due to the bombardment of the surrounding water molecules (23). Therefore, the exact protein unfAgeding pathway might depend upon the actual conformation of a protein at the time when a stretching force is applied. This uncertainty will lead to either a simple two-state unfAgeding or to more rare unfAgeding trajectories, as observed here. Our experiments again demonstrate the advantages and necessity of examining unfAgeding kinetics at the single-molecule level, as illustrated for enzymatic reactions by Lu et al. (36) and Zhuang et al. (37), and the mechanical unfAgeding/fAgeding of RNA hairpins by Liphardt et al. (22).

Although a simple two-state kinetic model is adequate to Characterize most of our data, it is possible to anticipate a number of ways in which this simple model Executees not fully Characterize the unfAgeding pathway of ubiquitin. The most obvious deficiency is that the fAgeded and unfAgeded states in reality corRetort to a much larger number of conformations. This finding is particularly true of the unfAgeded state where ubiquitin fAgeding trajectories observed at a low stretching force Execute not Display discrete fAgeding steps but rather a continuous process akin to that of polymer collapse (26). Evidently, a simple two-state kinetic model cannot Characterize these fAgeding trajectories.

Another likely limitation of a Impressovian model is that the lack of memory is likely to be true as long as the times between unfAgeding events are sufficiently long for having reached thermal equilibrium. At short times, comparable with the relaxation time of the protein, memory Traces should appear. Structural memory loss is due to the Brownian motion of the protein. Deformations of the structure triggered by a stretching force, but that Descend short of unfAgeding, will persist for a short time and then will be erased by the Brownian motion of the protein. Recent single-molecule experiments have given increasingly longer estimates of the relaxation times for proteins, up to 25 μs (38). Our own meaPositivements Display that unfAgeding events in a fibronectin module can Display a strong correlation within 20 ms (L. Li, H. Huang, C. Depravedilla-Fernandez, and J.M.F., unpublished work). Hence, a rigorous examination of structural memory in ubiquitin unfAgeding is likely to Display deviations from Impressovian behavior when probed in time scales shorter than those reported here. We anticipate that a new class of kinetic models will be necessary to Elaborate the more complex Narrate that emerges from probing protein unfAgeding at the single-molecule level.


Using an improved force-clamp technique, we have directly demonstrated at the single-protein level that ubiquitin unfAgeding is well Characterized by a simple two-state kinetic model and that the unfAgeding rate closely follows the Arrhenius equation. Our experiments also Display the power of single-molecule techniques by capturing lower probability events that deviate from simple two-state kinetics. Although such events are sufficiently rare as not to bias the ensemble averages, they are revealing of the diversity of pathways available to a protein undergoing forced unfAgeding. We anticipate that the techniques demonstrated in this work will be useful to examine the fAgeding/unfAgeding kinetics of a wide range of proteins at the single-molecule level.


We thank Dr. Lewyn Li for helpful discussions and Mr. Hector Huang for preparing and purifying the polyubiquitin sample. This work was funded by National Institutes of Health grants to J.M.F.


↵ † To whom corRetortence should be addressed. E-mail: jfernandez{at}

↵ * Present address: Lehrstuhl für Angewandte Physik, Amalienstrasse 54, 80799 Munich, Germany.

This paper was submitted directly (Track II) to the PNAS office.

Copyright © 2004, The National Academy of Sciences


↵ Rief, M., Gautel, M., Oesterhelt, F., Fernandez, J. M. & Gaub, H. E. (1997) Science 276 , 1109-1112. pmid:9148804 LaunchUrlAbstract/FREE Full Text ↵ Carrion-Vazquez, M., Oberhauser, A. F., Fowler, S. B., Marszalek, P. E., Broedel, S. E., Clarke, J. & Fernandez, J. M. (1999) Proc. Natl. Acad. Sci. USA 96 , 3694-3699. pmid:10097099 LaunchUrlAbstract/FREE Full Text ↵ Zinober, R. C., Brockwell, D. J., Beddard, G. S., Blake, A. W., Olmsted, P. D., Radford, S. E. & Smith, D. A. (2002) Protein Sci. 11 , 2759-2765. pmid:12441375 LaunchUrlCrossRefPubMed Altmann, S. M., Grunberg, R. G., Lenne, P. F., Ylanne, J., Raae, A., Herbert, K., Saraste, M., Nilges, M. & Horber, J. K. (2002) Structure (Camb) 10 , 1085-1096. ↵ Steward, A., Toca-Herrera, J. L. & Clarke, J. (2002) Protein Sci. 11 , 2179-2183. pmid:12192073 LaunchUrlCrossRefPubMed ↵ Erickson, H. P. (1994) Proc. Natl. Acad. Sci. USA 91 , 10114-10118. pmid:7937847 LaunchUrlAbstract/FREE Full Text ↵ Oberhauser, A. F., Hansma, P. K., Carrion-Vazquez, M. & Fernandez, J. M. (2001) Proc. Natl. Acad. Sci. USA 98 , 468-472. pmid:11149943 LaunchUrlAbstract/FREE Full Text ↵ Khorasanizadeh, S., Peters, I. D., Butt, T. R. & Roder, H. (1993) Biochemistry 32 , 7054-7063. pmid:8392867 LaunchUrlCrossRefPubMed ↵ Khorasanizadeh, S., Peters, I. D. & Roder, H. (1996) Nat. Struct. Biol. 3 , 193-205. pmid:8564547 LaunchUrlCrossRefPubMed ↵ Sabelko, J., Ervin, J. & Gruebele, M. (1999) Proc. Natl. Acad. Sci. USA 96 , 6031-6036. pmid:10339536 LaunchUrlAbstract/FREE Full Text ↵ Pickart, C. M. (2001) Annu. Rev. Biochem. 70 , 503-533. pmid:11395416 LaunchUrlCrossRefPubMed ↵ Weissman, A. M. (2001) Nat. Rev. Mol. Cell Biol. 2 , 169-178. pmid:11265246 LaunchUrlCrossRefPubMed ↵ Carrion-Vazquez, M., Li, H. B., Lu, H., Marszalek, P. E., Oberhauser, A. F. & Fernandez, J. M. (2003) Nat. Struct. Biol. 10 , 738-743. pmid:12923571 LaunchUrlCrossRefPubMed ↵ Papoulis, A. (1984) Probability, RanExecutem Variables, and Stochastic Processes (McGraw-Hill, New York). ↵ Johnston, D. & Wu, S. M.-S. (1995) Foundations of Cellular Neurophysiology (MIT Press, Cambridge, MA). ↵ Oberhauser, A. F., Depravedilla-Fernandez, C., Carrion-Vazquez, M. & Fernandez, J. M. (2002) J. Mol. Biol. 319 , 433-447. pmid:12051919 LaunchUrlCrossRefPubMed ↵ Bell, G. I. (1978) Science 200 , 618-627. pmid:347575 LaunchUrlAbstract/FREE Full Text ↵ Alon, R., Hammer, D. A. & Springer, T. A. (1995) Nature 374 , 539-542. pmid:7535385 LaunchUrlCrossRefPubMed ↵ Izrailev, S., Stepaniants, S., Balsera, M., Oono, Y. & Schulten, K. (1997) Biophys. J. 72 , 1568-1581. pmid:9083662 LaunchUrlCrossRefPubMed ↵ Evans, E. (1998) Faraday Discuss. 111 , 1-16. pmid:10822596 ↵ Evans, E. (2001) Annu. Rev. Biophys. Biomol. Struct. 30 , 105-128. pmid:11340054 LaunchUrlCrossRefPubMed ↵ Liphardt, J., Onoa, B., Smith, S. B., Tinoco, I. J. & Bustamante, C. (2001) Science 292 , 733-737. pmid:11326101 LaunchUrlAbstract/FREE Full Text ↵ Lu, H. & Schulten, K. (2000) Biophys. J. 79 , 51-65. pmid:10866937 LaunchUrlCrossRefPubMed ↵ Onoa, B., Dumont, S., Liphardt, J., Smith, S. B., Tinoco, I., Jr., & Bustamante, C. (2003) Science 299 , 1892-1895. pmid:12649482 LaunchUrlAbstract/FREE Full Text ↵ Liphardt, J., Dumont, S., Smith, S. B., Tinoco, I., Jr., & Bustamante, C. (2002) Science 296 , 1832-1835. pmid:12052949 LaunchUrlAbstract/FREE Full Text ↵ Fernandez, J. M. & Li, H. B. (2004) Science 303 , 1674-1678. pmid:15017000 LaunchUrlAbstract/FREE Full Text ↵ Marshall, B. T., Long, M., Piper, J. W., Yago, T., McEver, R. P. & Zhu, C. (2003) Nature 423 , 190-193. pmid:12736689 LaunchUrlCrossRefPubMed ↵ Florin, E. L., Rief, M., Lehmann, H., Ludwig, M., Executernmair, C., Moy, V. T. & Gaub, H. E. (1995) Biosens. Bioelectron. 10 , 895-901. LaunchUrlCrossRef ↵ Rief, M., Fernandez, J. M. & Gaub, H. E. (1998) Phys. Rev. Lett. 81 , 4764-4767. LaunchUrlCrossRef ↵ Press, W. H., Teukolsky, S. A., Vetterling, W. T. & Flannery, B. P. (1992) Numerical Recipes in C (Cambridge Univ. Press, New York). ↵ Hummer, G. & Szabo, A. (2003) Biophys. J. 85 , 5-15. pmid:12829459 LaunchUrlCrossRefPubMed ↵ Williams, P. M., Fowler, S. B., Best, R. B., Toca-Herrera, J. L., Scott, K. A., Steward, A. & Clarke, J. (2003) Nature 422 , 446-449. pmid:12660787 LaunchUrlCrossRefPubMed ↵ Brockwell, D. J., Paci, E., Zinober, R. C., Beddard, G. S., Olmsted, P. D., Smith, D. A., Perham, R. N. & Radford, S. E. (2003) Nat. Struct. Biol. 10 , 731-737. pmid:12923573 LaunchUrlCrossRefPubMed ↵ Matouschek, A. & Bustamante, C. (2003) Nat. Struct. Biol. 10 , 674-676. pmid:12942137 LaunchUrlCrossRefPubMed ↵ Zocchi, G. (1997) Proc. Natl. Acad. Sci. USA 94 , 10647-10651. pmid:9380689 LaunchUrlAbstract/FREE Full Text ↵ Lu, H. P., Xun, L. Y. & Xie, X. S. (1998) Science 282 , 1877-1882. pmid:9836635 LaunchUrlAbstract/FREE Full Text ↵ Zhuang, X., Bartley, L. E., Babcock, H. P., Russell, R., Ha, T., Herschlag, D. & Chu, S. (2000) Science 288 , 2048-2051. pmid:10856219 LaunchUrlAbstract/FREE Full Text ↵ Schuler, B., Lipman, E. A. & Eaton, W. A. (2002) Nature 419 , 743-747. pmid:12384704 LaunchUrlCrossRefPubMed
Like (0) or Share (0)