Multispherical shapes of vesicles with intramembrane domains

Abstract Phase separation of biomembranes into two fluid phases, a and b, leads to the formation of vesicles with intramembrane a- and b-domains. These vesicles can attain multispherical shapes consisting of several spheres connected by closed membrane necks. Here, we study the morphological complexity of these multispheres using the theory of curvature elasticity. Vesicles with two domains form two-sphere shapes, consisting of one a- and one b-sphere, connected by a closed ab-neck. The necks’ effective mean curvature is used to distinguish positive from negative necks. Two-sphere shapes of two-domain vesicles can attain four different morphologies that are governed by two different stability conditions. The closed ab-necks are compressed by constriction forces which induce neck fission and vesicle division for large line tensions and/or large spontaneous curvatures. Multispherical shapes with one ab-neck and additional aa- and bb-necks involve several stability conditions, which act to reduce the stability regimes of the multispheres. Furthermore, vesicles with more than two domains form multispheres with more than one ab-neck. The multispherical shapes described here represent generalized constant-mean-curvature surfaces with up to four constant mean curvatures. These shapes are accessible to experimental studies using available methods for giant vesicles prepared from ternary lipid mixtures. Graphic abstract


Introduction
Biological and biomimetic membranes represent twodimensional liquids.Biological membranes contain a large assortment of lipids and membrane proteins, whereas biomimetic membranes typically consist of a few lipid and protein components.These membranes should be able to undergo phase separation into two types of fluid domains, in close analogy to phase separation of liquid mixtures in three dimensions.This conclusion seems quite obvious from a theoretical point of view but, at the beginning of the 1990 s, it was rather difficult to find experimental evidence for it [1].
a e-mail: lipowsky@mpikg.mpg.de(corresponding author) branes have been determined using spectroscopic methods [16] as well as fluorescence microscopy of giant vesicles and X-ray diffraction of membrane stacks [9,[17][18][19].Fluid-fluid coexistence has even been found in giant plasma membrane vesicles that contain a wide assortment of different lipids and proteins [20,21].
Direct evidence for the formation of two types of fluid domains was provided by single particle tracking that showed that both phases exhibit relatively fast lateral diffusion [2].In addition, using GUVs, several theoretical predictions [22][23][24] could be directly confirmed: the growth and coalescence of small domains into larger ones; domain-induced budding; and small shifts of the domain boundary away from the waist-line of the membrane neck.
The three examples in Fig. 1 display vesicles with two intramembrane domains which are labeled by two different fluorophores.In all three examples, the boundary between the two domains forms an open membrane neck, which prefers to close when the volume of the vesicles is further reduced by osmotic deflation.Closed membrane necks can lead to a variety of multispherical shapes as observed for giant vesicles with laterally uniform membranes [25,26].In the latter case, the multispheres are built up from spheres with up to two different curvature radii, corresponding to large and small spheres, which are connected by closed membrane necks.Some examples for such multispherical shapes Fig. 1 Domain-induced budding of giant unilamellar vesicles (GUVs) as predicted by theory [22,23] and observed by fluorescence microscopy [4,7,10].The two intramembrane domains consist of liquid-disordered (Ld) and liquid-ordered (Lo) lipid phases: a Cross-section through a vesicle that formed two lipid phase domains after a decrease in temperature.Reprinted with permission from Ref. [4]  are displayed in Fig. 2. Each multisphere consists of large and small spheres but exhibits only two different radii, one for the large and one for the small spheres.More precisely, each large and small sphere is actually a punctured sphere that is connected to the punctures of neighboring spheres via closed membrane necks.
Here, the observed behavior of two-domain vesicles as illustrated in Fig. 1 and of multispherical vesicle shapes as found for uniform membranes, see Fig. 2, will be used to elucidate the morphological complexity of multispheres formed by vesicles with two or more intramembrane domains.The analysis is based on the theory of curvature elasticity.We consider different curvatureelastic properties of the two membrane domains apart from their Gaussian curvature moduli, which are taken to have identical values in the two domains.This simplifying assumption has several advantages.First, the domain boundaries between the intramembrane domains are located within the closed membrane necks [23].Second, the multispherical shapes can be obtained by elementary calculus, without the need to use numerical methods for their computation.Third, identical Gaussian curvature moduli do not affect the vesicle shapes, which then depend on five membrane-elastic parameters as provided by two spontaneous curvatures and two bending rigidities as well as the line tension.Because the two bending rigidities are usually of the same order of magnitude, the morphological complexity of the multispheres depends primarily on the two spontaneous curvatures and on the line tension.At the end, we look at the changes arising from different Gaussian curvature moduli for the two domains.The most important change is related to the constriction forces at closed membrane necks as discussed in Sect.10.2.
The paper is organized as follows: Section 2 provides a brief summary of multispherical shapes formed by uniform membranes and introduces the notion of positive and negative membrane necks.In Sect.3, we will look at two-domain vesicles that form two-sphere shapes with a single ab-neck.Depending on the sign of the ab-neck and on the relative size of the aand b-sphere, four two-sphere morphologies will be distinguished.The stability of these morphologies is governed by two stability relations, which apply to positive and negative abnecks, respectively.Each closed ab-neck is subject to a constriction force that acts to compress the neck as described in Sect. 4. This constriction force depends primarily on the line tension of the domain boundary and on the spontaneous curvatures of the aand b-domain.Large line tensions and/or large spontaneous curvatures generate constriction forces that drive the fission of the closed necks, whereas smaller line tensions and moderate spontaneous curvatures are unlikely to induce such a fission process.
The two-sphere shapes formed by two-domain vesicles are the simplest examples for multispheres of multi-domain vesicles.More complex morphologies are described in Sect.5, corresponding to two-domain vesicles with multispheres formed by individual domains and to vesicles with more than two domains, which can transform into multispheres with more than one ab-neck.The stability regimes for multispheres with one ab-neck are determined in Sect.6 and nested multispheres arising from nested domains are described in Sect.7. The last two Sects.8 and 9 interpret multispheres as generalized constant-mean-curvature surfaces and show how available methods for the experimental study of GUVs can be applied to multispheres Fig. 2 Multispherical shapes of GUVs with laterally uniform membranes: a Two-sphere shape consisting of one large and one small sphere, forming an out-bud; b Two-sphere shape with one large and one small sphere, forming an inbud; c multisphere with one large sphere and a linear chain of six outward-pointing small spheres; d Multisphere with one large sphere and a linear chain of six inward-pointing small spheres; e Multisphere consisting of two large and one small spheres; f Multisphere with two large and two small spheres; and g Multisphere consisting of 24 equally sized spheres.In (a), (c), (e)-(g), all spheres have a positive mean curvature which implies that the membranes have a positive spontaneous curvature.In (b) and (d), the mean curvature of the small spheres is negative caused by a negative spontaneous curvature.All scale bars are 10 μm [25] with intramembrane domains.The changes arising from different Gaussian curvature moduli are described in Sect.10.

Multispherical shapes of uniform membranes
This section contains a brief review of the multispherical shapes as formed by uniform membranes, which are characterized by a uniform molecular composition and thus by uniform membrane-elastic parameters.A more detailed discussion of these shapes can be found in Ref. [26].

Basic aspects of multispherical shapes
Each multisphere as displayed in Fig. 2 involves only a single fluid membrane, which encloses both the spherical compartments and the membrane necks connecting the spheres.Thus, each sphere is actually a punctured sphere, with its punctures being connected to the punctures of neighboring spheres via closed membrane necks.If we added a fluorescent probe to one spherical membrane segment, the probe would diffuse across the membrane necks and eventually spread over the whole multispherical,membrane.Likewise, the closed necks may undergo shape fluctuations, which lead to short-lived open necks, which transiently provide narrow water channels between the adjacent aqueous compartments.Inspection of Fig. 2 reveals that each multisphere involves large and small spheres with only up to two different curvature radii, R l and R s .These radii are intimately related to the mean curvatures M l and M s of the large and small spheres.As explained in the next subsection, the coexistence of two different sphere sizes on the same multispherical shape is a direct consequence of the shape equation for spherical membrane segments and implies that all spheres are subject to the same membrane tension and, thus, formed by a single membrane.
The multispherical shapes in Fig. 2 were experimentally observed to remain unchanged for many hours.This stability is primarily determined by the stability of the closed necks against neck opening.Stably closed necks require sufficiently large spontaneous curvatures of the vesicle membrane.As described in the next but one subsection, two stability conditions must be distinguished depending on the signs of the mean curvatures M l and M s .The mean curvature M l of the large spheres is always positive, but the mean curvature M s of the small spheres can be positive or negative.Examples for small spheres with negative mean curvature, corresponding to inverted spheres, are provided by the small spheres in Fig. 2b, d.

Local shape equation for uniform membranes
First, let us consider a membrane with uniform molecular composition that can be characterized by uniform spontaneous curvature m and uniform bending rigidity κ.When such a membrane forms a spherical segment, this segment attains a constant mean curvature M that satisfies the local shape (or Euler-Lagrange) equation [27] in which the pressure difference between the interior and exterior aqueous solution is balanced by the linear term proportional to the total membrane tension and by a second term, which is quadratic in the mean curvature M .Here, Σ is the mechanical tension acting within the membrane and 2κm 2 is the spontaneous tension [28] arising from the spontaneous curvature m.Alternatively, the two parameters ΔP and Σ can be viewed as two Lagrange multipliers used to minimize the bending energy for certain, prescribed values of the vesicle volume V and the membrane area A. For such a constrained minimization, ΔP and Σ represent auxiliary variables that are conjugate to the geometric variables V and A. When we consider vesicles with a certain volume V and a certain membrane area A, the shape functional for these vesicles depends on the bending rigidity κ and the spontaneous curvature m as well as on the two geometric parameters V and A. Using the bending rigidity as the basic energy scale and the vesicle size R ve = A/(4π) as the basic length scale, the vesicle shapes are found to depend only on two dimensionless shape parameters, the volume-to-area ratio (or reduced volume) v which is proportional to V /A 3/2 and the rescaled spontaneous curvature m = mR ve [29].

Case-by-case analysis of mean curvature
For zero spontaneous curvature, m = 0, the local shape equation in (1) reduces to ΔP = 2 ΣM , which has the same form as the classical Young-Laplace equation for liquid droplets.In this special case, the shape equation has the single solution or root for the mean curvature M of the spherical segment.For nonzero spontaneous curvature, m = 0, the local shape equation in (1) can be rewritten in the form with the two parameter combinations Inspection of Eq. ( 5) directly shows that this equation has no (real-valued) solution or root for δ − σ 2 > 0 ; (7) one degenerate (double) root as given by M = σ for δ − σ 2 = 0 ; (8) and two different roots and for the parameter range

Multispherical architectures for uniform membranes
In principle, the two parameter combinations σ and δ as defined in Eq. ( 6) can be positive or negative, depending in particular on the sign of the spontaneous curvature m.A detailed analysis as described in Ref. [26] reveals, however, that physically meaningful solutions M + and M − are only obtained for two cases, I and II.Case I is characterized by In this case, the multispheres consist of large and small spheres, both of which have positive mean curvature.Furthermore, the two radii R l and R s of the large and small spheres are given by Here and below, all radii are taken to be positive.Examples for Case I are provided by panels a, c, e, and f of Fig. 2. In the last panel g of this figure, we see an example for many equally sized spheres, corresponding to the doubly degenerate root in Eq. ( 8).
On the other hand, Case II is given by corresponding to one large sphere with positive mean curvature M + and inverted small spheres with negative mean curvature M − .For case II, the curvature radii R l and R s of the large and small spheres have the form Examples for case II are shown in panels b and d of Fig. 2. For both cases I and II, the formation of a multispherical shape provides direct evidence that all spherical membrane segments experience the same mechanical tension Σ and that the whole multisphere is formed by a single bilayer membrane.

Closure of open membrane necks
The second ingredient from curvature elasticity that is necessary to understand multispherical shapes is the formation of closed membrane necks.Two-sphere shapes of uniform vesicle membranes were originally obtained as limit shapes of smoothly curved shapes with open necks [29,30], using numerical methods applied to curvature models.These models describe the membranes as elastic surfaces, governed by certain curvature-elastic parameters.Uniform membranes as considered in the present section are characterized by curvature-elastic parameters, which are laterally uniform along the whole membrane, reflecting the uniform molecular composition of the membrane.For axisymmetric shapes, the minimization of the shape functional leads to a set of ordinary differential equations [29].The solutions of these equations form a discrete set of energy branches.Along each of these branches, the vesicle shape evolves smoothly as we vary one of the model parameters until we encounter a limit shape that can no longer be obtained by solving the differential equations.The two-sphere shapes considered here represent such limit shapes, which involve kinks of the membrane contours at the membrane necks and discontinuities of the mean curvature across this neck.
Furthermore, the geometry of a multispherical shape does not depend on the spontaneous curvature but only on the volume-to-area ratio v as well as on the number of large and small spheres [26].As a consequence, the limit shapes continue to exist when the energy branches are further continued, keeping the multispherical geometry fixed but changing a single curvatureelastic parameter such as the spontaneous curvature.Even though the vesicle shape remains unchanged along this continuation, the bending energy of the vesicle changes because this energy depends on the curvatureelastic parameters.
When the vesicle forms an axisymmetric shape with an open neck, this neck has a finite radius R ne , which represents the radius of the waist-line around the neck.When the neck closes, the radius R ne goes to zero which implies that the second principal curvature C 2,wl = 1/R ne parallel to the waist-line diverges.However, the mean curvature M remains finite on both sides of the neck.Therefore, the divergence of the second principal curvature must be canceled by another divergence arising from the first principal curvature C 1,wl , which is equal to the contour curvature.Furthermore, as the neck becomes closed, the mean curvature of the membrane attains two finite but different values on the two sides of the neck which implies that the mean curvature develops a discontinuity across the closed neck.

Stability of closed membrane necks
Each closed neck provides a connection between two spherical membrane segments i and j with mean curvatures M i and M j .The stability of such a closed neck is governed by a stability condition that involves the spontaneous curvature m of the adjacent membrane segments and the effective mean curvature of the closed neck as defined by [31] Note that the neck curvature M eff ij represents a purely geometric quantity.When the large and small spheres can be resolved by optical microscopy as in Fig. 2, the neck curvature M eff ij can be directly deduced from the optical images.Therefore, this curvature represents an observable quantity.
The form of the stability condition depends on the sign of the effective neck curvature M eff ij .For positive neck curvature M eff ij > 0, the stability condition is given by: which can only be fulfilled for a sufficiently large and positive spontaneous curvature m.For the multispheres displayed in panels a, c, e, and f of Fig. 2, all closed membrane necks have positive neck curvatures M eff ij > 0. Furthermore, if the multisphere consists of a chain of equally sized spheres as in panel g of Fig. 2, all membrane necks have the same neck curvature, which is positive as well.
For negative neck curvature M eff ij < 0, the stability condition has the form which requires a sufficiently large and negative spontaneous curvature m.For the multispheres displayed in panels b and d of Fig. 2, all closed membrane necks have negative neck curvatures M eff ij < 0. The stability conditions for a closed membrane neck as given by Eqs.(17) and (18) are local in the sense that they depend only on the geometry and on the spontaneous curvature of the two membrane segments adjacent to the membrane neck.In particular, these stability conditions do not depend on the global morphology of the vesicle as characterized by its volume and surface area or by the number of large and small spheres formed by the vesicle [26].

Positive and negative membrane necks
It will be convenient to characterize the membrane necks by the sign of their effective neck curvature and to distinguish positive from negative necks.By definition, a "positive neck" has a positive effective mean curvature M eff ij > 0, whereas a "negative neck" has a negative effective mean curvature M eff ij < 0. Thus, the multispheres shown in Fig. 2 involve only positive necks apart from those in panels b and d, which involve only negative necks.Using these definitions, we obtain an alternative characterization of the two cases I and II distinguished in Sect.2.2.2.Indeed, multispheres belonging to case I have only positive membrane necks, whereas multispheres belonging to case II have only negative necks.
Inspection of the different examples in Fig. 2 shows that positive membrane necks connect two interior subcompartments, whereas negative membrane necks connect two exterior subcompartments.Therefore, positive necks can be regarded as interior necks and negative necks as exterior necks [26].In the following, we will focus on the distinction between positive and negative necks and will only occasionally refer to the equivalent distinction between interior and exterior necks.

Two-sphere shapes of two-domain vesicles
In this section, we go back to Fig. 1, which displays several examples of giant vesicles with two intramembrane domains, visualized by different fluorophores.The two domains are now distinguished by the domain labels a and b.In Fig. 1, the budding process is incomplete in the sense that each two-domain vesicle assumes a dumbbell shape with an open neck.Furthermore, in each example, the domain boundary between the a and b domains is located within this open neck.In order to close the neck, we now imagine to reduce the vesicle volume, which can be achieved experimentally by osmotic deflation.As a result of this deflation process, we obtain a two-sphere shape consisting of an a-sphere and a b-sphere, which are connected by a closed ab-neck.
In the following subsections, we will first demonstrate that the geometry of two-sphere shapes formed by twodomain vesicles is completely determined by the area fractions of the two domains.Second, we will examine the stability of the closed ab-neck and determine the stability and instability regimes.These regimes will be visualized by morphology diagrams, which are defined in terms of the spontaneous curvatures of the aand b-domains.

Basic geometry of two-domain vesicles
The geometry of a single vesicle with two domains is determined by the vesicle volume V , the surface area A of its membrane, and the area fractions of the two domains.The vesicle size R ve is defined in terms of the membrane area A and given by R ve ≡ A/(4π) , (19) which represents the radius of a sphere with area A and is taken to provide the basic length scale of the vesicles.Likewise, the rescaled vesicle volume has the form with 0 ≤ v ≤ 1 where the limiting value v = 1 corresponds to a spherical shape of the vesicle.Now, consider a vesicle as in Fig. 1 with one a-domain and one b-domain with surface areas A a and A b .The total surface area A of the vesicle membrane is given by: and the area fractions Φ a and Φ b of the two domains are defined by with Φ a + Φ b = 1.

Geometry of two-sphere shapes with two domains
Two-sphere shapes consisting of one a-sphere and one b-sphere are the simplest multispherical shapes that can be formed by vesicles with two domains, with the domain boundary being located within the closed membrane neck between the two spheres.The geometry of such two-sphere shapes depends on the radius R a of the a-sphere and the radius R b of the b sphere.As before, all radii are taken to be positive.

Radii and mean curvatures of two-sphere shapes
In general, the a-sphere may be larger than the bsphere or vice versa as illustrated in Fig. 3.In addition, both spheres may have a positive mean curvature as in Fig. 3a, b or the smaller sphere may have a negative mean curvature as in Fig. 3c, d.Indeed, for the examples in Fig. 3a, b, the mean curvatures M a and M b of the aand b-sphere are both positive and given by: corresponding to out-budded two-sphere vesicles.On the other hand, for the example in Fig. 3c, which represents a two-sphere vesicle with an in-bud formed by the b-domain, these mean curvatures have the values (24) whereas they are equal to for the example in Fig. 3d, which displays a two-sphere vesicle with an in-bud formed by the a-domain.

Positive and negative ab-necks
Generalizing the definition for uniform membranes as given by Eq. ( 16), the effective mean curvature of the ab-necks is taken to be where M a and M b are the mean curvature of the aand b-sphere adjacent to the neck.Using this definition, the out-budded two-sphere shapes in Fig. 3a, b have a positive ab-neck with neck curvature as follows from Eq. (23).In contrast, the in-budded twosphere shape in Fig. 3c with R a > R b has a negative ab-neck with Likewise, the shape in Fig. 3d with R b > R a involves a negative ab-neck as well with effective neck curvature These effective neck curvatures will be useful to classify the different patterns of multispherical shapes as discussed further below.
For the out-budded two-spheres, the positive ab-neck provides a closed channel between two interior subcompartments.For the in-budded two-spheres, the negative ab-neck represents a closed channel between two exterior subcompartments.Thus, positive and negative ab-necks can again be regarded as interior and exterior necks in the sense, that interior ab-necks provide a connection between two interior subcompartments whereas exterior ab-necks connect two exterior subcompartments.

Two-sphere geometry determined by area fractions
In terms of the surface areas A a and A b of the two domains, the radii of the aand b-sphere are given by: R a = A a /(4π) and R b = A b /(4π) .(30) To simplify the mathematical formula, it will be convenient to define the rescaled radii with the vesicle size R ve = A/(4π).For the twosphere shapes formed by a two-domain vesicle as considered here, the rescaled radii become and Furthermore, the rescaled and dimensionless mean curvatures Ma ≡ M a R ve and Mb ≡ M b R ve (34) are now given by where the plus and minus signs are determined by Eqs. ( 23)-( 25), corresponding to the different twosphere morphologies in Fig. 3.The area decomposition in Eq. ( 21) now attains the simple form which applies to both out-budded and in-budded twosphere shapes.As far as the rescaled volume v is concerned, we have to distinguish three cases.For outbudded two-sphere vesicles as in Fig. 3a, b, the rescaled volume is given by with Φ a = 1 − Φ b .For in-budded two-sphere vesicles with the in-bud formed by the b-domain (Fig. 3c), the rescaled volume is Finally, when the in-bud is formed by the a-domain (Fig. 3d), the two-sphere vesicle has the rescaled volume: Thus, all geometric properties of the two-sphere vesicles with one a-sphere and one b-sphere can be expressed in terms of the area fractions Φ b and Φ a = 1 − Φ b .
In order to illustrate the formation and characterization of two-sphere vesicles, we consider the examples in Figs. 4 and 5.We start from spherical vesicles with rescaled volume v = 1 and different area fractions Φ b .The vesicles are then exposed to an increased osmotic pressure in the exterior compartment, which acts to reduce the vesicle volume by osmotic deflation, a standard experimental procedure.Likewise, osmotic inflation can be applied to increase the vesicle volume.As a result of the deflation, the spherical vesicles may transform into out-budded two-sphere vesicles, for which both mean curvatures M a and M b are positive as in Fig. 4, or into two-sphere vesicles with an in-budded b-domain as in Fig. 5.The different cases of two-sphere vesicles with (i) M a > 0 and M b > 0, (ii) M a > 0 and M b < 0, as well as (iii) M a < 0 and M b > 0, see Fig. 3, can be distinguished by different stability conditions for the closed ab-necks as described after the next subsection.

Curvature discontinuities at domain boundary
The second ingredient from curvature elasticity that is necessary to understand the formation and stability of two-sphere vesicles is the stability of the closed ab-necks.The corresponding stability conditions for these necks are more involved than for uniform membranes.In fact, even for open necks, axisymmetric vesicles with two domains exhibit curvature discontinuities at the domain boundaries.These discontinuities can be computed explicitly for axisymmetric shapes parametrized by arc length s.The curvature discontinuities follow from the matching conditions for the mean curvatures M a (s db ) and M b (s db ) at the aand b-sides of the domain boundary, which is located at arc length s = s db .
These matching conditions are obtained from the first variation of the shape functional F as given by Eq. (A7) in Appendix A. The shape functional depends on the bending rigidity κ a and the spontaneous curvature m a of the a-domain as well as on the bending rigidity κ b and the spontaneous curvature m b of the b-domain.If we allowed the aand b-domains to have different Gaussian curvature moduli κ Ga and κ Gb , the first variation of the shape functional would lead to the matching condition [31] where C 2 (s db ) is the second principal curvature parallel to the domain boundary, which is continuous across this boundary.
For an axisymmetric dumbbell shape with an open neck, the second principal curvature C 2 (s db ) is directly related to the neck radius R ne via C 2 (s db ) = 1/R ne .Therefore, this second principal curvature diverges if the neck radius vanishes.In order to avoid this divergence, the domain boundary moves away from the waist-line of the open neck during the neck closure process as shown by numerical calculations [23].On the other hand, when the two Gaussian curvature moduli κ Ga and κ Gb have the same value, the matching condition in Eq. ( 40) simplifies and becomes which is equivalent to the mean curvature discontinuity at the domain boundary.Thus, in contrast to the smoothly curved dumbbells formed by uniform membranes, the dumbbell shape of a two-domain vesicle exhibits a mean curvature discontinuity at the domain boundary as given by Eq. ( 42), even for κ Ga = κ Gb , that is, when both domains have the same Gaussian curvature modulus.Therefore, one should expect that the stability condition for a closed ab-neck is more complex than the corresponding condition for uniform membranes as shown in the next subsection.

Stability of closed ab necks
The stability of closed ab-necks with respect to neck opening depends on the curvature-elastic parameters of the two membrane domains as provided by the spontaneous curvatures m a and m b as well as the bending rigidities κ a and κ b of the two domains.In addition, the stability of a closed ab-neck also depends on the line tension λ of the domain boundary between the aand b-domain.
In this subsection, we describe the stability conditions for the ab-necks of the different types of twosphere shapes displayed in Fig. 3.The form of these conditions is somewhat different for the out-budded two-spheres in Fig. 3a, b, for the in-budded b-domains in Fig. 3c, and for the in-budded a-domains in Fig. 3d.These conditions can be visualized in terms of morphology diagrams that depend on the spontaneous curvatures m a and m b of the two membrane domains.

Neck stability for out-budded domains
The stability conditions for closed ab-necks can be obtained by looking at dumbbell shapes with slightly open necks and parametrizing these shapes by piecewise constant-mean-curvature surfaces.For out-budded two-sphere vesicles, such a parametrization was first considered in Ref. [23] generalizing an analogous parametrization for uniform membranes in Ref. [32].In this parametrization, one considers two hemispheres connected by an intermediate unduloid segment with neck radius R ne .In the limit of small neck radius, the bending energy of the out-budded dumbbell shape behaves as: up to first order in the neck radius R ne , with the mean curvatures M a and M b of the aand b-sphere. 1The closed neck with R ne = 0 is stable if the bending energy This closed neck condition applies to both Fig. 3a, b, that is, to both a larger a-sphere with R a > R b and to a larger b-sphere with R b > R a .A simple cross-check of the closed neck condition in Eq. ( 44) is obtained when we look at the limiting case of two identical domains with κ a = κ b , m a = m b , and λ = 0.In this limit, Eq. ( 44) reduces to 2m ≥ M a + M b , the correct closed neck condition for uniform membranes as in Eq. ( 17)

Stability condition in terms of rescaled variables
The stability condition in Eq. ( 44 In terms of these rescaled variables, the stability condition in Eq. ( 44) becomes This closed neck condition applies to both panels a and b of Fig. 3, that is, to 0 < Φ b ≤ 1/2 as in Fig. 3a and to 1/2 ≤ Φ b < 1 as in Fig. 3b.The line of limit shapes L ab is now described by the equality For the two-sphere vesicles discussed in the present section, the rescaled radii r a and r b can be expressed in terms of the area fractions Φ a and Φ b which leads to r a = √ Φ a and r b = √ Φ b , see Eqs. (32) and (33).
To visualize the stability regime for the closed abnecks, it is convenient to rename the rescaled spontaneous curvatures and to define the coordinates x ≡ ma and y ≡ mb (48) for the two-dimensional morphology diagrams in Fig. 6.When Eq. ( 47) is solved for mb = y, the line of limit shapes L ab is described by the linear relation and the intercept value The asymptotic equality (≈) in Eq. ( 50) applies to giant vesicles with a large vesicle size R ve κ b /λ.In the (x, y)-plane, the line of limit shapes L ab as given by Eq. ( 49) is a straight line with negative slope dy/dx = d mb /d ma = −κ a /κ b , which intersects the yaxis at the intercept value y ab , see Fig. 6a.Likewise, the L ab -line intersects the x-axis at the intercept value where the asymptotic equality again applies to giant vesicles with large size R ve , which is implicitly assumed in Fig. 6a.For a given value of the area fraction Φ b = r 2 b , the line of limit shapes L ab divides the (x, y)-plane into two parameter regimes corresponding to two-sphere vesicles with closed and with open ab-necks.As shown in Fig. 6a, the positive ab-neck is stably closed for y ≥ h out (x) or mb ≥ h out ( ma ) , (52) but opens up for y < h out (x) or mb < h out ( ma ), with the linear function h out (x) defined by Eqs. ( 49) and (50).The neck opens up in a continuous manner, that is, the neck radius R ne increases continuously from R ne = 0 in the yellow stability regime above the L ab -line in Fig. 6a to a nonzero value below this line.

Neck stability for in-budded domains
For two-sphere vesicles with in-budded b-domains as in Fig. 5b, the membrane shapes can again be parametrized by smoothly curved surface segments with piece-wise constant mean curvatures.In the limit of small neck radius R ne , the bending energy of the in-budded shape then behaves as [31] 27); and b Stability regime for negative ab-necks of two-sphere vesicles with in-budded b-domains (blue), corresponding to positive mean curvature Ma > 0 of the a-sphere, negative mean curvature M b < 0 of the b-sphere, and negative effec-tive mean curvature M eff ab < 0 of the ab-neck as in Eq. ( 28).In (a), all vesicles have constant volume v = 0.728 as in Fig. 4b.In (b), all vesicles have constant volume v = 0.296 as in Fig. 5b.In both panels, the line of limit shapes L ab (purple) separates vesicle shapes with stably closed necks from those with open necks.The intercepts of these L ablines with the coordinate axes are denoted by y ab = mb,ab and x ab = ma,ab .For out-budded shapes as in (a), these intercepts are given by Eqs. ( 50) and ( 51); for in-budded shapes as in (b), they are provided by Eqs. ( 59) and (60) further below up to first order in the neck radius R ne .For two-sphere vesicles with an in-budded b-domain, the mean curvatures M a and M b of the aand b-sphere are equal to M a = 1/R a and M b = −1/R b .Compared to the bending energy of the out-budded shape, see Eq. ( 43), the bending energy of the two-sphere vesicle with an in-budded b-domain as given by Eq. ( 53) involves two changes of sign.First, the mean curvature M b = −1/R b is now negative, whereas M b = +1/R b for out-budded shapes.In addition, the whole curvature-elastic term, which depends on the bending rigidities κ a and κ b , is negative in Eq. ( 53), whereas it is positive in Eq. (43).The form of the closed neck condition in Eq. ( 53) does not change when we swap the domain labels a and b which implies that this closed neck condition also applies for in-budded a-domains.
The ab-neck of an in-budded b-domain is stably closed, if the bending energy E be (R ne ) as given by Eq. ( 53) increases with increasing neck radius R ne .Therefore, the closed ab-neck of an in-budded b-domain is stable if A simple cross-check of this stability criterion is obtained for two identical domains with κ a = κ b , m a = m b , and λ = 0.In this case, Eq. ( 55) reduces to 2m ≤ M a + M b , the correct stability condition for uniform membranes as given by Eq. ( 18) with M i = M a and M j = M b .Neck stability for in-budded b-domains For an inbudded b-domain, the rescaled mean curvatures are given by Ma = M a R ve = +1/r a and Mb = M b R ve = −1/r b .In terms of these rescaled curvatures, the stability condition in Eq. ( 55) becomes In addition, an in-budded b-domain is only possible if the radius r b of the b-sphere does not exceed the radius r a of the a-sphere.
When the inequality in Eq. ( 56) becomes an equality, we obtain the line of limit shapes L ab for in-budded bdomains.Thus, for such b-domains, the line of limit shapes L ab is now given by For the two-sphere vesicles discussed in this section, the rescaled radii r a and r b are related to the area fractions Φ a and Φ b via r a = √ Φ a and r b = √ Φ b .Using the previously introduced coordinates x = ma and y = mb , the L ab -line is described by with the intercept value for the intersection of the L ab -line with the y-axis, see the morphology diagram in Fig. 6b.Likewise, the L abline intersects the x-axis at the intercept value For giant vesicles with a large value of R ve , both intercepts x ab = mab a and y ab = mab b become large and positive as described by the asymptotic equalities in Eqs. ( 59) and (60), which is implicitly assumed in Fig. 6b.As shown in this figure, the ab-neck of an inbudded b-domain is stably closed for but opens up for y > h b-in (x) or mb > h b-in ( ma ) with the linear function h b-in (x) defined in Eq. ( 58).The neck opens up in a continuous manner, that is, the neck radius R ne increases continuously from R ne = 0 in the yellow stability regime below the L ab -line in Fig. 6b to a nonzero value above this line.
For two-sphere vesicles with in-budded a-domains as in Fig. 3d, the rescaled mean curvatures are equal to Ma = −1/r a and Mb = +1/r b .Therefore, the relationships for in-budded a-domains can be obtained from Eqs. (56), (57), (59), and (60), which have been derived for in-budded b-domains, by replacing +1/r a by −1/r a as well as −1/r b by +1/r b in all of these equations.

Out-budded versus in-budded domains
The line of limit shapes L ab , which separates the stability regimes of the closed ab-necks from their instability regimes, corresponds to the purple lines in Fig. 6a, b.For two-sphere shapes with out-budded b-domains and positive ab-necks, the necks are stable above the purple L ab -line in Fig. 6a.For in-budded b-domains and negative ab-necks, the necks are stable below the purple L ab -line in Fig. 6b.For a fixed value of the area fraction Φ b , these two stability regimes exhibit a substantial overlap region, which is located between the purple line in Fig. 6a and the purple line in Fig. 6b, which are parallel to each other.Within this overlap region, the closed ab-necks are stable both for out-budded and for in-budded b-domains which implies the stability of both positive and negative ab-necks.
The two purple lines in Fig. 6a, b cross the x-axes at the two intercept values x ab as given by Eqs. ( 51) and (60).The difference between these two intercept values is: where the asymptotic equality applies to giant vesicles with large R ve -values.In such a situation, λR ve /κ a is of the order of 10 2 , which implies that the separation of the two lines of limit shapes L ab is quite large, leading to a broad overlap region.The overlap region includes the parameter values close to the origin of the (x, y)-plane, corresponding to small spontaneous curvatures ma and mb .Therefore, for small spontaneous curvatures, the ab-neck is stably closed both for out-budded and for in-budded bdomains.Consider, for instance, the two-domain vesicles in Figs.4b and 5b, corresponding to area fraction Φ b = 0.36, which display an out-budded and in-budded b-domain, respectively.Thus, we predict that both twosphere vesicles are stable for small spontaneous curvatures.Comparison of Figs.4b with 5b also shows that these two-sphere vesicles have a rather different volume as given by v = 0.728 for the out-budded b-domain in Fig. 4b and by v = 0.296 for the in-budded b-domain in Fig. 5b.Therefore, reducing the volume of a spherical two-domain vesicle as in Fig. 4a will first lead to a two-sphere vesicle as in Fig. 4b with an out-budded b-domain.Further reduction of the volume may then transform the out-budded b-domain into an in-budded one as shown in Fig. 5b.

Constriction forces and neck fission
The yellow stability regimes in Fig. 6 describe twosphere shapes for a fixed value Φ b = 0.36 of the bdomain's area fraction.As emphasized in Sect.3.2, such a fixed value of the area fraction completely determines the geometry of the two-sphere vesicle, provided we distinguish out-budded from in-budded shapes.Thus, when we move across the yellow stability regimes in Fig. 6 by varying the spontaneous curvatures x = ma and y = mb , we will always encounter the same twosphere shape.However, such variations in the spontaneous curvatures have another important consequence: they change the constriction force acting against the closed neck.This constriction force is defined by and represents the force acting against the closed neck.

Constriction force for out-budded two-spheres
The yellow stability regime for out-budded two-sphere vesicles with fixed area fraction Φ b = 0.36 and rescaled volume v = 0.728 is displayed in Fig. 6a.Thus, when we move within this stability regime by changing the spontaneous curvatures ma = x and mb = y of the two membrane domains, the shape of the two-sphere vesicle remains unchanged.On the other hand, using the form of the bending energy E be as given by Eq. ( 43), the constriction force f as defined by Eq. ( 63) becomes with M a = +1/R a and M b = +1/R b .This constriction force vanishes along the line of limit shapes L ab , as described by Eq. ( 47) and illustrated in Fig. 6a.The force is positive within the yellow stability regime above the L ab -line in Fig. 6a and increases with increasing line tension λ as well as with increasing excess curvatures m a −M a and m b −M b .The constriction force in Eq. ( 64) reduces to the particularly simple form when the mean curvatures of the two spheres are equal to the spontaneous curvatures, that is, for corresponding to out-budded two-sphere vesicles with vanishing bending energy E be as follows from Eqs. (A2) and (A3).

Constriction force for in-budded two-spheres
For in-budded b-domains, the limiting behavior of the bending energy E be for small neck radius R ne is given by Eq. (53).Using the definition ot the constriction force f in Equ (63), this force now becomes For in-budded a-domains, the constriction force f has the same form as in Eq. ( 67) but with M a = −1/R a and M b = +1/R b .For both types of in-budded two-sphere shapes, the constriction force attains the simple form f = 2πλ when the mean curvatures are equal to the spontaneous curvatures, that is, for or In both cases, the in-budded two-sphere vesicles have vanishing bending energy.

Fission of closed membrane necks
Sufficiently large constriction forces lead to the fission of closed membrane necks as observed experimentally for GUVs with uniform membranes [33].More precisely, the closed necks of the GUVs were cleaved when the constriction forces exceeded about 20 pN.
A similar threshold value for the constriction force is expected to apply to the two-domain vesicles considered here.Indeed, both for uniform and for two-domain membranes, the constriction force has to overcome an energy barrier provided by the formation of two ringlike bilayer edges across the closed membrane neck [27,31].
For both out-budded and in-budded two-sphere vesicles, the constriction force includes the line tension term 2πλ, see Eqs. ( 64)-(67).The line tension λ is equal to the excess free energy of the domain boundary per unit length.When the domain extends across both leaflets of the lipid bilayer, the domain boundary represents a cut through the whole bilayer.The crosssection of such a cut consists of three distinct regions: two hydrophilic headgroup regions with a combined thickness of about 1 nm and an intermediate hydrophobic tail region with a thickness of about 3 nm.For 3dimensional fluid phases, a typical value for the interfacial free energy is of the order of 10 mN/m.If one assumes that this value is also applicable to the headgroup region of the lipid bilayer and that the latter region gives the main contribution to the line tension, one obtains the rough estimate λ 10 pN.The latter value would lead to a contribution of about 63 pN to the constriction force f , which is equal to about three times the observed threshold value of 20 pN and, thus, sufficient to cleave the neck.This simple estimate ignores the possible vicinity of a critical demixing point, at which the line tension must vanish.Therefore, close to such a critical point, the line tension can be reduced by orders of magnitude [22].For the ternary mixture DOPC, sphingomyelin (SM), and cholesterol (CHOL), different compositions have been studied using giant vesicles.A detailed comparison of the experimentally observed two-domain shapes with the shapes computed in the framework of curvature elasticity [23] led to line tension values between 1 pN and 0.01 pN [4,6,8].The line tension contribution 2πλ to the constriction force then varies between 6.3 pN and 0.063 pN.
The other contributions to the constriction forces f in Eqs. ( 64) and (67) are proportional to the bending rigidities κ a and κ b as well as to the excess cur- ).The bending rigidities are of the order of 10 −19 J.For giant vesicles, the excess curvatures are dominated by the spontaneous curvatures m a and m b .Moderate spontaneous curvatures as generated by sugar asymmetries [25] are of the order of 1/(μm).Larger spontaneous curvatures up to about 10/(μm) can be obtained by the binding of His-tagged GFP to the outer membrane leaflet [33].Therefore, the excess curvature terms contribute about 1 pN for m a m b 1/(μm) and about 10 pN for m a m b 10/(μm) to the constriction forces.
Combining the line tension contribution with the excess curvature contributions, we conclude that the constriction forces are sufficiently large to cleave the closed membrane neck when the line tension λ 1 pN and the spontaneous curvatures are of the order of 10/(μm).On the other hand, line tensions below 0.1 pN and moderate spontaneous curvatures of the order of 1/(μm) are unlikely to induce neck fission.

Morphological complexity of multispheres
The two-sphere shapes formed by two-domain vesicles as discussed in the previous section represent the simplest examples for multispheres that can be formed by vesicles with several membrane domains.In general, more complex shapes are also possible.First, each domain of a two-domain vesicle can form a multispherical shape itself.Second, vesicles with several aand/or b-domains can attain multispheres with several ab-necks.

Multispheres of two-domain vesicles
For a two-domain vesicle, both the aand the b-domains can attain a multispherical shape.When the a-domain transforms into a multispherical shape, this shape consists of two or more (punctured) a-spheres, which are connected by closed aa-necks.Likewise, when the bdomain forms a multispherical shape, this shape consists of two or more (punctured) b-spheres, which are connected by closed bb-necks.Thus, each multispherical shape formed by a two-domain vesicle involves both a single ab-neck as discussed in the previous Sect.3 and additional membrane necks between two a-spheres or between two b-spheres, which are governed by the stability conditions for uniform membranes as described in Sect.2.4.Some examples for such multispheres are displayed in Fig. 7

Multispheres of multi-domain vesicles
Next, let us consider vesicle membranes with several aand b-domains and, thus, with more than one domain boundary.All domains are taken to be in chemical equilibrium as described in the next subsection.

Chemical equilibrium between all domains
In order to distinguish the different aand b-domains, we label them by the integers k and n, respectively.The membrane areas of the a k -and b n -domains are denoted by A ak and A bn .In chemical equilibrium, the coexisting aand b-phases are characterized by two different molecular compositions, where one composition applies to all a-domains and the other composition to all b-domains.In order to allow the domains to form different domain patterns, we introduce Lagrange multipliers Σ a and Σ b , which are conjugate to the total surface area of all a-domains and to the total surface area of all b-domains, respectively.The corresponding shape functional F >2Do is obtained by generalizing the shape functional F 2Do for a two-domain vesicle as given by Eq. (A7) in Appendix A. Indeed, apart from the pressure term, each term of the shape functional F 2Do in Eq. (A7) is replaced by a sum over the different a kand b n -domains.In particular, the two tension terms in Eq. (A7) are substituted according to [34] In general, each a k -domain and each b n -domain can now form a multisphere, in close analogy to the multispheres formed by uniform membranes as described in Sect. 2. When the a k -domain forms a multisphere, the individual a k -spheres are labeled by the index i and have the area A aki as well as the mean curvature M aki .It then follows from the first variation of the generalized shape functional F >2Do that the mean curvature M aki satisfies the local shape equation with the total membrane tension which is completely analogous to the local shape equation for uniform membranes as given by Eq. ( 1).One should note that the Lagrange multiplier Σ a and the curvature-elastic parameters κ a and m a are independent of the domain-index k and of the individual sphere index i. 2hen the b n -domain forms a multisphere, the individual b n -spheres are labeled by the index j.The mean curvature M bnj of an individual b nj -sphere formed by the b n -domain fulfills the local shape equation In the examples displayed here, the aand b-sphere connected by the abneck have positive mean curvatures which implies that each ab-neck is positive.Multispheres with negative ab-necks are discussed further below, see Sect.6.2 and Fig. 11 with the total membrane tension which is again completely analogous to Eq. ( 1).The quadratic form of Eq. ( 71) for the mean curvature M aki implies that each a k -domain forms a kispheres with up to two different radii, provided by large a ki -spheres with radius R al and by small a kispheres with radius R as .Likewise, the quadratic form of Eq. ( 73) for M bnj has the consequence that the b ndomain forms b nj -spheres with up to two different radii, R bl and R bs , corresponding to large and small b njspheres.Some examples for multispheres arising from three-domain vesicles in chemical equilibrium are displayed in Fig. 8.In these examples, the vesicle membrane consists of one a-domain and two b-domains, forming different clusters of aand b-spheres, which are connected by two ab-necks.

Geometry of multispheres
As before, the individual a-spheres are labeled by a ki , where the integer k labels the a k -domain and the integer i the individual a ki -spheres formed by the a k -domain.Likewise, the individual b-spheres are labeled by b nj where the integer n is the index of the b n -domain and the integer j labels a certain b nj -sphere formed by the b n -domain.The a ki -sphere has the radius R aki and the rescaled radius r aki = R aki /R ve ; the b nj -sphere has the radius R bnj and the rescaled radius r bnj = R bnj /R ve .
The a k -domain with area A ak and the b n -domain with area A bn are now characterized by the area fractions and with the total area fractions of the aand b-domains as given by and Therefore, in contrast to the two-sphere vesicles in Sect.3.3, all rescaled radii now fulfill the inequalities where the equality signs apply to a single a-domain forming a single a-sphere and to a single b-domain forming a single b-sphere, respectively.Furthermore, the total volume enclosed by all aspheres and all b-spheres is denoted by V a and V b , respectively, which leads to the rescaled volumes If all spheres have a positive mean curvature as in Figs.7a, c and 8a, b, the rescaled volumes are given by If some aor b-spheres have a negative mean curvature and enclose some part of the exterior compartment as in Figs.7b, d and 8c, d, we have to substract the subvolumes of these spheres from the combined volume of the other spheres with a positive mean curvature.

Stability of ab-necks
To discuss the stability of an ab-neck between the a kand the b n -domain, we label the a-sphere adjacent to which has the same form as Eq. ( 46) but with the radii r a and r b replaced by the radii r ak1 ≤ r a and r bn1 ≤ r b , leading to the mean curvatures M ak1 = 1/r ak1 ≥ 1/r a = 1/ √ Φ a and M bn1 = 1/r bn1 ≥ 1/r b = 1/ √ Φ b , as follows from Eq. (79).The stability condition in Eq. ( 82) can be rewritten in the form mb ≥ h out ( ma ) or y ≥ h out (x) (83) with x = ma and y = mb as before.The limiting case y = h out (x) describes the limit shapes L ab for a positive ab-neck.These limit shapes define a straight line in the (x, y)-plane, which is quite similar to the purple L abline in Fig. 6a.The stability condition for a negative ab-neck with an in-budded b-sphere has the form: which has the same form as Eq. ( 56) but with the radii r a and r b replaced by the radii r ak1 ≤ r a and r bn1 ≤ r b , leading to the mean curvatures M ak1 = +1/r ak1 ≥ 1/r a and M bnj = −1/r bn1 ≤ −1/r a .The stability condition in Eq. ( 84) can be rewritten in the form: The limiting case mb = h b-in ( ma ) describes the limit shapes L ab for a negative ab-neck.These limit shapes are located along a straight line in the (x, y) plane, which is quite similar to the purple L ab -line in Fig. 6b.

Stability of aa-and bb-necks
In general, the multispheres consist of ab-necks as well as bband aa-necks, see the examples in Figs.7 and 8.Each ab-neck can be positive or negative as described in Sect.3.2.2.In addition, each a-domain forms one asphere or a cluster of several a-spheres connected by aa-necks.Likewise, each b-domain forms one b-sphere or a cluster of b-spheres connected by bb-necks.All aanecks are either positive or negative and likewise for the bb-necks.Indeed, the stability regimes for positive and negative aaor bb-necks have no overlap in the morphology diagrams defined by the two spontaneous curvatures x = ma and y = mb , see Fig. 9.
First, consider two a-spheres with rescaled radii r ak1 = R ak1 /R ve and r ak2 = R ak2 /R ve , which are connected by a closed aa-neck.It follows from Eq. ( 17) that a positive aa-neck is stable if the rescaled curvatures fulfill the inequality which defines the right stability regime in Fig. 9a.On the other hand, a negative aa-neck with r ak1 > r ak2 is stable provided as follows from Eq. ( 18), leading to the left stability regime in Fig. 9a.Next, consider two b-spheres with rescaled radii r bn1 = R bn1 /R ve and r bn2 = R bn2 /R ve , which are connected by a closed bb-neck.Equation ( 17) now implies Fig. 9 Morphology diagrams with coordinates x = ma and y = mb : a Stability regimes (yellow) for positive and negative aa-necks.For positive aa-necks, the spontaneous curvature ma ≥ M eff aa > 0 as in Eq. ( 86), which defines the right stability regime.For negative aa-necks, the spontaneous curvature ma ≤ M eff aa < 0 as in Eq. ( 87), leading to the left stability regime.Because the left and the right stability regime have no overlap, all aa-necks must be either positive or negative; and b Stability regimes (yellow) for positive and negative bb-necks.For positive bb-necks, the spontaneous curvature mb ≥ M eff bb > 0 as in Eq. ( 88), corresponding to the upper stability regime.For negative bbnecks, the spontaneous curvature ma ≤ M eff aa < 0 as in Eq. ( 89), which defines the lower stability regime.Because the lower and the upper stability regime have no overlap, all bb-necks must be either positive or negative that a positive bb-neck is stable if corresponding to the upper stability regime in Fig. 9b.On the other hand, a negative bb-neck with r bn1 > r bn2 is stable for as in Eq. ( 18), which defines the lower stability regime in Fig. 9b.

Multispherical shapes of two-domain vesicles
In this section, the stability of multispherical shapes formed by vesicles with one aand one b-domain will be examined in more detail.These multispheres involve both a single ab-neck and additional aaand bb-necks between two a-spheres and two b-spheres.The stability of each neck is governed by its own stability condition as described in Sects.5.2.3 and 5.2.4.In order to identify the parameter regimes of stable multispheres, we need to impose and combine the stability conditions for all membrane necks, which are present in the multisphere.

Four-sphere shapes with positive ab-neck
Representative examples for four-sphere shapes with one positive ab-neck are displayed in Fig. 10.These multispheres consist of two a-spheres connected by a single aa-neck and of two b-spheres connected by a single bbneck.Each of the four spheres can have a different mean curvature, in accordance with the shape equations for the aand b-spheres.Both the aa-neck and the bb-neck can be positive or negative, which implies four different types of four-sphere shapes with positive ab-necks as in Fig. 10.In each panel of this figure, the top and bottom subpanels display one of the four-sphere shapes together with the corresponding stability regime within the morphology diagram defined by the rescaled spontaneous curvatures x = ma and y = mb .Each four-sphere shape in the upper row of Fig. 10 involves three closed necks, each of which is governed by a different stability condition.These three stability conditions determine three half-planes, which represent the stability regimes for the three individual necks.The intersection of these three half-planes determines the stability regime of the four-sphere shape under consideration.Each half-plane is bounded by a line of limit shapes, denoted by L ab , L aa , and L bb .Thus, the lower subpanels of Fig. 10 display three lines of limit shapes, purple L ab -lines as in Fig. 6a as well as red L aa -lines and blue L bb -lines as in Fig. 9.
In Fig. 10a, all four spheres have positive mean curvatures which implies that all three necks are positive.Four-sphere shapes with stable ab-necks must be located to the right of the purple L ab -line as in Fig. 6a.Furthermore, the positive aa-neck confines the stability regime of the four-sphere shape to positive values of x = ma as in Fig. 9a, and the positive bb-neck is only (yellow) within the morphology diagrams defined by the rescaled spontaneous curvatures x = ma and y = mb .The purple, red, and blue lines represent the lines of limit shapes L ab , Laa, and L bb , respectively, which provide the boundaries for the stability regimes.In (d), the three lines lead to a small, triangular stability regime, which requires finetuning of both spontaneous curvatures stable for sufficiently positive values of y = mb as in Fig. 9b.As a consequence, the yellow stability regime in the bottom row of Fig. 10a is confined to the upper right quadrant of the (x, y)-plane, which implies that the stability of the four-sphere shape in Fig. 10a requires sufficiently large spontaneous curvatures ma = x and mb = y.
In Fig. 10b, the bb-neck is negative which moves the yellow stability regime to negative values of y = mb and thus to the lower right quadrant of the morphology diagram On the other hand, Fig. 10c involves a negative aa-neck, which moves the stability regime to negative values of x = ma and thus to the upper left quadrant of the morphology diagram.Finally, the four-sphere shape in Fig. 10d involves both a negative aaand a negative bb-neck.In the latter case, the stability regime is confined to the small triangle formed by the three lines of limit shapes.Therefore, the formation of the multisphere in Fig. 10d requires fine-tuning of the two spontaneous curvatures ma and mb .

Four-sphere shapes with negative ab-neck
Four-sphere shapes with one negative ab-neck are displayed in Fig. 11.In these examples, the negative abneck arises from the in-budded b-domain.Each foursphere shape consists of two a-spheres and two bspheres.Each of these four spheres can have a different radius, in accordance with the two shape equations for the aand b-domain.In addition to the negative ab-neck, the multispheres in Fig. 11 again involve a single aa-neck and a single bb-neck, both of which can be positive or negative, generating four different foursphere shapes with a negative ab-neck.In each panel of Fig. 11, the top and bottom subpanels display one of these four-sphere shapes and the corresponding stability regime within the morphology diagram as defined by the rescaled spontaneous curvatures x = ma and y = mb .
All four-sphere shapes in the upper row of Fig. 11 involve one aa-neck and one bb-neck in addition to the negative ab-neck.Each neck is stably closed when it fulfills the associated stability condition.Each of these conditions again defines a half-plane in the morphology diagram.The intersection of these three halfplanes determines the stability regime of the foursphere shape.Furthermore, each half-plane in Fig. 11 is bounded by a line of limit shapes as displayed in the lower subpanels of Fig. 11: purple L ab -lines as in Fig. 6b as well as red L aa -lines and blue L bb -lines as in Fig. 9.
In Fig. 11a, the a-domain forms an in-bud which implies a negative aa-neck.Furthermore, the in-budded b-domain consists of two b-spheres with negative mean curvature, which leads to a negative bb-neck.The negative aa-neck and the negative bb-neck are stably closed for sufficiently large negative values of the spontaneous curvatures x = ma and y = mb .Therefore, the yellow stability regime of this shape is confined to the lower left quadrant of the (x, y)-plane.
In Fig. 11b, the aa-neck is positive whereas the bbneck is negative.The positive aa-neck shifts the stability regime to positive values of x = ma and thus to the lower right quadrant of the (x, y)-plane.In Fig. 11c, the aa-neck is negative whereas the bb-neck is positive.The positive bb-neck shifts the stability regime to positive values of y = mb and thus to the upper left quadrant of the (x, y)-plane.Finally, the four-sphere shape in Fig. 11d involves both a positive aaand a positive bbneck.In the latter case, the stability regime is confined to the small triangle enclosed by the three lines of limit Fig. 11 Four-sphere shapes and their stability regimes for negative ab-necks.The multispheres consist of one a-domain (red) and one in-budded b-domain (blue), forming two aspheres and two b-spheres: a Negative aa-neck and negative bb-neck; b Positive aa-neck and negative bb-neck; c Negative aa-neck and positive bb-neck; and d Positive aa-neck and positive bb-neck.The subpanels in the bottom row display the corresponding stability regimes (yellow) within the morphology diagram defined by the rescaled spontaneous curvatures x = ma and y = mb .The purple, red, and blue lines represent the lines of limit shapes L ab , Laa, and L bb , respectively, which provide the boundaries for the stability regimes of the four-sphere shapes.In (d), the three lines lead to a small, triangular stability regime, which requires fine-tuning of both spontaneous curvatures shapes.Therefore, the formation of the shape displayed in Fig. 11d requires fine-tuning of the two spontaneous curvatures ma and mb .

Multispheres with multiple aaand bb-necks
In general, a two-domain vesicle can form multispheres that consist of an a-cluster with more than two aspheres and a b-cluster with more than two b-spheres.The a-cluster is built up from large a-spheres with radius R al and small spheres with radius R as as follows from the local shape equation in Eq. (71).Furthermore, Sect.2.2.2 and 2.5 for uniform membranes imply that all aa-necks are either positive or negative.Thus, the a-cluster can attain two global architectures, corresponding to cases I and II for uniform membranes.
For case I, the large and small a-spheres have positive mean curvature and are connected by positive aa-necks.For case II, the a-cluster is provided by one large a-sphere with positive mean curvature and multiple small a-spheres with negative mean curvature, with all a-spheres being connected by negative aa-necks.The same two cases can be distinguished for the b-cluster.For case I, the large and small b-spheres have positive mean curvature and are connected by positive bbnecks.For case II, the b-cluster consists of one large b-sphere with positive mean curvature and one or several small b-spheres with negative mean curvature, with all b-spheres being connected by negative bb-necks.
In general, both the aand the b-clusters can involve different types of necks: ss-necks between two small spheres; ls-necks between a large and a small sphere; and ll-necks between two large spheres.The aand b-cluster of the seven-sphere shapes in Fig. 7c, d, for example, involve both ls-necks and ss-necks.The stability of the multisphere is then determined by the least stable necks which impose the strongest closed neck condition on the spontaneous curvatures.
If the cluster of a-spheres belongs to case I with positive mean curvatures of the large and small aspheres, the cluster consists, in general, of large and small spheres, which can be connected by ss-, ls-, or ll-necks.The effective mean curvatures of these necks are ordered according to Therefore, all necks of the a-cluster are stable for sufficiently large and positive spontaneous curvature A special case I is obtained if all spheres of the acluster have the same rescaled radius r a * .Such a multisphere consisting of equally sized a-spheres has the smallest rescaled volume v a of all multispheres with the same total number of a-spheres [25,26].In the latter case, all necks have the same effective mean curvature M eff * * = 1/r a * .These necks are stable if the spontaneous curvature is large and positive with ma ≥ M eff * * (Case I, equally sized spheres).( 92 One example for a multisphere consisting of equally sized spheres as formed by a uniform membrane is displayed in Fig. 2g.
If the cluster of a-spheres belongs to case II, it consists of one large a-sphere with positive mean curvature and one or several small a-spheres with negative mean curvature.Such an a-cluster involves only lsand ss-necks with negative neck curvatures M eff ls < 0 and M eff ss < 0. These necks are stable if The neck stability of the b-cluster is obtained by replacing the domain label a in Eqs. ( 91), (92), and (93) by the domain label b.It follows from these stability conditions for the aand b-cluster that the qualitative features of the morphology diagrams as shown in Figs. 10 and 11 for four-sphere shapes also apply to twodomain vesicles with more than two a-spheres and/or more than two b-spheres.

Nested multispheres from nested domains
A special case of multispheres with several ab-necks is obtained starting from nested domains.For a spherical vesicle, the simplest example for such a domain pattern is displayed by the three-domain vesicle in Fig. 12a.In this example, the southern hemisphere of the vesicle together with a small fraction of the northern hemisphere is covered by a large a 1 -domain, while the northern hemisphere contains a ring-like b 1 -domain, which encloses an even smaller a 2 -domain at the north pole.Note that the a 1 -and the a 2 -domain have the same molecular composition and thus possess the same spontaneous curvature ma .Likewise, the b 1 -and the b 2domain in Fig. 12c are characterized by the same spontaneous curvature mb .Furthermore, all domain boundaries have the same line tension λ.

Nested multispheres with two ab-necks
Deflation of the vesicle in Fig. 12a can lead to the nested multisphere displayed in Fig. 12b.The latter multisphere consists of three spheres that are nested into each other.The a 1 -domain forms the largest sphere with radius R a1 and positive mean curvature M a1 = +1/R a1 > 0, whereas the ring-like b 1 -domain has transformed into the b 1 -sphere with radius R b1 and negative mean curvature M b1 = −1/R b1 < 0. In addition, the small a 2 -domain, which was located close to the north pole of the spherical vesicle in Fig. 12a, now forms the smallest sphere with radius R a2 and positive mean curvature M a2 = +1/R a2 > 0. The requirement that the membrane of the nested multisphere in Fig. 12b should not intersect itself implies that the three spherical radii must satisfy the inequalities R a1 > R b1 > R a2 and that the areas of the three domains are thus ordered according to A a1 > A b1 > A a2 .
The nested multisphere in Fig. 12b involves one negative and one positive ab-neck.The negative ab-neck connects the outer a 1 -sphere with the in-budded b 1 -sphere.The corresponding stability condition is described by Eq. ( 84) and leads to a line of limit shapes L ab as displayed in Fig. 6b.The negative ab-neck is stable for xand y-values below the L ab -lines in Fig. 6b, with the coordinates x and y of the morphology diagram provided by the spontaneous curvatures ma and mb , respectively.The positive ab-neck of the nested multisphere in Fig. 12b is located between the in-budded b 1 -sphere and the out-budded a 2 -sphere.This positive ab-neck is stable for xand y-values above the L ab -line in Fig. 6a.As previously discussed in Sect.3.5.3, the stability regimes for negative and positive ab-necks exhibit a large overlap region in the morphology diagram.The nested multisphere in Fig. 12b is stable for values of the spontaneous curvatures ma and mb within this overlap region.
The formation of the a 1 -and a 2 -spheres with positive mean curvatures M a1 > 0 and M a2 > 0 will be facilitated by positive spontaneous curvature x = ma > 0. Likewise, the formation of the b 1 -sphere with negative mean curvature M b1 < 0 will be supported by negative spontaneous curvature y = mb < 0. However, the absolute values of ma and mb must be sufficiently small so that the two spontaneous curvatures x = ma and y = mb define a point (x, y) of the morphology diagram that is located within the overlap region of the yellow stability regimes in Fig. 6a, b.

Nested multispheres with three ab-necks
Another example for a nested domain pattern is displayed by the four-domain vesicle in Fig. 12c.This domain pattern involves the same a 1 -and b 1 -domains as the pattern in Fig. 12a but the a 2 -domain now forms another ring-like domain, enclosing a fourth b 2 -domain at the north pole.Deflation of such a vesicle can lead to four nested spheres as displayed in Fig. 12d.The a 1 -domain again forms the largest sphere with radius R a1 and positive mean curvature M a1 = +1/R a1 > 0, whereas the ring-like b 1 -domain again forms the b 1sphere with radius R b1 and negative mean curvature M b1 = −1/R b1 < 0. Furthermore, the ring-like a 2domain now turns into the a 2 -sphere with positive mean curvature M a2 = +1/R a2 > 0. Finally, the b 2 -domain close to the north pole becomes the in-budded b 2 -sphere with negative mean curvature M b2 = −1/R b2 < 0. The requirement that the membrane of the nested multisphere in Fig. 12d should not intersect itself implies that the four spherical radii satisfy the inequalities The nested multisphere in Fig. 12d involves two negative and one positive ab-necks.One negative ab-neck connects the outer a 1 -sphere with the in-budded b 1sphere.The positive ab-neck is located between the inbudded b 1 -sphere and the out-budded a 2 -sphere.The Fig. 12 Trasformation of nested multi-domain vesicles into nested multispheres: a a spherical three-domain vesicle consisting of a large a1-domain (red) on the southern hemisphere, a smaller ring-like b1-domain (blue) on the northern hemisphere, and an even smaller a2-domain (red) close to the north pole.Deflation of such a three-domain vesicle can lead to the nested multisphere in (b) with an outer a1-sphere, an in-budded b1-sphere, and an out-budded a2sphere; and c a spherical four-domain vesicle consisting of a large a1-domain (red), a ring-like b1-domain (blue), a ringlike a2-domain (red), and a b2-domain (blue) close to the north pole.For simplicity, the three-and four-domain vesicles in (a) and (c) are taken to be axisymmetric, which implies that the nested domains are concentric as well.However, shifting the concentric domains against each other will again lead to nested multispheres, which are quite similar to those in (b) and (d) stability regimes for these two ab-necks are very similar to the stability regimes for the two ab-necks in Fig. 12b.In addition, the nested multisphere in Fig. 12d contains another negative ab-neck, which connects the out-budded a 2 -sphere with the in-budded b 2 -sphere.
The nested multisphere displayed in Fig. 12d involves two a-spheres with positive mean curvatures M a1 > 0 and M a2 > M a1 and two b-spheres with negative mean curvatures M b1 < 0 and M b2 < M b1 .We could try to add another level of nesting, by adding an a 3 -domain within the b 2 -domain close to the north pole.Such a domain pattern is possible but cannot lead, in chemical equilibrium, to a nested multisphere because such a multisphere would involve three a-spheres with three different mean curvatures, which is inconsistent with the local shape equation as given by Eq. (71).On the other hand, we can connect another a-sphere with positive mean curvature M a2 to the existing a 1 -sphere, thereby creating a positive aa-neck or another b-sphere with negative mean curvature M b2 to the existing b 1sphere via a negative bb-neck.These additional aaor bb-necks will introduce L aa -or L bb -lines as in Fig. 9, which further restrict the stability regimes as discussed in Sect.5.2.4.

Constant-mean-curvature (CMC) surfaces 8.1 Conventional CMC surfaces
In the differential geometry of surfaces [35], multispherical shapes consisting of equally-sized spheres have been studied in the context of constant-mean-curvature (CMC) surfaces, generalizing the concept of minimal surfaces with zero mean curvature M = 0.For a long time, the only examples for freely suspended CMC surfaces with M = 0 were provided by the unduloids of Delaunay [36], which provide a one-parameter family of tubular shapes that interpolate smoothly between multispherical tubes consisting of equally sized (and punctured) spheres and cylindrical tubes.More recently, additional CMC surfaces have been constructed by pertubing a cluster of identical spheres that touch each other [37][38][39][40][41].One example are triunduloids [40,41] that consist of three unduloidal arms connected by a central core as displayed in Fig. 13.
The physical system typically used to motivate CMC surfaces are the shapes of soap films and liquid droplets.However, when the initial cluster of identical and touching spheres is viewed as a cluster of liquid droplets, the resulting CMC surface is not stable.Indeed, the cluster will either fall apart and then form many small droplets or it will coalesce into one large droplet that will eventually attain the shape of a single sphere.However, when the cluster of droplets is enclosed by a membrane, this membrane can lead to stable multispherical shapes as described in Sects.3-7.

Multispheres as generalized CMC surfaces
The conventional CMC surfaces considered in differential geometry have a constant mean curvature that is uniform along the whole surface.For the vesicle surfaces as considered here, examples for such conventional CMC surfaces are provided by multispheres consisting of equally sized spheres as in Fig. 2g.The latter vesicles are bounded by uniform membranes which have a laterally uniform composition as well as laterally uniform elastic properties.However, even uniform membranes can form multispheres with two different piece-wise constant mean curvatures as displayed in most panels of Fig. 2. Such multispheres with different piece-wise constant mean curvature should be regarded as generalized CMC surfaces.
As discussed in the previous sections, vesicles with intramembrane domains can form multispheres with up

Experiments on multispherical vesicles
To experimentally study the multispherical shapes of multi-domain vesicles as determined here theoretically, it will be useful to combine and extend several experimental protocols.The first protocol corresponds to the same procedure as recently used for uniform membranes [25], see Fig. 2, but now applied to ternary lipid mixtures that undergo phase separation into two fluid phases.The second protocol, also developed quite recently [33], allows to control and fine-tune the membrane's spontaneous curvature by the binding of Histagged fluorophores to anchor lipids in the membranes.The third protocol, which has been introduced already some time ago [7,9], generates multi-domain vesicles by electrofusion of membranes.

Deflation combined with solution asymmetry
The multispheres displayed in Fig. 2 have been obtained for giant vesicles by a combination of osmotic deflation and solution asymmetry between the interior and exterior compartments of the vesicles, which were first prepared in a symmetric sucrose solution.Subsequently, a small aliquot of the prepared sucrose-vesicle solution was transferred into the observation chamber where they were added to a larger aqueous droplet that contained primarily glucose, with a glucose concentration that exceeded the sucrose concentration in the aliquot.This transfer or dilution step led to the reduction of the vesicle volume by fast osmotic deflation and, at the same time, to the generation of bilayer asymmetry and spontaneous curvature.
The lipid bilayers of the giant vesicles in Fig. 2 contained binary mixtures of the phospholipid POPC and cholesterol.In order to obtain a lipid bilayer that forms two coexisting fluid phases, a third lipid component such as another phospholipid or sphingomyelin should be added.Thus, the simplest experimental approach to generate multispheres with two or more intramembrane domains will be obtained when the experimental protocol developed in [25] is applied to such ternary lipid mixtures.
For the binary mixtures of POPC and cholesterol, the spontaneous curvatures generated by the sugar asymmetry between sucrose and glucose within the uniform bilayer membranes was of the order of 1/μm.

Fine-tuning of spontaneous curvatures
In order to fine-tune the spontaneous curvatures m a and m b of the aand b-domains, it will be useful to dope the ternary lipid mixtures with some anchor lipids that bind His-tagged fluorophores from the exterior compartment.Such a method has been successfully applied to ternary lipid mixtures of POPC, POPG, and cholesterol exposed to His-tagged GFP [33].In general, the density of the membrane-bound His-tagged fluorophores depends on the density of the anchor lipids as well as on the solution concentration.In Ref. [33], rather large spontaneous curvatures of the order of 1/(100 nm) have been achieved by exposing the giant vesicles to nanomolar concentrations of Histagged GFP.
As explained in Sects.3 and 5 and illustrated in Figs. 6, 10, and 11, the morphology of the theoretically predicted multispheres depends strongly on the rescaled spontaneous curvatures x = ma and y = mb of the aand b-domains.In general, the partitioning of the anchor lipids is expected to lead to different anchor lipid densities within the two types of domains.In addition, one might use two types of anchor lipids, each of which becomes enriched in one of the two membrane domains.
As a result, one should be able to obtain significantly different densities of membrane-bound fluorophores in the aand b-domains and, thus, significantly different spontaneous curvatures ma and mb .The bilayer asymmetries can be further enhanced by asymmetric sugar or ion solutions.In this way, it will become possible to explore large parameter regions of the morphology diagrams in Figs. 6, 10, and 11.To experimentally obtain nested multispherical shapes as displayed in Fig. 12 will be particularly challenging.First, the a and b domains need to possess finetuned values of the spontaneous curvatures x = ma and y = mb , which belong to the overlap region of the stability regimes in Figs.6a and b.Second, to create such nested multispheres from a multidomain vesicle with a spherical shape, one should start from nested aand b-domains as displayed in Fig. 12a, c.

Multi-domain vesicles via electrofusion
In principle, one can create multi-domain vesicles by fusing two GUVs that are aspirated by two micropipettes.When the compositions of the two GUVs correspond to coexisting aand b-phases, the fused GUVs should be composed of stable aand b-domains.Fusion of two aspirated GUVs is, however, difficult to achieve because the aspirated vesicles tend to rupture.A more robust method is provided by electrofusion which has indeed been used to create multi-domain vesicles [7,9].
In these studies, two populations of GUVs have been prepared from different lipid compositions.In Ref. [9], for example, the membranes of one vesicle population consisted of DOPC and CHOL whereas the membranes of the other population were composed of SM and CHOL.Electrofusing a vesicle from the DOPC/CHOL population with a vesicle from the SM/CHOL population then leads to a GUV membrane that contains all three lipid components.
The standard protocol for electrofusion of GUVs is similar to the electrofusion of cells [42,43].This protocol consists of two steps.First, the vesicles are aligned by an alternating electrical field, which brings them into close contact at their poles.Second, a short pulse of a high electric field is applied to the aligned vesicles, thereby electroporating the two vesicle membranes within their contact area.

Different Gaussian curvature moduli
So far, the Gaussian curvature moduli κ Ga and κ Gb of the aand b-domains were taken to have the same value.As mentioned in the introduction, it then follows that the domain boundaries are located within the closed ab-necks.At the end, let us look at the changes arising from different Gaussian curvature moduli, κ Ga = κ Gb in the two domains.In this case, the matching condition in Eq. (40) implies that the domain boundary moves out of the waist-line of the closing neck.The closed neck is then formed by the membrane domain with the larger κ G -value.Indeed, the Gaussian curvature G is negative around the closed neck and makes a more negative contribution to the Gaussian curvature energy in Eq. (A4) when this neck is formed by the domain with the larger κ G -value. 3 This conclusion also applies to negative values of the Gaussian curvatures moduli, that is, to κ Ga < 0 and κ Gb < 0..

Shift of domain boundary
A simple estimate for the displacement of the domain boundary away from the closed neck can be obtained as follows.Such a displacement leads to the change ΔE G in the Gaussian curvature energy as given by Eq. (A5).The largest possible energy gain arising from ΔE G is given by min This energy gain must overcompensate the line energy ΔE λ of the domain boundary with radius R db , which is equal to with positive line tension λ.Therefore, the radius R db of the domain boundary satisfies the inequality For a uniform membrane, the Gaussian curvature modulus κ G is expected to be negative with a magnitude that is comparable to the bending rigidity κ [44][45][46].Therefore, the term |κ Ga − κ Gb | should be comparable to |κ a − κ b |, which is of the order of 10 −19 J. Based on experimental studies of GUVs, the line tension λ was estimated to lie within the range 1 pN and 0.01 pN, depending on the lipid composition of the GUV membranes [4,6,8].It then follows from Eq. (96) that the radius R db of the domain boundary satisfies R db < 100 nm for λ = 1 pN and R db < 1 μm for λ = 0.1 pN.Therefore, the displacement of the domain boundary away from the closed membrane neck will not be detectable by conventional fluorescence microscopy if the line tension is of the order of 1 pN but should become visible for line tension values below 0.1 pN, in accordance with experimental observations [6,8].

Reduction of constriction forces
Now, assume that the a-domain represents the Ld lipid phase and the b-domain the Lo lipid phase, com-pare Fig. 1.The Lo phase is more rigid than the Ld phase which implies that the bending rigidity κ b of the b-domain exceeds the bending rigidity κ a of the a-domain, that is κ b > κ a .Using the previously mentioned estimates κ Ga −κ a and κ Gb −κ b , the inequality κ a < κ b implies that κ Ga > κ Gb and that the neck is formed by the a-domain, that is, by the more flexible domain with the lower bending rigidity.
When the domain boundary moves out of the neck for κ Ga > κ Gb , the closed neck is located within the adomain.More precisely, this neck provides a connection between the complete a-sphere and the narrow a-strip between the neck and the b-domain.The effective mean curvature of this neck is given by where M a−sp is the mean curvature of the complete asphere and M a−st is the mean curvature of the narrow a-strip on the other side of the neck.A positive neck with neck curvature M eff a|a > 0 then experiences the constriction force f = 8πκ a (m a − M eff a|a ) for m a ≥ M eff a|a > 0 (98) as follows from Eq. ( 64) with λ = 0 and m b = m a .This constriction force has the same form as for the positive neck of a uniform GUV membrane with bending rigidity κ a and spontaneous curvature m a [33].
On the other hand, a negative neck with neck curvature M eff a|a < 0 is subject to the constriction force f = 8πκ a (M eff a|a − m a ) for m a ≤ M eff a|a < 0 (99) as follows from Eq. (67) with λ = 0 and m b = m a .The constriction force as given by Eq. (99) has the same form as for the negative neck of a uniform GUV membrane with bending rigidity κ a and spontaneous curvature m a [47].Therefore, if the ab domain boundary moves out of the neck during neck closure, the constriction force f as given by Eq. (98) for positive necks and by Eq. (99) for negative necks contains no contribution from the line tension λ, in contrast to Eqs. (64) and (67), which contain the term 2πλ for both positive and negative ab-necks.Because the line tension is necessarily positive, the constriction forces as given by Eqs.(98) and (99) are reduced compared to the constriction forces in Eqs.(64) and (67).Nevertheless, the constriction forces in Eqs. ( 98) and (99), which have the same form as the forces experienced by the closed necks of a uniform GUV membrane, can be sufficiently large to induce neck fission as demonstrated experimentally in Ref. [33].

Conclusion
In this paper, multispherical shapes of vesicles were studied using the theory of curvature elasticity.We started with a brief review of multispheres formed by uniform membranes and introduced the distinction between positive and negative membrane necks based on the sign of the necks' effective mean curvature (Sect.2.5).We then described multispheres formed by vesicles with two intramembrane domains, one aand one b-domain, which arise from membrane phase separation into two fluid phases.These two-domain vesicles can form two-sphere shapes consisting of one aand one b-sphere, connected by a single closed ab-neck.
Depending on the mean curvatures M a and M b of these two spheres, four different two-sphere morphologies can be distinguished as shown in Fig. 3.The morphologies with out-budded domains have positive abnecks, those with in-budded domains have negative abnecks.The stability of the four two-sphere morphologies as formed by two-domain vesicles depends on the stability of their closed ab-necks.The corresponding stability relations are given by Eq. ( 44) for positive ab-necks and by Eq. ( 55) for negative ab-necks.The resulting morphology diagrams are displayed in panels a and b of Fig. 6.
The closed ab-necks experience constriction forces as defined by Eq. ( 63), which act to compress these necks.The form of the constriction forces is provided by Eq. (64) for out-budded domains with positive abnecks and by Eq. (67) for in-budded domains with negative ab-necks.These constriction forces must exceed about 20 pN in order to induce membrane fission across the closed neck [33], thereby dividing the budded vesicle into two daughter vesicles.It is argued in Sect.4.3 that the membrane necks undergo fission for large line tensions of the domain boundaries and/or large spontaneous curvatures but remain stable against fission for smaller line tensions and moderate spontaneous curvatures.If the Gaussian curvature moduli of the aand b-domains are different, the constriction forces are given by Eqs.(98) and (99), which contain no contribution from the line tenion of the domain boundary.As a consequence, different Gaussian curvature moduli κ Gb = κ Ga act to reduce the constriction forces at closed membrane necks, see Sect.10.2.
The morphological complexity of multispherical shapes formed by multi-domain vesicles arises from two different mechanisms.First, each domain of a two-domain vesicle with a single ab-neck can form a multispherical shape by itself.Second, vesicles with more than two domains can form multispheres with more than one abneck.Examples for multispherical shapes with one and several ab-necks are displayed in Figs.7 and 8, respectively.In addition to the ab-necks, these shapes involve closed aaand bb-necks.The stability regimes for the latter necks are displayed in Fig. 9.In general, all necks of a multisphere must be stably closed, a condition that acts to reduce the stability regime of the respective multisphere as illustrated in Figs. 10 and 11.Partic-ularly interesting multispheres are formed by vesicles with nested aand b-domains as shown in Fig. 12.
From a mathematical point of view, multispheres represent generalized CMC surfaces, which exhibit up to four different piece-wise constant mean curvatures as discussed in Sect.8. Examples for conventional CMC surfaces with one constant mean curvature are provided by multispheres consisting of equally sized spheres as in Fig. 2g and by the multispherical triunduloid in Fig. 13c.The latter shape provides a model for the three-way junctions of membrane nanotubes as observed in the endoplasmic reticulum [48].As explained in Sect.9, the multispherical shapes obtained here from the theory of curvature elasticity can be studied experimentally, generalizing available protocols for the multisphere formation of uniform membranes, for the fine-tuning of the spontaneous curvatures, and for the preparation of multidomain vesicles by electrofusion.
Two-domain vesicles have also been studied on the nanoscale by simulations using dissipative particle dynamics [49][50][51][52].One of these studies provided a series of simulation snapshots for the closure of the ab-neck [52].The snapshots indicate that the ab domain boundary stayed in the membrane neck during the whole neck closure process of the nanovesicle.It then follows from Eq. ( 40) and Sect. 10 that the two Gaussian curvature moduli κ Ga and κ Gb were identical for the aand bdomains of the vesicles studied in Ref. [52].Additional simulation studies are required in order to determine the location of the domain boundary during the neck closure process for other two-domain vesicles.Based on recent simulation results for nanovesicles [53], one would expect that the difference in the two Gaussian curvature moduli will depend on the stress asymmetry between the leaflet tensions of the lipid bilayers.

Fig. 3
Fig. 3 Different types of two-sphere shapes formed by vesicles with one a-domain (red) and one b-domain (blue).The radius of the a-sphere is denoted by Ra, the radius of the b-sphere by R b : a, b Out-budded two-sphere shapes with Ra > R b in (a) and R b > Ra in (b); and (c, d) In-budded two-sphere shapes with Ra > R b in (c) and R b > Ra in (d)

Fig. 4 Fig. 5
Fig. 4 Formation of two-sphere shapes with positive abnecks by osmotic deflation (→), which reduces the vesicle volume whereas osmotic inflation (←) increases this volume: a, b When deflated, a spherical two-domain vesicle with rescaled volume v = 1 and area fraction Φ b = 9/25 = 0.36 forms a two-sphere shape with one out-budded a-sphere of radius ra = 4/5 and one out-budded b-sphere of radius r b = 3/5, thereby reducing the rescaled volume from v = 1 ) becomes more transparent when we use the rescaled and dimensionless mean curvatures Ma = M a R ve = 1/r a and Mb = M b R ve = 1/r b as well as the rescaled and dimensionless spontaneous curvatures defined by ma ≡ R ve m a and mb ≡ R ve m b .(

Fig. 6
Fig.6 Stability regimes (yellow) for closed ab-necks of two-sphere vesicles within the morphology diagrams defined by the coordinates x ≡ ma and y ≡ mb .All two-domain vesicles have the same area fraction Φ b = 0.36: a Stability regime for positive ab-necks of two-sphere vesicles with positive mean curvature Ma > 0 of the a-sphere, positive mean curvature M b > 0 of the b-sphere, and positive effective mean curvature M eff ab > 0 of the ab-neck as in Eq. (27); and b Stability regime for negative ab-necks of two-sphere vesicles with in-budded b-domains (blue), corresponding to positive mean curvature Ma > 0 of the a-sphere, negative mean curvature M b < 0 of the b-sphere, and negative effec- 67) for closed necks of in-budded b-domains with M a = +1/R a and M b = −1/R b .The constriction force f in Eq. (67) vanishes along the line of limit shapes L ab as described by Eq. (57) and illustrated in Fig. 6b.The force is positive within the yellow stability regime below the L ab -line in Fig. 6b and increases with increasing line tension λ as well as with increasing excess curvatures M a − m a and M b − m b .

Fig. 7
Fig. 7 Some examples for multspherical shapes formed by two-domain vesicles.Each multisphere has a single ab-neck that contains the domain boundary between the aand bdomain: a, b Four-sphere shapes with two aand two bspheres connected by one aa-neck and one bb-neck; and c, d Seven-sphere shapes with three aand four b-spheres con-

4 4 Fig. 8
Fig.8 Some examples for multispherical shapes with two ab-necks that can be formed by vesicles with three domains in chemical equilibrium: a Three-sphere vesicles with one a-sphere (red) and two b-spheres (blue); b-d Multisphere vesicles with one cluster of three a-spheres (red) and two

Fig. 10
Fig. 10 Four-sphere shapes and their stability regimes for positive ab-necks.The multispheres consist of one a-domain (red) and one b-domain (blue) forming two a-spheres and two b-spheres: a Positive aa-neck and positive bb-neck; b positive aa-neck and negative bb-neck; c negative aa-neck and positive bb-neck; and d Negative aa-neck and negative bb-neck.The bottom row displays the stability regimes

4 4 Fig. 13
Fig. 13 Multispherical limit of triunduloids: a, b Triunduloids with open membrane necks as studied in differential geometry.Reprinted with permission from Ref. [40] (Copyright 1997, Springer-Verlag); and c Triunduloid built Eur. Phys.J. E (2024) 47 :4 Page 27 of 28 4 ψ(s db ) Tilt angle of surface normal at domain boundary R Curvature radius of spherical segment, which is always positive r Rescaled curvature radius of spherical segment, r ≡ R/Rve Ra Radius of a-sphere for two-sphere vesicles R aki Radius of a-sphere formed by a k -domain R ak1 Radius of a-sphere formed by a k -domain and connected to b-sphere R b Radius of b-sphere for two-sphere vesicles R bnj Radius of b-sphere formed by bn-domain R bn1 Radius of b-sphere formed by bn-domain and connected to a-sphere R db Radius of circular domain boundary Rne Radius of circular membrane neck Rve Basic length scale provided by vesicle size, Rve = Coordinate of morphology diagrams with x = ma y Coordinate of morphology diagrams with y = mb