Arxiv:1303.1793V2 [Cond-Mat.Soft]

arXiv:1303.1793v2 [cond-mat.soft] 9 Jul 2013 h xeietlsd aee om r on in freez- found the temperature, are below ing forms vesicles faceted phosphatidylcholine small side experimental the ag iue,tepoenecoe atra or- bacterial many carboxysomes protein-enclosed of as capsids known the ex- ob- protein ganelles Prominent viruses, the often are large one well. nature in as smooth amples structures are faceted systems serves biological in ueo noi n aincamphiphiles, cationic and anionic of ture vesicles, polymerized hpdbcei on nstrtdbiepools. brine saturated in found bacteria shaped odnanocages. gold e ihavreyo ebaeadtransmem- and membrane of variety proteins. brane a with embed- of ded structures made bilayer phospholipid vesicles typically of are membranes aid Biological compartmentaliza- the membranes. The with molecularly-thin achieved to organism. was able con- an tion be of to to parts that necessary ent differ- separate clear was and compartmentalize it became hierarchically vesicles. life it smooth complex of evolution struct physics in the Early understand to tial vertices. singular and edges sharp structures re- their develop- faceted with of in Only understanding interest theoretical growing a vesicles. a ing been smooth there under- of has sub- to cently shapes decades devoted four the been past has standing the effort Over research form- stantial by facets. energy pos- stable total of the ing feature lower common to pathways the sessing share systems liquid- these all with copolymers chains, block side from crystalline formed cles ically hslreapc ai uhmmrnscnbe can membranes mem- Vesicular such Given quasi-two-dimensional. ratio as microns. aspect treated to large nanometers this of hundreds from lhuhms ftelpdblyrvsce found vesicles bilayer lipid the of most Although oudrtn h rgno aeigi sessen- is it faceting of origin the understand To ∼ ooeeu ytm pn o-rva ahast h form the to vesi pathways al elastic structures. non-trivial or homogeneous opens properties, system, a elastic homogeneous the in in icosahedra heterogeneities co of troducing lattice I formation crystalline topology. the any spherical to in with present po ela I membranes necessarily irregular of elastic defects, and faceting of surface. regular the states of curved in ergy assortment a result an of can conditions, context structure right this and in shape frustration the affects foundly 5nm h nepa ewe emty oooyadodrcnlead can order and topology geometry, between interplay The .INTRODUCTION I. 14 hl h aea iesosrange dimensions lateral the while 13 3 h ebaetikesi typ- is thickness membrane The hl ies nterstructure, their in diverse While oectoi vesicles, cationic some 6 eilsasmldfo mix- a from assembled vesicles 9–11 yaueUiest,Srcs,NwYr 34,USA 13244, York New Syracuse, University, Syracuse in fullerenes, giant eateto hsc n otMte Program, Matter Soft and Physics of Department ahast aeigo vesicles of faceting to Pathways akJ Bowick J. Mark 1 n h square- the and 5,7,8 4,5 12 cross- 2 vesi- and On ∗ n atoSknepnek Rastko and asincurvature, Gaussian htrltstebnigeeg oteitga of curvature, integral mean the the functional to of energy energy square free bending the a the shape constructed relates changing He that to corresponding space. in bending, deforma- is dominant their tion nature that quasi-two-dimensional implies membranes the of that argued per, where omn hr ed scmaal omlclren- molecular at to off for comparable is cut penalty bends energy is sharp the divergence forming scale, the length molecular systems some diverging real or in high While very structures with Sharp curvature regions edges. to or for- the correspond corners by sharp generated of be mation would re- that to stress themselves any move rearrange quickly molecules since rn tef sse,freape nlpdrafts. lipid in example, for mem- seen, also the as It itself, within quickly. brane reorganizations closed be structural can allows from kind, or some of transport, aper- damage for transmembrane small by important that formed ensures is tures it cost Fluidity as energy membranes no biological deformations. is there shear as for membranes diffusion lateral fluid free describes details. with energy molecular free the Helfrich on The rigidities depend Gaussian themselves and which bending the respectively are clpoete a egie ycntutn sim- constructing Helfrich, by models. gained mesoscopic be ple phys- can their into properties insight still much ical understood, and fully complex be exceedingly to is membranes logical compounds. inorganic nanocages, gold of case the r a ecntutdi h aoaoyuiglipid using laboratory the in contain- molecules, constructed sys- micron-sized be biological can and ers to Nanometer-sized however, limited, tems. not are branes ooeeu udvsce r yial smooth typically are vesicles fluid Homogeneous hl h ealdsrcueadfnto fbio- of function and structure detailed the While fie oshrcltplg,ntrlylead naturally topology, spherical to nfined oigfrnnlna edn epneo a of response bending non-linear for lowing E H yerlsrcue a etelwen- low the be may structures lyhedral tcvsce.W hwta,udrthe under that, show We vesicles. stic bend atclr eso o topological how show we particular, n 0 to ffctd e non-icosahedral, yet faceted, of ation 20 l.Frhroe eso htin- that show we Furthermore, cle. stesotnoscraueand curvature spontaneous the is 15,16 hscmetr eso how show we commentary this n ogoercfutainta pro- that frustration geometric to = n r hrfr nreial costly. energetically therefore are and ˆ block-copolymers dS † h K 2 i.e. , κ ( H − H 0 18 17 ) 2 nasmnlpa- seminal a in + n vn sin as even, and H κ n othe to and , G K κ i , and 19 (1) κ G 2

ergies and far exceeds the energy cost ( kBT ) of briefly review the complex question of the existence long-wavelength undulations. To form sharp∼ struc- of long-range order on curved surfaces. tures, therefore, there must be a mechanism that ei- Constructing a regular lattice in the plane (or any ther reduces the energy cost associated with the for- other flat surface) is straightforward. The simplest mation of bends or that inhibits the flow of molecules example is the ground state of identical point-like or imposes additional “in-vesicle” order (e.g. liquid particles in an external confining potential interact- crystalline) that would select a non-smooth underly- ing via long range repulsive interactions, e.g. identi- ing geometry. The first mechanism can be achieved cal charges confined to a planar region. The ground by introducing asymmetry in the composition of the state is a triangular lattice with each particle having two leaflets of the bilayer induced by a phase segre- exactly six equidistant nearest neighbors and a lat- gation of the two molecular species, thus resulting 21 tice spacing that is determined by the density. Such in local non-zero spontaneous curvature or by seg- a state is not frustrated and is stress free. This is regating excess amphiphiles along long ridges and no longer the case when one attempts a similar con- inducing a spontaneous curvature that is commen- 7,22 struction on the 2-sphere (the surface of a ball in surate with the dihedral angle between two faces. R3). Any attempt to map a regular lattice on to the In this paper we discuss mechanisms that lead to surface of a 2-sphere will result in excess particles. faceting of elastic vesicles. An interesting account In order to build a smooth covering of the 2-sphere of faceting of giant vesicles in the rippled gel phase it is necessary, in other words, to remove parts of Pβ′ by osmotic deflation has recently been given by 23 the original planar lattice. As a result, even in the Quemeneur, et al. ground state any crystalline lattice on a 2-sphere This paper is organized as follows. In Section II A will necessarily have a finite number of sites that we discus pathways to the faceting of elastic mem- have coordination different from six, i.e. so called branes that describe viral capsids or bilayer systems topological defects. Starting from Euler’s polyhe- cooled below the gelation transition or for which dral formula, which states that for any polyhedron the molecules are cross-polymerized (i.e. tethered V E +F = χ, where V is the number of vertices, E to each other). In particular, we show, in Section is− the number of edges, F is the number of faces, and IIA, how topological defects can act as the seeds χ is the Euler characteristic of the polyhedron, one of buckling to icosahedra, thus providing a simple can show that the total number of vertices of a crys- explanation for the observation that viruses with talline lattice on a sphere has to satisfy the following icosahedral symmetry are spherical for small sizes 25 relation: (6 z) Nz =6χ = 12. Here Nz is the but well-faceted icosahedra for larger sizes. In Sec- z − total numberP of sites with coordination z and we tion II B we show how the presence of elastic inho- have explicitly used that χ =2 for a 2-sphere. If we mogeneities also can lead to faceting to regular and restrict ourselves to coordination numbers z = 5, 6 irregular polyhedra quite different from icosahedra. or z = 7 (physically coordination numbers are near Finally, in Section II C, we argue that faceting can 6), the configuration with the lowest number of de- also occur in homogeneous elastic vesicles if one al- fects on a sphere has N5 = 12, i.e. twelve five-fold lows for a non-linear bending response in terms of a coordinated sites embedded in a regular triangular critical curvature. lattice. We note that, somewhat counterintuitively, this configuration is not always the ground state of the system as shown by Bowick, et al.26 For suffi- II. FACETING OF ELASTIC MEMBRANES ciently large total number of particles, and reasonable defect core energies, the total energy can be A. Buckling into icosahedra lowered by forming pairs of five- and seven-fold defects that emanate from each of the twelve five-fold One of the hallmarks of viruses is the regular coordinated sites. In the following discussion we will structure of their capsids. The unit building blocks ignore this complication and assume that only twelve of capsids are capsomeres, robust protein complexes, five-fold defects are present. 10nm in diameter, that tile the entire capsid in Despite their appearance five-fold disclination de- a≈ regular (crystalline) array. The lattice structures fects are not local since they owe their existence of icosahedral viruses can be described in terms of to removing whole sections of a planar lattice so pairs of non-negative integers, (p, q) that form a that it conforms to a 2-sphere. As such they are T number, T = p2 + q2 + pq.24 The existence of endowed with considerable elastic energy. To es- a− lattice combined with the spherical topology is at timate how the elastic energy of a disclination de- the center of the argument for size-driven buckling pends on its size we assume that it is surrounded by into icosahedra. To properly explain the mechanism a Hookean elastic medium. The elastic energy of a 27 2 ijkl behind this buckling phenomenon we need first to two-dimensional crystal is Eel = d ~rA uij ukl, ´ 3

1 where uij (~r) = 2 (∂iuj + ∂j ui + ∂kui∂kuj ) is the strain tensor field, with ui (~r) (i = 1, 2) being the components of the displacement field, and Aijkl is the material-dependent elastic tensor. For a hexag- onal lattice, the elastic tensor has only two independent components,27 i.e. two Lamé coefficients µ and λ and the elastic energy simplifies to27 1 E = d2~r λ (u )2 +2µu u , (2) el 2 ˆ ii ij ij (a) (b) and we assume summation over pairs of repeated Figure 1: The transition is triggered if the dimensionless indices. FvK number γ exceeds a critical value γ ∼ 102.30 Simu- A planar five-fold disclination defect can be con- lated structures are constructed using the (p, q) = (12, 3) ◦ structed by removing a 60 wedge of a triangular icosadeltahedral triangulation with ≈ 2 × 103 vertices lattice and sewing the two cut edges together. A de- and radius R ≈ 11.5l0, with l0 = 1 being the aver- tailed calculation28 finds that the stretching energy age lattice spacing. (a) γ ≈ 170 and (b) γ ≈ 3500. 5−fold 2 Red dots designate five-fold defects. All snapshots are of a five-fold disclination is Eel = (Y/32π) R , generated with the Visual Molecular Dynamics (VMD) where Y is the two-dimensional Young’s modulus 31 32 and R is the system size. package and rendered with the Tachyon ray-tracer. If the disclination is allowed to buckle out of the Disclination-driven buckling into an icosahedron is discussed in Ref. 30. Results of that study are reproduced plane it can lower its energy by forming a conical here for completeness. structure with the apex at the disclination. In this case the stretching energy is relieved at the expense of a bending penalty. The bending penalty is given 28,29 numerical prefactors due to partial screening of the by strain induced by the non-zero Gaussian curvature of the sphere.26 As in the planar case, therefore, E = dS 2κH2 + κ K . (3) bend ˆ G the onset of buckling depends solely on the relative strengths of the stretching and bending en- It can be shown28 that the bending energy of such ergies. Since the stretching energy grows rapidly a conical structure is E(5) (π/3) κ log (R/R )+ with the radius of the sphere, all defects buckle bend ≈ b simultaneously into cones once a critical radius is 1 Y R2, where R 154κ/Y . By comparing the 32π b b reached and an icosahedral vesicle is formed.30 In expressions for the≈ stretching and bending energies p Fig. 1 we show the shape of the shell below and it is clear that for small radii stretching energy is above the buckling transition. The transition is not smaller than bending energy and the disclination sharp but rather a rounded smooth crossover, with will remain flat at the expense of stretching its sur- vesicles having FvK number γ . 150 being spheri- rounding. As R increases, the logarithm grows less cal, for 200 . γ . 1500 being noticeably buckled, rapidly than R2 and it becomes energetically favor- and for γ & 2000 being very sharp. For extremely able to buckle into a cone. A detailed calculation28 large FvK numbers (γ & 107) the vesicle is non- shows that the buckling transition occurs when the extensible and one observes characteristic ridge-like dimensionless ratio, the so-called Föppl-von Kármán 2 YR structures with very interesting scaling properties, or FvK number, γ = κ & 154. 7/3 2 1/6 e.g. Eridge/κ α YL /κ , where α is the Another consequence of their non-local character ∼ 33,34 is that disclination defects, much like charges, ex- angle of the ridge and L is its length. Finally we note that imposing internal pressure,35 ert long-range forces on each other. This long-range 35,36 force is transmitted via the elastic deformation of or a volume constraint, shifts the buckling tran- the medium with the speed of sound playing the sition to higher values of γ but leaves the basic role of the speed of light in electrostatics. Discli- physics unchanged. nations with the same topological charge (s =1 for a five-fold defect) repel each other. On a 2-sphere, therefore, disclinations adopt a configuration that B. Buckling of multicomponent vesicles maximizes their separation; this is the icosahedron for twelve defects.26 An implicit assumption in the discussion thus far Moving from planar to spherical geometry does has been that the vesicle is made of an elastically not change the functional forms of the stretching homogenous material. Introducing heterogeneities and bending energies but it does renormalize the by allowing for spatially varying Young’s modulus 4

and bending rigidity can lead to alternative buckling tween two sites at ~r and ~r , respectively, ε i j ij ∈ mechanisms. The presence of heterogeneities, for ex- εhard,εsoft , and the sum is carried out over all ample, is believed to be responsible for the struc- pairs{ of nearest} neighbors. It is convenient to set ture of the recently synthesized assemblies of op- l0 =1 and use it as the unit of length. positely charged amphiphiles.6,8 The structures ob- A proper discretization of the bending en- served are 100nm bilayer vesicles that are faceted ergy is more complicated with a number of dif- but not icosahedral.∼ SAXS/WAXS measurements ferent approaches having been proposed in the combined with atomistic and detailed coarse-grained physics, applied mathematics and computer graph- molecular dynamics simulations have revealed the ics literature.28,39–42 Different approaches vary in existence of strongly correlated crystalline domains the level of computational complexity and accuracy. within the bilayer.8 A typical domain size was found For the present discussion it is sufficient to adopt to be 25nm. Simulations show that the crystal- the simple discretization suggested by Seung and lographic≈ axes of two neighboring domains are not Nelson.28 The discrete bending energy is computed mutually aligned and that the domains are separated as by narrow regions of disordered, or even liquid, am- 1 phiphiles. It is reasonable to expect that the bilayer E˜ = κ˜ (~n ~n )2 bend 2 IJ I − J will be softer along those disordered boundary re- hXI,Ji gions. These regions represent grain boundary de- = κ˜ (1 cos(θ )) , (5) fects with very different physical properties from the IJ − IJ topological defects discussed in the previous section. hXI,Ji To model a vesicle with inhomogeneous elastic pa- where θ is the angle between the unit-length nor- rameters we generalize the expressions for stretching IJ mals ~n and ~n to the neighboring triangles I and J, (Eq. (2)) and bending (Eq. (3)) energies and allow I J and κ˜ is the position dependent discrete bending the two Lamé coefficients and the bending rigidity IJ rigidity also chosen from a binary set κ ,κ . to depend on the spatial position, i.e. λ = λ (~r), hard soft The spring constant and the discrete{ bending rigid-} µ = µ (~r), and κ = κ (~r). We assume for simplicity are related to the Young’s modulus and the bend- ity that each of these coefficients can take only two ing rigidity of the continuum theory by Y =2ε/√328 values, denoted as “hard” and “soft” and loosely cor- √ 28,43 responding to strongly correlated facets and grain (with Poisson ratio ν =1/3) and κ = 3˜κ/2. boundary regions, respectively. We also need to specify which element of the The elasticity theory of thin objects is non-linear37 discrete mesh carries which type of information. and if one allows for spatially varying elastic pa- We note that discrete version of both stretching rameters (even if they have a simple binary distri- and bending energies are most naturally defined on bution) it quickly becomes analytically intractable. edges, i.e. stretching energy is modeled by assigning We therefore turn to numerical simulations. The a linear elastic spring to each edge and bending en- vesicle is represented as a discrete triangulated sur- ergy is proportional to the dihedral angle between face. In the case of modeling viral capsids, the ver- two neighboring triangles sharing an edge. One tices of the triangulation correspond to individual is therefore tempted to assign component types to capsomeres. If the constitutive units are smaller, edges. While possible, this is not the most conve- as in the case of amphiphilic vesicles, the triangula- nient approach as either a vertex or a triangle, and tion can be thought of as a coarse-graining scheme not an edge, is a natural discrete representation of an that assigns discrete elements to patches of actual infinitesimal area element in the continuum descrip- molecules. Each patch is large enough to contain tion. Defining a discrete form of line tension energy, a sufficient number of microscopic degrees of free- as discussed below, is also exceedingly hard for com- dom that molecular details are of no importance, ponent types assigned to edges. We opt rather for but small enough to be considered homogeneous. defining component types on vertices. In this case On a triangular lattice the simplest form of the a vertex i is assumed to carry a half-spring with a discrete stretching energy for a multicomponent sys- spring constant εi which is connected in series to tem is obtained if we assume that the neighbor- the half-spring of its neighbor j with a spring constant εj . The spring constant of the edge (ij) is then ing sites are connected to each other with harmonic −1 −1 −1 springs with spatially-dependent spring constants given by εij = εi + εj . For the bending energy 28 we choose to define the bending rigidity of an edge εij and a uniform rest length l0 as the arithmetic mean of the bending rigidities as- 1 2 signed to its vertices, i.e. κij = (κi + κj ) /2. There E˜el = εij (lij l0) , (4) 2 − are two side effects of defining components on ver- hXi,ji tices. First, the spring constant of an edge shared where l = ~r ~r is the Euclidean distance be- by two vertices of the same type is ε = ε /2; and ij | i − j | ij i 5

a simulation, various cooling protocols can be applied such as linear, exponential or power-law. We note that the annealing temperature is not the actual physical temperature but rather a parameter that controls the acceptance rate of the Monte Carlo moves which increase the energy. During the component swap stage the component-type of a vertex (either soft or hard) is preserved so that the total fraction of each component is also constant. Finally we note that the component swap move is purely a (a) (b) convenient simulation tool that allows sampling of the space of component arrangements – it does not correspond to an actual rearrangement of the material within the vesicle which is assumed solid without diffusion.Throughout the simulation the triangulation is preserved – no edge-flip moves44 that would create or annihilate defects are performed. The main drawback of the simulated annealing approach is that it typically converges very slowly once a minimum has been approached (for the studied system sizes a minimum is typically reached after (c) (d) 105 106 Monte Carlo moves) and there is no guar- antee− that the obtained structure is a true global Figure 2: Snapshots of faceted structures with 3% (a), minimum, i.e. the ground state configuration. We 4% (b), 4.5% (c), and 5% (d) of the soft component (yel- thus refer to the shapes we obtain as typical. On low). Twelve five-fold disclination defects are present, as the other hand the experimentally observed shapes required by the spherical topology, but appear to have no influence on the position of the edges and the size of the may not be true equilibrium shapes either – they are facets. A detailed account of faceting of two-component likely to depend in some part on the details of the elastic vesicles is presented in Ref. 38 and reproduced assembly process. here for completeness. A detailed study38 of the shapes of two-component elastic vesicles has found a wide variety of regular and irregular polyhedral shapes, some of which are second, while the total number of vertices of a given shown in Fig. 2. If a small amount of soft compo- type is conserved, the vertex-type swaps during the nent is added to an otherwise hard vesicle it arranges Monte Carlo optimization do not preserve the total into branching lines. As the amount of soft compo- number of edge spring constants and bending rigidi- nent is increased, these lines start to merge and form ties. Extensive tests, however, did not reveal any facets. Facets are nearly flat with very sharp bends qualitative difference between systems with compo- along the soft ridges. The total number of facets in- nents defined on edges and systems with components creases with the amount of the soft component un- defined on vertices or triangles. til a non-universal, size-dependent concentration is For multicomponent vesicles we require a numeri- reached. Beyond this size-dependent concentration cal minimization technique that simultaneously op- of the soft component it is no longer favorable to timizes the position of the individual components generate new facets but instead the soft component as well as the relative population of components. is randomly distributed inside the existing ones. Al- This can be achieved by using stochastic minimiza- though each configuration still contains twelve five- tion techniques, e.g. the simulated annealing Monte fold defects, it appears that the size of the facets Carlo method. Two components are assigned at ran- and the position of the soft ridges is not sensitive to dom to vertices (or equivalently to edges or triangles) their presence (Fig. 2). We point out that it is possi- of the initial sphere such that the total fraction of the ble to entirely remove the residual strain induced by soft component is f. A Monte Carlo move has two the defects either by appropriately tuning the spring 45 steps: 1) a vertex is displaced by a vector ∆~r with rest lengths or by directly working with reference 46,47 48 components chosen at random from a uniform distri- and realized metrics. Such an approach shows bution in an interval [ ζl0,ζl0] (typically, ζ 0.05) that faceting of the soft-hard two-component vesicles followed by 2) swap of− types (soft hard) of≈ a pair occurs even if the defect contribution to the elastic of randomly selected vertices. In both↔ steps moves energy is removed. are accepted according to Metropolis rules. During Thus far we have not assumed any mixing penalty 6

2 nI

qIJ

k 1

bend

E/ (a) (b)

0 0 qc=p/6 qc=p/3 qc=p/2 p q

Figure 4: Bending energy (in units of the bending rigidity κ) as a function of the angle θ between unit-length (c) (d) normals to two neighboring triangles without a critical π π angle (black), with θc = 6 (red), with θc = 3 (green), Figure 3: Snapshots of faceted elastic vesicles with var- π and with θc = (blue). Inset: Side view of two neigh- ious mixing penalties between soft (yellow) and hard 2 boring triangles I and J with their corresponding unit (gray) components for a system with an equal number length normals ~nI and ~nJ . of hard and soft vertices. The mixing penalty increases from (a) to (d). A detailed discussion of the possible phases in two-component elastic vesicles with line tension is given in Ref. 49. The key results of that study gle between neighboring facets. Once the radius are reproduced here for completeness. of curvature becomes comparable to the molecular length scales the deformation is no longer locally small and microscopic details of the vesicle start to 28 for the two components. The components segregated play a role. Eq. (5) was derived with the assump- solely as a result of the bending preference. If a mix- tion that the curvature is smooth and slowly varying ing penalty, in the form of an effective line tension, Γ, between neighboring points, enabling a perturbative is introduced it competes with the preferred elastic expansion of the bending energy in terms of small organization of the components and leads to an even deviations of the normal vectors to a planar surface. richer variety of shapes.49 In an experimental real- For a locally large deformation this is no longer jus- ization of this system one can envision using the re- tifiable and the expression for the bending energy cently discovered linactant molecules,50–52 which are needs to be modified. This is clearly a very hard two-dimensional analogues of surfactants. For weak problem that requires detailed knowledge of the mi- line tension elasticity dominates and the soft com- croscopic structure of the vesicle shell. Rather than ponent forms highly bent ridges separating flat hard constructing such a detailed model we make the sim- faces. The width of these soft ridges grows as the line ple assumption that the bending energy saturates tension increases. If Γ is increased even further the beyond a critical angle. In other words, Eq. (5) is faces start to coarsen and merge together until the modified to components fully separate at Γl /κ 10−2,49 0 hard E˜crit. =κ ˜ (1 F (θ )) , (6) suggesting that even a small mixing penalty≈ can bend IJ X − lead to full phase segregation. A subset of shapes hI,Ji of two-component vesicles with line tension is shown where in Fig. 3. We note that line tension can also lead to very in- cos(θ) for θ θc F (θ)= ≤ (7) teresting shapes in two-53 and three-component liq- cos(θc) for θ>θc uid vesicles.54 and, as in Eq. (5), θIJ measures the angle between normals to a pair of neighboring triangles and θc is the critical value of that angle beyond which the C. Faceting due to critical curvature ˜crit. bending penalty saturates. A plot of Ebend vs. θ is shown in Fig. 4. In this model θc is treated as an Some of the faceted shapes observed in experi- input parameter. In essence, θc controls the defor- ments, as well as those obtained in simulations (e.g. mation angle beyond which the vesicle shell is suf- Fig. 2), show sharp ridges with a large dihedral an- ficiently deformed that atomistic details start to be 7

this system mapped with 200 independent pairs ≈ of γ and θc. Five distinct regimes are identified. For low values of γ and moderate to large θc (green region, denoted as “sphere” in Fig. 5) the bending penalty is too large and vesicles remain spherical. In this region, the angles between normals on any pair of neighboring triangles are smaller than θc and the existence of an angle cutoff plays no role in determin- ing the shape: we recover the situation discussed in Section II A. If γ is increased, while keeping θc large, one crosses over into the regime where icosahedral shapes are favorable. Since no deformation exceeds θc the transition occurs at the same value of γ as in the case without the critical bending angle (yellow region, denoted as “smooth icosahedra” in Fig. 5). If θc is reduced the existence of an angle cutoff is seen to affect the low energy shapes. In the region where γ is small (red region denoted as “polyhedra” in Fig. 5) one observes irregular polyhedra with sizes and positions of facets and edges that are insensitive Figure 5: Phase diagram for a vesicle with critical bend- to the presence of the five-fold defects. This is to be ing angle as a function of the FvK number (γ) and the expected since for small γ the bending penalty com- value of the critical angle (θc, shown in degrees). We pletely dominates over stretching and defect-driven identify five distinct regimes. In the “sphere” (green) buckling into a cone is suppressed. In this region and “smooth icosahedra” (yellow) regimes θc is large and the faceting mechanism is similar to buckling of the γ does not affect the shape. For small the vesicle is soft component as discussed in the previous section. spherical and crosses over into a smooth icosahedron as As γ is increased topological defects start to play a γ is increased. At lower values of θc the presence of a critical angle starts to play a role and the vesicle shape role and the system transitions to the regime where changes substantially. For small γ the elastic energy is defects drive buckling. With a relatively small value completely dominated by the bending contribution and of θc, however, the system can further lower its en- vesicles take irregular polyhedral shapes (red region). ergy by forming sharp edges that emanate from the For larger values of γ the stretching energy caused by defects. This region is denoted as “sharp icosahedra” topological defects cannot be ignored and the icosahe- (blue in Fig. 5). We note that the direction of the dral symmetry reemerges, but this time accompanied by sharp edges appears to depend on the details of the sharp ridges (blue region). Finally, for very small values triangulation and thus is not universal. Finally, if θc of γ the vesicle is unstable to crumpling (gray). is too small the vesicle becomes unstable and crumples (gray region denoted as “crumpled” in Fig. 5). In this region the simulations often fail to converge important. We assume that the vesicle is homoge- and instead get locked in high-energy states that of- nous sphere with uniform values for κ˜ and ε, and a ten occur in the early stages of a simulation. In standard triangular lattice with twelve five-fold de- Fig. 6 we show typical low energy configurations in fects positioned at the corners of an inscribed icosa- each of the five regimes. hedron. While one could argue that the stretching energy should also be modified by introducing a crit- We point out that this is a toy model with very ical edge length beyond which the energy is constant, little ingredients of an actual system. Therefore, all we use, for simplicity, the expression in Eq. (4) and conclusions should be taken as qualitative with very note that test simulations show that the results are little to no quantitative bearing for real experimental not qualitatively affected by this simplification. systems such as lipid vesicles. However, we believe As before, the low energy states were found using that some general qualitative conclusions can still simulated annealing Monte Carlo simulations. No be deduced. For example, it is natural to ask which constraints on volume or area were imposed. Mini- experimental systems could follow bending mecha- mum energy configurations were found after 106 nism described by this model. We note that for the Monte Carlo moves using a linear cooling protocol.≈ continuum model to be applicable the deformation We studied vesicles with approximately 6 103 ver- has to be small at molecular length scales. In a lipid tices ((p, q) = (18, 10)-chiral; (25, 0)-achiral)· over a vesicle, on the other hand, the angular difference range of values of the FvK number (γ) and critical between the long axes of two neighboring molecules angle (θc). In Fig. 5 we display a phase diagram of should be sufficient (a few degrees) for the non-linear 8

tween crystalline order and spherical topology implies topological defects and residual stress even in the ground state. The stress can be relieved, in part, by buckling into cones seeded at the defects. The buckled structure typically inherits the icosahedral symmetry associated with the defects. This mechanism can qualitatively account for the observation that smaller icosahedral viruses (e.g. Poly- (a) (b) omaviruses55) are smooth (spherical) while larger viruses (e.g. Mimivirus56) are buckled. For viruses the size of each capsomere is sufficiently large ( 10nm) that the capsid is a nearly perfect lattice∼ with 102 103 units and no additional defects beyond∼ the minimum− twelve five-fold sites required by the topology. In the case of lipid bilayer vesicles lipid molecules are much smaller (occupying an area / 0.5nm2) and each vesicle is assembled from hundreds of thousands of lipids. If such a vesicle (c) (d) develops crystalline order as a result of cooling or strong electrostatic correlations it is much harder to form a perfect lattice and one instead expects grain boundaries to form. This effect is enhanced by the spherical topology which hinders formation of a uniform crystal to begin with. In the grain boundary region molecules are disordered and may even be in the liquid state. As a result, the shell is expected to be softer. We have shown that, if one allows (e) for the presence of a softer component in an otherwise crystalline vesicle, the low energy configura- Figure 6: Snapshots of the low energy configurations of tion can be regular and irregular polyhedra that are a vesicle constructed using (p, q) = (18, 10) T −number triangulation with ≈ 6 × 103 vertices and radius R ≈ not icosahedra. The soft component forms ridges 20.5l0. (a) For low a FvK number γ ≈ 50 and a large that are highly bent and separate hard flat facets. ◦ critical angle θc = 20 the vesicle remains spherical. (b) The faceting mechanism is quite different from the 3 ◦ For γ ≈ 10 and θc = 20 one observes a smooth buckled buckling into an icosahedron seeded at the five-fold structure with icosahedral symmetry akin to the shapes sites. Furthermore, it is reasonable to expect that discussed in Section II A. (c) For small γ ≈ 60 and small the bending response of the vesicle shell is not lin- ◦ θc = 5 we find an irregular polyhedral vesicle with facets ear and significantly changes if the local deforma- being insensitive to the position of the defects. (d) As tion becomes large. We introduced a simple model FvK number if increased to γ ≈ 103 and the critical ◦ that takes into account such effects at the most ba- angle set to θc = 14 the structure is icoshahedral but sic level. Even with a such simplified model we were with sharp ridges emanating for each defect. (e) Finally, 3 ◦ able to show that a non-linearity in the bending re- at γ ≈ 2.5 · 10 and very small critical angle, θc = 2 the vesicle is no longer stable and crumples. Red points sponse can drastically affect the low energy shapes indicate positions of twelve five-fold defects. and lead to interesting faceted structures. Understanding the mechanisms that lead to the faceting of vesicles goes beyond the realm of an aca- effects to be appreciable. It is reasonable to ex- demic exercise and could have far-reaching impact pect such conditions to be met for vesicles in the on nano-technology. Flat faces, for example, could sub 100nm size range as are most of those for which be easier to functionalize or otherwise chemically faceting has been reported.3,6,8 treat. Biochemical reactions are expected to have different rates near flat compared to curved surfaces, which could be utilized to construct highly sensitive III. SUMMARY AND CONCLUSIONS transport agents. With the ability to accurately con- trol the shape of vesicles it could be possible to engi- We have shown that faceted structures can be low neer vesicles that would, e.g. selectively target only energy configurations of elastic vesicles as a result particular agents. We hope this work will broaden of several different mechanisms. Incompatibility be- the interest in this emergent subject of modern nano- 9 science. multicomponent vesicles as well as to R. Everaers This work was supported by NSF grant DMR- who suggested to analyze the effects of a critical 0808812. RS would like to thank M. Demers, angle bending. RS and MJB also thank the Soft C. Funkhouser, C. Leung, M. Olvera de la Cruz, Matter Program at Syracuse University for financial L. Palmer, B. Qiao, and G. Vernizzi for collabora- support. tions on various aspects of the study of faceting in

∗ Electronic address: [email protected] Chemical Society 129, 756 (2007). † Electronic address: [email protected] 22 C. A. Haselwandter and R. Phillips, Phys. Rev. Lett. 1 C. V. Iancu, D. M. Morris, Z. Dou, S. Heinhorst, G. C. 105, 228101 (2010). Cannon, and G. J. Jensen, Journal Of Molecular Bi- 23 F. Quemeneur, C. Quilliet, M. Faivre, A. Viallat, and ology 396, 105 (2010). B. Pépin-Donat, Physical Review Letters 108, 108303 2 A. Walsby, Nature 283, 69 (1980). (2012). 3 A. Blaurock and R. Gamble, Journal of Membrane 24 D. L. D. Caspar and A. Klug, Cold Spring Harbor Biology 50, 187 (1979). Symposia on Quantitative Biology 27, 1 (1962). 4 E. Marques, O. Regev, A. Khan, M. da Graca Miguel, 25 S. Sachdev and D. R. Nelson, J. Phys. C: Solid State and B. Lindman, Macromolecules 32, 6626 (1999). Phys. 17, 5473 (1984). 5 F. Antunes, R. Brito, E. Marques, B. Lindman, and 26 M. Bowick, D. Nelson, and A. Travesset, Physical Re- M. Miguel, The Journal of Physical Chemistry B 111, view B 62, 8738 (2000), ISSN 0163-1829. 116 (2007). 27 L. Landau and E. Lifshitz, Theory of Elasticity 6 M. A. Greenfield, L. C. Palmer, G. Vernizzi, (Butterworth-Heinemann, London, 1995), 3rd ed. M. Olvera de la Cruz, and S. I. Stupp, Journal of 28 H. S. Seung and D. R. Nelson, Physical Review A 38, the American Chemical Society 131, 12030 (2009). 1005 (1988). 7 M. Dubois, B. Deme, T. Gulik-Krzywicki, J.-C. 29 W. Koiter, A consistent first approximation in the Dedieu, C. Vautrin, S. Desert, E. Perez, and T. Zemb, general theory of thin elastic shells: Foundations and Nature 411, 672 (2001). linear theory, pt. 1 (Laboratorium voor Toegepaste 8 C. Leung, L. Palmer, B. Qiao, S. Kewalra- Mechanica der Technische Hogeschool, 1959). mani, R. Sknepnek, C. Newcomb, M. Greenfield, 30 J. Lidmar, L. Mirny, and D. R. Nelson, Physical Re- G. Vernizzi, S. Stupp, M. Bedzyk, et al., ACS Nano view E 68, 051910 (2003). 6, 10901 (2012). 31 W. Humphrey, A. Dalke, and K. Schulten, Journal of 9 L. Jia, A. Cao, D. Lévy, B. Xu, P. Albouy, X. Xing, Molecular Graphics 14, 33 (1996). M. Bowick, and M. Li, Soft Matter 5, 3446 (2009). 32 J. Stone, Master’s thesis, Computer Science Depart- 10 B. Xu, R. Piñol, M. Nono-Djamen, S. Pensec, ment, University of Missouri-Rolla (1998). P. Keller, P. Albouy, D. Lévy, and M. Li, Faraday 33 A. E. Lobkovsky, Phys. Rev. E 53, 3750 (1996). discussions 143, 235 (2009). 34 A. E. Lobkovsky and T. Witten, Physical Review E 11 L. Jia, D. Lévy, D. Durand, M. Impéror-Clerc, 55, 1577 (1997). A. Cao, and M. Li, Soft Matter 7, 7395 (2011). 35 A. Siber, Physical Review E 73, 061915 (2006). 12 L. Lamb, D. Huffman, R. Workman, S. Howells, 36 C. M. Funkhouser, R. Sknepnek, and M. Olvera de la T. Chen, D. Sarid, and R. Ziolo, Science 255, 1413 Cruz, Soft Matter 9, 60 (2012). (1992). 37 B. Audoly and Y. Pomeau, Elasticity and Geometry - 13 M. Hu, J. Chen, M. Marquez, Y. Xia, and G. V. Hart- From Hair Curls to the Non-linear Response of Shells land, Journal of Physical Chemistry C 111, 12558 (Oxford University Press, Oxford, UK, 2010). (2007). 38 G. Vernizzi, R. Sknepnek, and M. de la Cruz, Pro- 14 B. Alberts, A. Johnson, J. Lewis, M. Raff, K. Roberts, ceedings of the National Academy of Sciences 108, and P. Walter, Molecular Biology of the Cell (Garland 4292 (2011). Science, New York, 2002). 39 K. Brakke, Experimental Mathematics 1, 141 (1992). 15 A. Bangham and R. Horne, Journal of molecular bi- 40 G. Gompper and D. Kroll, Journal De Physique I 6, ology 8, 660 (1964). 1305 (1996), ISSN 1155-4304. 16 S. Veatch and S. Keller, Biophysical Journal 85, 3074 41 M. Meyer, M. Desbrun, P. Schröder, and A. H. Barr, (2003). Visualization and mathematics 3, 34 (2002). 17 C. Nardin, T. Hirt, J. Leukel, and W. Meier, Lang- 42 E. Magid, O. Soldea, and E. Rivlin, Computer Vision muir 16, 1035 (2000). and Image Understanding 107, 139 (2007). 18 W. Helfrich, Zeitschrift Fur Naturforschung C-a Jour- 43 B. Schmidt and F. Fraternali, Journal of the Mechan- nal Of Biosciences C 28, 693 (1973), ISSN 0939-5075. ics and Physics of Solids 60, 172 (2012). 19 S. Munro, Cell 115, 377 (2003). 44 D. Nelson, T. Piran, and S. Weinberg, Statistical me- 20 T. A. Witten, Rev. Mod. Phys. 79, 643 (2007). chanics of membranes and surfaces (World Scientific 21 A. González-Pérez, M. Schmutz, G. Waton, Singapore, 2004), 2nd ed. M. Romero, and M. Krafft, Journal of the American 45 E. Katifori, S. Alben, E. Cerda, D. R. Nelson, and 10

J. Dumais, Proceedings of the National Academy of D. K. Schwartz, The Journal of Physical Chemistry Sciences 107, 7635 (2010). B 114, 8616 (2010). 46 A. Green and W. Zerna, Theoretical elasticity (Dover 53 J. Hu, T. Weikl, and R. Lipowsky, Soft Matter 7, 6092 Pubns, 2002). (2011). 47 E. Efrati, E. Sharon, and R. Kupferman, Journal of 54 M. F. Demers, R. Sknepnek, and M. O. de la Cruz, the Mechanics and Physics of Solids 57, 762 (2009). Physical Review E 86, 021504 (2012). 48 R. Sknepnek and M. Olvera de la Cruz, Physical Re- 55 B. Sweet and M. Hilleman, in Proceedings of the So- view E 85, 050501 (2012). ciety for Experimental Biology and Medicine. Society 49 R. Sknepnek, G. Vernizzi, and M. Olvera de la Cruz, for Experimental Biology and Medicine (New York, Soft Matter 8, 636 (2012). NY) (Royal Society of Medicine, 1960), vol. 105, pp. 50 S. Trabelsi, S. Zhang, T. Lee, and D. Schwartz, Phys- 420–427. ical review letters 100, 37802 (2008). 56 B. La Scola, S. Audic, C. Robert, L. Jungang, 51 S. Trabelsi, Z. Zhang, S. Zhang, T. R. Lee, and D. K. X. de Lamballerie, M. Drancourt, R. Birtles, J.-M. Schwartz, Langmuir 25, 8056 (2009). Claverie, and D. Raoult, Science 299, 2033 (2003). 52 S. M. Malone, S. Trabelsi, S. Zhang, T. R. Lee, and