Shape deformation of lipid membranes by banana-shaped protein rods: Comparison with isotropic inclusions and membrane rupture

Shape deformation of lipid membranes by banana-shaped protein rods: Comparison with isotropic inclusions and membrane rupture


The assembly of curved protein rods on fluid membranes is studied using implicit-solvent meshless membrane simulations. As the rod curvature increases, the rods on a membrane tube assemble along the azimuthal direction first and subsequently along the longitudinal direction. Here, we show that both transition curvatures decrease with increasing rod stiffness. For comparison, curvature-inducing isotropic inclusions are also simulated. When the isotropic inclusions have the same bending rigidity as the other membrane regions, the inclusions are uniformly distributed on the membrane tubes and vesicles even for large spontaneous curvature of the inclusions. However, the isotropic inclusions with much larger bending rigidity induce shape deformation and are concentrated on the region of a preferred curvature. For high rod density, high rod stiffness, and/or low line tension of the membrane edge, the rod assembly induces vesicle rupture, resulting in the formation of a high-genus vesicle. A gradual change in the curvature suppresses this rupture. Hence, large stress, compared to the edge tension, induced by the rod assembly is the key factor determining rupture. For rod curvature with the opposite sign to the vesicle curvature, membrane rupture induces inversion of the membrane, leading to division into multiple vesicles as well as formation of a high-genus vesicle.


I Introduction

In living cells, biomembranes often dynamically change their shapes to carry out their biological functions. Membrane budding, fusion, and fission occur during endo/exocytosis and vesicle transports. Cell organelles have specific shapes depending on their functions. Various types of proteins participate in regulation of these dynamic and static membrane shapes McMahon and Gallop (2005); Shibata et al. (2009); Drin and Antonny (2010); Baumgart et al. (2011); Johannes et al. (2014); McMahon and Boucrot (2015). These proteins mainly control local membrane shapes in two ways: hydrophobic insertions (wedging) and scaffolding. In the former mechanism, a part of the protein, such as an amphipathic -helix, is inserted into the lipid bilayer membrane. In the latter mechanism, the protein domain has a strong affinity for the lipid polar head groups and adsorbs onto the lipid membrane. A BAR (Bin/Amphiphysin/Rvs) domain, which consists of a banana-shaped dimer, mainly bends the membrane along the domain axis via scaffolding Itoh and De Camilli (2006); Masuda and Mochizuki (2010); Zhao et al. (2011); Mim and Unger (2012); Simunovic et al. (2015). Some of the BAR superfamily proteins, such as N-BAR proteins, also have hydrophobic insertions. Experimentally, the formation of membrane tubes and curvature-sensing by various types of BAR superfamily proteins have been observed Itoh and De Camilli (2006); Masuda and Mochizuki (2010); Zhao et al. (2011); Mim and Unger (2012); Simunovic et al. (2015); Peter et al. (2004); Mattila et al. (2007); Frost et al. (2008); Wang et al. (2009); Zhu et al. (2012); Tanaka-Takiguchi et al. (2013); Shi and Baumbart (2015); Prévost et al. (2015); Isas et al. (2015); Adam et al. (2015). Dysfunctional BAR proteins are considered to be implicated in neurodegenerative, cardiovascular, and neoplastic diseases. Thus, it is important to understand the mechanism by which membrane shaping is regulated by proteins, not only from a basic science perspective but also for medical applications.

Homogeneous lipid membranes in a fluid phase are laterally isotropic and have zero spontaneous curvature. A local non-zero isotropic spontaneous curvature can be induced by adhesion of spherical colloids and polymer-anchoring, as well as by transmembrane and other proteins Lipowsky (2013). Here, we call these objects that induce isotropic spontaneous curvature an isotropic inclusion. Their assembly into preferred curvature regions Góźdź (2006); Phillips et al. (2009); Greenall and Gompper (2011); Noguchi (2012); Aimon et al. (2014) and membrane-mediated interactions between the colloids Reynwar et al. (2007); Auth and Gompper (2009); Šarić and Cacciuto (2012) have been previously explored.

In contrast, BAR domains, which are banana-shaped, generate an anisotropic curvature. Amphipathic -helices can also yield an anisotropic curvature Gómez-Llobregat et al. (2016). Recently, the anisotropic nature of curvature has received increasing attention theoretically. The classical Canham–Helfrich curvature free energy Canham (1970); Helfrich (1973) has been extended to anisotropic curvatures Fournier (1996); Kabaso et al. (2011); Iglič et al. (2006). To simplify the interactions, the protein and membrane underneath it have been often modeled together as an undeformable object with a fixed curved shape such as a point-like object with an anisotropic curvature Dommersnes and Fournier (1999, 2002) and a bent elliptical surface Schweitzer and Kozlov (2015). Furthermore, it has also been clarified that two undeformable parallel rods have an attractive interaction but the interaction is repulsive for a perpendicular orientation.

Atomic and coarse-grained molecular simulations Blood and Voth (2006); Arkhipov et al. (2008); Khelashvili et al. (2009); Yu and Schulten (2013); Simunovic et al. (2013a); Simunovic and Voth (2015) have been employed to investigate molecular-scale interactions between BAR proteins and lipids. The scaffold formation Yu and Schulten (2013) and linear assembly Simunovic et al. (2013a) of BAR domains have been demonstrated. To investigate large-scale membrane deformations, a dynamically triangulated membrane model Ramakrishnan et al. (2012, 2013) and meshless membrane models Noguchi (2014, 2015, 2016); Ayton et al. (2009) have been employed; consequently, various (meta)stable vesicle shapes Ramakrishnan et al. (2012, 2013); Noguchi (2014, 2015) and the tubule formation dynamics Noguchi (2016) have been reported. Using meshless membrane and molecular simulations, vesicle rupture into high-genus vesicles has also been investigated Ayton et al. (2009); Simunovic et al. (2013b). The high-genus vesicles obtained in this way resemble electron microscopic images of high-genus liposomes induced by N-BAR proteins well Ayton et al. (2009); Simunovic et al. (2013b). Despite these numerous advancements, many questions related to the coupling between membrane shape deformation and the assembly of the protein rods remain.

In this paper, we focus on three questions: (i) How does protein elasticity modify protein assembly? (ii) How is rod assembly different from the isotropic inclusions? (iii) How is membrane rupture induced by protein rods? In the previous rupture simulations Ayton et al. (2009); Simunovic et al. (2013b), the effects of bending rigidity and rod curvature have been investigated but the other mechanical properties were not varied. We show here that the line tension of the membrane edge and the annealing speed, as well as the rod stiffness and density, are important parameters that determine the condition of membrane rupture.

In Sec. II, the simulation model and method are described. We simulate membrane tubes and vesicles using an implicit-solvent meshless membrane model Noguchi (2009); Noguchi and Gompper (2006); Shiba and Noguchi (2011); Noguchi (2013, 2014, 2015, 2016). A banana-shaped protein rod is assumed to be strongly adsorbed onto the membrane and the protein and membrane region below it are modeled as a linear string of particles with a bending stiffness and preferred curvature. In order to investigate the membrane-mediated interactions, no direct attractive interaction is considered between the rods.

In Sec. III, the coupling between the assembly of the isotropic inclusions and shape deformation of membrane tubes and vesicles are presented. In Sec. IV, the assembly of the protein rods in the membrane tubes is shown. The dependence on the rod stiffness is investigated and the results are compared with those of the isotropic inclusions. In Sec. V, the vesicle rupture into high-genus vesicles and vesicle division are presented. The summary is given in Sec. VI.

Ii Simulation Model and Method

ii.1 Membrane Model

We employ a spin meshless membrane model Shiba and Noguchi (2011); Noguchi (2014, 2015, 2016). The details of this meshless membrane model are described in Ref. Shiba and Noguchi, 2011. The position and orientational vectors of the -th particle are and , respectively. The membrane particles interact with each other via a potential,


where , , and . Each particle has an excluded volume with a diameter that results from the repulsive potential, , with a cutoff at .

The second term in Eq. (1) represents the attractive interaction between particles. An attractive multibody potential is employed to allow the formation of a fluid membrane over wide parameter ranges. The potential is given by


Here, and , where is a cutoff function Noguchi and Gompper (2006)


with , , and the cutoff radius . Here, is employed, as described in Ref. Noguchi (2014), instead of , as described in Ref. Shiba and Noguchi (2011), to use a less steep function. The density in is the characteristic density at which the attraction is smoothly truncated. For , acts as a pairwise attractive potential, while it approaches a constant value for .

The third and fourth terms in Eq. (1) are discretized versions of the tilt and bending potentials, respectively. A smoothly truncated Gaussian function Noguchi and Gompper (2006) is employed as the weight function


where , , and .

The membranes are in a fluid phase over a wide range of these parameters, and the properties of the fluid membrane can be widely varied. The spontaneous curvature of the membrane is given by Shiba and Noguchi (2011). For the membrane particles, which do not consist of proteins, and are used. The bending rigidity and edge tension can be controlled by and , respectively. is a linear function of : at , while increases only % from to . is a monotonic increasing function of : , , and at , , and , respectively. We fix except for the vesicle-rupture simulations.

ii.2 Protein Model

A laterally isotropic membrane inclusion is modeled as a membrane particle with times larger bending rigidity and isotropic spontaneous curvature . For the neighbor pair of inclusions, and are employed in Eq. (1). For the pair of an inclusion and a membrane particle, averaged values are used as and .

Note that, if the same values are used for a pair of an inclusion and a membrane particle, an additional attraction between the inclusions is induced by depletion, since the inclusion assembly reduces the area of the large bending rigidity. A similar attraction has previously been obtained for the binding sites of two membranes, when the membranes around the binding sites are hardened Noguchi (2013).

The protein rod is modeled as a linear chain of membrane particles. We use and a rod length of , which corresponds to the typical aspect ratio of BAR domains. The BAR domain width is approximately nm, and its length ranges from to nm Masuda and Mochizuki (2010). The protein particles in each protein rod are connected by a bond potential . The bending potential is given by , where . We use , , and . The membrane potential parameters between neighboring protein particles in each rod are modified as and in order to ensure bending of the rod along the normal to the membrane surface.

ii.3 Simulation Method

The motion of the particle position and the orientation are given by underdamped Langevin equations:


where and are the mass and moment of inertia of the particle, respectively. The forces are given by and with the perpendicular component and a Lagrange multiplier to keep . According to the fluctuation–dissipation theorem, the friction coefficients and and the Gaussian white noises and obey the following relations of their averages and variances:


where and . The Langevin equations are integrated by the leapfrog algorithm Allen and Tildesley (1987); Noguchi (2011). In this study, we use , , , and where . The simulation results are displayed with a time unit of , where is the diffusion coefficient of the membrane particles in the tensionless membranes; is calculated from the mean square displacement of the particles: , so that . This time unit is estimated as from nm and m/s of transmembrane proteins Ramadurai et al. (2009).

The assemblies of the isotropic inclusions and protein rods on membrane tubes and vesicles with were investigated at the inclusion density and the rod density , where and are the numbers of the isotropic inclusions and rods, respectively. The replica exchange molecular dynamics Hukushima and Nemoto (1996); Okamoto (2004) for or Noguchi (2014, 2015) is used to obtain the thermal equilibrium states. Membrane rupture was investigated for vesicles with for . The error bars are estimated from four and independent runs for the replica exchange simulations and vesicle rupture, respectively.

The tube length is fixed in the longitudinal () direction and periodic boundary conditions are employed. The radius of the tube is for . This tube radius is used, if not otherwise specified. In Sec. IV, is also varied in order to investigate the tube radius dependence. The radii of the vesicles are and at and , respectively, in the absence of the rods.

Figure 1: (Color online) Deformation of membrane tubes induced by isotropic inclusions at and . (a)–(d) Snapshots for (a) , (b) , (c) , and (d) at . The front and side views are shown. The inclusion is displayed as a sphere whose halves are colored in red (dark gray) and in yellow (light gray). The orientation vector lies along the direction from the yellow (light gray) to red (dark gray) hemispheres. Transparent gray particles represent membrane particles. (e)–(h) Fourier amplitudes of (e), (f) membrane shape and (g), (h) the inclusion densities as functions of the spontaneous curvature . The amplitudes of the lowest Fourier mode along the azimuthal () and longitudinal () directions are calculated for the membrane shape ( and ) and densities ( and ). The mean axial force of the membrane tube is shown in the inset of (e). Error bars are displayed at several data points.
Figure 2: (Color online) Vesicle deformation induced by isotropic inclusions at and . (a)–(c) Snapshots for (a) , (b) , and (c) at . (d) Three mean eigenvalues of the gyration tensor of the vesicle at , , and . (e) Probability distribution of the inclusions along the axis for , , , and at . Error bars are displayed at several data points.

Iii Isotropic inclusions

First, we investigated the assembly of the isotropic inclusions in order to compare it with the rod assembly. When the bending rigidity difference between the inclusions and membrane is small, the inclusions are isotropically distributed on a membrane tube and vesicle. However, the inclusion assembly occurs at (see Figs. 1 and 2). The assembly regions have curvatures closer to the preferred curvature of the inclusions. Hence, for hard inclusions, assembly reduces the bending energy, which is greater than the loss of the mixing entropy.

iii.1 Membrane Tube

For at , the membrane tube deforms into an elliptic cylinder and the inclusions assemble into two longitudinal domains in two flatter regions [see Fig. 1(a)]. The tips of the ellipse have a large curvature, but the inclusions in the flatter regions can have a low bending energy; hence, the total bending energy is reduced. This phase separation in the azimuthal direction is captured by the amplitudes of the Fourier modes [see Figs. 1(e)–(h)]. The lowest Fourier modes of the membrane shape and inclusion density along the azimuthal () direction are given by and , respectively, where . In the axial () direction, and . With increasing , the amplitudes of and along the direction increase concurrently at .

As increases, the membrane also deforms sigmoidally along the axial () direction and the inclusions become concentrated in the convex region [see Fig. 1(b)]. As a result, the Fourier amplitudes of and decrease and increase, respectively. At , the Fourier amplitudes of reach the values of the homogeneous membranes; the membrane deforms into a circular sigmoidal shape and the inclusions are concentrated along a convex ring [see Fig. 1(c)]. With a further increase in , the inclusions become uniformly distributed and the membrane is less undulated in the axial direction [see Figs. 1(d), (g), and (h)]. The spontaneous curvatures between the inclusion pair and between the membrane particle and inclusion are and , so that for , local membrane spontaneous curvature can match the curvature of the membrane tube.

An even further increase in induces large membrane fluctuations, which induce the contact of the membranes at an hourglass-like neck of the membrane tube. This contact results in rupture of the membrane and formation of a spherical vesicle. At , the rupture occurs at . With decreasing , larger values of are required for the rupture.

In a cylindrical tube of a homogeneous membrane, the bending energy yields an axial force


since an increase in the axial length results in a decrease in the cylindrical radius, i.e., an increase in the membrane mean curvature Shiba and Noguchi (2011). At , the force linearly decreases with increasing as shown in the inset of Fig. 1(e). This indicates that the inclusions are homogeneously mixed in the membrane. A similar linear dependence has previously been obtained for the density of anchored ideal-polymer chains Wu et al. (2013). In contrast, for , the curve deviates from a straight line due to the inclusion assembly. A larger shape change induces greater deviation [compare the inset of Fig. 1(e) with Figs. 1(g) and (h)].

iii.2 Vesicle

Similarly to the membrane tube, at and , the isotropic inclusions deform a vesicle into an oblate shape and the inclusions are concentrated in the two flatter regions [see Fig. 2(a)]. As increases, the vesicle becomes more spherical and the inclusions are distributed more uniformly [see Figs. 2(b) and (c)]. For at , vesicle division occurs and two spherical vesicles are formed.

These changes in the vesicle shape and inclusion assembly can be captured by the changes of the three principal lengths and inclusion distribution as shown in Figs. 2(d) and (e), respectively. The squared principal lengths , and are the eigenvalues of the gyration tensor, , where and is the center of mass. At , they are almost independent of so that the vesicle maintains its spherical shape. Note that the differences among three eigenvalues at are due to the thermal fluctuations. At , small deviations are recognized at and at , a clear decrease in is obtained. The inclusion distribution along the eigenvector of the smallest eigenvalue () is calculated as , where and are probabilities of finding the inclusions and all particles at each bin, respectively. The peaks of at both ends indicate the inclusion assembly on the flatter regions of the oblate vesicle [see the (red) solid line in Fig. 2(e)]. For , the inclusions are uniformly distributed.

Figure 3: (Color online) Snapshots of membrane tubes for (a) , (b) , (c) , and (d) at and . The front and side views are shown. The protein rod is displayed as a chain of spheres whose halves are colored in red (dark gray) and in yellow (light gray).
Figure 4: (Color online) Rod orientation and axial force of the membrane tube at . (a), (b) Rod curvature, , dependence of (a) the orientation degree and (b) axial force at , , , , and . The (black) dashed line in (a) shows the relation of for undeformable rods without thermal fluctuations. (c) Effective bending rigidity of the rods estimated by Eq. (10). Error bars are displayed at several or all data points in (a) and (b) or (c), respectively.
Figure 5: (Color online) Rod curvature, , dependence of Fourier amplitudes of (a), (c) membrane shape and (b), (d) the rod densities at and , , , , and . The Fourier amplitudes are normalized by the values at (denoted by the superscript ). Error bars are displayed at several data points.
Figure 6: (Color online) Phase diagrams of the membrane tube. (a) diagram at and . (b) diagram at and . (c) diagram at and .

Compared to the assembly of the anisotropic rods, as described in the next section, the inclusions assemble weakly and the shapes of the membrane tubes and vesicles fluctuate greatly, and for small values of , a single domain, instead of two domains, is occasionally formed. High bending rigidity of the inclusion is required for assembly in the absence of direct attraction between the inclusions.

Iv Membrane tube with protein rods

The protein rods exhibit a two-step assembly with increasing rod curvature (see Figs. 36). In our previous papers Noguchi (2014, 2015), we reported the assembly at . Here, we show the dependence and summarize the assembly processes.

At , the protein rods are oriented along the axial () direction and uniformly distributed in the membrane tube [see Fig. 3(a)]. As increases, the rod orientation changes to the azimuthal () direction and the mean orientational order parameter decreases from to , where the orientational order parameter is defined as [see Fig. 4(a)]. At larger , decreases more rapidly. For the limit of , i.e., undeformable rods, the curvature along the rod axis is exactly . If the thermal fluctuations are neglected, a linear relation is obtained for . As increases, approaches this linear relation, but relatively large deviations remain at . These deviations are likely due to membrane undulations, which allow orientation fluctuations even for undeformable rods.

The axial force behaves differently from the close-to-linear dependence of the isotropic inclusions [compare Fig. 4(b) with the inset of Fig. 1(e)]. During orientation changes, is almost constant for small , while it increases slightly for large . This increase may be due to entropy reduction by tilted rods, since the tilted rods suppress membrane undulation in both the and directions. In this region (), changes in the Fourier amplitudes are very small in both directions for all values of (see Fig. 5). Therefore, the membrane shapes and axial stress are modified only a little by the rods in this region.

With a further increase in , decreases linearly until the azimuthal assembly commences. When the rods are assumed to be completely oriented in the azimuthal direction, the axial force is given by Noguchi (2015),


The first term is the force in the absence of rods. The second term is the force generated by the rods, which is proportional to and linear with respect to . The effective bending rigidity of the rods is estimated from the slope of curves in the linear-decrease regions as shown in Fig. 4(c). increases with but is not linear with . This is because the orientation is not completely in the azimuthal direction and the interactions between rods and neighboring membrane particles are also involved in the bending deformation along the rod axis.

With an even further increase in , the rods assemble along the azimuthal direction and the membrane deforms into an elliptic tube [see Fig. 3(c)]. As increases more, the rod assembly also occurs along the axial direction [see Fig. 3(d)]. The increases in the Fourier amplitudes of the azimuthal () and axial () modes indicate the azimuthal and axial assemblies, respectively (see Fig. 5). With increasing , both assemblies occur at smaller [see Figs. 5 and 6(a)]. The curvatures of the azimuthal and axial assembly points are determined by the inflection points of and , respectively. The assembly is enhanced by the large rod stiffness in a manner similar to that of the isotropic inclusions. Thus, rod elasticity is one of the important factors that determine the assembly curvatures.

The phase diagram for the rod density and the tube radius are shown in Figs. 6(b) and (c), respectively. For the azimuthal assembly, gives an effect that is very similar to that of the rod stiffness . However, the axial assembly is different. At a large density (), axial assembly does not occur, since the elliptic edges are filled by the rods Noguchi (2015). As increases further, the membrane deforms into a triangular or other polygonal tube, instead of the elliptic tube, and the rods assemble at the edges of the polygonal tube Noguchi (2015). As increases, a slightly larger is needed for assembly in both directions. At a large tube radius (), axial assembly does not occur for the same reason as for a large . The length decreases as for the constant membrane area, so that the elliptic edges can be filled by a smaller number of rods.

Figure 7: (Color online) Membrane rupture at and . (a)–(d) Snapshots of vesicles. (a) Genus-11 vesicle at , , and . (b) Genus-1 vesicle at , , and . (c) Genus-5 vesicle at , , and . (d) Genus-0 vesicle at , , and . Membrane particles are displayed by nontransparent gray particles. (e)–(f) Mean number of the genus of vesicles as functions of the line tension . (e) , , , , and for . (f) , , and for .
Figure 8: (Color online) Membrane rupture at and . (a)–(c) Snapshots of vesicles. (a) Genus-11 vesicle at , , and . (b) Invaginated genus-0 vesicle at , , and . (c) Ruptured membrane at , , and . Membrane particles are displayed by nontransparent or transparent gray particles in (a) and (b) or (c), respectively, for clarity. (d)–(e) Mean number of the genus of vesicles as functions of the line tension . (d) , , and for . (e) , , , and for .
Figure 9: (Color online) Annealing time, , dependence of membrane rupture at , , , and . (a) Sequential snapshots of a vesicle for , , , , , and at and . Membrane particles are displayed by nontransparent gray particles. (b) Mean number of the genus as functions of for and .

Protein rods with large rod curvatures form a tight assembly in the membrane tubes as well as in the vesicles Noguchi (2014, 2015). It differs from the assembly of isotropic inclusions into flat membranes, which requires a larger bending elasticity. The rod assembly forms a saddle shape in the membrane tube. In contrast, such a saddle membrane is not stabilized by isotropic inclusions with a large positive spontaneous curvature since its mean curvature is small.

V Vesicle rupture

Membranes can rupture under a large stress. It is determined by the competition between the membrane deforming force (by the protein rods in this study) and the line tension of the membrane edge. For example, under a large positive surface tension , a membrane pore can expand until the pore radius reaches the stable value Tolpekina et al. (2004); Noguchi and Gompper (2006). A vesicle can spontaneously transform into a disk-shaped flat membrane when Fromherz (1983); Hu et al. (2012); Nakagawa and Noguchi (2015). Here, we use sufficiently large to prevent membrane rupture in the absence of rods.

Figures 79 show vesicle shapes resulting from membrane rupture by rod bending. Rupture changes vesicle topology. Figure 7 shows the vesicle shapes as the rod curvature changed suddenly from to . As the edge tension decreases, the membrane is ruptured more frequently and vesicles with higher genus are formed [see Figs. 7(e) and (f)]. As the rod stiffness or rod density increases, the rods deform the membrane more rapidly and membrane rupture is enhanced. For , depends linearly on while it is saturated at . A greater number of indicates the occurrence of more rupture. The obtained values of have a narrow distribution as indicated in small error bars that represent standard errors for samples. The saturation of is likely caused by the upper limit of genus, , which is determined by the vesicle size. A higher genus can be obtained for larger vesicles. The obtained high-genus vesicles, shown in Fig. 7(a), agree with the previously reported shapes in simulations and experiments Ayton et al. (2009); Simunovic et al. (2013b). Here, we have clarified that the edge tension and rod stiffness and density are important factors for vesicle rupture. At low , rod-free membranes are phase-separated and form flat regions while the rod assemblies form tubes and semicylindrical edges of the flat membranes [see Figs. 7(c) and (d)]. A similar coexistence of tubular networks and flat membranes is seen in endoplasmic reticulum (ER) Shibata et al. (2009). When changes slowly rather than rapidly, vesicle rupture is suppressed. The lower (red) line in Fig. 9(b) shows as a function of the annealing time , for which is changed linearly from to . No vesicles are ruptured at ( ms). Thus, the rapid adsorption of proteins onto the membrane is also important for obtaining vesicle rupture in experiments.

We also investigated vesicle rupture by protein rods with a curvature opposite to that of the vesicle curvature (see Figs. 8 and 9). At first, the rods induce many tubular invaginations into the inside of the vesicle [see Fig. 8(b)]. For large with relatively small , the membrane is not ruptured and densely-packed invaginations remain. Similar invaginations were observed by electron microscopy for liposomes with I-BAR proteins Mattila et al. (2007). For small or large , the membrane is ruptured and subsequently the inside surface of the vesicle turns to the outside [see Fig. 9(a)]. At high , this inversion (inside out) occurs completely and a high-genus vesicle is formed [see Fig. 8(a)]. However, at low , the inversion is only partial and rod-free membrane regions are not inverted, leading to division and partial connection of the membranes. In the membrane shown in Fig. 8(c), the rod-free membrane forms a pored vesicle and inverted membrane tubes partially remain in the vesicle. The ends of the tubes can be connected to the vesicle, but some of them are eventually pinched off. For genus estimation, complete membrane fission and partial connection are counted as and , respectively [ in the case of Fig. 8(c)]. Since the bending energy of an initially spherical vesicle is larger than the positive rod curvature with the same amplitude, membrane rupture occurs at smaller , , or larger . Even when is annealed slowly, high-genus vesicles are still formed [see Fig. 9(b)].

Vi Summary

We have investigated the shape transformation of vesicles and membrane tubes induced by protein rod assembly. As the rod curvature increases, the protein rods in the membrane tube assemble by two steps; first in the azimuthal direction and next in the longitudinal direction. These assemblies occur at lower for stiffer rods and/or higher rod density. Compared to the anisotropic rods, laterally isotropic inclusions assemble weakly. The isotropic inclusions assemble only when they induce very large bending rigidity locally; The inclusions with small spontaneous curvatures assemble in flatter regions of an elliptic membrane tube and oblate vesicle, while inclusions with large spontaneous curvatures induce membrane fission into vesicles. The rod–rod excluded-volume interactions make protein rods that are closer than align in the same orientation. Thus, the elongated shapes of membrane-reshaping proteins also assist protein assembly.

When the protein rods induce a large membrane stress, the membrane can be ruptured and high-genus vesicles form. We have clarified that membrane rupture is induced by the large bending stiffness of the rods, high density, rapid protein adhesion, and/or low line tension of the membrane edge. Thus, the choice of lipids is also important. When is negative with respect to the initial vesicle curvature, the membrane inversion also results in vesicle division as well as in the formation of high-genus vesicles. Our simulation results suggest that rapid exposure of liposomes to a protein solution is a key factor for observing high-genus liposomes due to protein-induced membrane rupture. Vesicle inversion was previously observed during lysis of a liposome by detergents Nomura et al. (2001). Here, we have demonstrated that the protein adhesion can also induce the vesicle inversion.

Here, we considered that the spontaneous (side) curvature and bending rigidity of the rods perpendicular to the rod axis are the same as the other membrane regions. Excluded-volume and other interactions between the proteins or between the protein and membrane can generate effective side spontaneous curvatures. When the rod and side curvatures are in opposite directions, saddle-shaped membranes, such as egg-carton Dommersnes and Fournier (1999, 2002), ring, and network structures Noguchi (2016), can be stabilized. Thus, anisotropic inclusions can induce much more variety in membrane structures than isotropic inclusions.

The replica exchange simulations were carried out on SGI Altix ICE 8400EX and ICE XA at ISSP Supercomputer Center, University of Tokyo. This work was partially supported by a Grant-in-Aid for Scientific Research on Innovative Areas “Fluctuation & Structure” (No. 25103010) from the Ministry of Education, Culture, Sports, Science, and Technology of Japan.


  1. H. T. McMahon and J. L. Gallop, Nature 438, 590 (2005).
  2. Y. Shibata, J. Hu, M. M. Kozlov, and T. A. Rapoport, Annu. Rev. Cell Dev. Biol. 25, 329 (2009).
  3. G. Drin and B. Antonny, FEBS Lett. 584, 1840 (2010).
  4. T. Baumgart, B. R. Capraro, C. Zhu, and S. L. Das, Annu. Rev. Phys. Chem. 62, 483 (2011).
  5. L. Johannes, C. Wunder, and P. Bassereau, Cold Spring Harbor Perspect. Biol. 6, a016741 (2014).
  6. H. T. McMahon and E. Boucrot, J. Cell Sci. 128, 1065 (2015).
  7. T. Itoh and P. De Camilli, Biochim. Biophys. Acta 1761, 897 (2006).
  8. M. Masuda and N. Mochizuki, Semin. Cell Dev. Biol. 21, 391 (2010).
  9. H. Zhao, A. Pykäläinen, and P. Lappalainen, Curr. Opin. Cell Biol. 23, 14 (2011).
  10. C. Mim and V. M. Unger, Trends Biochem. Sci. 37, 526 (2012).
  11. M. Simunovic, G. A. Voth, A. Callan-Jones, and P. Bassereau, Trends Cell Biol. 25, 780 (2015).
  12. B. J. Peter, H. M. Kent, I. G. Mills, Y. Vallis, P. J. G. Butler, P. R. Evans, and H. T. McMahon, Science 303, 495 (2004).
  13. P. K. Mattila, A. Pykäläinen, J. Saarikangas, V. O. Paavilainen, H. Vihinen, E. Jokitalo, and P. Lappalainen, J. Cell Biol. 176, 953 (2007).
  14. A. Frost, R. Perera, A. Roux, K. Spasov, O. Destaing, E. H. Egelman, P. De Camilli, and V. M. Unger, Cell 132, 807 (2008).
  15. Q. Wang, M. V. A. S. Navarro, G. Peng, E. Molinelli, S. L. Goh, B. L. Judson, K. R. Rajashankarc, and H. Sondermann, Proc. Natl. Acad. Sci. USA 106, 12700 (2009).
  16. C. Zhu, S. L. Das, and T. Baumgart, Biophys. J. 102, 1837 (2012).
  17. Y. Tanaka-Takiguchi, T. Itoh, K. Tsujita, S. Yamada, M. Yanagisawa, K. Fujiwara, A. Yamamoto, M. Ichikawa, and K. Takiguchi, Langmuir 29, 328 (2013).
  18. Z. Shi and T. Baumbart, Nature Comm. 6, 5974 (2015).
  19. C. Prévost, H. Zhao, J. Manzi, E. Lemichez, P. Lappalainen, A. Callan-Jones, and P. Bassereau, Nature Comm. 6, 8529 (2015).
  20. J. M. Isas, M. R. Ambroso, P. B. Hegde, J. Langen, and R. Langen, Structure 23, 873 (2015).
  21. J. Adam, N. Basnet, and N. Mizuno, Sci. Rep. 5, 15452 (2015).
  22. R. Lipowsky, Faraday Discuss. 161, 305 (2013).
  23. W. T. Góźdź, J. Phys. Chem. B 110, 21981 (2006).
  24. R. Phillips, T. Ursell, P. Wiggins, and P. Sens, Nature 459, 379 (2009).
  25. M. J. Greenall and G. Gompper, Langmuir 27, 3416 (2011).
  26. H. Noguchi, Soft Matter 8, 8926 (2012).
  27. S. Aimon, A. Callan-Jones, A. Berthaud, M. Pinot, G. E. Toombes, and P. Bassereau, Dev. Cell 28, 212 (2014).
  28. B. J. Reynwar, G. Ilya, V. A. Harmandaris, M. M. Müller, K. Kremer, and M. Deserno, Nature 447, 461 (2007).
  29. T. Auth and G. Gompper, Phys. Rev. E 80, 031901 (2009).
  30. A. Šarić and A. Cacciuto, Phys. Rev. Lett. 108, 118101 (2012).
  31. J. Gómez-Llobregat, F. Elías-Wolff, and M. Lindén, Biophys. J. 110, 197 (2016).
  32. P. B. Canham, J. Theor. Biol. 26, 61 (1970).
  33. W. Helfrich, Z. Naturforsch 28c, 693 (1973).
  34. J.-B. Fournier, Phys. Rev. Lett. 76, 4436 (1996).
  35. D. Kabaso, E. Gongadze, P. Elter, U. van Rienen, J. Gimsa, V. Kralj-Iglič, and A. Iglič, Mini Rev. Med. Chem. 11, 272 (2011).
  36. A. Iglič, H. Hägerstrand, P. Veranič, A. Plemenitaš, and V. Kralj-Iglič, J. Theor. Biol. 240, 368 (2006).
  37. P. G. Dommersnes and J. B. Fournier, Eur. Phys. J. B 12, 9 (1999).
  38. P. G. Dommersnes and J. B. Fournier, Biophys. J. 83, 2898 (2002).
  39. Y. Schweitzer and M. M. Kozlov, PLoS Comput. Biol. 11, e1004054 (2015).
  40. P. D. Blood and G. A. Voth, Proc. Natl. Acad. Sci. USA 103, 15068 (2006).
  41. A. Arkhipov, Y. Yin, and K. Schulten, Biophys. J. 95, 2806 (2008).
  42. G. Khelashvili, D. Harries, and H. Weinstein, Biophys. J. 97, 1626 (2009).
  43. H. Yu and K. Schulten, PLoS Comput. Biol. 9, e1002892 (2013).
  44. M. Simunovic, A. Srivastava, and G. A. Voth, Proc. Natl. Acad. Sci. USA 110, 20396 (2013a).
  45. M. Simunovic and G. A. Voth, Nature Comm. 6, 7219 (2015).
  46. N. Ramakrishnan, J. H. Ipsen, and P. B. Sunil Kumar, Soft Matter 8, 3058 (2012).
  47. N. Ramakrishnan, P. B. Sunil Kumar, and J. H. Ipsen, Biophys. J. 104, 1018 (2013).
  48. H. Noguchi, EPL 108, 48001 (2014).
  49. H. Noguchi, J. Chem. Phys. 143, 243109 (2015).
  50. H. Noguchi, Sci. Rep. 6, 20935 (2016).
  51. G. S. Ayton, E. Lyman, V. Krishna, R. D. Swenson, C. Mim, V. M. Unger, and G. A. Voth, Biophys. J. 97, 1616 (2009).
  52. M. Simunovic, C. Mim, T. C. Marlovits, G. Resch, V. M. Unger, and G. A. Voth, Biophys. J. 105, 711 (2013b).
  53. H. Noguchi, J. Phys. Soc. Jpn. 78, 041007 (2009).
  54. H. Noguchi and G. Gompper, Phys. Rev. E 73, 021903 (2006).
  55. H. Shiba and H. Noguchi, Phys. Rev. E 84, 031926 (2011).
  56. H. Noguchi, EPL 102, 68001 (2013).
  57. M. P. Allen and D. J. Tildesley, Computer Simulation of Liquids (Clarendon Press, Oxford, 1987).
  58. H. Noguchi, J. Chem. Phys. 134, 055101 (2011).
  59. S. Ramadurai, A. Holt, V. Krasnikov, G. van den Bogaart, J. A. Killian, and B. Poolman, J. Am. Chem. Soc. 131, 12650 (2009).
  60. K. Hukushima and K. Nemoto, J. Phys. Soc. Jpn. 65, 1604 (1996).
  61. Y. Okamoto, J. Mol. Graph. Model. 22, 425 (2004).
  62. H. Wu, H. Shiba, and H. Noguchi, Soft matter 9, 9907 (2013).
  63. T. V. Tolpekina, W. K. den Otter, and W. J. Briels, J. Chem. Phys. 121, 8014 (2004).
  64. P. Fromherz, Chem. Phys. Lett. 94, 259 (1983).
  65. M. Hu, J. J. Briguglio, and M. Deserno, Biophys. J. 102, 1403 (2012).
  66. K. M. Nakagawa and H. Noguchi, Soft matter 11, 1403 (2015).
  67. F. Nomura, M. Nagata, T. Inaba, H. Hiramatsu, H. Hotani, and K. Takiguchi, Proc. Natl. Acad. Sci. USA 98, 2340 (2001).
Comments 0
Request Comment
You are adding the first comment!
How to quickly get a good reply:
  • Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
  • Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
  • Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
The feedback must be of minumum 40 characters
Add comment
Loading ...
This is a comment super asjknd jkasnjk adsnkj
The feedback must be of minumum 40 characters
The feedback must be of minumum 40 characters

You are asking your first question!
How to quickly get a good answer:
  • Keep your question short and to the point
  • Check for grammar or spelling errors.
  • Phrase it like a question
Test description