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

Originally published as Biophys J. BioFAST on February 24, 2006.
doi:10.1529/biophysj.105.076224
This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow All Versions of this Article:
biophysj.105.076224v1
90/10/3806    most recent
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Walther, K. A.
Right arrow Articles by Fernández, J. M.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Walther, K. A.
Right arrow Articles by Fernández, J. M.
Biophysical Journal 90:3806-3812 (2006)
© 2006 The Biophysical Society

Sub-Angstrom Conformational Changes of a Single Molecule Captured by AFM Variance Analysis

Kirstin A. Walther * {dagger}, Jasna Brujic {dagger}, Hongbin Li {ddagger} and Julio M. Fernández {dagger}

* Department of Physics, {dagger} Department of Biological Sciences, Columbia University, New York, New York; and {ddagger} 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
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS AND DISCUSSION
 CONCLUSION
 APPENDIX A: MACROSCOPIC...
 APPENDIX B: VARIANCE FOR...
 ACKNOWLEDGEMENTS
 REFERENCES
 
A system's equilibrium variance can be analyzed to probe its underlying dynamics at higher resolution. Here, using single-molecule atomic-force microscope techniques, we show how the variance in the length of a single dextran molecule can be used to establish thermodynamic equilibrium and to detect conformational changes not directly observable with other methods. Dextran is comprised of a chain of pyranose rings that each undergoes an Angstrom-scale transition from a chair to boat conformation under a stretching force. Our analysis of the variance of the molecule's fluctuations verifies equilibrium throughout the force-extension curve, consistent with the expected thermodynamic ensemble. This validates further analysis of the variance in the transition region, which reveals an intermediate conformation between the chair and the boat on the sub-Angstrom scale. Our test of thermal equilibrium as well as our variance analysis can be readily extended to a wide variety of molecules, including proteins.


    INTRODUCTION
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS AND DISCUSSION
 CONCLUSION
 APPENDIX A: MACROSCOPIC...
 APPENDIX B: VARIANCE FOR...
 ACKNOWLEDGEMENTS
 REFERENCES
 
Recent single-molecule force-spectroscopy experiments have provided insight into the mechanical stability and conformational changes of molecules such as proteins, polysaccharides, DNA, and RNA (1Go–4Go). However, the resolution of these experiments has been limited to several nanometers. Sub-Angstrom resolution can be reached with nuclear magnetic resonance (NMR) spectroscopy (5Go), x-ray crystallography (6Go), and scanning transmission electron microscopy (7Go). However, these techniques require large ensembles of molecules that are all in the same configuration. It is therefore difficult to observe conformations of a molecule that are transient, dynamic, and only occupied by a fraction of the ensemble. Einstein's 1905 work on Brownian motion famously demonstrated how a system's equilibrium variance can reveal its underlying dynamics at higher resolution (8Go). A prominent example of this principle used in biology was the estimation of single ion-channel conductance from the variance of the macroscopic membrane current (9Go) before the era of patch-clamp methods. This type of analysis has never before been used in force-spectroscopy experiments even though fluctuations have previously been quantified (10Go–14Go). Here, we demonstrate a fundamental thermodynamic test of equilibrium using variance. This is a prerequisite for further variance analysis of a single dextran molecule stretched with an atomic force microscope (AFM) to observe the dynamic conformational changes of its ring subunits with sub-Angstrom resolution. The method we describe here represents the most sensitive approach for capturing dynamic conformations in a single molecule to date.

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 (15Go–21Go), 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 (22Go,23Go).

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 (24Go,25Go). These conformational changes occur on length changes smaller than one Angstrom.


    MATERIALS AND METHODS
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS AND DISCUSSION
 CONCLUSION
 APPENDIX A: MACROSCOPIC...
 APPENDIX B: VARIANCE FOR...
 ACKNOWLEDGEMENTS
 REFERENCES
 
Single-molecule atomic force microscopy
We used a custom-made AFM under constant-velocity and force-clamp conditions (26Go,27Go). Each cantilever (Si3N4 from Digital Instruments, Santa Barbara, CA) was calibrated in solution by applying the equipartition theorem (28Go). The spring constant was typically found to be ~50 pN nm–1.

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 s–1 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 (29Go). Furthermore, this pulling rate results in smaller piezo velocities than the maximum velocity for near-equilibrium conditions (34Go).

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 s–1. To rule out any effect of the pulling rate on the variance, one molecule was pulled at both 150 nm s–1 and 250 nm s–1. 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 (15Go).

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 (30Go). 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 Å (25Go). 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 {Delta}lbi, {Delta}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 {alpha}-(1 -> 6) D-glucopyranose rings.

The sample was prepared by dissolving dextran in water at a concentration of 0.1–1% (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
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS AND DISCUSSION
 CONCLUSION
 APPENDIX A: MACROSCOPIC...
 APPENDIX B: VARIANCE FOR...
 ACKNOWLEDGEMENTS
 REFERENCES
 
Force-clamp experiments
In typical single-molecule stretching experiments, either the distance between the cantilever and coverslip is externally controlled (constant-velocity) (25Go), or the force on the molecule (force-clamp) (27Go). In thermal equilibrium, these methods respectively yield the Helmholtz or the Gibbs ensemble. Below, we discuss two sets of stretching experiments that employ force-clamp and constant-velocity conditions.

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 (27Go,29Go). 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 (31Go). 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% (25Go).


Figure 1
View larger version (26K):
[in this window]
[in a new window]
 
FIGURE 1  Variance reveals that force-clamp experiments are not in thermal equilibrium in the Gibbs ensemble. (A) A single dextran molecule was stretched (blue) and relaxed (red) consecutively by linearly varying the force. The figure shows the first three out of 108 cycles. In the transition region, between 800 and 1000 pN, the pyranose rings in dextran are thought to flip from a chair to a boat conformation, corresponding to a 20% increase in length (25Go). The noise in the length increases markedly in the transition region. (B) The experimental (green) and predicted (pink) length variances are highly discrepant, indicating that the molecule is not in thermal equilibrium in the Gibbs ensemble.

 
The fluctuations in the length of the molecule are given by the difference between the measured extension and its mean. The variance is then calculated by squaring the fluctuations and averaging over all traces for one molecule. The resulting variance for the molecule from Fig. 1 A is shown as the green curve in Fig. 1 B.

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 (15Go),

Formula 1(1)
where <...> 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 (4Go). 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.


Figure 2
View larger version (25K):
[in this window]
[in a new window]
 
FIGURE 2  Variance reveals that constant-velocity experiments are in thermal equilibrium in the Helmholtz ensemble. (A) A single dextran molecule was stretched (blue) and relaxed (red) consecutively at a rate of 250 nm s–1. The figure shows the first four out of 667 cycles. (B) The experimental (green) and predicted (pink) length variances are in very good agreement, indicating that the molecule is indeed in thermal equilibrium in the Helmholtz ensemble.

 
Applying a similar analysis as for the Gibbs ensemble above, the predicted variance of the coupled cantilever-molecule system in equilibrium in the Helmholtz ensemble is (15Go)

Formula 2(2)
with kc as the spring constant of the cantilever and km = d<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 (25Go,29Go). We now examine whether this assumption is indeed appropriate. The mean extension of the molecule is <L>=N(1–pc){Delta}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 {Delta}lbc = lboatlchair is the difference in length between the boat and the chair conformation. From the binomial distribution, the variance is Formula 2 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 (9Go):

Formula 3(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

Formula 4(4)
where <{delta}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 Å ({Delta}lbc, ai = 1.284 Å) from ab initio calculations (25Go) (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, {Delta}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.


Figure 3
View larger version (19K):
[in this window]
[in a new window]
 
FIGURE 3  Two-state versus three-state model. (A) The force-extension curve (blue) from one molecule is fit in the low-force (chair) and high-force (boat) regime simultaneously with the modified FJC. In the fit (black), the Kuhn segment lengths of the chair and the boat, lk,chair and lk,boat, were fixed to 4.412 and 5.696 Å, as given by ab initio calculations (25Go). The number of rings, N = 142, and the segment elasticities, Se, chair = 11, 519 pN nm–1 and Se, boat = 90, 230 pN nm–1 (all values rounded to closest integer), are the results of the fit. (B) The magnitude of the predicted variance from the two-state model (red) clearly disagrees with the experimental variance (green) (fixing N = 142 from the FJC fit). A simple three-state model (black), however, results in good agreement with the experimental variance (again fixing N = 142). The results of the fit are {Delta}lbi = 0.488 ± 0.005 Å, {Delta}lbc = 0.994 ± 0.004 Å, and B = 0.617 ± 0.010.

 
Variance reveals sub-Angstrom intermediate
We find, however, that a three-state model fits the experimental variance very well. This model introduces an intermediate state between the chair and the boat. The variance in the Gibbs ensemble using the trinomial distribution is (see Appendix B):

Formula 5(5)
Here, pi is the probability that a given ring is in the intermediate conformation, {Delta}lbc = lboat lchair is again the difference in length between boat and chair, and {Delta}lbi = lboatlintermediate 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),

Formula 6(6)
where B is the maximum value of the probability.

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): {Delta}lbi = 0.488 ± 0.005 Å, {Delta}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 {Delta}lbc, the slope of the variance at the end of the transition region determines {Delta}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 = 74–289. Despite their different lengths, we obtained fairly narrow distributions for the three fitting parameters (see Table 1) with average values <{Delta}lbi> = 0.518 ± 0.018 Å, <{Delta}lbc> = 0.988 ± 0.020 Å, and <B> = 0.626 ± 0.037.


Figure 4
View larger version (22K):
[in this window]
[in a new window]
 
FIGURE 4  Fits using our three-state model (black) result in good agreements with the variances (green) of molecules of very different lengths. For each molecule, N is fixed to the value obtained from fitting the force-extension curves using the FJC model, while {Delta}lbi, {Delta}lbc, and B are left free to vary. The values obtained for the fitting parameters are very similar even though the molecules vary in length by more than a factor of 3.

 

View this table:
[in this window]
[in a new window]
 
TABLE 1  Results from fits with three-state model

 
Note that <{Delta}lbc> = 0.988 Å is significantly smaller than the value from ab initio calculations, {Delta}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 (32Go).


    CONCLUSION
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS AND DISCUSSION
 CONCLUSION
 APPENDIX A: MACROSCOPIC...
 APPENDIX B: VARIANCE FOR...
 ACKNOWLEDGEMENTS
 REFERENCES
 
We have demonstrated for the first time the use of variance analysis to capture dynamic molecular conformations of a single molecule mechanically stretched by an AFM. The fluctuations that we have detected occur at the sub-Angstrom level and are invisible to any other type of spectroscopy, including NMR and x-ray crystallography, as those techniques require large ensembles of molecules that are all in the same configuration. It is therefore difficult to observe conformations of a molecule that are transient, rapidly fluctuating, and only occupied by a fraction of the ensemble. By contrast, our single molecule force spectroscopy analysis of the variance in the length of a stretched dextran molecule readily captures transient conformations on the sub-Angstrom scale. Furthermore, we developed a rigorous test for the applicability of our variance analysis based on establishing that the molecule under study was in thermodynamic equilibrium.

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 (2Go,23Go).

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 (1Go) as well as more complex molecules like DNA (33Go), RNA (2Go), and proteins (26Go) that undergo discrete conformational changes under a stretching force.


    APPENDIX A: MACROSCOPIC ELASTICITY OF THE MOLECULE
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS AND DISCUSSION
 CONCLUSION
 APPENDIX A: MACROSCOPIC...
 APPENDIX B: VARIANCE FOR...
 ACKNOWLEDGEMENTS
 REFERENCES
 
In Eq. 2 we assume km is constant over the fluctuation range at each point in the force-extension curve, which is correct to first order. The exact expression for km can be written as a Taylor expansion, as

Formula 7(7)
where <...> indicates an average. Substituting Formula 7 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
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS AND DISCUSSION
 CONCLUSION
 APPENDIX A: MACROSCOPIC...
 APPENDIX B: VARIANCE FOR...
 ACKNOWLEDGEMENTS
 REFERENCES
 
Our three-state model introduces an intermediate between the chair and the boat conformation. The probability density function is given by the trinomial distribution

Formula 8(8)
with a variance, in the Gibbs ensemble, of

Formula 9(9)
Here, ni and nc are the number of monomers in the intermediate and chair conformation, pi and pc are the probabilities that a given ring is in the intermediate or chair, and {Delta}lbi and {Delta}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. Its magnitude has to be ≤1.
  2. It has to be 0 outside the transition region.
  3. Its maximum has to occur at N{Delta}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
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS AND DISCUSSION
 CONCLUSION
 APPENDIX A: MACROSCOPIC...
 APPENDIX B: VARIANCE FOR...
 ACKNOWLEDGEMENTS
 REFERENCES
 
K.A.W. thanks Ali Kinkhabwala for useful discussions and editing.

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
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS AND DISCUSSION
 CONCLUSION
 APPENDIX A: MACROSCOPIC...
 APPENDIX B: VARIANCE FOR...
 ACKNOWLEDGEMENTS
 REFERENCES
 
1. Marszalek, P. E., H. Li, and J. M. Fernandez. 2001. Fingerprinting polysaccharides with single-molecule atomic force microscopy. Nat. Biotechnol. 19:258–262.[CrossRef][Medline]

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:733–737.[Abstract/Free Full Text]

3. Bustamante, C., Z. Bryant, and S. B. Smith. 2003. Ten years of tension: single-molecule DNA mechanics. Nature. 421:423–427.[CrossRef][Medline]

4. Fernandez, J. M., and H. Li. 2004. Force-clamp spectroscopy monitors the folding trajectory of a single protein. Science. 303:1674–1678.[Abstract/Free Full Text]

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:8618–8626.[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:3171–3176.[Abstract/Free Full Text]

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.[Abstract/Free Full Text]

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:549–560.

9. Sigworth, F. J. 1980. The variance of sodium current fluctuations at the node of Ranvier. J. Physiol. (Lond.). 307:97–129.[Abstract/Free Full Text]

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:7927–7931.[Abstract/Free Full Text]

11. Liu, Y. Z., S. H. Leuba, and S. M. Lindsay. 1999. Relationship between stiffness and force in single molecule pulling experiments. Langmuir. 15:8547–8548.[CrossRef]

12. Humphris, A. D. L., J. Tamayo, and M. J. Miles. 2000. Active quality factor control in liquids for force spectroscopy. Langmuir. 16:7891–7894.[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:9299–9303.[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:2505–2514.[Abstract/Free Full Text]

16. Titantah, J. T., C. Pierleoni, and J. P. Ryckaert. 1999. Different statistical mechanical ensembles for a stretched polymer. Phys. Rev. E. 60:7010–7021.[CrossRef]

17. Gerland, U., R. Bundschuh, and T. Hwa. 2001. Force-induced denaturation of RNA. Biophys. J. 81:1324–1332.[Abstract/Free Full Text]

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:13544–13548.[Abstract/Free Full Text]

19. Gerland, U., R. Bundschuh, and T. Hwa. 2003. Mechanically probing the folding pathway of single RNA molecules. Biophys. J. 84:2831–2840.[Abstract/Free Full Text]

20. Keller, D., D. Swigon, and C. Bustamante. 2003. Relating single-molecule measurements to thermodynamics. Biophys. J. 84:733–738.[Abstract/Free Full Text]

21. Tinoco, I. 2004. Force as a useful variable in reactions: unfolding RNA. Annu. Rev. Biophys. Biomol. 33:363–385.[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:1832–1835.[Abstract/Free Full Text]

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:539–544.[Free Full Text]

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:661–664.[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:1109–1112.[Abstract/Free Full Text]

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:7299–7304.[Abstract/Free Full Text]

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:895–901.[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:4278–4283.[Abstract/Free Full Text]

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:1295–1297.[Abstract/Free Full Text]

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:1456–1465.[Abstract/Free Full Text]

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:795–799.[Abstract]

34. Braun, O., and U. Seifert. 2005. Force spectroscopy of single multidomain biopolymers: a master equation approach. Eur. Phys. J. E. 18:1–13.[CrossRef][Medline]




This article has been cited by other articles:


Home page
Biophys. JHome page
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]


Home page
Biophys. JHome page
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]


Home page
ScienceHome page
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]


Home page
Biophys. JHome page
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]


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


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