Home | About | Journals | Submit | Contact Us | Français |

**|**HHS Author Manuscripts**|**PMC3169433

Formats

Article sections

Authors

Related links

Phys Rev Lett. Author manuscript; available in PMC 2011 September 8.

Published in final edited form as:

Published online 2010 November 24. doi: 10.1103/PhysRevLett.105.228101

PMCID: PMC3169433

NIHMSID: NIHMS319326

Department of Applied Physics, California Institute of Technology, Pasadena, California 91125, USA

The publisher's final edited version of this article is available at Phys Rev Lett

See other articles in PMC that cite the published article.

Motivated by recent experiments on bilayer polyhedra composed of amphiphilic molecules, we study the elastic bending energies of bilayer vesicles forming polyhedral shapes. Allowing for segregation of excess amphiphiles along the ridges of polyhedra, we find that bilayer polyhedra can indeed have lower bending energies than spherical bilayer vesicles. However, our analysis also implies that, contrary to what has been suggested on the basis of experiments, the snub dodecahedron, rather than the icosahedron, generally represents the energetically favorable shape of bilayer polyhedra.

In an aqueous environment, amphiphilic molecules such as lipids are observed to self-organize into bilayer vesicles [1,2], thus forming the physical basis for cell membranes. Bilayer vesicles generally exhibit shapes with constant or smoothly varying curvature [3,4]. But in recent experiments [5,6], bilayer vesicles with polyhedral shape, consisting of flat faces connected by ridges and vertices with high local curvature, have been observed. In these experiments, two types of oppositely charged, single-tailed amphiphiles were used, with a slight excess of one amphiphile species over the other. At high temperatures, the amphiphiles were found to form spherical bilayer vesicles. However, provided that the number of excess, unpaired amphiphiles was tuned to some optimal range, cooling the system below the chain melting temperature yielded polyhedral bilayer vesicles. It was reported that the bilayer polyhedra were stable over weeks, and that their shape was consistently reproduced upon thermal cycling. Furthermore, it was suggested [5,6] that the observed polyhedra had icosahedral symmetry, although some uncertainty regarding the polyhedral symmetry remained.

What is the mechanism governing the formation and symmetry of bilayer polyhedra? It was argued [5,6] on the basis of the experimental phenomenology that elastic contributions to the polyhedron free energy dominate over entropic or electrostatic contributions, and that minimization of elastic bending energy alone determines the shape of bilayer polyhedra. In this Letter, we take these intriguing observations as our starting point and address the general problem of finding polyhedral shapes with minimal bending energy. Questions regarding the minimal energy shape of bilayer vesicles are commonly answered using a variational approach [1–3]. However, a given polyhedral shape is defined by the geometric parameters characterizing its vertices and ridges. Thus, polyhedra are inherently of singular nature, which severely restricts the applicability of variational calculus. We therefore employ a complementary method, in which we allow for ridges and vertices with arbitrary geometric properties and, on this basis, calculate polyhedron bending energies as a function of polyhedron symmetry. In the remainder of this Letter, we first consider the most straightforward case of bilayers with *uniform* composition, and then turn to the richer case in which there is segregation of excess amphiphiles.

The solution of the two-dimensional equations of elasticity is a formidable challenge, and has only been achieved for the vertices and ridges of polyhedra in certain limiting cases [7–11]. Thus, in order to determine the elastic energies of arbitrary polyhedral shapes, we mainly employ simple expressions based on the Helfrich-Canham-Evans free energy of bending [1–3],

$$G=\frac{{K}_{b}}{2}\int dS{\left(\frac{1}{{R}_{1}}+\frac{1}{{R}_{2}}-{H}_{0}\right)}^{2},$$

(1)

where *K _{b}* is the bilayer bending rigidity,

Figure 1(a) shows schematic illustrations of a bilayer bending gradually (left panel) and sharply (right panel) along a ridge with dihedral angle *α _{i}*. The first model is inspired by the electron micrographs of bilayer polyhedra in Refs. [5,6,12], while the second model provides a more faithful representation of the polyhedral geometry. Based on the picture presented in the left panel of Fig. 1(a), we approximate ridges by a bilayer bending partially around a cylinder of radius

Illustration of the contributions to the elastic bending energies of polyhedra: (a) Side view of a ridge with dihedral angle *α*_{i}, (b) vertex with face angle *β*_{j}, (c) cross section of half of a pore around the tip of a cone (see inset) with **...**

$${G}_{r}=\frac{{\overline{K}}_{b}}{2}{(\pi -{\alpha}_{i})}^{2}l,$$

(2)

where *K̄ _{b}* =

$${G}_{v}=\frac{{K}_{b}}{2}\sum _{j}{(\pi -{\beta}_{j})}^{2},$$

(3)

where *β _{j}* denotes the face angle subtended by two ridges meeting at a given polyhedron vertex.

Closed bilayer vertices may break up to form pores [5,6,12], which was suggested [5] as a mechanism for avoiding the curvature singularity associated with closed vertices. Figures 1(c) and 1(d) show two models for pores at the vertices of polyhedra which, similarly as before, are inspired by the experimental images in Refs. [5,6,12] [Fig. 1(c)] and a stricter interpretation of the polyhedral geometry of bilayer vesicles [Fig. 1(d)]. In our first model [see Fig. 1(c)] we approximate the vertex of a given polyhedron by a cone with apex angle *π* − 2*θ*, where *θ* = *π*/2 − arccos(1 − Ω/2*π*) for a solid angle Ω subtended by the polyhedron vertex. In our second model we assume that, along each face, the pore consists of a straight cylindrical edge [see Fig. 1(d), left panel], which bends through an angle *γ _{j}* across a ridge from one face to a neighboring face [see Fig. 1(d), right panel]. In both cases, the elastic pore energy

We now turn to the presence of excess amphiphiles in bilayer polyhedra. The expulsion of excess amphiphiles from flat bilayers and the resulting molecular segregation, together with the high spontaneous curvature of single-tailed excess amphiphiles, are thought [5,6] to have two principal effects on the bending energies of polyhedra. On the one hand, excess amphiphiles can seed pores into bilayers [5,6] and, thus, pores may have a role beyond reducing the elastic energy of polyhedron vertices. In particular, it has been suggested [6] that excess amphiphiles produce pores in the spherical bilayer vesicles from which bilayer polyhedra originate upon cooling. On the other hand, it has been found [6] that excess amphiphiles preferentially accumulate along the ridges of polyhedra. As a result, molecular segregation can decrease the bending energy of the outer monolayer at ridges.

The above observations suggest a simple description of how vertex and ridge energies are modified by the presence of excess amphiphiles. Ideally, excess amphiphiles are arranged along ridges such that they induce an anisotropic spontaneous curvature commensurate with the dihedral angle. Assuming such “perfect segregation,” only the bending energy of the inner layer must be considered when computing ridge energies. We therefore obtain a lower bound on the modified ridge energy which takes a similar form as Eq. (2), but with the rescaled *bilayer* bending modulus *K̄ _{b}* replaced by the rescaled

From the simple model for perfectly segregated bilayer polyhedra described above one finds that the optimal ratio of the amphiphile species in excess to the total amphiphile content is given by *r _{I}* ≈ 0.51 for the polyhedron sizes observed in experiments [5,6]. This optimal value for

For certain limits of the equations of elasticity, approximate solutions corresponding to polyhedron vertex [7–9] and ridge [10,11] energies have been obtained. In particular, it has been found [7] that fivefold disclinations in hexagonal lattices are accommodated for small lattice sizes through a stretching of lattice vectors. However, for large enough lattice sizes it becomes energetically favorable to buckle out of the plane [7], in which case the energetics of the system are dominated by bending. This behavior is characterized by a dimensionless quantity known as the Föppl-von Kármán number Γ = *YR*^{2}*/K _{b}*, where

The lowest energy states of icosadeltahedral triangulations of the sphere are found to resemble icosahedra for Γ ≥ 10^{7} [8,9], which corresponds to a vertex energy greater than 8*K _{b}* with, for instance, a value 12

Ridges impose an energetic cost and, hence, one expects that for a fixed area and dihedral angle the faces of polyhedra relax to form regular polygons. We therefore focus here on the convex polyhedra with regular polygons as faces, but we have also considered polyhedra with irregular faces [13]. The class of convex polyhedra with regular faces encompasses the five Platonic solids, the 13 Archimedean solids, the two (infinitely large) families of prisms and antiprisms, and the 92 Johnson solids [15,16]. It has been shown [17] that this list exhausts all convex polyhedra with regular faces. Thus, counting prisms and antiprisms as one solid each, there are exactly 112 convex polyhedra with regular faces. With each of these polyhedra a specific set of parameters characterizing ridges, vertices, and faces is associated, leading to distinct contributions to the total elastic bending energy.

We have evaluated the total elastic bending energies of all convex polyhedra with regular faces. As illustrated in Fig. 2, polyhedron energies are compared by plotting elastic energy as a function of the polyhedron radius *R _{p}* [8], which is related to the polyhedron area

Total elastic bending energies of the convex polyhedra with regular faces, obtained for the case of perfect segregation of excess amphiphiles, and bending energy of the sphere, normalized by the total bending energy of the icosahedron, *G*_{i}, with (a) pores **...**

Image representations of (a) the icosahedron, (b) the snub dodecahedron, (c) the snub cube, and (d) the great rhombicosidodecahedron. The polyhedra in (b) and (c) are chiral.

Allowing for an optimal number of excess amphiphiles, we find (see Fig. 2) that polyhedra can have lower bending energies than the sphere, but only if we permit molecular segregation along *ridges* as observed in Ref. [6]. Segregation at pores, which was originally suggested in Ref. [5] as a potential mechanism stabilizing polyhedral shapes, is *not* sufficient to produce polyhedra with bending energies which are favorable compared to the sphere [13]. Indeed, if we assume that pores are closed, bilayer polyhedra are energetically favorable for the experimentally observed polyhedron radius *R _{p}* ≈ 500 nm [see Fig. 2(a)], with the snub dodecahedron as the minimum energy shape among the convex polyhedra with regular faces. If we allow pores of a finite size, a sequence of polyhedral shapes is obtained which minimize bending energy for smaller polyhedron radii [see Fig. 2(b)]. The most notable of these polyhedral shapes is the great rhombicosidodecahedron [see Fig. 3(d)], which surpasses the snub dodecahedron in bending energy at

In summary, we have used Eq. (1) and a variety of other expressions [7–11] to systematically evaluate the bending energies of bilayer polyhedra. We find that, contrary to what has been suggested on the basis of experiments [5,6], the snub dodecahedron and the snub cube generally have lower total elastic bending energies than the icosahedron. This result is consistent with several complementary theoretical studies which suggest that the elastic energies of chiral shapes such as the snub dodecahedron and the snub cube can be favorable compared to the icosahedron [18–20] and that, even if the icosahedral shape is imposed, the minimum energy structure may still be chiral [21]. While we followed here the experimental phenomenology [5,6] and assumed that minimization of bending energy governs the shape of bilayer polyhedra, other contributions to the free energy, as well as kinetic effects [14,22], could, in principle, modify the preferred polyhedral symmetry. In light of our results, we suggest revisiting the symmetry of bilayer polyhedra, and the thermodynamic or kinetic mechanisms potentially governing their formation and stability, in greater experimental detail.

This work was supported by a Collaborative Innovation Award of the Howard Hughes Medical Institute, and the National Institutes of Health through NIH Award number R01 GM084211 and the Director’s Pioneer Award. We thank A. Agrawal, M. B. Jackson, W. S. Klug, R. W. Pastor, T. R. Powers, D. C. Rees, M. H. B. Stowell, D. P. Tieleman, T. S. Ursell, and H. Yin for helpful comments.

1. Boal D. Mechanics of The Cell. Cambridge University Press; Cambridge: 2002.

2. Safran SA. Statistical Thermodynamics of Surfaces, Interfaces, and Membranes. Westview Press; Boulder: 2003.

3. Seifert U. Adv Phys. 1997;46:13.

4. Lim GHW, Wortis M, Mukhopadhyay R. Proc Natl Acad Sci USA. 2002;99:16 766 . [PubMed]

5. Dubois M, et al. Nature (London) 2001;411:672. [PubMed]

6. Dubois M, et al. Proc Natl Acad Sci USA. 2004;101:15 082 .

7. Seung HS, Nelson DR. Phys Rev A. 1988;38:1005. [PubMed]

8. Lidmar J, Mirny L, Nelson DR. Phys Rev E. 2003;68:051910. [PubMed]

9. Nguyen TT, Bruinsma RF, Gelbart WM. Phys Rev E. 2005;72:051923. [PubMed]

10. Lobkovsky AE. Phys Rev E. 1996;53:3750. [PubMed]

11. Lobkovsky AE, Witten TA. Phys Rev E. 1997;55:1577.

12. Glinel K, et al. Langmuir. 2004;20:8546. [PubMed]

13. Haselwandter CA, Phillips R. to be published

14. Hartmann MA, et al. Phys Rev Lett. 2006;97:018106. [PubMed]

15. Cromwell PR. Polyhedra. Cambridge University Press; Cambridge: 1997.

16. Berman M. J Franklin Inst. 1971;291:329.

17. Johnson NW. Can J Math. 1966;18:169.Zalgaller VA. Sem in Math. Vol. 2. Steklov Math. Institute; Leningrad: Consultants Bureau; New York: 1967. 1969, (English translation)

18. Bruinsma RF, et al. Phys Rev Lett. 2003;90:248101. [PubMed]

19. Zandi R, et al. Proc Natl Acad Sci USA. 2004;101:15 556 . [PubMed]

20. Dodgson MJW. J Phys A. 1996;29:2499.

21. Vernizzi G, Olvera de la Cruz M. Proc Natl Acad Sci USA. 2007;104:18382. [PubMed]

22. Noguchi H. Phys Rev E. 2003;67:041901. [PubMed]

PubMed Central Canada is a service of the Canadian Institutes of Health Research (CIHR) working in partnership with the National Research Council's national science library in cooperation with the National Center for Biotechnology Information at the U.S. National Library of Medicine(NCBI/NLM). It includes content provided to the PubMed Central International archive by participating publishers. |