Shape instabilities in vesicles: A phase-field model

A phase field model for dealing with shape instabilities in fluid membrane vesicles is presented. This model takes into account the Canham-Helfrich bending energy with spontaneous curvature. A dynamic equation for the phase-field is also derived. With this model it is possible to see the vesicle shape deformation dynamically, when some external agent instabilizes the membrane, for instance, inducing an inhomogeneous spontaneous curvature. The numerical scheme used is detailed and some stationary shapes are shown together with a shape diagram for vesicles of spherical topology and no spontaneous curvature, in agreement with known results.     

Simulating the dynamics of inextensible vesicles by the penalty immersed boundary method

In this paper, we develop an immersed boundary (IB) method to simulate the dynamics of inextensible vesicles interacting with an incompressible fluid. In order to take into account the inextensibility constraint of the vesicle, the penalty immersed boundary (pIB) method is used to virtually decouple the fluid and vesicle dynamics. As numerical tests of our current pIB method, the dynamics of single and multiple inextensible vesicles under shear flows have been extensively explored, and compared with the previous literature. The method is also validated by a series of convergence study, which confirms its consistent first-order accuracy on the velocity field, the vesicle configuration, the vesicle area and the perimeter errors. In addition, the method is also applied to study a binary-component vesicle problem.     

Domain-induced budding in buckling membranes

We present a phase field model on buckling membranes to analyze phase separation and budding on soft membranes. By numerically integrating dynamic equations, it turns out that the formation of caps is greatly influenced by the presence of a little excess area due to the surface area constraint. When cap-shaped domains are created, domain coalescence is mainly observed not between domains with same budding directions, but between domains with opposite budding directions, because the bending energy between two domains is larger in the former case. Although we do not introduce spontaneous curvature like Helfrich model, we obtain some suggestions related to the slow dynamics of the phase separation on vesicles.     

Adhesion of binary giant vesicles containing negative spontaneous curvature lipids induced by phase separation

We report the adhesion of binary giant vesicles composed of two types of phospholipids, one has negative spontaneous curvature which tends to bend toward the head group and the other has zero spontaneous curvature. In a homogeneous one-phase region, the giant vesicles do not adhere to each other, whereas in a coexisting two-phase region, the giant vesicles show adhesion. A fluorescence microscope observation reveals that the adhesion takes place through the domains rich in phospholipids having negative spontaneous curvature. We propose a phase separation induced hemifusion model where two apposed monolayers of adjacent vesicles are hemifused in order to reduce the bending energy of monolayers with negative spontaneous curvature and the boundary energy between the domains and matrix. We provide a strong evidence for the hemifusion model by lipid transfer experiments.     

Dynamic simulation of locally inextensible vesicles suspended in an arbitrary two-dimensional domain,a boundary integral method

We consider numerical algorithms for the simulation of hydrodynamics of two-dimensional vesicles suspended in a viscous Stokesian fluid. The motion of vesicles is governed by the interplay between hydrodynamic and elastic forces. Continuum models of vesicles use a two-phase fluid system with interfacial forces that include tension (to maintain local “surface” inextensibility) and bending. Vesicle flows are challenging to simulate. On the one hand, explicit time-stepping schemes suffer from a severe stability constraint due to the stiffness related to high-order spatial derivatives in the bending term. On the other hand, implicit time-stepping schemes can be expensive because they require the solution of a set of nonlinear equations at each time step.     

Advanced flicker spectroscopy of fluid membranes

The bending elasticity of a fluid membrane is characterized by its modulus and spontaneous curvature. We present a new method, advanced flicker spectroscopy of giant nonspherical vesicles, which makes it possible to simultaneously measure both parameters for the first time. Our analysis is based on the generation of a large set of reference data from Monte Carlo simulations of randomly triangulated surfaces. As an example of the potential of the procedure, we monitor thermal trajectories of vesicle shapes and discuss the elastic response of zwitterionic membranes to transmembrane pH gradients. Our technique makes it possible to easily characterize membrane curvature as a function of environmental conditions.     

Quadrupolar ordering of phospholipid molecules in narrow necks of phospholipid vesicles

Shapes of phospholipid vesicles that involve narrow neck(s) were studied theoretically. It is taken into account that phospholipid molecules are intrinsically anisotropic with respect to the membrane normal and that they exhibit quadrupolar orientational ordering according to the difference between the local principal membrane curvatures. Direct interactions between oriented molecules were considered within a linear approximation of the energy coupling with the deviatoric field. The equilibrium shapes of axisymmetric closed vesicles were studied by minimization of the free energy of the phospholipid bilayer membrane under relevant geometrical constraints. The variational problem was stated by a system of Euler-Lagrange differential equations that revealed a singularity in the derivative of the meridian curvature at points where the effect of the orientational ordering exactly counterbalances the effect of the isotropic bending. The system of Euler-Lagrange differential equations was solved numerically to yield consistently related equilibrium orientational distribution of the phospholipid molecules and vesicle shape. According to our estimation of the model constants the formation of the neck is promoted if direct interactions between the oriented molecules are taken into account. It was shown that the energy of the equilibrium shapes is considerably affected by the quadrupolar ordering of phospholipid molecules.     

Morphological phase diagram for lipid membrane domains with entropic tension

Circular domains in phase-separated lipid vesicles with symmetric leaflet composition commonly exhibit three stable morphologies: flat, dimpled, and budded. However, stable dimples (i.e., partially budded domains) present a puzzle since simple elastic theories of domain shape predict that only flat and spherical budded domains are mechanically stable in the absence of spontaneous curvature. We argue that this inconsistency arises from the failure of the constant surface tension ensemble to properly account for the effect of entropic bending fluctuations. Formulating membrane elasticity within an entropic tension ensemble, wherein tension represents the free energy cost of extracting membrane area from thermal bending of the membrane, we calculate a morphological phase diagram that contains regions of mechanical stability for each of the flat, dimpled, and budded domain morphologies.     

Accurate determination of elastic parameters for multicomponent membranes

Heterogeneities in the cell membrane due to coexisting lipid phases have been conjectured to play a major functional role in cell signaling and membrane trafficking. Thereby the material properties of multiphase systems, such as the line tension and the bending moduli, are crucially involved in the kinetics and the asymptotic behavior of phase separation. In this Letter we present a combined analytical and experimental approach to determine the properties of phase-separated vesicle systems. First we develop an analytical model for the vesicle shape of weakly budded biphasic vesicles. Subsequently experimental data on vesicle shape and membrane fluctuations are taken and compared to the model. The parameters obtained set limits for the size and stability of nanodomains in the plasma membrane of living cells.     

Configurations of fluid membranes and vesicles

Vesicles consisting of a bilayer membrane of amphiphilic lipid molecules are remarkably flexible surfaces that show an amazing variety of shapes of different symmetry and topology. Owing to the fluidity of the membrane, shape transitions such as budding can be induced by temperature changes or the action of optical tweezers. Thermally excited shape fluctuations are both strong and slow enough to be visible by video microscopy. Depending on the physical conditions, vesicles adhere to and unbind from each other or a substrate.

This article describes the systematic physical theory developed to understand the static and dynamic aspects of membrane and vesicle configurations. The preferred shapes arise from a competition between curvature energy, which derives from the bending elasticity of the membrane, geometrical constraints such as fixed surface area and fixed enclosed volume, and a signature of the bilayer aspect. These shapes of lowest energy are arranged into phase diagrams, which separate regions of different symmetry by continuous or discontinuous transitions. The geometrical constraints affect the fluctuations around these shapes by creating an effective tension.

For vesicles of non-spherical topology, the conformal invariance of the curvature energy leads to conformal diffusion, which signifies a one-fold degeneracy of the ground state. Unbinding and adhesion transitions arise from the balance between attractive interactions and entropic repulsion or a cost in bending energy, respectively. Both the dynamics of equilibrium fluctuations and the dynamics of shape transformations are governed not only by viscous damping in the surrounding liquid but also by internal friction if the two monolayers slip over each other. More complex membranes such as that of the red blood cell exhibit a variety of new phenomena because of coupling between internal degrees of freedom and external geometry.     

The role of Gauss curvature in a membrane phase separation problem

Consider a two-phase lipid vesicle. Below the transition temperature, the phases separate into non-connecting domains that coarsen into larger areas. The free energy of phase properties determines the length of the boundaries separating the regions. The two phases correspond to different lipid compositions, and in cells, this fluctuation in composition is a dynamic process vital to its function. We prove that a small patch of the minority lipids forms at a point of the membrane where the Gauss curvature attains a maximum. This patch has a round shape approximately and its boundary has a constant geodesic curvature. The proof consists of three steps. The construction of a family of good approximate solutions, an improvement of the approximate solutions so that their geodesic curvature is a constant modulo translation, and the identification of an exact solution from the family of the improved approximate solutions. Our theoretical results are supported by vesicle experiments.     

Simulations of flexible membranes using a coarse-grained particle-based model with spontaneous curvature variables

We propose a simple mathematical model for lipid bilayer membranes of flexible fluid or solid state. The model consists of interacting coarse-grained particles with two extra variables. One indicates the spontaneous curvature at the particle position, and the other indicates the vector representing the direction normal to the membrane. When the spontaneous curvature variable is allowed to fluctuate significantly, the fluctuation causes a softening of the membrane and growth of large undulations as the amplitude of the fluctuation is increased. By changing the amplitude of the fluctuation in simulations, the bending rigidity of the membrane can be easily controlled. Because the proposed model includes anisotropic interactions between the particles, multilayered vesicles can be obtained through a reversible transition by weakening the strength of the anisotropic interactions.     

Cellular uptake of elastic nanoparticles

A fundamental understanding of cell-nanomaterial interaction is of essential importance to nanomedicine and safe applications of nanotechnology. Here we investigate the adhesive wrapping of a soft elastic vesicle by a lipid membrane. We show that there exist a maximum of five distinct wrapping phases based on the stability of full wrapping, partial wrapping, and no wrapping states. The wrapping phases depend on the vesicle size, adhesion energy, surface tension of membrane, and bending rigidity ratio between vesicle and membrane. These results are of immediate interest to the study of vesicular transport and endocytosis or phagocytosis of elastic particles into cells.     

Elastic energy of tilt and bending of fluid membranes

Tilt of hydrocarbon chains of lipid molecules with respect to membrane plane is commonly considered to characterize the internal structure of a membrane in the crystalline state. However, membranes in the liquid state can also exhibit tilt resulting from packing constraints imposed on the lipid molecules in diverse biologically relevant structures such as intermediates of membrane fusion, pores in lipid bilayers and others. We analyze the energetics of tilt in liquid membranes and its coupling with membrane bending. We consider three contributions to the elastic energy: constant tilt, variation of tilt along the membrane surface and membrane bending. The major assumption of the model is that the core of a liquid membrane has the common properties of an elastic continuum. We show that the variation of tilt and membrane bending are additive and that their energy contributions are determined by the same elastic coefficient: the Helfrich bending modulus, the modulus of Gaussian curvature and the spontaneous curvature known from previous studies of pure bending. The energy of a combined deformation of bending and varying tilt is determined by an effective tensor accounting for the two factors. In contrast, the deformation of constant tilt does not couple with bending and its contribution to the elastic energy is determined by an independent elastic constant. While accurate determination of this constant requires additional measurements, we estimate its value using a simplified approach. We discuss the relationships between the obtained elastic Hamiltonian of a membrane and the previous models of membrane elasticity. Received 10 February 2000 and Received in final form 19 June 2000     

Adhesion of polymer vesicles

The adhesion and bending modulus of polybutadiene-poly(ethylene oxide) block copolymer vesicles made from a bidisperse mixture of polymers is measured using micropipette aspiration. The adhesion energy between biotinylated vesicles and avidin beads is modeled by incorporating the extension of the adhesive ligands above the surface brush of the vesicle according to the blob model of bidisperse polymer mixtures of Komura and Safran assuming the polymer brush at the surface of the vesicle is compact. The same model accurately reproduces the scaling of the bending modulus with polymer composition.     

A boundary integral method for simulating the dynamics of inextensible vesicles suspended in a viscous fluid in 2D

We present a new method for the evolution of inextensible vesicles immersed in a Stokesian fluid. We use a boundary integral formulation for the fluid that results in a set of nonlinear integro-differential equations for the vesicle dynamics. The motion of the vesicles is determined by balancing the non-local hydrodynamic forces with the elastic forces due to bending and tension. Numerical simulations of such vesicle motions are quite challenging. On one hand, explicit time-stepping schemes suffer from a severe stability constraint due to the stiffness related to high-order spatial derivatives and a milder constraint due to a transport-like stability condition. On the other hand, an implicit scheme can be expensive because it requires the solution of a set of nonlinear equations at each time step. We present two semi-implicit schemes that circumvent the severe stability constraints on the time step and whose computational cost per time step is comparable to that of an explicit scheme. We discretize the equations by using a spectral method in space, and a multistep third-order accurate scheme in time. We use the fast multipole method (FMM) to efficiently compute vesicle–vesicle interaction forces in a suspension with a large number of vesicles. We report results from numerical experiments that demonstrate the convergence and algorithmic complexity properties of our scheme.     

Adhesion of fluid vesicles at chemically structured substrates

The adhesion of fluid vesicles at chemically structured substrates is studied theoretically via Monte Carlo simulations. The substrate surface is planar and repels the vesicle membrane apart from a single surface domain γ , which strongly attracts this membrane. If the vesicle is larger than the attractive γ domain, the spreading of the vesicle onto the substrate is restricted by the size of this surface domain. Once the contact line of the adhering vesicle has reached the boundaries of the γ domain, further deflation of the vesicle leads to a regime of low membrane tension with pronounced shape fluctuations, which are now governed by the bending rigidity. For a circular γ domain and a small bending rigidity, the membrane oscillates strongly around an average spherical cap shape. If such a vesicle is deflated, the contact area increases or decreases with increasing osmotic pressure, depending on the relative size of the vesicle and the circular γ domain. The lateral localization of the vesicle's center of mass by such a domain is optimal for a certain domain radius, which is found to be rather independent of adhesion strength and bending rigidity. For vesicles adhering to stripe-shaped surface domains, the width of the contact area perpendicular to the stripe varies nonmonotonically with the adhesion strength.     

Aggregates of two-dimensional vesicles: rouleaux, sheets, and convergent extension

Using both numerical and variational minimization of the bending and adhesion energy of two-dimensional lipid vesicles, we study their aggregation, and we find that the stable aggregates include an infinite number of vesicles and that they arrange either in a columnar or in a sheetlike structure. We calculate the stability diagram and we show that the sheetlike aggregate can be transformed into the columnar aggregate via vesicle intercalation, which makes the transformation reminiscent of the process of convergent extension observed in some biological tissues.     

Dynamic model and stationary shapes of fluid vesicles

A phase-field model that takes into account the bending energy of fluid vesicles is presented. The Canham-Helfrich model is derived in the sharp-interface limit. A dynamic equation for the phase-field has been solved numerically to find stationary shapes of vesicles with different topologies and the dynamic evolution towards them. The results are in agreement with those found by minimization of the Canham-Helfrich free energy. This fact shows that our phase-field model could be applied to more complex problems of instabilities.