ABSTRACT
The introduction of disulfide bonds into proteins creates additional mechanical barriers and limits the unfolded contour length (i.e., the maximal extension) measured by single-molecule force spectroscopy. Here, we engineer single disulfide bonds into four different locations of the human cardiac titin module (127) to control the contour length while keeping the distance to the transition state unchanged. This enables the study of several biologically important parameters. First, we are able to precisely determine the end-to-end length of the transition state before unfolding (53 [Angstrom]), which is longer than the end-to-end length of the protein obtained from NMR spectroscopy (43 A). Second, the measured contour length per amino acid from five different methods (4.0 � 0.2 [Angstrom]) is longer than the end-to-end length obtained from the crystal structure (3.6 [Angstrom]). Our measurement of the contour length takes into account all the internal degrees of freedom of the polypeptide chain, whereas crystallography measures the end-to-end length within the "frozen" protein structure. Furthermore, the control of contour length and therefore the number of amino acids unraveled before reaching the disulfide bond (n) facilitates the test of the chain length dependence on the folding time (τ^sub F^). We find that both a power law scaling τ^sub F^ ∝ n^sup λ^ with λ = 4.4, and an exponential scaling with n^sup 0.6^ fit the data range, in support of different protein-folding scenarios.
Abbreviations used: 127, 27th immunoglobulin-like domain of human cardiac titin;WLC, worm-like chain model; DTT, 1,4-dithiothreitol; SMD, steered molecular dynamics; AFM, atomic force microscope; PBS, phosphate-buffered saline.
� 2007 by the Biophysical Society
0006-3495/07/01/225/09 $2.00
INTRODUCTION
Many important predictions of the physical models for protein folding andl unfolding show a clear dependence on the number of amino acids, n, in a protein chain. These include the WLC theory of polymer elasticity (1) to measure the contour length of the molecule (2-7) and the scaling of protein-folding times with n (8-11). There is a wide discrepancy in the contour length contribution per amino acid for different proteins in the literature and, surprisingly, no apparent trend for the folding times as a function of the length of the chain. The diversity of the proteins used in previous studies and the variety of their respective transition state positions have precluded the investigation of the effect of n alone on contour length and folding rates.
With the aim of retaining the same transition state, yet unfolding varying numbers of amino acids, we engineer single disulfide bonds (S-S) into various locations of the human cardiac titin 127 module (12). We then study the folding and unfolding properties of the different constructs using the force-extension technique with an AFM (13). When a protein is placed under force, it first extends up to its mechanically stable transition state, where the protein's tertiary structure is still intact and able to store elastic energy. Stretching beyond this point disrupts the mechanical architecture and leads to the unraveling of the protein (14). These experiments are usually interpreted using the WLC theory of entropic elasticity (1), which predicts the total contour length. ΔL^sub u^, of the protein, proportional to the number of unraveled amino acids in the chain. The contour length measured by force-extension experiments is therefore defined as the maximum distance between the mechanically stable transition state and the fully extended state of the linear protein chain with no degrees of freedom. All previous studies using force-extension experiments with an AFM have calculated the contour length as the distance between the native slate and the extended state along the reaction coordinate, neglecting the important contribution of the distance up to the transition state (2-6,15,16). For this reason, the cited values for the contour length per amino acid, l, vary from protein to protein, spanning from 3.4 [Angstrom] (2) to 3.6 [Angstrom] (3.15.16). 3.8 [Angstrom] (4-6), and further up to 4.0 [Angstrom] (7).
The position of a disulfide bond was previously shown to uniquely determine the contour length associated with protein unfolding, ΔL^sub u^, in a number of different proteins (17-20). In our study, varying ΔL^sub u^ within the same protein allows for the most accurate determination of the contour length per amino acid to date as well as the characterization of the end-to-end length of the transition state from the extrapolated value at zero contour length. We employ four independent measurements, including force-extension and force-clamp techniques as well as SMD simulations, to obtain the contour length of an amino acid.
The same disulfide-honded constructs that allow for the controlled variation in the unfolded number of amino acids also provide a unique platform for studying the effects of chain length on protein folding. The architecture of the protein free energy landscape can be inferred from the scaling of the protein folding time, τ^sub F^, with n. Several such theoretical predictions have been pustulated in the literature, satisfying different folding scenarios (21-24). Furthermore, a variety of numerical simulations have predicted a strong dependence of folding limes on n in model proteins ( 10,25,26) as well as its dependence on the interaction energies between the native contacts, the presence of side chains, and the topology of the native-state fold. However, previous ensemble refolding experiments have argued that the folding times are solely determined by the native stale topology, as quantilied by contact order (CO), exhibiting little or no correlation with the chain length (8-11,27). These experiments were performed on a variety of proteins, where the distinct amino acid sequences and the topologies of the native structures may have played a dominant role in the time scales of protein folding.
Because we precisely control the contour length with engineered disultide bonds in our experiments, we can also measure the refolding time dependence on chain length within the same protein molecule and thus directly test the existing theoretical models for folding. The data indicate that an energy landscape with a power law scaling with n, as well as glassy models with a heterogeneous distribution of energy barriers successfully describe the experimentally observed trend.
MATERIALS AND METHODS
Protein engineering
We used I27^sub Cys-free^ protein, a cysteine-tree version (Cys^sup 47^ to Ala, Cys^sup 63^ to Ala) of the 27th module of the I band of human cardiac titin (I27) as a platform to engineer the 127 mutants used in this study. Disulfide bonds at specific positions were engineered by introducing additional mutations according to methods described elsewhere (13): Gly^sup 32^ to Cys, Ala^sup 75^ to Cys (I27^sub G32C-A75C^); Pro^sup 28^ to Cys, Lys^sup 54^ to Cys (127^sub P28C-K54C^); Glu^sup 24^ to Cys, Lys^sup 55^ to Cys (I27^sub E24C-K55C^): Asp^sup 46^ to Cys. His^sup 61^ to Cys (I27^sub D46C-H61C^). The residue pair for cysteine mutations was chosen such that their C^sub α^ atoms are separated by 3.5-5.5 [Angstrom] in the NMR structure of win 127 (Protein Data Bank file: 1TIT). This simple selection criterion proved in be very useful in comparison to the computationally intensive approach based on the disulfide bond prediction algorithm (28). In our experiments we used polyproteins made of eight identical repeats of a single type of mutant i27 protein. The (I27^sub P28C-K54C^)^sub 8^ polyprotein was engineered with two C-terminal cysteines for attachment as described previously (13). However, because these cysteines proved unnecessary, they were not used in the other constructs. The synthetic polyproteins were cloned into pT7Blue vector-XL1 Blue E. coli system and expressed in pQE80L vector-BLR(DE3) E. coli cells. The proteins were puritied using Ni^sup 2+^ affinity chrommography followed by size-exclusion chromatography.
Single-molecule force spectroscopy
The details of the AFM and its mode of operation have been described elsewhere (13,29). The spring constant of the cantilevers used in our experiment was measured using the equipartition theorem (30). We used silicon nitride cantilevers ( Veeco, Santa Barbara, CA) with an average spring conslant of ~45 pN/nm in force-extension experiments and those with ~15 pN/ nm in force-damp experiments. An aliquot of 5-10 �l of the polyprotein (0.1 mg/ml) solution was added to the PBS solulion, pH 7.0, on a gold-coated cover slide and immobilized for 10 min. The presence of disulfide bond in the disulfide -engineered 127 modules was confirmed by recoding 2-10 sawtooth-patterned forte-extension traces containing unsequestered unfolding events only. The disulfide reducing agent. DTT (5-100 mM), was then added to the solution from a stock solution of 1 M in PBS and thoroughly mixed. The pulling of the polyproteins was resumed to record the unsequestered unfolding and disulfide bond reduction events under a stretching force. However, in protein-refolding experiments and force-clamp experiments, the poly protein is directly added to the PBS solution containing the specitied amount of DTT. In all our force-extension experiments the pulling rate was 400 nm/s unless otherwise mentioned. All our experiments were performed at room temperature.
RESULTS AND DISCUSSION
The fingerprint of an engineered disulfide bond
Fig. 1 demonstrates the engineering of a single disulfide bond into the I27^sub Cys-free^ protein. We engineered two Cys residues into positions 32 and 75 (G32C-A75C). which arc distant in the sequence (Fig. 1 A) but are vicinal in thc structure (Fig. 1 B). They are situated in the β-strands, which are Hexible enough to allow tor their spontaneous oxidation, forming a disulfide bond. In Fig. 1 B, we show the 46 amino acids that remain "unsequestered" by the disulfide bond (red) as well as the 43 amino acids that are "trapped" behind the mechanically strong disulfide bond (green). The reduction of the disulfide bond, which is a prerequisite tor the unfolding of the trapped amino acids, is never observed in the absence of a reducing agent, such as DTT, independent of the position of the disulfide bond. This result is in agreement with the observation that covalent bonds are mechanically too strong for rupture below 1 nN of force (31). The kinetics of reduction of the disulfide bond under a stretching force has been established previously (32). The experiments are therefore performed in the presence of DTT, such that both the unfolding and reduction under a stretching force can be observed within the same force-extension trace. In all our experiments, we are concerned with measuring the contour length increase after the disulfide bond reduction but not the kinetics of reduction itself. Thus, we vary the concentration of DTT (5-100 mM) to maximize the observed disulfide reduction events in each experiment while minimizing the reduction of the disulfide bonds before mechanical unfolding. Force-extension curves obtained from stretching the (I27^sub G32C-A75C^)^sub 8^ polyprotein in the presence of DTT (100 mM) reveal two consecutive sawtooth patterns, as shown in Fig. 1 C. The first series of equally spaced peaks corresponds to the unsequestered unfolding depicted in red, and the second sawtooth pattern at a much higher force corresponds to disullide bond reduction events and the release of the trapped amino acids, completing the full unraveling of the 127 protein. We can therefore distinguish between unsequestered unfolding events and disulfide bond reduction events from the difference in the measured contour lengths, ΔL^sub u^ and ΔL^sub r^, respectively. The unfolding force peaks are titled with the WLC theory (1), yielding the persistence length and the contour length of each module (14). The unsequestered portion of the molecule unfolds at an unfolding force of F^sub u^ = 182 � 33 pN with a contour length of ΔL^sub u^ = 12.7 � 0.5 nm, whereas the disulfide bond reduction takes place at forces of F^sub r^ = 303 � 135 pN with a contour length of ΔL^sub r^ = 16.5 � 0.7 nm. The sum of the two contour lengths observed in the double sawtooth patterns, ΔL^sub u^ + ΔL^sub r^ = 29.2 nm, is comparable to the contour length observed from a fully unraveled 127 module (ΔL^sub u^ = 28.3 � 1.1 nm).
Contour lengths, ΔL^sub u^ and ΔL^sub r^, depend on the position of the disulfide bond
Fig. 2 A shows force-extension traces for all four of our polyprotein constructs: (I27^sub G32C-A75C^)^sub 8^, (I27^sub E24C-K55C^)^sub 8^, (I27^sub P28C-K54C^)^sub 8^, and (I27^sub D46C-H61C^)^sub 8^. All of these experiments were performed in the presence of DTT. The contour lengths for the unfolding, ΔL^sub u^, and the disulfide bond reduction, ΔL^sub r^, are derived from the WLC theory, shown by the red and green curves. The characteristics of each construct in terms of the number of unsequestered amino acids, the number of trapped amino acids, their respective contour lengths, and force distributions are presented in Table 1. Fig. 2, B and C, shows histograms of the contour length increments measured from all tour constructs (ΔL^sub u^ in Fig. 2 B and ΔL^sub r^ in Fig. 2 C). In all four cases, the sums of ΔL^sub u^ + ΔL^sub r^ (~28.8 nm) are within the error of the contour length measured for the fully reduced proteins (~28.4 nm), where the disulfide bond is reduced before mechanical stretching (Table 1). In addition, all the fully reduced proteins (shown in blue in Fig. 2 A) have the same contour length (ΔL^sub f^) as the wild-type I27 (13) to within the error in the experiment. Because the contour length is measured from the transition state up to the fully extended molecule, the equivalence of these results suggests that the transition state remains the same between the different constructs. Next, we further confirm this assumption for one of the protein constructs.
The distance to the transition state is unaltered by the disulfide bond
We vary the pulling speed from 80 nm/s to 4000 nm/s and measure the corresponding force distribution of the unsequestered unfolding (F^sub u^) for the (I27^sub E24C-K55C^)^sub 8^ polyprotein. The average unfolding force versus the pulling speed is shown in Fig. 3. Monte Carlo simulations using a two-state model for unfolding were performed to reproduce the experimental data, as described by Carrion-Vazquez et al. (13). The fit of the data in Fig. 3 has an unfolding distance of Δx^sub u^ = 0.25 nm and a spontaneous unfolding rate constant (i.e., at zero force) of k^sup 0^^sub u^ = 18.5�10^sup -4^ s^sup -1^ as its only adjustable parameters. The unfolding distance (Δx^sub u^) is identical to that of the wild-type protein, 0.25 nm (13), which indicates that the engineered disulfide bond has not altered the distance to the transition state. However, the spontaneous unfolding rate constant is increased 5.6 times, suggesting a reduced stability of the mutant protein.
The contour length contribution of a single amino acid
Each polyprotein construct provides two contour length increments (ΔL^sub u^ and ΔL^sub r^), and the number of amino acids contributing to each of them is determined by the position of the disulfide bond (Table 1). Fig. 4 A shows a plol of ΔL^sub u^ (triangles) and ΔL^sub r^ (circles) as a function of the number of contributing amino acids. Linear fits to the data. ΔL^sub u^ and ΔL^sub r^, give the contour length contribution per amino acid, l^sub u^ and l^sub r^, respectively (solid lines in Fig. 4 A). The measured slopes give l^sub u^ = 3.9 �/aa and l^sub r^ = 4.3 �/aa for the different regions of the protein. Because of the very broad mechanical stability of the disulfide bond (spanning from 100 pN to 1 nN), the typical sawtooth pattern is often jagged and leads to very narrow force ranges over which the WLC is fit to the data, as seen in Fig. 2 A. Fitting these force curves leads to a much higher uncertainty in the obtained values for the contour length, such that the value of 4.3 �/aa is prone to error. Although it seems to signify a dependence of the measured contour length on the force range, force-clamp experiments and SMD simulations shown below confirm the absence of any force dependence.
Both of these values are significantly longer than the end-to-end length between neighboring amino acids obtained from x-ray crystallography data on extended antiparallel β-sheets (3.6 �/aa) (33). This is because the contour length takes into account the configurational entropy of the poly-peptide chain arising from the internal degrees of freedom of the constituent bonds, discussed in detail below. They also deviate from the l values assumed by many groups employing AFM force spectroscopy (2-6,15,16). These experiments have relied on the folded length of the protein (obtained from NMR or x-ray crystallography) to determine l in the absence of an independent measurement of the end-to-end length of the transition state. However, our results are in good agreement with the only other experiment in which the transition state end-to-end length was accounted for by the measurement (7). In that experiment, five glycine amino acids were introduced into the I27 wild-type protein, such that the transition state remained unperturbed by the insertion. The difference in contour length between the wild-type and the modified I27 corresponded to l = 4.0 �/aa (7).
In Fig. 4 A, both of me fitted straight lines to ΔL^sub u^ and ΔL^sub r^ have nonzero intercepts with the abscissa: 13.4 amino acids for the ΔL^sub u^ data and 4.3 amino acids for the ΔL^sub r^ data, respectively. The zero contour length measurement for ΔL^sub u^ corresponds to the end-to-end length of the protein (13.4 aa � l^sub u^ = 53 � (L^sub NC^)) at the transition state before unfolding, in agreement with previous studies (7). Earlier studies on wild-type I27 revealed the structure of the transition state (34,35) and its end-to-end length, which corresponds to 17 amino acids when u contour length per amino acid of ~3.8 �/aa is assumed (7). The measurement presented here is more accurate because it does not assume the contour length per amino acid. It should be noted that the transition state is longer than the end-to-end length of the native state of I27, measured to be 43 � by NMR spectroscopy (12).
The assumption is that the intercept for the ΔL^sub r^ data in Fig. 4 A corresponds to the distance between the α-carbon atoms of the cysteines, L^sub SS^, through the disulfide bond (C^sub α^-C^sub β^-S-S-C^sub β^-C^sub α^). However, the measured value of 4.3 aa � l^sub r^ = 18 � is in excess of the range expected from the average crystallographic measurement of disulfide bonds (5.5 �) in proteins (36), suggesting that further elements are present in the protein structure.
Independent confirmation of the contour length per amino acid
The preceding experiments measure the widely distributed force-versus-extension behavior as the protein is stretched under a constant velocity. To ensure a constant unfolding force to measure the corresponding end-to-end length as a function of time, we now stretch the protein using force-damp spectroscopy (29). Such experiments on the (I27^sub G32C-A75C^)^sub 8^ poly-protein have been done previously to investigate the kinetics of the disulfide bond reduction under a stretching force (32). The protein is first placed under a force of 130 pN for 1 s, during which the protein unfolds much faster than the disultide bonds are reduced. Alter 1 s, the force is changed to a value between 50 pN and 800 pN for 5 s, during which the disulfide bond reduction takes place. The 43 amino acids trapped behind the disulfide bond are then unraveled under the stretching force in a two-state, stepwise manner for each protein module. This results in staircases, with each step corresponding to the reduction of a single disulfide bond (inset of Fig. 4 B). The constant stretching force is plotted against the average step length in Fig. 4 B. As the force is increased, the step size for the unraveling of me 43 amino acids increases according to the WLC theory of entropic elasticity. A Levenberg-Marquardt fit of the WLC results in a persistence length of 3 � and a contour length of L^sub r^ = 16.4 nm. This contour length corresponds to the maximum distance between the transition state for disulfide bond reduction and the fully extended stale of the 43 trapped amino acids. The transition-state length, which corresponds to the end-to-end length of the disulfide bond, is measured to be 0.55 nm by x-ray crystallography (36). Therefore, the total contour length obtained by extending the 43 residues is 17 nm (= 16.4 nm + 0.55 nm). The contour length per amino acid is then calculated to be l^sub r^ = 4.0 �/aa.
All experimentally measured contour lengths per amino acid are much higher than the end-to-end length values of 3.6 �/aa from x-ray crystallography (33). To examine the origin of this discrepancy, we next perform a series of SMD simulations to mimic our force-clamp experiments (34,37). We stretch the G strand of I27 protein (consisting of 12 amino acids) at four constant forces. In each simulation, the peptide segment is equilibrated under force for a period of 1 ns in explicit solvent with periodic boundary conditions. The equilibrium end-to-end distances of this peptide are then plotted in Fig. 4 B (cross symbols). They are in good agreement with the experimental data and are fitted well with the WLC theory, resulting in a contour length per amino acid of 3.9 �/aa. This is an important result in that it shows that no enthalpic stretching elasticity component is necessary to explain the data in the range of 50-800 pN. Furthermore, studies on vibrational spectroscopy of bonds measure an elastic constant of >5 N/cm for C-C or C-N bonds (38), which translates to a stretching of <0.05 �/aa at the highest force measured in our experiments (1 nN). Such a small elongation is not detectable from our data. As a result, the contour length is independent of the force.
Surprisingly, we find that within the polypeptide backbone, the dihedral angles and the bond angles allow the end-to-end length of the individual chemical bonds to fluctuate in the range of 0.1 �, as compared to the 0.02 � fluctuations usually observed in x-ray crystallography data. These fluctuations reveal the internal degrees of freedom of the chain that account for the difference between the measured contour length in our experiments and the end-to-end lengths of the individual bonds, as measured by the crystal structure. We thus demonstrate that AFM force spectroscopy can measure the degrees of freedom of chemical bonds that other structural characterization techniques overlook.
In summary, the disulfide bond variation first allows us to precisely measure the transition-state length under force, with no a priori assumptions. We then determine l by force-extension experiments (3.9 �/aa and 4.3 �/aa), force-clamp experiments (4.0 �/aa), as well as SMD simulations (3.9 �/aa), as shown in Fig. 4 C. Therefore, the average value for l from the rive well-controlled measurements (dashed line), including the previous study using inserted glycines (4.0 �/aa) (7), is found to be 4.02 � 0.16 �/aa. This value is a measure of the con figurational entropy of the amino acids in a polypeptide chain arising from the internal degrees of freedom of the constituent bonds. If all the bonds could rotate freely, the sum of the lengths of the skeletal bonds of a single amino acid backbone, which is 4.3 �/aa, would correspond to the accessible contour length in the experiment. However, the reduced degrees of freedom within the peptide backbone limit the measured contour length of an amino acid to l = 4.0 �/aa, a value between the end-to-end length within the protein crystal of 3.6 �/aa and the extended length of a purely linear molecule of 4.3 �/aa. This is the first study that systematically varies the number of amino acids contributing to the contour length, takes the transition state for unfolding into account, and utilizes different experimental approaches to measure the contour length per amino acid.
Dependence of folding rates on protein length
We have demonstrated how the introduction of a disulfide bond precisely determines the unfolding contour length of the protein chain as the maximum distance between the transition state position and the rigid disulfide bond. This capability presents a unique opportunity to study the folding time (τ^sub F^), defined as the time it lakes to reform the transition state and therefore the protein's mechanical stability, solely as a function of the unfolded contour length, using the different disulfide bonded constructs. Previous attempts have been made to relate protein length with the experimental rate of folding between different proteins, showing a wide range of correlation (8-11,39). This is because other factors, such as the native state stability and topology, also play important roles in the folding mechanism, camouflaging the effect of protein length alone. Moreover, these ensemble refolding experiments lack a well-defined reaction coordinate such that the initial conditions, which influence the folding rate (40), vary for different proteins as well as within the ensemble of the same protein molecules. By contrast, limiting the unfolded length beyond the transition state within the same I27 molcule, with the same initial conditions, for the first time allows for a rigorous test of the theories of length dependence on folding time.
In force-extension mode, we first unfold the amino acids within each I27 module that are not arrested by the covalent disulfide bond (up to ΔL^sub u^), as shown by the force peaks depicted in red in Fig. 1 C, such that the initial conditions are always kept the same. This is achieved by pulling the molecule a known distance away from the surface at a constant velocity of 400 nm/s. We then relax the molecule at very high velocity (4000 nm/s) and wait for the individual modules to refold during a waiting time Δ^sub t^. The molecule is then extended tor a second time to measure the fraction of protein modules (N^sub fold^/N^sub total^) in the single protein chain that refolded, that is, that regained mechanical stability. This experimental scheme is shown in Fig. 5 A for the disulfide-bonded construct (I27^sub E24C-K55C^)^sub 8^, in which n = 58 amino acids refold in all the modules after a waiting time of 0.1 s, as seen by the two consecutive unfolding traces.
As a signature of the disulfide-bonded I27 polyprotein, in the same figure we show that on addition of the reducing agent DTT (20 mM), the disulfide bond in each module is reduced, leading to an increase in ΔL^sub f^ (28.4 nm) that is consistent with the complete unraveling of all the 89 amino acids in the protein. The refolding of the 89 amino acids drastically slows down the refolding process to a rate comparable to that of the wild-type. After the same waiting time of Δ^sub t^ = 0.1 s, none of the modules refold, as seen by the absence of peaks in the second unfolding trace.
By varying the waiting time Δ^sub t^ from 20 ms to 10 s and measuring the ratio of folded modules, N^sub fold^/N^sub total^, we deduce the rate constants of refolding for the polyprotein constructs with varying n. The plot in Fig. 5 B shows the time evolution of the fraction of refolded modules, in other words, the probability of refolding. The refolding rate constants (k^sub F^) are deduced from the single exponential tits of a two-state model to the data. It is clear that the introduction of disulfide bonds that change the length of the folding protein chain drastically alters the rate of folding. The folding rate constants for the disulfide mutants are given in Table 2 and indicate a steep decrease in the folding rates with increasing length of the chain. This result is expected because of the increased conformational degrees of freedom for longer molecules. The effect of the disulfide bond restriction is observed between the rate of 46.7 s^sup -1^ to fold 58 amino acids of the I27^sub E24C-K55C^ protein, as compared to the rule of 0.3 s^sup -1^ for refolding the reduced form of the same molecule, corresponding to 89 amino acids, thus marking a > 150-fold increase in the rate of folding.
The dependence of the protein length n on folding times is shown in Fig. 6. where we plot the average refolding lime (τ^sub F^ = 1/k^sub F^ determined in Fig. 5 B) against the number of amino acids in the chain. A strong trend is immediately apparent, in contrast to previous experimental studies on diverse proteins that show a wide range of correlation. Instead, this trend can be rationalized using the polymer physics models that have been applied to the study of proteins by means of numerical simulations. Simplified models of proteins have predicted a variety of scaling laws for folding times with size, ranging from a power law (23,25,26,41,42) in an exponential dependence (23,24), depending on the roughness of the underlying free energy landscape. Because of the limited size range of the protein molecules, between 46 and 89 amino acids, it is difficult to demonstrate a sharp contrast between the distinct theories. We show that the data points are fitted equally well with a power law scaling, τ^sub F^ ~ n^sup λ^ (solid line), with the scaling exponent λ = 4.4, and with a barrier-activated process, giving τ^sub F^ [asymptotically =] exp(0.6n^sup 0.6^) (dotted line), in remarkable agreement with the different stages in protein folding in the model developed by Thirumalai (23). This does not preclude a fit with other rapidly rising functions, but we explain the data in terms of these existing models, which have been tested extensively using numerical simulations. The obtained power law exponent agrees well with the value predicted by Shakhanovich and co-workers (25) for a Go-model simulation of a designed peptide sequence but differs from other numerical simulation models (23,26,41). The scaling exponent is not universal because it depends on the energy landscape roughness of the models and therefore provides insight into the actual distribution of contact energies within the protein molecule (42).
We also investigated the influence of the proportion of native contacts, qualified by CO, on the folding rates (8). We found that CO calculated according to Baker et al. (8) from the 127 structure does not vary between the disulfide mutants (0.17-0.21). Yet, the folding times span almost three orders of magnitude (0.02-3.3 s). Therefore, we find that CO does not correlate with our observed folding times.
After the entropic collapse, the protein globule reptates, giving rise to the power law (21,22), followed by the search for a native state within a glassy energy landscape in a barrier-activated process (24). The obtained scaling with n^sup 0.6^ in the latter stage corresponds to the assumption of a Gaussian distribution of barrier heights. Most likely the two mechanisms are competing during folding, such that both play a role in the measured folding time scales. Further experiments are currently being performed to determine the rate-determining step in protein folding.
Most importantly, we have presented the first set of experiments to show that chain length plays a significant role in the rate at which proteins fold. We cannot, however, exclude the role of the polypeptide chain sequestered behind the disulfide bond during the folding process, but it clearly does not dominate over the observed length dependence. Moreover, the trend supports the view of protein folding in a complex energy landscape, in which the conformational degrees of freedom are affected by the size of the protein chain. These results emphasize the importance of single-molecule experiments in understanding the physical picture underlying the folding process.
Here we have shown how the introduction of disulfide bonds into the same protein molecule serves as the ideal control experiment to test several important biological questions: the length of the contour length per amino acid (l), the end-to-end length of the transition state of the protein, as well as the folding time (τ^sub F^) dependence on the length of the protein chain (n). This technique allows for the characterization of the unfolding transition slate in other proteins, as well as the study of their underlying protein folding dynamics.
This work has been supported by National Institutes of Health grants HL66030 and HL61228 to J.M.F. L.L. is a Damon Runyon Fellow (DRG-No. 1792-03), and J.B. holds a Career Award at the Scientific Interface from the Burroughs Wellcome Fund.
[Reference]
REFERENCES
1. Bustamante, C., J. F. Marko, E. D. Siggia, and S. Smith. 1994. Entropic elasticity of lambda-phage DNA. Science. 265:1599-1600.
2. Yang, G., C. Cecconi, W. A. Baase, I. R. Vetter, W. A. Breyer, J. A. Haack, B. W. Matthews, F. W. Dahlquist, and C. Bustamante. 2000. Solid-state synthesis and mechanical unfolding of polymers of T4 lysozyme. Proc. Natl. Acad. Sci. USA. 97:139-144.
3. Dietz, H., and M. Rief. 2006. Protein structure by mechanical triangulation. Proc. Natl. Acad. Sci. USA. 103:1244-1247.
4. Erickson, H. P. 1994. Reversible unfolding of fibronectin type III and immunoglobulin domains provides the structural basis for stretch and elasticity of titin and fibronectin. Proc. Nalt. Acad. Sci. USA. 91:10114-10118.
5. Trombitas, K., M. Greaser, S. Labeit, J. P. Jin, M. Kellermayer, M. Helmes, and H. Granzier. 1998. Titin extensibility in situ: entropic elasticity of permanently folded and permanently unfolded molecular segments. J. Cell Biol. 140:853-859.
6. Sarkar, A., S. Caamano, and J. M. Fernandez. 2005. The elasticity of individual titin PEVK exons measured by single molecule atomic force microscopy. J. Biol. Chem. 280:6261-6264.
7. Carrion-Vazquez, M., P. E. Marszulek, A. F. Oberhauser, and J. M. Fernandez. 1999. Atomic force microscopy captures length phenotypes in single proteins. Proc. Natl. Acad. Sci. USA. 96:11288-11292.
8. Plaxco, K. W., K. T. Simons, and D. Baker. 1998. Contact order, transition state placement and the refolding rates of single domain proteins. J. Mol. Biol. 277:985-994.
9. Plaxco, K. W., K. T. Simons, I. Ruczinski, and D. Baker. 2000. Topology, stability, sequence, and length: defining the determinants of two-state protein folding kinetics. Biochemistry. 39:11177-11183.
10. Li, M. S., D. K. Klimov, and D. Thirumalai. 2004. Thermal denaturation and folding rates of single domain proteins: size matters. Polym. 45:573-579.
11. Naganathan, A. N., and V. Munoz. 2005. Scaling of folding times with protein size. J. Am. Chem. Soc. 127:480-481.
12. Improta, S., A. S. Politou, and A. Pastore. 1996. Immunoglubulin-like modules from titin I-hand: extensible components of muscle elasticity. Structure. 4:323-337.
13. Carrion-Vazquez, M., A. F. Oberhauser, S. B. Fowler, P. E. Marszalek, S. E. Broedel, J. Clarke, and J. M. Fernandez. 1999. Mechanical and chemical unfolding of a single protein: a comparison. Proc. Natl. Acad. Sci. USA. 96:3694-3699.
14. Fisher, T. E., A. F. Oberhauser, M. Carrion-Vazquez, P. E. Marszalek, and J. M. Fernandez. 1999. The study of protein mechanics with the atomic force microscope. Trends Biochem. Sci. 24:379-384.
15. Carrion-Vazquez, M., H. 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.
16. Oesterhelt, F., D. Oesterhelt, M. Pfeiffer, A. Engel, H. E. Gaub, and D. J. Muller. 2000. Unfolding pathways of individual bacteriorhodopsins. Science. 288:143-146.
17. Li, H., and J. M. Fernandez. 2003. Mechanical design of the first proximal Ig domain of human cardiac titin revealed by single molecule force spectroscopy. J. Mol. Biol. 334:75-86.
18. Bhasin, N., P. Carl, S. Harper, G. Feng, H. Lu. D. W. Speicher, and D. E. Discher. 2004. Chemistry on a single protein, vascular cell adhesion molecule-1, during forced unfolding. J. Biol. Chem. 279:45865-45874.
19. Bustanji, Y., and B. Samori. 2002. The mechanical properties of human angiostatin can be modulated by means of its disulfide bonds: a single-molecule force-spectroscopy study. Angew. Chem. Int. Ed. Engl. 41:1546-1548.
20. Carl, P., C. H. Kwok, G. Manderson, D. W. Speicher, and D. E. Discher. 2001. Forced unfolding modulated by disulfide bonds in the Ig domains of a cell adhesion molecule. Proc. Natl. Acad. Sci. USA. 98:1565-1570.
21. de Gennes, P. G. 1985. Kinetics of collapse for a flexible coil. J. Phys. Lett. 46:L639-L642.
22. Grosberg, A. Y., S. K. Nechaev, and E. I. Shakhanovich. 1988. The role of topolopical constraints in the kinetics of collapse of macromolecules. J. Phys. France. 49:2095-2100.
23. Thirumalai, D. 1995. From minimal models to real proteins: time scales for protein folding kinetics. J. Phys. I (France). 5:1457-1467.
24. Bryngelson, J. D., and P. G. Wolynes. 1987. Spin glasses and the statistical mechanics of protein folding. Proc. Natl. Acad. Sci. USA. 84:7524-7528.
25. Gutin, A. M., V. V. Abkevich, and E. I. Shakhanovich. 1996. Chain length scaling of protein folding time. Phys. Rev. Lett. 77:5433-5436.
26. Cieplak, M., and T. X. Hoang. 2003. Universality classes in folding times of proteins. Biophys. J. 84:475-488.
27. Cieplak, M., T. X. Hoang, and M. O. Robbins. 2004. Stretching of proteins in the entropic limit. Phys. Rev. E. 69:011912.
28. Hazes, B., and B. W. Dijkstra. 1988. Model building of disulfide bonds in proteins with known three-dimensional structure. Protein Eng. 2:119-125.
29. Oberhauser, A. F., P. K. Hansma, M. Carrion-Vazquez, and J. M. Fernandez. 2001. Stepwise unfolding of titin under force-clamp atomic force microscopy. Proc. Natl. Acad. Sci. USA. 98:468-472.
30. Florin, E. L., M. Rief, H. Lehmann, M. Ludwig, K. Dornmair, V. T. Moy, and H. E. Gaub. 1995. Sensing specific molecular interactions with the atomic force microscope. Biosens. Bioelectron. 10:895-901.
31. Grandbois, M., M. Beyer, M. Rief, H. Clausen-Schaumann, and H. E. Gaub. 1999. How strong is a covalent bond? Science. 283:1727-1730.
32. Wiita, A. P., S. R. K. Ainavarapu, H. H. Huang, and J. M. Fernandez. 2006. Force-dependent chemical kinetics of disulfide bond reduction observed with single-molecule techniques. Proc. Natl. Acad. Sci. USA. 103:7222-7227.
33. Pauling, L., and R. B. Corey. 1951. The pleated sheet, a new layer configuration of polypeptide chains. Proc. Natl. Acad. Sci. USA. 37:251-256.
34. Lu, H., B. Isralewitz, A. Krammer, V. Vogel, and K. Schulten. 1998. Unfolding of titin immunoglobulin domains by steered molecular dynamics simulation. Biophys. J. 75:662-671.
35. Li, H., M Carrion-Vazquez, A. F. Oberhauser, P. E. Marszalek, and J. M. Fernandez. 2000. Point mutations after the mechanical stability of immunoglobulin modules. Nat Struct. Biol. 7:1117-1120.
36. Katz, B. A., and A. Kossiakoff. 1986. The crystallographically determined structures of atypical strained disulfides engineered into subtilisin. J. Biol. Chem. 261:15480-15485.
37. Lu, H., and K. Schulten. 2000. The key event in force-induced unfolding of Titin's immunoglobulin domains. Biophys. J. 79:51-65.
38. Lide, D. R. 1995. CRC Handbook of Chemistry and Physics. CRC Press, Boca Raton, FL. 9-74.
39. Ivankov, D. N., S. O. Garbuzynskiy, E. Alm, K. W. Plaxco, D. Baker, and A. V. Finkelstein. 2003. Contact order revisited: influence of protein size on the folding rate. Protien Sci. 12:2057-2062.
40. Li, M. S., C. K. Hu, D. K. Klimov, and D. Thirumalai. 2006. Multiple stepwise refolding of immunoglobulin domain 127 upon force quench depends on initial conditions. Proc. Natl. Acad. Sci. USA. 103:93-98.
41. Zhdanov, V. P. 1998. Folding time of ideal β-sheet vs. chain length. Europhys. Lett. 42:577-581.
42. Cieplak, M., T. X. Hoang, and M. S. Li. 1999. Scaling of folding properties in simple models of proteins. Phys. Rev. Lett. 83:1684-1687.
[Author Affiliation]
Sri Rama Koti Ainavarapu,* Jasna Brujic,* Hector H. Huang,* Arun P. Wiita.* Hui Lu,[double dagger] Lewyn Li,* Kirstin A. Walther,*[dagger] Mariano Carrion-Vazquez,� Hongbin Li,� and Julio M. Fernandez*
* Department of Biological Sciences and [dagger] Department of Physics, Columbia University, New York, New York 10027; [double dagger] Department of Bioengineering, University of Illinois, Chicago, Illinois 60607; � Institute Cajal, Consejo Superior de Investigaciones Cientificas, 28002 Madrid, Spain; and � University of British Columbia, Vancouver, British Columbia, Canada
[Author Affiliation]
Submitted June 17, 2006, and accepted for publication September 13, 2006.
Sri Rama Koti Ainavarapu, Jasna Brujic, and Hector H. Huang contributed equally to this work.
Address reprint requests to J. M. Fernandez, Dept. of Biological Sciences, Columbia University, New York, NY 10027.E-mail: jfernandez@Columbia. edu.
