Elasticity of cell membranes

Last updated

A cell membrane defines a boundary between a cell and its environment. The primary constituent of a membrane is a phospholipid bilayer that forms in a water-based environment due to the hydrophilic nature of the lipid head and the hydrophobic nature of the two tails. In addition there are other lipids and proteins in the membrane, the latter typically in the form of isolated rafts.

Contents

Of the numerous models that have been developed to describe the deformation of cell membranes, a widely accepted model is the fluid mosaic model proposed by Singer and Nicolson in 1972. [1] In this model, the cell membrane surface is modeled as a two-dimensional fluid-like lipid bilayer where the lipid molecules can move freely. The proteins are partially or fully embedded in the lipid bilayer. Fully embedded proteins are called integral membrane proteins because they traverse the entire thickness of the lipid bilayer. These communicate information and matter between the interior and the exterior of the cell. Proteins that are only partially embedded in the bilayer are called peripheral membrane proteins. The membrane skeleton is a network of proteins below the bilayer that links with the proteins in the lipid membrane.

Elasticity of closed lipid vesicles

The simplest component of a membrane is the lipid bilayer which has a thickness that is much smaller than the length scale of the cell. Therefore, the lipid bilayer can be represented by a two-dimensional mathematical surface. In 1973, based on similarities between lipid bilayers and nematic liquid crystals, Helfrich [2] proposed the following expression for the curvature energy per unit area of the closed lipid bilayer

 

 

 

 

(1)

where are bending rigidities, is the spontaneous curvature of the membrane, and and are the mean and Gaussian curvature of the membrane surface, respectively.

The free energy of a closed bilayer under the osmotic pressure (the outer pressure minus the inner one) as:

 

 

 

 

(2)

where dA and dV are the area element of the membrane and the volume element enclosed by the closed bilayer, respectively, and λ is the Lagrange multiplier for area inextensibility of the membrane, which has the same dimension as surface tension. By taking the first order variation of above free energy, Ou-Yang and Helfrich [3] derived an equation to describe the equilibrium shape of the bilayer as:

 

 

 

 

(3)

They also obtained that the threshold pressure for the instability of a spherical bilayer was

 

 

 

 

(4)

where being the radius of the spherical bilayer.

Using the shape equation (3) of closed vesicles, Ou-Yang predicted that there was a lipid torus with the ratio of two generated radii being exactly . [4] His prediction was soon confirmed by the experiment [5] Additionally, researchers obtained an analytical solution [6] to (3) which explained the classical problem, the biconcave discoidal shape of normal red blood cells. In the last decades, the Helfrich model has been extensively used in computer simulations of vesicles, red blood cells and related systems. From a numerical point-of-view bending forces stemming from the Helfrich model are very difficult to compute as they require the numerical evaluation of fourth-order derivatives and, accordingly, a large variety of numerical methods have been proposed for this task. [7]

Elasticity of open lipid membranes

The opening-up process of lipid bilayers by talin was observed by Saitoh et al. [8] arose the interest of studying the equilibrium shape equation and boundary conditions of lipid bilayers with free exposed edges. Capovilla et al., [9] Tu and Ou-Yang [10] carefully studied this problem. The free energy of a lipid membrane with an edge is written as

 

 

 

 

(5)

where and represent the arclength element and the line tension of the edge, respectively. This line tension is a function of dimension and distribution of molecules comprising the edge, and their interaction strength and range. [11] The first order variation gives the shape equation and boundary conditions of the lipid membrane: [12]

 

 

 

 

(6)

 

 

 

 

(7)

 

 

 

 

(8)

 

 

 

 

(9)

where , , and are normal curvature, geodesic curvature, and geodesic torsion of the boundary curve, respectively. is the unit vector perpendicular to the tangent vector of the curve and the surface normal vector of the membrane.

Elasticity of cell membranes

A cell membrane is simplified as lipid bilayer plus membrane skeleton. The skeleton is a cross-linking protein network and joints to the bilayer at some points. Assume that each proteins in the membrane skeleton have similar length which is much smaller than the whole size of the cell membrane, and that the membrane is locally 2-dimensional uniform and homogenous. Thus the free energy density can be expressed as the invariant form of , , and :

 

 

 

 

(10)

where is the in-plane strain of the membrane skeleton. Under the assumption of small deformations, and invariant between and , (10) can be expanded up to second order terms as:

 

 

 

 

(11)

where and are two elastic constants. In fact, the first two terms in (11) are the bending energy of the cell membrane which contributes mainly from the lipid bilayer. The last two terms come from the entropic elasticity of the membrane skeleton.

Related Research Articles

<span class="mw-page-title-main">Biological membrane</span> Enclosing or separating membrane in organisms acting as selective semi-permeable barrier

A biological membrane, biomembrane or cell membrane is a selectively permeable membrane that separates the interior of a cell from the external environment or creates intracellular compartments by serving as a boundary between one part of the cell and another. Biological membranes, in the form of eukaryotic cell membranes, consist of a phospholipid bilayer with embedded, integral and peripheral proteins used in communication and transportation of chemicals and ions. The bulk of lipids in a cell membrane provides a fluid matrix for proteins to rotate and laterally diffuse for physiological functioning. Proteins are adapted to high membrane fluidity environment of the lipid bilayer with the presence of an annular lipid shell, consisting of lipid molecules bound tightly to the surface of integral membrane proteins. The cell membranes are different from the isolating tissues formed by layers of cells, such as mucous membranes, basement membranes, and serous membranes.

<span class="mw-page-title-main">Vesicle (biology and chemistry)</span> Any small, fluid-filled, spherical organelle enclosed by a membrane

In cell biology, a vesicle is a structure within or outside a cell, consisting of liquid or cytoplasm enclosed by a lipid bilayer. Vesicles form naturally during the processes of secretion (exocytosis), uptake (endocytosis), and the transport of materials within the plasma membrane. Alternatively, they may be prepared artificially, in which case they are called liposomes. If there is only one phospholipid bilayer, the vesicles are called unilamellar liposomes; otherwise they are called multilamellar liposomes. The membrane enclosing the vesicle is also a lamellar phase, similar to that of the plasma membrane, and intracellular vesicles can fuse with the plasma membrane to release their contents outside the cell. Vesicles can also fuse with other organelles within the cell. A vesicle released from the cell is known as an extracellular vesicle.

<span class="mw-page-title-main">Lipid bilayer</span> Membrane of two layers of lipid molecules

The lipid bilayer is a thin polar membrane made of two layers of lipid molecules. These membranes are flat sheets that form a continuous barrier around all cells. The cell membranes of almost all organisms and many viruses are made of a lipid bilayer, as are the nuclear membrane surrounding the cell nucleus, and membranes of the membrane-bound organelles in the cell. The lipid bilayer is the barrier that keeps ions, proteins and other molecules where they are needed and prevents them from diffusing into areas where they should not be. Lipid bilayers are ideally suited to this role, even though they are only a few nanometers in width, because they are impermeable to most water-soluble (hydrophilic) molecules. Bilayers are particularly impermeable to ions, which allows cells to regulate salt concentrations and pH by transporting ions across their membranes using proteins called ion pumps.

<span class="mw-page-title-main">Thylakoid</span> Membrane enclosed compartments in chloroplasts and cyanobacteria

Thylakoids are membrane-bound compartments inside chloroplasts and cyanobacteria. They are the site of the light-dependent reactions of photosynthesis. Thylakoids consist of a thylakoid membrane surrounding a thylakoid lumen. Chloroplast thylakoids frequently form stacks of disks referred to as grana. Grana are connected by intergranal or stromal thylakoids, which join granum stacks together as a single functional compartment.

<span class="mw-page-title-main">Fluid mosaic model</span> Describe the fluid mosaic model of plasma membrane

The fluid mosaic model explains various characteristics regarding the structure of functional cell membranes. According to this biological model, there is a lipid bilayer in which protein molecules are embedded. The phospholipid bilayer gives fluidity and elasticity to the membrane. Small amounts of carbohydrates are also found in the cell membrane. The biological model, which was devised by Seymour Jonathan Singer and Garth L. Nicolson in 1972, describes the cell membrane as a two-dimensional liquid that restricts the lateral diffusion of membrane components. Such domains are defined by the existence of regions within the membrane with special lipid and protein cocoon that promote the formation of lipid rafts or protein and glycoprotein complexes. Another way to define membrane domains is the association of the lipid membrane with the cytoskeleton filaments and the extracellular matrix through membrane proteins. The current model describes important features relevant to many cellular processes, including: cell-cell signaling, apoptosis, cell division, membrane budding, and cell fusion. The fluid mosaic model is the most acceptable model of the plasma membrane. In this definition of the cell membrane, its main function is to act as a barrier between the contents inside the cell and the extracellular environment.

Phosphatidic acids are anionic phospholipids important to cell signaling and direct activation of lipid-gated ion channels. Hydrolysis of phosphatidic acid gives rise to one molecule each of glycerol and phosphoric acid and two molecules of fatty acids. They constitute about 0.25% of phospholipids in the bilayer.

<span class="mw-page-title-main">SNARE protein</span> Protein family

SNARE proteins – "SNAPREceptors" – are a large protein family consisting of at least 24 members in yeasts, more than 60 members in mammalian cells, and some numbers in plants. The primary role of SNARE proteins is to mediate the fusion of vesicles with the target membrane; this notably mediates exocytosis, but can also mediate the fusion of vesicles with membrane-bound compartments. The best studied SNAREs are those that mediate the release of synaptic vesicles containing neurotransmitters in neurons. These neuronal SNAREs are the targets of the neurotoxins responsible for botulism and tetanus produced by certain bacteria.

<span class="mw-page-title-main">ARFGAP1</span> Protein-coding gene in the species Homo sapiens

ADP-ribosylation factor GTPase-activating protein 1 is an enzyme that in humans is encoded by the ARFGAP1 gene. Two transcript variants encoding different isoforms have been found for this gene.

<span class="mw-page-title-main">Erich Sackmann</span>

Erich Sackmann is a German experimental physicist and a pioneer of biophysics in Europe.

Lipid bilayer mechanics is the study of the physical material properties of lipid bilayers, classifying bilayer behavior with stress and strain rather than biochemical interactions. Local point deformations such as membrane protein interactions are typically modelled with the complex theory of biological liquid crystals but the mechanical properties of a homogeneous bilayer are often characterized in terms of only three mechanical elastic moduli: the area expansion modulus Ka, a bending modulus Kb and an edge energy . For fluid bilayers the shear modulus is by definition zero, as the free rearrangement of molecules within plane means that the structure will not support shear stresses. These mechanical properties affect several membrane-mediated biological processes. In particular, the values of Ka and Kb affect the ability of proteins and small molecules to insert into the bilayer. Bilayer mechanical properties have also been shown to alter the function of mechanically activated ion channels.

One property of a lipid bilayer is the relative mobility (fluidity) of the individual lipid molecules and how this mobility changes with temperature. This response is known as the phase behavior of the bilayer. Broadly, at a given temperature a lipid bilayer can exist in either a liquid or a solid phase. The solid phase is commonly referred to as a “gel” phase. All lipids have a characteristic temperature at which they undergo a transition (melt) from the gel to liquid phase. In both phases the lipid molecules are constrained to the two dimensional plane of the membrane, but in liquid phase bilayers the molecules diffuse freely within this plane. Thus, in a liquid bilayer a given lipid will rapidly exchange locations with its neighbor millions of times a second and will, through the process of a random walk, migrate over long distances.

<span class="mw-page-title-main">Lipid bilayer fusion</span>

In membrane biology, fusion is the process by which two initially distinct lipid bilayers merge their hydrophobic cores, resulting in one interconnected structure. If this fusion proceeds completely through both leaflets of both bilayers, an aqueous bridge is formed and the internal contents of the two structures can mix. Alternatively, if only one leaflet from each bilayer is involved in the fusion process, the bilayers are said to be hemifused. In hemifusion, the lipid constituents of the outer leaflet of the two bilayers can mix, but the inner leaflets remain distinct. The aqueous contents enclosed by each bilayer also remain separated.

A model lipid bilayer is any bilayer assembled in vitro, as opposed to the bilayer of natural cell membranes or covering various sub-cellular structures like the nucleus. They are used to study the fundamental properties of biological membranes in a simplified and well-controlled environment, and increasingly in bottom-up synthetic biology for the construction of artificial cells. A model bilayer can be made with either synthetic or natural lipids. The simplest model systems contain only a single pure synthetic lipid. More physiologically relevant model bilayers can be made with mixtures of several synthetic or natural lipids.

Membrane curvature is the geometrical measure or characterization of the curvature of membranes. The membranes can be naturally occurring or man-made (synthetic). An example of naturally occurring membrane is the lipid bilayer of cells, also known as cellular membranes. Synthetic membranes can be obtained by preparing aqueous solutions of certain lipids. The lipids will then "aggregate" and form various phases and structures. According to the conditions and the chemical structures of the lipid, different phases will be observed. For instance, the lipid POPC tends to form lamellar vesicles in solution, whereas smaller lipids, such as detergents, will form micelles if the CMC was reached. There are five commonly proposed mechanisms by which membrane curvature is created, maintained, or controlled: lipid composition, shaped transmembrane proteins, protein motif insertion/BAR domains, protein scaffolding, and cytoskeleton scaffolding.

Hydrophobic mismatch is the difference between the thicknesses of hydrophobic regions of a transmembrane protein and of the biological membrane it spans. In order to avoid unfavorable exposure of hydrophobic surfaces to water, the hydrophobic regions of transmembrane proteins are expected to have approximately the same thickness as the hydrophobic region of the surrounding lipid bilayer. Nevertheless, the same membrane protein can be encountered in bilayers of different thickness. In eukaryotic cells, the plasma membrane is thicker than the membranes of the endoplasmic reticulum. Yet all proteins that are abundant in the plasma membrane are initially integrated into the endoplasmic reticulum upon synthesis on ribosomes. Transmembrane peptides or proteins and surrounding lipids can adapt to the hydrophobic mismatch by different means.

Mechanosensitive channels (MSCs), mechanosensitive ion channels or stretch-gated ion channels are membrane proteins capable of responding to mechanical stress over a wide dynamic range of external mechanical stimuli. They are present in the membranes of organisms from the three domains of life: bacteria, archaea, and eukarya. They are the sensors for a number of systems including the senses of touch, hearing and balance, as well as participating in cardiovascular regulation and osmotic homeostasis (e.g. thirst). The channels vary in selectivity for the permeating ions from nonselective between anions and cations in bacteria, to cation selective allowing passage Ca2+, K+ and Na+ in eukaryotes, and highly selective K+ channels in bacteria and eukaryotes.

Vesicle fusion is the merging of a vesicle with other vesicles or a part of a cell membrane. In the latter case, it is the end stage of secretion from secretory vesicles, where their contents are expelled from the cell through exocytosis. Vesicles can also fuse with other target cell compartments, such as a lysosome. Exocytosis occurs when secretory vesicles transiently dock and fuse at the base of cup-shaped structures at the cell plasma membrane called porosome, the universal secretory machinery in cells. Vesicle fusion may depend on SNARE proteins in the presence of increased intracellular calcium (Ca2+) concentration.

<span class="mw-page-title-main">Saffman–Delbrück model</span> Mathematical model of lipid membranes

The Saffman–Delbrück model describes a lipid membrane as a thin layer of viscous fluid, surrounded by a less viscous bulk liquid. This picture was originally proposed to determine the diffusion coefficient of membrane proteins, but has also been used to describe the dynamics of fluid domains within lipid membranes. The Saffman–Delbrück formula is often applied to determine the size of an object embedded in a membrane from its observed diffusion coefficient, and is characterized by the weak logarithmic dependence of diffusion constant on object radius.

<span class="mw-page-title-main">Cell membrane</span> Biological membrane that separates the interior of a cell from its outside environment

The cell membrane is a biological membrane that separates and protects the interior of a cell from the outside environment. The cell membrane consists of a lipid bilayer, made up of two layers of phospholipids with cholesterols interspersed between them, maintaining appropriate membrane fluidity at various temperatures. The membrane also contains membrane proteins, including integral proteins that span the membrane and serve as membrane transporters, and peripheral proteins that loosely attach to the outer (peripheral) side of the cell membrane, acting as enzymes to facilitate interaction with the cell's environment. Glycolipids embedded in the outer lipid layer serve a similar purpose. The cell membrane controls the movement of substances in and out of a cell, being selectively permeable to ions and organic molecules. In addition, cell membranes are involved in a variety of cellular processes such as cell adhesion, ion conductivity, and cell signalling and serve as the attachment surface for several extracellular structures, including the cell wall and the carbohydrate layer called the glycocalyx, as well as the intracellular network of protein fibers called the cytoskeleton. In the field of synthetic biology, cell membranes can be artificially reassembled.

A unilamellar liposome is a spherical liposome, a vesicle, bounded by a single bilayer of an amphiphilic lipid or a mixture of such lipids, containing aqueous solution inside the chamber. Unilamellar liposomes are used to study biological systems and to mimic cell membranes, and are classified into three groups based on their size: small unilamellar liposomes/vesicles (SUVs) that with a size range of 20–100 nm, large unilamellar liposomes/vesicles (LUVs) with a size range of 100–1000 nm and giant unilamellar liposomes/vesicles (GUVs) with a size range of 1-200 µm. GUVs are mostly used as models for biological membranes in research work. Animal cells are 10–30 µm and plant cells are typically 10–100 µm. Even smaller cell organelles such as mitochondria are typically 1-2 µm. Therefore, a proper model should account for the size of the specimen being studied. In addition, the size of vesicles dictates their membrane curvature which is an important factor in studying fusion proteins. SUVs have a higher membrane curvature and vesicles with high membrane curvature can promote membrane fusion faster than vesicles with lower membrane curvature such as GUVs.

References

  1. Singer, S. Jonathan; Nicolson, Garth L. (1972), "The fluid mosaic model of the structure of cell membranes", Science, 175 (23): 720–731, Bibcode:1972Sci...175..720S, doi:10.1126/science.175.4023.720, PMID   4333397, S2CID   83851531
  2. Helfrich, Wolfgang (1973), "Elastic properties of lipid bilayers: theory and possible experiments", Zeitschrift für Naturforschung C, 28 (11): 693–703, doi: 10.1515/znc-1973-11-1209 , PMID   4273690, S2CID   24949930
  3. Zhong-Can, Ou-Yang; Helfrich, Wolfgang (1987), "Instability and deformation of a spherical vesicle by pressure", Physical Review Letters, 59 (21): 2486–2488, Bibcode:1987PhRvL..59.2486Z, doi:10.1103/physrevlett.59.2486, PMID   10035563
  4. Zhong-Can, Ou-Yang (1990), "Anchor ring-vesicle membranes", Physical Review A, 41 (8): 4517–4520, Bibcode:1990PhRvA..41.4517O, doi:10.1103/physreva.41.4517, PMID   9903652
  5. Mutz, M.; Bensimon, D. (1991), "Observation of toroidal vesicles", Physical Review A, 43 (8): 4525–4527, Bibcode:1991PhRvA..43.4525M, doi:10.1103/physreva.43.4525, PMID   9905557
  6. Naito, Hiroyoshi; Okuda, Masahiro; Zhong-Can, Ou-Yang (1993), "Counterexample to some shape equations for axisymmetric vesicles", Physical Review E, 48 (3): 2304–2307, Bibcode:1993PhRvE..48.2304N, doi:10.1103/physreve.48.2304, PMID   9960853
  7. Guckenberger, Achim; Gekle, Stephan (2017), "Theory and algorithms to compute Helfrich bending forces: a review", J. Phys. Condens. Matter, 29 (20): 203001, Bibcode:2017JPCM...29t3001G, doi:10.1088/1361-648X/aa6313, PMID   28240220
  8. Saitoh, Akihiko; Takiguchi, Kingo; Tanaka, Yohko; Hotani, Hirokazu (1998), "Opening-up of liposomal membranes by talin", Proceedings of the National Academy of Sciences, 95 (3): 1026–1031, Bibcode:1998PNAS...95.1026S, doi: 10.1073/pnas.95.3.1026 , PMC   18660 , PMID   9448279
  9. Capovilla, R.; Guven, J.; Santiago, J. A. (2002), "Lipid membranes with an edge", Physical Review E, 66 (2): 021607, arXiv: cond-mat/0203335 , Bibcode:2002PhRvE..66b1607C, doi:10.1103/physreve.66.021607, PMID   12241189, S2CID   8529667
  10. Tu, Z. C.; Z. C., Ou-Yang (2003), "Lipid membranes with free edges", Physical Review E, 68 (6): 061915, arXiv: cond-mat/0305700 , Bibcode:2003PhRvE..68f1915T, doi:10.1103/physreve.68.061915, PMID   14754242, S2CID   30907597
  11. Asgari, M.; Biria, A. (2015), "Free energy of the edge of an open lipid bilayer based on the interactions of its constituent molecules", International Journal of Non-linear Mechanics, 76: 135–143, arXiv: 1502.05036 , Bibcode:2015IJNLM..76..135A, doi:10.1016/j.ijnonlinmec.2015.06.001, PMC   4509687 , PMID   26213414
  12. Biria, A.; Maleki, M.; Fried, E (2013), "Continuum theory for the edge of an open lipid bilayer", Advances in Applied Mechanics, 46: 1–68, doi:10.1016/B978-0-12-396522-6.00001-3, ISBN   9780123965226

Bibliography

Reviews on configurations of lipid vesicles

[1] R. Lipowsky, The Conformation of Membranes, Nature 349 (1991) 475-481.

[2] U. Seifert, Configurations of Fluid Membranes and Vesicles, Adv. Phys. 46 (1997) 13-137.

[3] Z. C. Ou-Yang, J. X. Liu and Y. Z. Xie, Geometric Methods in the Elastic Theory of Membranes in Liquid Crystal Phases (World Scientific, Singapore, 1999).

[4] A. Biria, M. Maleki and E. Fried, (2013). Continuum theory for the edge of an open lipid bilayer, Advances in Applied Mechanics 46 (2013) 1-68.

Research papers on closed vesicles

[1] W. Helfrich, Elastic Properties of Lipid Bilayers—Theory and Possible Experiments, Z. Naturforsch. C 28 (1973) 693-703.

[2] O.-Y. Zhong-Can and W. Helfrich, Instability and Deformation of a Spherical Vesicle by Pressure, Phys. Rev. Lett. 59 (1987) 2486-2488.

[3] O.-Y. Zhong-Can, Anchor Ring-Vesicle Membranes, Phys. Rev. A 41 (1990) 4517-4520.

[4] H. Naito, M. Okuda, and O.-Y. Zhong-Can, Counterexample to Some Shape Equations for Axisymmetric Vesicles, Phys. Rev. E 48 (1993) 2304-2307.

[5] U. Seifert, Vesicles of toroidal topology, Phys. Rev. Lett. 66 (1991) 2404-2407.

[6] U. Seifert, K. Berndl, and R. Lipowsky, Shape transformations of vesicles: Phase diagram for spontaneous- curvature and bilayer-coupling models, Phys. Rev. A 44 (1991) 1182-1202.

[7] L. Miao, et al., Budding transitions of fluid-bilayer vesicles: The effect of area-difference elasticity, Phys. Rev. E 49 (1994) 5389-5407.

Research papers on open membranes

[1] A. Saitoh, K. Takiguchi, Y. Tanaka, and H. Hotani, Opening-up of liposomal membranes by talin, Proc. Natl. Acad. Sci. 95 (1998) 1026-1031.

[2] R. Capovilla, J. Guven, and J.A. Santiago, Lipid membranes with an edge, Phys. Rev. E 66 (2002) 021607.

[3] R. Capovilla and J. Guven, Stresses in lipid membranes, J. Phys. A 35 (2002) 6233-6247.

[4] Z. C. Tu and Z. C. Ou-Yang, Lipid membranes with free edges, Phys. Rev. E 68, (2003) 061915.

[5] T. Umeda, Y. Suezaki, K. Takiguchi, and H. Hotani, Theoretical analysis of opening-up vesicles with single and two holes, Phys. Rev. E 71 (2005) 011913.

[6] A. Biria, M. Maleki and E. Fried, (2013). Continuum theory for the edge of an open lipid bilayer, Advances in Applied Mechanics 46 (2013) 1-68.

Numerical solutions on lipid membranes

[1] J. Yan, Q. H. Liu, J. X. Liu and Z. C. Ou-Yang, Numerical observation of nonaxisymmetric vesicles in fluid membranes, Phys. Rev. E 58 (1998) 4730-4736.

[2] J. J. Zhou, Y. Zhang, X. Zhou, Z. C. Ou-Yang, Large Deformation of Spherical Vesicle Studied by Perturbation Theory and Surface Evolver, Int J Mod Phys B 15 (2001) 2977-2991.

[3] Y. Zhang, X. Zhou, J. J. Zhou and Z. C. Ou-Yang, Triconcave Solution to the Helfrich Variation Problem for the Shape of Lipid Bilayer Vesicles is Found by Surface Evolver, In. J. Mod. Phys. B 16 (2002) 511-517.

[4] Q. Du, C. Liu and X. Wang, Simulating the deformation of vesicle membranes under elastic bending energy in three dimensions, J. Comput. Phys. 212 (2006) 757.

[5] X. Wang and Q. Du, physics/0605095.

Selected papers on cell membranes

[1] Y. C. Fung and P. Tong, Theory of the Sphering of Red Blood Cells, Biophys. J. 8 (1968) 175-198.

[2] S. K. Boey, D. H. Boal, and D. E. Discher, Simulations of the Erythrocyte Cytoskeleton at Large Deformation. I. Microscopic Models, Biophys. J. 75 (1998) 1573-1583.

[3] D. E. Discher, D. H. Boal, and S. K. Boey, Simulations of the Erythrocyte Cytoskeleton at Large Deformation. II. Micropipette Aspiration, Biophys. J. 75 (1998) 1584-1597.

[4] E. Sackmann, A.R. Bausch and L. Vonna, Physics of Composite Cell Membrane and Actin Based Cytoskeleton, in Physics of bio-molecules and cells, Edited by H. Flyvbjerg, F. Julicher, P. Ormos And F. David (Springer, Berlin, 2002).

[5] G. Lim, M. Wortis, and R. Mukhopadhyay, Stomatocyte–discocyte–echinocyte sequence of the human red blood cell: Evidence for the bilayer–couple hypothesis from membrane mechanics, Proc. Natl. Acad. Sci. 99 (2002) 16766-16769.

[6] Z. C. Tu and Z. C. Ou-Yang, A Geometric Theory on the Elasticity of Bio-membranes, J. Phys. A: Math. Gen. 37 (2004) 11407-11429.

[7] Z. C. Tu and Z. C. Ou-Yang, Elastic theory of low-dimensional continua and its applications in bio- and nano-structures,arxiv:0706.0001.