| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Department of Bioengineering, Rice University, Houston, Texas
Correspondence: Address reprint requests to Dr. Robert M. Raphael, Rice University, Dept. of Bioengineering, MS-142, PO Box 1892, Houston, TX 77251-1892.E-mail: rraphael{at}rice.edu.
| ABSTRACT |
|---|
|
|
|---|
| INTRODUCTION |
|---|
|
|
|---|
Though many basic membrane properties can be characterized by pulling tethers from aspirated bilayers, tether experiments are often conducted on unaspirated membranes. In fact, tethers were first observed when red blood cells and neutrophils were exposed to fluid flow (12
). Following the development of more recent tether formation methods, such as optical tweezers and fiber-deflection techniques, tether studies have helped characterize the membrane-cytoskeleton adhesion energy (13
16
). Protrusive tether-like structures have been observed in several experiments which examine the interaction between cytoskeletal components and lipid membranes. As tubulin polymerizes inside of a vesicle, the biopolymers induce an axial tubular membrane protrusion with a geometry similar to a membrane tether (17
,18
). In addition, microtubule motor proteins can form tethers from lipid vesicles in vitro (19
,20
).
To understand the forces these molecules exert on lipid membranes, optical tweezers have been used to measure the force required to form tethers from both unaspirated lipid bilayers as well as from reconstituted Golgi and endoplasmic reticulum networks (20
,21
). Distinct from the aspirated case, a lateral tension is not imposed a priori on the membrane. As a result, tether studies can be conducted on membranes with tensions much lower than those induced by an aspiration pipette and more similar to those of many cellular membranes.
For a membrane under low tension, the thermally induced fluctuations of the vesicle surface are more pronounced than those of a membrane under moderate tension (22
24
). Since the bending energy is on the order of kT, the fluctuations manifest as out-of-plane undulations of the vesicle surface (25
). As a result, the total surface area of the membrane is greater than the apparent or measured surface area. This "excess" surface area influences the response of a lipid membrane to a mechanical deformation. Expansion of the membrane area, for example, smoothes the surface undulations (22
). Because flattening of the surface reduces the entropy, a force is needed to increase the apparent surface area. Thus, smoothing the surface undulations leads to an increase in membrane tension. The relationship between applied tension and area expansion has been demonstrated by micropipette aspiration experiments. When a vesicle is first aspirated into a micropipette, an increase in the observable surface area is accompanied by a small increase in tension. Upon further aspiration, the tension increases linearly with the area dilation (22
,26
,27
). The observed behavior correlates with theoretical predictions in which the first phase is dominated by smoothing of the surface undulations and the second by expansion of the area per lipid molecule (28
). According to this theory, the mean-square amplitude of these fluctuations depends upon the effective vesicle tension.
The effective tension of a vesicle can be manipulated by altering the relative concentration of internal and external impermeable solutes (11
). Because the osmotic pressure is typically greater than the pressures created by mechanical manipulation, the volume of the vesicle is assumed to remain constant during tether formation (29
). Under the constraint of constant volume, excess surface area is required to form a tether. To provide this area for micropipette aspiration tether experiments, vesicles are usually transferred to a hypertonic medium before tether extraction (2
,8
,30
). Theoretical studies have considered the shape and stability of tethers formed from both axially strained vesicles and from lipid reservoirs (1
,29
,31
,32
). In another model, the unbinding of a strongly adhered, tense vesicle from the surface provides area for tether formation (33
). Approaches applied for cellular membranes and unaspirated vesicles model the membrane from which the tether is formed as a disc with a far-field tension
(34
). By approximating the tether geometry as a simple cylinder, the tether force is expressed in terms of
and the tether length Lt: F = 2
(2kc
)1/2 + 4
2krLt/Ao where Ao is the area of the vesicle, and kc and kr are the local and nonlocal bending moduli, respectively (3
,34
). The tether radius Rt is determined by the membrane tension
, where Rt = (kc/2
)1/2. Thus, if the tension does not change as a tether is lengthened, then the tether radius should remain constant. Recent experiments, however, have demonstrated a decrease in tether radius with an increase in tether length, suggesting membrane tension increases as a tether is elongated (20
).
Here we develop a model of tether formation from an unaspirated lipid vesicle which provides a quantitative relation between membrane tension and tether length. In this model, we consider how entropically driven tensions may influence tether conformation; therefore, this model allows characterization of tethers formed from membranes that are not elastically tense. Following the approach first presented by Bo
i
et al. (4
), the vesicle shape is parameterized by a simple geometrical model, and only significant contributions to the thermodynamic potential are considered (4
). Energy minimization provides an analytical relation between the equilibrium pulling force and the tether radius. Since tether radii often cannot be measured directly, an analytical relation between tether force and length is obtained by excluding the contributions of the nonlocal bending energy. Given many analyses of tether experiments do not consider nonlocal bending, we have analyzed the appropriate bounds of this approximation.
| MODEL |
|---|
|
|
|---|
|
![]() | (1) |
![]() | (2) |
The first term of Eq. 2 is the local bending energy (7
,35
). The local bending modulus kc represents both monolayers' resistance to bending;
is the sum of the two principal curvatures, and co is the initial curvature of the bilayer. Because the bending energy of the tether is at least several orders of magnitude greater than the bending energy of other vesicle regions, only the bending energy of the tether is included in the energy expression (4
).
The second term of Eq. 2, the nonlocal bending energy, accounts for the relative expansion and compression of each monolayer. The constant kr is the nonlocal bending modulus and h is the distance between the neutral surfaces of each monolayer (4
,6
,7
,36
). Here
A represents the relative area difference between the membrane leaflets;
Ao is included for completeness to account for the initial area difference between the monolayers. Using the expression
Eq. 2 may be rearranged so that the contributions of
Ao and co are grouped together into the reduced area difference
Ao* and a constant wo (5
,37
):
![]() | (3) |
![]() |
![]() |
For this analysis,
Ao* is taken to be zero; the consequences of a nonzero
Ao* are considered in the Discussion.
A arising from tether formation is estimated from the difference between the areas of the inner tether leaflet, 2
Lt(Rt h/2), and the outer leaflet, 2
Lt(Rt + h/2) (4
).
The third term of Eq. 2 is the apparent area expansion energy, where
is the effective membrane tension (25
,38
). The percent increase in measured area due to the flattening of the thermal fluctuations
is (Am Am,o)/Ao, where Am is the measured area, Am,o is the initial measured area, and Ao is the total vesicle area (22
). Because expansion of molecular area is not considered, the total area is held constant. Using a planar wave approximation for thermal fluctuations (22
,25
), the percent change in measured area is estimated as:
![]() | (4) |
o is the initial effective membrane tension. This approximation is valid for
<<
2kc/a2; the constant a is the molecular spacing (a
5 Å). For the geometry of this model, the change in measured area is
![]() | (5) |
For a fixed volume (V = Vo; Rv = Rv,o),
reduces to 2
RtLt/Ao. When the geometric and energetic approximations are applied, Eq. 1 becomes:
![]() | (6) |
Ao*/2
h. The constant xo is wo
kT/8. By selecting a characteristic length scale as the radius of a sphere of area Ao and energy scale as the bending energy of a sphere (8
kc), Eq. 6 can be written in nondimensional form in terms of the variables lt and rt (lt = Lt/R; rt = Rt/R):
![]() | (7) |
Here, the nondimensional force
is FR/(8
kc), and the initial tension
o is
oR2/kc. At equilibrium, the partial derivatives of Eq. 7 are zero [
=
(rt, lt) = 0]. This provides two equations which can be solved to determine the minimum energy tether conformations:
![]() | (8) |
![]() | (9) |
The tether force equation, obtained by combining the dimensional forms of Eqs. 8 and 9, mirrors the equations derived by previous energy variational and force balance methods (3
,4
):
![]() | (10) |
The constant Lt* contributes to the total force. The difference between this model and those for aspirated vesicles, however, lies in the relation between Lr and Rt. This difference can be examined by solving Eq. 9 for Lt:
![]() | (11) |
2kc and b = kc/(2
o + 2
2kc/Ao). Note that the relationship between Lt and Rt is independent of the nonlocal bending stiffness kr. | RESULTS AND DISCUSSION |
|---|
|
|
|---|
5 pN. The membrane tension for a specified tether length depends upon the value of the initial tension (Fig. 2 B). If the initial vesicle tension is changed by an order of magnitude, the change in the force versus tether length curve is within experimental resolution (Fig. 2 A). These trends are significant as they motivate a method to determine the initial tension of a vesicle and indicate that even relatively small tensions influence tether behavior.
|
|
kc/Rt and substituted into Eq. 11 to provide a relationship between the tether force and length:
![]() | (12) |
and
A fit of Eq. 12 to the F versus Lt tether behavior for the SLO DOPC vesicle (20
o obtained from parameters c and d are 1.1 x 108 m2 and 7.2 x 105 mN/m, respectively. These values are similar to those calculated from a fit of the Rt versus Lt data to Eq. 11: 1.2 x 108 m2 and 8.2 x 105 mN/m (Fig. 3 A).
Since Eq. 12 was obtained by excluding the contributions of the nonlocal bending energy, we have analyzed the appropriateness of this assumption. The extent to which the energy affects equilibrium tether conformation depends upon the magnitude of the nonlocal bending coefficient kr, the total vesicle area and the initial tension difference between the two membrane leaflets. The percent contribution of the nonlocal bending energy to the tether force is plotted for a range of tether lengths and ratios of kr/kc (Fig. 4). For a typical ratio of kr/kc (3.0 (2
,8
)), the percent contribution to the force for a tether in which lt = 4 is <5%. The reported bending stiffness (5.5 x 1020 J) used to obtain the fits in Fig. 3 was obtained from the slope of the experimentally measured plot F vs. 2
/Rt (20
). Provided the initial area difference between the membrane leaflets is negligible, nonlocal bending increases the slope of the plot of F vs. 1/4r (Fig. 4 B). Although neglecting the nonlocal bending energy may lead to an overestimation of the bending stiffness kc, for a typical ratio of kr/kc of 3, the change in slope is <5% (
4%).
|
Ao*/2
h. For a stable shape of the area difference elasticity model,
Ao* can be related to the effective curvature, Co*:
Ao* = Co*Aohkc/kr, where the effective curvature is given by Co* = Co + kr/kc(
Ao/Aoh) (41
kr. Including this expression for Lt* into Eq. 10 shifts the mechanical tether force by 2
Co*kc (and F by co*/4, where co* = RCo*). Note that the effective spontaneous curvature depends not only upon the vesicle composition, formation, and history but also upon its aqueous bathing environment. Values of co* for vesicles suspended in an asymmetric glucose solution, for example, range in value from
2 to 6 (42
The results obtained from this analysis can help to probe the microarchitecture of bilayer membranes. Experimental evidence indicates some vesicles contain surface microvesiculations that are not fully resolvable optically (41
). The propensity to form these vesiculations may in part be due to the lipid composition of the vesicle. For example, under thermal expansion, the mean-square fluctuation amplitude of DMPC vesicles increases, whereas for SOPC vesicles microvesiculations form (41
). Many cellular membranes also exhibit a rich spectrum of microarchitectural contours. Since a typical tether (Lt = 10 µm and Rt = 100 nm) formed from a vesicle in which R = 10 µm encompasses
0.5% of the total vesicle area, microvesiculations could conceivably serve as an area reservoir for tether formation. To estimate the resulting tether behavior, we consider the simplest theoretical case in which area extraction from this reservoir occurs under constant membrane tension. For this case, the tether force equation provided in the introduction is applicable. As the tether is lengthened, the radius would remain constant and the force would increase slightly with the slope of 4
2kr/Ao. This might continue until a critical area Ac is extracted from the surface vesiculations. This behavior would be in contrast to the force versus length trends predicted for a tether formed by extracting area from the thermal undulations.
The tether force may reflect the microarchitectural properties in ways more subtle through the constant Lt*. Since the differential area between the membrane leaflets is believed to contribute in part to the development of surface microvesiculations, those membranes which exhibit significant budding may also have a nontrivial Lt* (43
). Consequently, care should be taken when interpreting the response of tethers formed from vesicles with marked surface microvesiculations.
According to the analysis, as area is drawn from the thermal fluctuations to form a tether, the effective vesicle tension increases (Fig. 2 B). Due to the constant area constraint, the tension arising from expansion of molecular area is not included. For the small tensions considered in this model, the energy due to area expansion is much less than the thermal undulation term. As the tension increases, though, the area expansivity energy becomes larger than the undulation energy. The transition tension
Trans can be estimated by determining the tension for which the entropic and elastic stretching energies are equal (22
,25
):
Trans
kT/kc. For a pure DOPC vesicle (
= 0.265 N/m and kc = 8.5 x 1020 J (28
)),
Trans is 0.84 mN/m. Because the initial tension influences the relation between tether force and tether length (Fig. 2 B), the exact tether length at which this transition tension is reached depends upon the initial effective tension.
The thermally driven surface undulations of the vesicle are approximated using a simple planar approach. This approach is valid for a spherical surface provided the amplitudes of the undulations are much smaller than the radius of the vesicle. Alternate models which employ either spherical harmonics or a coarse-grained approach may provide a more accurate description of the surface undulations, especially those for vesicles with significant excess area (23
,24
). Of interest, theoretical studies have suggested if the vesicle has enough excess area, the effective tension may be <0 (24
,44
). This would have important implications for the physical actuality of the effective tension, a term often added as a Lagrange multiplier to enforce a constant area constraint (44
). A "negative tension" might influence for example the vesicle microarchitecture. We also note that the presence of the tether will alter the overall characteristics of the surface undulations. The simple planar approximation employed, however, provides a basic description of the undulations of the vesicle surface as well as analytical relations for the dependence of the tether geometry on the tether force.
Previous work has considered how fluctuation-induced differential area dilation suppresses the mean-square amplitude of the bending undulations (11
,45
). Tether formation, however, imposes an additional area dilation field on the bilayer that could further suppress the surface undulations. Tether studies, especially at longer lengths, may provide a means to examine how differential area dilation affects the amplitude of surface undulations.
In summary, we have developed a model that provides a quantitative relation between force and tether conformation for a tether formed from an unaspirated lipid vesicle. According to the model, the thermal undulations of the vesicle surface can provide the excess area required for tether formation. The model captures both the increase in tether force with length and the decrease in tether radius observed experimentally (Fig. 3). Since the initial tension influences the extent of the undulations, the analysis indicates that the tether force is sensitive to these extraordinarily small tensions. Because these tensions are on the order of the resting tensions of cellular membranes, analyses of tether experiments can further our understanding of the role that thermally driven tensions play in biological processes.
| ACKNOWLEDGEMENTS |
|---|
|
|
|---|
This work was supported by the National Science Foundation under grant No. 0114624.
Submitted on June 9, 2005; accepted for publication April 5, 2006.
| REFERENCES |
|---|
|
|
|---|
2. Waugh, R. E., J. Song, S. Svetina, and B. Zeks. 1992. Local and nonlocal curvature elasticity in bilayer membranes by tether formation from lecithin vesicles. Biophys. J. 61:974982.
3. Waugh, R. E., and R. M. Hochmuth. 1987. Mechanical equilibrium of thick, hollow, liquid membrane cylinders. Biophys. J. 52:391400.
4. Bo
i
, B., S. Svetina, B.
ek
, and R. E. Waugh. 1992. Role of lamellar membrane structure in tether formation from bilayer vesicles. Biophys. J. 61:963973.
5. Miao, L., U. Seifert, M. Wortis, and H. G. Dobereiner. 1994. Budding transitions of fluid-bilayer vesicles: The effect of area-difference elasticity. Phys. Rev. E Stat. Phys. Plasmas Fluids Relat. Interdiscip. Topics. 49:53895407.[Medline]
6. Helfrich, W. 1974. Blocked lipid exchange in bilayers and its possible influence on the shape of vesicles. Z. Naturforsch. [C]. 29C:510515.
7. Evans, E. A. 1974. Bending resistance and chemically induced moments in membrane bilayers. Biophys. J. 14:923931.
8. Raphael, R. M., and R. E. Waugh. 1996. Accelerated interleaflet transport of phosphatidylcholine molecules in membranes under deformation. Biophys. J. 71:13741388.
9. Svetina, S., B. Zeks, R. E. Waugh, and R. M. Raphael. 1998. Theoretical analysis of the effect of the transbilayer movement of phospholipid molecules on the dynamic behavior of a microtube pulled out of an aspirated vesicle. Eur. Biophys. J. 27:197209.[CrossRef][Medline]
10. Raphael, R. M., R. E. Waugh, S. Svetina, and B. Zeks. 2001. Fractional occurrence of defects in membranes and mechanically driven interleaflet phospholipid transport. Phys. Rev. E Stat. Nonlin. Soft Matter Phys. 64:051913.[Medline]
11. Yeung, A., and E. Evans. 1995. Unexpected dynamics in shape fluctuations of bilayer vesicles. J. Phys. II. 5:15011523.[CrossRef]
12. Hochmuth, R. M., N. Mohandas, E. E. Spaeth, J. R. Williamson, P. L. Blackshear Jr., and D. W. Johnson. 1972. Surface adhesion, deformation and detachment at low shear of red cells and white cells. Trans. Am. Soc. Artif. Intern. Organs. 18:325334.[Medline]
13. Sheetz, M. P. 2001. Cell control by membrane-cytoskeleton adhesion. Nat. Rev. Mol. Cell Biol. 2:392396.[CrossRef][Medline]
14. Hochmuth, R. M., and W. D. Marcus. 2002. Membrane tethers formed from blood cells with available area and determination of their adhesion energy. Biophys. J. 82:29642969.
15. Waugh, R. E., and R. G. Bauserman. 1995. Physical measurements of bilayer-skeletal separation forces. Ann. Biomed. Eng. 23:308321.[Medline]
16. Hwang, W. C., and R. E. Waugh. 1997. Energy of dissociation of lipid bilayer from the membrane skeleton of red blood cells. Biophys. J. 72:26692678.
17. Hotani, H., and H. Miyamoto. 1990. Dynamic features of microtubules as visualized by dark-field microscopy. Adv. Biophys. 26:135156.[CrossRef][Medline]
18. Fygenson, D. K., J. F. Marko, and A. Libchaber. 1997. Mechanics of microtubule-based membrane extension. Phys. Rev. Lett. 79:44974500.[CrossRef]
19. Roux, A., G. Cappello, J. Cartaud, J. Prost, B. Goud, and P. Bassereau. 2002. A minimal system allowing tubulation with molecular motors pulling on giant liposomes. Proc. Natl. Acad. Sci. USA. 99:53945399.
20. Koster, G., M. VanDuijn, B. Hofs, and M. Dogterom. 2003. Membrane tube formation from giant vesicles by dynamic association of motor proteins. Proc. Natl. Acad. Sci. USA. 100:1558315588.
21. Upadhyaya, A., J. R. Chabot, A. Andreeva, A. Samadani, and A. van Oudenaarden. 2003. Probing polymerization forces by using actin-propelled lipid vesicles. Proc. Natl. Acad. Sci. USA. 100:45214526.
22. Evans, E., and W. Rawicz. 1990. Entropy-driven tension and bending elasticity in condensed-fluid membranes. Phys. Rev. Lett. 64:20942097.[CrossRef][Medline]
23. Milner, S. T., and S. A. Safran. 1987. Dynamical fluctuations of droplet microemulsions and vesicles. Phys. Rev. A. 36:43714379.[CrossRef][Medline]
24. Fournier, J. B., A. Ajdari, and L. Peliti. 2001. Effective-area elasticity and tension of micromanipulated membranes. Phys. Rev. Lett. 86:49704973.[CrossRef][Medline]
25. Helfrich, W. S. R.M. 1984. Undulations, steric interactions and cohesion of fluid membranes. Il Nuovo Cimento. 3D:137151.[CrossRef]
26. Discher, B. M., Y. Y. Won, D. S. Ege, J. C. Lee, F. S. Bates, D. E. Discher, and D. A. Hammer. 1999. Polymersomes: tough vesicles made from diblock copolymers. Science. 284:11431146.
27. Zhou, Y., and R. M. Raphael. 2005. Effect of salicylate on the elasticity, bending stiffness and strength of SOPC membranes. Biophys. J. 89:17891801.
28. Rawicz, W., K. C. Olbrich, T. McIntosh, D. Needham, and E. Evans. 2000. Effect of chain length and unsaturation on elasticity of lipid bilayers. Biophys. J. 79:328339.
29. Heinrich, V., B. Bozic, S. Svetina, and B. Zeks. 1999. Vesicle deformation by an axial load: from elongated shapes to tethered vesicles. Biophys. J. 76:20562071.
30. Heinrich, V., and R. E. Waugh. 1996. A piconewton force transducer and its application to measurement of the bending stiffness of phospholipid membranes. Ann. Biomed. Eng. 24:595605.[Medline]
31. Derenyi, I., F. Julicher, and J. Prost. 2002. Formation and interaction of membrane tubes. Phys. Rev. Lett. 88:238101.[CrossRef][Medline]
32. Powers, T. R., G. Huber, and R. E. Goldstein. 2002. Fluid-membrane tethers: minimal surfaces and elastic boundary layers. Phys. Rev. E. 65:041901.[CrossRef]
33. Smith, A., E. Sackmann, and U. Seifert. 2004. Pulling tethers from adhered vesicles. Phys. Rev. Lett. 92:208101.[Medline]
34. Hochmuth, R. M., J. Y. Shao, J. Dai, and M. P. Sheetz. 1996. Deformation and flow of membrane into tethers extracted from neuronal growth cones. Biophys. J. 70:358369.
35. Helfrich, W. 1973. Elastic properties of lipid bilayers: theory and possible experiments. Z. Naturforsch. [C]. 28:693703.[Medline]
36. Svetina, S., M. Brumen, and B. Zeks. 1985. Lipid bilayer elasticity and the bilayer couple interpretation of red-cell shape transformations and lysis. Stud. Biophys. 110:177184.
37. Döbereiner, H. G., E. Evans, M. Kraus, U. Seifert, and M. Wortis. 1997. Mapping vesicle shapes into the phase diagram: a comparison of experiment and theory. Phys. Rev. E Stat. Nonlin. Soft Matter Phys. 55:44584474.
38. Helfrich, W. 1975. Out-of-plane fluctuations of lipid bilayers. Z. Naturforsch. [C]. 30:841842.[Medline]
39. Palmer, M., R. Harris, C. Freytag, M. Kehoe, J. Tranum-Jensen, and S. Bhakdi. 1998. Assembly mechanism of the oligomeric streptolysin O pore: the early membrane lesion is lined by a free edge of the lipid membrane and is extended gradually during oligomerization. EMBO J. 17:15981605.[CrossRef][Medline]
40. Upadhyaya, A., and M. P. Sheetz. 2004. Tension in tubulovesicular networks of Golgi and endoplasmic reticulum membranes. Biophys. J. 86:29232928.
41. Häckl, W., U. Seifert, and E. Sackmann. 1997. Effects of fully and partially solubilized amphiphiles on bilayer bending stiffness and temperature dependence of the effective tension of giant vesicles. J. Phys. II. 7:11411157.[CrossRef]
42. Döbereiner, H. G., O. Selchow, and R. Lipowsky. 1999. Spontaneous curvature of fluid vesicles induced by trans-bilayer sugar asymmetry. Eur. Biophys. J. 28:174178.[CrossRef]
43. Iglic, A., and H. Hagerstrand. 1999. Amphiphile-induced spherical microexovesicle corresponds to an extreme local area difference between two monolayers of the membrane bilayer. Med. Biol. Eng. Comput. 37:125129.[Medline]
44. Seifert, U. 1995. The concept of effective tension for fluctuating vesicles. Z. Phys. B: Condens. Matter. 97:299309.[CrossRef]
45. Seifert, U. 1993. Curvature-induced lateral phase segregation in two-component vesicles. Phys. Rev. Lett. 70:13351338.[CrossRef][Medline]
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |