| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||




* Department of Physics,
Department of Biological Sciences, Columbia University, New York, New York; and
Department of Chemistry, University of British Columbia, Vancouver, British Columbia, Canada
Correspondence: Address reprint requests to J. M. Fernández, E-mail: jfernandez{at}columbia.edu.
| ABSTRACT |
|---|
|
|
|---|
| INTRODUCTION |
|---|
|
|
|---|
For the measured variance to be a useful quantity, equilibrium must first be established. Although in the past few years, numerous researchers have investigated the appropriate statistical mechanics relevant to force-spectroscopy experiments (15
21
), their theories concerning the establishment of thermal equilibrium have not yet been applied in experiments. In previous experiments, the existence of thermal equilibrium has been primarily determined by the absence or presence of hysteresis in the force-extension relation. However, the absence of hysteresis does not necessarily imply equilibrium with the bath (e.g., the parameter range may not be large enough to observe hysteresis). Additionally, interpretation of the measured molecular variance must also proceed with caution. This is especially true for measurement systems employing active feedback. The response time of the feedback must be fast enough to track the molecule's thermal motion. These considerations have not been fully appreciated in previous experiments, where it was crucial to know the thermodynamic state of the system (22
,23
).
Here, we report a variance analysis of the stretching of single dextran molecules with an AFM. First, the variance in length of dextran is used to prove thermal equilibrium of the molecule with the bath. Because the system is in equilibrium, we can then confidently use the variance to probe the conformational changes that occur when dextran's subunits, single pyranose rings, flip from the chair to the boat conformation under force (24
,25
). These conformational changes occur on length changes smaller than one Angstrom.
| MATERIALS AND METHODS |
|---|
|
|
|---|
50 pN nm1. Dextran molecules were manually picked up by pushing the cantilever onto the coverslip for up to several minutes. The piezo was then retracted while force and extension were observed in real-time on the oscilloscope. Once the user was confident that a single molecule was attached, a program was started to stretch and relax the molecule up to several hundred times before it detached. All experiments were performed at room temperature.
For the force-clamp experiments, the force was linearly changed at a rate of 1000 pN s1 from a minimum to a maximum force and back, encompassing the transition region. We verified that this pulling rate results in the same force-extension relation as slower pulling rates (29
). Furthermore, this pulling rate results in smaller piezo velocities than the maximum velocity for near-equilibrium conditions (34
).
For the constant-velocity experiments, the extension of the piezo was manually adjusted depending on the length of the molecule to capture the transition region. The pulling rate was typically 250 nm s1. To rule out any effect of the pulling rate on the variance, one molecule was pulled at both 150 nm s1 and 250 nm s1. The sampling rate was set around 12.5 kHz to acquire 10,000 data points over a period of
0.8 s, depending on the length of the dextran molecule. The anti-aliasing filter (eight-pole Bessel) was set to a fraction of the Nyquist frequency (6.25 kHz). Varying this cutoff frequency between 50 and 90% of the Nyquist frequency had no effect on the measured variance. However, at 7% of the Nyquist frequency, the variance was reduced by a factor of 2 and no increase in the transition region was observed.
Constant-velocity experiments do not involve any feedback that may interfere with the molecule's thermal equilibrium; however, instrumental noise could potentially still be a problem. The presence of electronic noise in the piezoelectric actuator should create fluctuations independent of the conformational state of the molecule, adding a roughly constant value to the variance independent of the molecule's extension. Based on the actual variance of the pulls alone, this value is small enough that we can still observe a change in the variance with respect to the extension of the molecule. More significantly, the striking accord between the actual variance and the predicted variance assuming thermodynamic equilibrium implies that any instrumental noise is completely negligible. We are therefore confident that there was no impact of the instrumentation on our measurements of the variance.
Data analysis
Calculation of the variance from all traces for a particular molecule was carried out in the following manner. Inherent drift in the piezo position and force over the time of the experiment required shifting all traces such that the transition regions overlapped. Each individual trace was smoothed to obtain the mean, which was then subtracted from the original trace, leaving only the fluctuations. These fluctuations were then squared, organized into bins of an appropriate size, and averaged over all traces for one molecule to obtain the variance in the length of the molecule (force-clamp experiments) or in the force on the molecule (constant-velocity experiments). For the latter case, the variance in the force was converted to the length variance using the spring constant of the cantilever (15
).
Force-extension curves were fit using the modified freely jointed chain (FJC) model, which is based on the Langevin equation but includes segment elasticity. The fits were performed simultaneously to the low-force (chair conformation) and the high-force (boat conformation) regime using the Levenberg-Marquardt method (30
). In the fits, the length of the pyranose ring in each conformation is fixed to the value obtained from ab initio calculations: lk, chair = 4.412 Å and lk, boat = 5.696 Å (25
). This leaves N, Se, chair, and Se, boat to be determined by the fit. The resulting N is rounded to the closest integer because it represents an integer quantity (number of rings in the dextran molecule). The standard deviation for N is at most five rings, and depends on how accurately the segment elasticity of the boat could be determined from the force-extension curve. This uncertainty in N does not significantly affect the results of the fits to the variance.
The variance was fit using our two- and three-state models. Contributions to the variance from the FJC elasticity are taken into account by a straight line connecting the value of the variance at the beginning and at the end of the transition region. Based on this straight line, we estimate the elasticity due to the FJC and add this to the elasticity due to the transitions between the chair and the boat conformations. In the fit, the number of rings, N, is fixed to the value obtained from the fit to the force-extension curve. This leaves
lbi,
lbc, and B to be determined by the fit to the variance.
Dextran
We used dextran T500 (Pharmacia Biotech, Uppsala, Sweden), which is a >95% linear homopolymer consisting of
-(1
6) D-glucopyranose rings.
The sample was prepared by dissolving dextran in water at a concentration of 0.11% (wt/vol). It was then applied onto a clean glass coverslip and dried overnight. By rinsing extensively with water, excessive dextran was washed off so that only a thin layer of molecules tightly adsorbed to the surface remained. All experiments were carried out in phosphate-buffered saline solution or in water.
| RESULTS AND DISCUSSION |
|---|
|
|
|---|
We first describe our force-clamp experiments, in which the force on the molecule is externally controlled. Fig. 1 A shows a typical recording of a single dextran molecule (comprised of many pyranose rings) as it is stretched (blue) and relaxed (red) throughout the chair-boat transition region by ramping the force while applying an electronic feedback (27
,29
). At low forces, the pyranose rings are in their ground-state chair conformation, and the molecule's elasticity is well described by a modified FJC model (31
). The accelerated increase in length between 800 and 1000 pN characterizes the transition region as the pyranose rings flip from the chair to the boat conformation. This change results in an elongation of the pyranose rings by 20% (25
).
|
If the system is in thermal equilibrium in the Gibbs ensemble, the variance can also be obtained from the equipartition theorem, using the slope of the force versus extension curve (15
),
![]() | (1) |
...
indicates an average, kB is the Boltzmann constant, and T is the temperature of the bath, and where an effective spring constant, km = d
F
/d
L
, describes the macroscopic elasticity of the molecule (see Appendix A). This predicted variance is shown as the pink curve in Fig. 1 B. The experimental and predicted variances are highly discrepant. In addition, the force variance increases in the transition region (black curve in Fig. 1 B). In the derivation of the variance, we assume perfect force feedback where the force is held constant with infinite time resolution. However, the molecule is clearly fluctuating on timescales faster than the response time of the electronic feedback, which has a bandwidth at
600 Hz. To fully account for the effect of the feedback we would need to know its complete transfer function, including the cantilever, the electronic feedback, the piezoelectric actuator, and the molecule, in order to predict the measured variance. This evaluation is most difficult and would have to be done for each individual cantilever-molecule combination. Therefore, we cannot confirm that the molecule is in equilibrium despite the reversibility of its force-extension relation. For these reasons, further information cannot be obtained from the variance under these conditions.
We note that the finite response time of the force-feedback has been discussed in previous publications (4
). However, for proteins, typical unfolding and folding time-series extend over several seconds. In this case, the millisecond response time of the feedback is trivial in comparison and therefore does not interfere with those previous measurements.
Constant-velocity experiments
Since thermal equilibrium cannot be confirmed in our force-clamp experiments, we consequently repeat the experiment under constant-velocity conditions where feedback is not employed. Fig. 2 A shows a typical recording of a single dextran molecule. The variance is calculated from the fluctuations in the measured force from all pulls of the molecule. This force variance is then converted into the variance in the length of the molecule using the spring constant of the cantilever. The resulting variance is shown as the green curve in Fig. 2 B.
|
![]() | (2) |
F
/d
L
(see Appendix A). This predicted variance is shown as the pink curve in Fig. 2 B. The exceptional agreement between the experimental and the predicted variance confirms that the system is indeed in equilibrium in the Helmholtz ensemble.
Two-state model fails to describe chair-boat transition
Because equilibrium has been established, the molecule's variance enables us to further investigate its conformational changes. Previously, the chair-boat transition in dextran was assumed to be a two-state process (25
,29
). We now examine whether this assumption is indeed appropriate. The mean extension of the molecule is
L
=N(1pc)
lbc, where L is measured relative to the beginning of the transition region, N is the total number of pyranose rings in the dextran molecule, pc is the probability that a ring is in the chair conformation, and
lbc = lboatlchair is the difference in length between the boat and the chair conformation. From the binomial distribution, the variance is
in the Gibbs ensemble. Eliminating pc from these two equations provides a functional relationship between the variance and the mean extension of the molecule, in analogy to standard ion channel variance analysis (9
):
![]() | (3) |
The theory assumes that each ring is statistically independent of all other rings.
The expression for the variance given by Eq. 3 is valid for the Gibbs ensemble, but our experiments were carried out in the Helmholtz ensemble. Fortunately, there exists a simple relationship between the variances in both ensembles. The stiffness of the molecule, km, is the same in each case, allowing us to combine Eqs. 1 and 2 to obtain
![]() | (4) |

L2
G is given by Eq. 3.
The resulting variance for a two-state system in the Helmholtz ensemble is shown as the red curve in Fig. 3 B. The number of rings for this molecule was determined to be N = 142 (rounded to an integer) from simultaneously fitting the low-force (chair) and the high-force (boat) regime of the force-extension curve with the FJC model (Fig. 3 A) using respective chair and boat lengths of 4.412 and 5.696 Å (
lbc, ai = 1.284 Å) from ab initio calculations (25
) (see Materials and Methods). Fixing N at this value, the fit of the two-state variance without (red dashed curve) and with (red solid curve) taking into account the FJC elasticity (see Materials and Methods) to the experimental variance (green curve) resulted in a value of 0.99 Å for the length difference between boat and chair,
lbc, when the widths of the two variances corresponded to each other. The magnitude of the variance of the two-state model in either case clearly disagrees with the experimental variance. This deviation was similar for all nine molecules that we observed. We therefore conclude that the chair-boat transition in dextran is not a two-state process.
|
![]() | (5) |
lbc = lboat lchair is again the difference in length between boat and chair, and
lbi = lboat lintermediate is the difference in length between boat and intermediate. We then assumed the following functional form for pi based on known constraints about its shape (see Appendix B),
![]() | (6) |
Our fit of this model to the molecule shown in Fig. 3 is executed similarly to the previous fit of the two-state model (see Materials and Methods). We again take the number of rings to be N = 142 from fits of the FJC to the force-extension curve. From the variance, we obtain the following values for the three fitting parameters (black curve):
lbi = 0.488 ± 0.005 Å,
lbc = 0.994 ± 0.004 Å, and B = 0.617 ± 0.010. Note that the fitted value for B is indeed <1 as required for a probability. In general, we find that the width of the transition region and the number of rings, N, essentially determine
lbc, the slope of the variance at the end of the transition region determines
lbi, and the maximum value of the variance determines B.
We analyzed a total of nine molecules, including the one shown in Fig. 3 and the three shown in Fig. 4. From the FJC fits to their force-extension curves, we find that these molecules had widely different lengths, ranging from N = 74289. Despite their different lengths, we obtained fairly narrow distributions for the three fitting parameters (see Table 1) with average values 
lbi
= 0.518 ± 0.018 Å, 
lbc
= 0.988 ± 0.020 Å, and
B
= 0.626 ± 0.037.
|
|

lbc
= 0.988 Å is significantly smaller than the value from ab initio calculations,
lcb, ai = 1.284 Å, used to fit the force-extension curves (see above). This is likely because the ab initio calculations considered the pyranose ring in equilibrium under no external force, whereas our molecule was under a stretching force in the transition region. From the slope of the force-extension curves, we know that the molecule's elasticity coefficient is smaller in the chair than in the boat conformation (i.e., the boat is stiffer than the chair). Therefore, the length difference between boat and chair should decrease under force, as observed.
Our results indicate an intermediate state that lies approximately half-way between the chair and the boat. Interestingly, recent SMD simulations of dextran suggest the existence of an intermediate conformation corresponding to simultaneous rotation of the protruding carbon-oxygen bond and transition of the pyranose ring to a boat structure distinct from the final one (32
).
| CONCLUSION |
|---|
|
|
|---|
This analysis should be equally applicable to any stretching experiment in order to probe the thermodynamic state of the system as well as equilibration timescales. In particular, future tests of fundamental thermodynamic theories such as the Jarzynski's equality should always start with a robust verification of nonequilibrium conditions, as suggested previously (2
,23
).
Finally, our analysis should allow for the investigation of conformational changes in a wide variety of molecules at significantly higher resolution than is now possible. These include simple molecules with a ring structure similar to dextran (1
) as well as more complex molecules like DNA (33
), RNA (2
), and proteins (26
) that undergo discrete conformational changes under a stretching force.
| APPENDIX A: MACROSCOPIC ELASTICITY OF THE MOLECULE |
|---|
|
|
|---|
![]() | (7) |
...
indicates an average. Substituting
we find that the second term in the series is at most 7% of the leading order term for all molecules, thereby validating our assumption. | APPENDIX B: VARIANCE FOR THREE-STATE MODEL |
|---|
|
|
|---|
![]() | (8) |
![]() | (9) |
lbi and
lbc, respectively, are the differences in length between boat-intermediate and boat-chair. In contrast to the two-state theory, it is not possible to eliminate both probabilities from the variance. We can, however, choose which probability to keep in Eq. 9. We decided to keep the probability of the intermediate, pi, which yields the relation given by Eq. 5. We chose pi because there are some constraints on its shape as a function of the average length of the molecule:
1.
lbi. A simple function fulfilling these conditions is Axc(1 x) with x in the interval [0, 1]. In our particular case, this becomes Eq. 6, given in the main text.
| ACKNOWLEDGEMENTS |
|---|
|
|
|---|
This work has been supported by National Institutes of Health grant No. R01 HL66030 to J.M.F.
Submitted on October 18, 2005; accepted for publication January 12, 2006.
| REFERENCES |
|---|
|
|
|---|
2. Liphardt, J., B. Onoa, S. B. Smith, I. Tinoco, and C. Bustamante. 2001. Reversible unfolding of single RNA molecules by mechanical force. Science. 292:733737.
3. Bustamante, C., Z. Bryant, and S. B. Smith. 2003. Ten years of tension: single-molecule DNA mechanics. Nature. 421:423427.[CrossRef][Medline]
4. Fernandez, J. M., and H. Li. 2004. Force-clamp spectroscopy monitors the folding trajectory of a single protein. Science. 303:16741678.
5. Zech, S. G., A. J. Wand, and A. E. McDermott. 2005. Protein structure determination by high-resolution solid-state NMR spectroscopy: application to microcrystalline ubiquitin. J. Am. Chem. Soc. 127:86188626.[CrossRef][Medline]
6. Jelsch, C., M. M. Teeter, V. Lamzin, V. Pichon-Pesme, R. H. Blessing, and C. Lecomte. 2000. Accurate protein crystallography at ultra-high resolution: valence electron distribution in crambin. Proc. Natl. Acad. Sci. USA. 97:31713176.
7. Nellist, P. D., M. F. Chisholm, N. Dellby, O. L. Krivanek, M. F. Murfitt, Z. S. Szilagyi, A. R. Lupini, A. Borisevich, W. H. Sides, Jr., and S. J. Pennycook. 2004. Direct sub-Angstrom imaging of a crystal lattice. Science. 305:1741.
8. Einstein, A. 1905. On the movement of small particles suspended in a stationary liquid demanded by the molecular-kinetic theory of heat. Annalen der Physik. 17:549560.
9. Sigworth, F. J. 1980. The variance of sodium current fluctuations at the node of Ranvier. J. Physiol. (Lond.). 307:97129.
10. Mehta, A. D., J. T. Finer, and J. A. Spudich. 1997. Detection of single-molecule interactions using correlated thermal diffusion. Proc. Natl. Acad. Sci. USA. 94:79277931.
11. Liu, Y. Z., S. H. Leuba, and S. M. Lindsay. 1999. Relationship between stiffness and force in single molecule pulling experiments. Langmuir. 15:85478548.[CrossRef]
12. Humphris, A. D. L., J. Tamayo, and M. J. Miles. 2000. Active quality factor control in liquids for force spectroscopy. Langmuir. 16:78917894.[CrossRef]
13. Helfer, E., S. Harlepp, L. Bourdieu, J. Robert, F. C. MacKintosh, and D. Chatenay. 2001. Viscoelastic properties of actin-coated membranes. Phys. Rev. E. 63:021904.[CrossRef]
14. Kawakami, M., K. Byrne, B. Khatri, T. C. B. Mcleish, S. E. Radford, and D. A. Smith. 2004. Viscoelastic properties of single polysaccharide molecules determined by analysis of thermally driven oscillations of an atomic force microscope cantilever. Langmuir. 20:92999303.[CrossRef][Medline]
15. Kreuzer, H. J., S. H. Payne, and L. Livadaru. 2001. Stretching a macromolecule in an atomic force microscope: statistical mechanical analysis. Biophys. J. 80:25052514.
16. Titantah, J. T., C. Pierleoni, and J. P. Ryckaert. 1999. Different statistical mechanical ensembles for a stretched polymer. Phys. Rev. E. 60:70107021.[CrossRef]
17. Gerland, U., R. Bundschuh, and T. Hwa. 2001. Force-induced denaturation of RNA. Biophys. J. 81:13241332.
18. Ritort, F., C. Bustamante, and I. Tinoco. 2002. A two-state kinetic model for the unfolding of single molecules by mechanical force. Proc. Natl. Acad. Sci. USA. 99:1354413548.
19. Gerland, U., R. Bundschuh, and T. Hwa. 2003. Mechanically probing the folding pathway of single RNA molecules. Biophys. J. 84:28312840.
20. Keller, D., D. Swigon, and C. Bustamante. 2003. Relating single-molecule measurements to thermodynamics. Biophys. J. 84:733738.
21. Tinoco, I. 2004. Force as a useful variable in reactions: unfolding RNA. Annu. Rev. Biophys. Biomol. 33:363385.[CrossRef]
22. Liphardt, J., S. Dumont, S. B. Smith, I. Tinoco, and C. Bustamante. 2002. Equilibrium information from nonequilibrium measurements in an experimental test of Jarzynski's equality. Science. 296:18321835.
23. Cohen, E. G. D., and D. Mauzerall. 2004. A note on the Jarzynski equality. J. Stat. Mech. Theory E. P07006.
24. Barton, D. H. R. 1970. Principles of conformational analysis. Science. 169:539544.
25. Marszalek, P. E., A. F. Oberhauser, Y.-P. Pang, and J. M. Fernandez. 1998. Polysaccharide elasticity governed by chair-boat transitions of the glucopyranose ring. Nature. 396:661664.[CrossRef][Medline]
26. Rief, M., M. Gautel, F. Oesterhelt, J. M. Fernandez, and H. E. Gaub. 1997. Reversible unfolding of individual titin immunoglobulin domains by AFM. Science. 276:11091112.
27. Schlierf, M., H. Li, and J. M. Fernandez. 2004. The unfolding kinetics of ubiquitin captured with single-molecule force-clamp techniques. Proc. Natl. Acad. Sci. USA. 101:72997304.
28. Florin, E. L., M. Rief, H. Lehmann, M. Ludwig, C. Dornmair, V. T. Moy, and H. E. Gaub. 1995. Sensing specific molecular-interactions with the atomic-force microscope. Biosens. Bioelectron. 10:895901.[CrossRef]
29. Marszalek, P. E., H. Li, A. F. Oberhauser, and J. M. Fernandez. 2002. Chair-boat transitions in single polysaccharide molecules observed with force-ramp AFM. Proc. Natl. Acad. Sci. USA. 99:42784283.
30. Press, W. H., S. A. Teukolsky, W. T. Vetterling, and B. P. Flannery. 1992. Numerical Recipes in C. Cambridge University Press, Cambridge, UK.
31. Rief, M., F. Oesterhelt, B. Heymann, and H. E. Gaub. 1997. Single molecule force spectroscopy on polysaccharides by atomic force microscopy. Science. 275:12951297.
32. Lee, G., W. Nowak, J. Jaroniec, Q. M. Zhang, and P. E. Marszalek. 2004. Molecular dynamics simulations of forced conformational transitions in 1,6-linked polysaccharides. Biophys. J. 87:14561465.
33. Smith, S. B., Y. J. Cui, and C. Bustamante. 1996. Overstretching B-DNA: the elastic response of individual double-stranded and single-stranded DNA molecules. Science. 271:795799.[Abstract]
34. Braun, O., and U. Seifert. 2005. Force spectroscopy of single multidomain biopolymers: a master equation approach. Eur. Phys. J. E. 18:113.[CrossRef][Medline]
This article has been cited by other articles:
![]() |
D. B. Staple, F. Hanke, and H. J. Kreuzer Comment on "Sub-Angstrom Conformational Changes of a Single Molecule Captured by AFM Variance Analysis" Biophys. J., July 15, 2008; 95(2): 1001 - 1002. [Full Text] [PDF] |
||||
![]() |
S. Garcia-Manyes, J. Brujic, C. L. Badilla, and J. M. Fernandez Force-Clamp Spectroscopy of Single-Protein Monomers Reveals the Individual Unfolding and Folding Pathways of I27 and Ubiquitin Biophys. J., October 1, 2007; 93(7): 2436 - 2446. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. T. Woodside, P. C. Anthony, W. M. Behnke-Parks, K. Larizadeh, D. Herschlag, and S. M. Block Direct measurement of the full, sequence-dependent folding landscape of a nucleic Acid. Science, November 10, 2006; 314(5801): 1001 - 1004. [Abstract] [Full Text] [PDF] |
||||
![]() |
Z. Lu, H. Hu, W. Yang, and P. E. Marszalek Simulating Force-Induced Conformational Transitions in Polysaccharides with the SMD Replica Exchange Method Biophys. J., September 15, 2006; 91(6): L57 - L59. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |