Shapes of bilayer vesicles with membrane embedded molecules

The interdependence of the lateral distribution of molecules which are embedded in a membrane (such as integral membrane proteins) and the shape of a cell with no internal structure (such as phospholipid vesicles or mammalian erythrocytes) has been studied. The coupling of the lateral distribution of the molecules and the cell shape is introduced by considering that the energy of the membrane embedded molecule at a given site of the membrane depends on the curvature of the membrane at that site. Direct interactions between embedded molecules are not considered. A simple expression for the interaction of the membrane embedded molecule with the local membrane curvature is proposed. Starting from this interaction, the consistently related expressions for the free energy and for the distribution function of the embedded molecules are derived. The equilibrium cell shape and the corresponding lateral distribution of the membrane embedded molecules are determined by minimization of the membrane free energy which includes the free energy of the membrane embedded molecules and the membrane elastic energy. The resulting inhomogeneous distribution of the membrane embedded molecules affects the cell shape in a nontrivial manner. In particular, it is shown that the shape corresponding to the absolute energy minimum at given cell volume and membrane area may be elliptically non-axisymmetric, in contrast to the case of a laterally homogeneous membrane where it is axisymmetric.     

Evidence of cholesterol accumulated in high curvature regions: implication to the curvature elastic energy for lipid mixtures

Recent experiments suggested that cholesterol and other lipid components of high negative spontaneous curvature facilitate membrane fusion. This is taken as evidence supporting the stalk-pore model of membrane fusion in which the lipid bilayers go through intermediate structures of high curvature. How do the high-curvature lipid components lower the free energy of the curved structure? Do the high-curvature lipid components modify the average spontaneous curvature of the relevant monolayer, thereby facilitate its bending, or do the lipid components redistribute in the curved structure so as to lower the free energy? This question is fundamental to the curvature elastic energy for lipid mixtures. Here we investigate the lipid distribution in a monolayer of a binary lipid mixture before and after bending, or more precisely in the lamellar, hexagonal, and distorted hexagonal phases. The lipid mixture is composed of 2:1 ratio of brominated di18:0PC and cholesterol. Using a newly developed procedure for the multiwavelength anomalous diffraction method, we are able to isolate the bromine distribution and reconstruct the electron density distribution of the lipid mixture in the three phases. We found that the lipid distribution is homogenous and uniform in the lamellar and hexagonal phases. But in the distorted hexagonal phase, the lipid monolayer has nonuniform curvature, and cholesterol almost entirely concentrates in the high curvature region. This finding demonstrates that the association energies between lipid molecules vary with the curvature of membrane. Thus, lipid components in a mixture may redistribute under conditions of nonuniform curvature, such as in the stalk structure. In such cases, the spontaneous curvature depends on the local lipid composition and the free energy minimum is determined by lipid distribution as well as curvature.     

Electrostatic free energy and spontaneous curvature of spherical charged layered membrane

A biophysical model for the equilibrium curvature of a composite membrane element is derived taking into account the mechanical bilayer properties and the adjacent charged protein layers. The minimum of the total free energy density with respect to the curvature of such a membrane curved was estimated from the sum of the electrostatic free energy density of the charges of the membrane and the elastic surface energy density due to bending the lipid bilayer membrane. It was shown that the equilibrium curvature, i.e. the spontaneous curvature, of such a charged composite sandwich-like membrane depends inversely on the bending stiffness of the lipid membrane itself and directly on the charge amount inside and outside the membrane to the second power. Furthermore the geometric and electrostatic structure of the protein layers and the physico-chemical environment conditions are involved. Corresponding to the model developed a "standard RBC" membrane element has a negative spontaneous curvature, accounting for a discocyte RBC shape. The shape change from a discocyte to a more stomatocytic shape (increase in the negative spontaneous curvature) after reducing the charges in the glycocalyx is also explained within this model.     

Exploring the binding dynamics of BAR proteins

We used a continuum model based on the Helfrich free energy to investigate the binding dynamics of a lipid bilayer to a BAR domain surface of a crescent-like shape of positive (e.g. I-BAR shape) or negative (e.g. F-BAR shape) intrinsic curvature. According to structural data, it has been suggested that negatively charged membrane lipids are bound to positively charged amino acids at the binding interface of BAR proteins, contributing a negative binding energy to the system free energy. In addition, the cone-like shape of negatively charged lipids on the inner side of a cell membrane might contribute a positive intrinsic curvature, facilitating the initial bending towards the crescent-like shape of the BAR domain. In the present study, we hypothesize that in the limit of a rigid BAR domain shape, the negative binding energy and the coupling between the intrinsic curvature of negatively charged lipids and the membrane curvature drive the bending of the membrane. To estimate the binding energy, the electric potential at the charged surface of a BAR domain was calculated using the Langevin-Bikerman equation. Results of numerical simulations reveal that the binding energy is important for the initial instability (i.e. bending of a membrane), while the coupling between the intrinsic shapes of lipids and membrane curvature could be crucial for the curvature-dependent aggregation of negatively charged lipids near the surface of the BAR domain. In the discussion, we suggest novel experiments using patch clamp techniques to analyze the binding dynamics of BAR proteins, as well as the possible role of BAR proteins in the fusion pore stability of exovesicles.     

Protein-Induced Membrane Curvature Alters Local Membrane Tension

Adsorption of proteins onto membranes can alter the local membrane curvature. This phenomenon has been observed in biological processes such as endocytosis, tubulation, and vesiculation. However, it is not clear how the local surface properties of the membrane, such as membrane tension, change in response to protein adsorption. In this article, we show that the partial differential equations arising from classical elastic model of lipid membranes, which account for simultaneous changes in shape and membrane tension due to protein adsorption in a local region, cannot be solved for nonaxisymmetric geometries using straightforward numerical techniques; instead, a viscous-elastic formulation is necessary to fully describe the system. Therefore, we develop a viscous-elastic model for inhomogeneous membranes of the Helfrich type. Using the newly available viscous-elastic model, we find that the lipids flow to accommodate changes in membrane curvature during protein adsorption. We show that, at the end of protein adsorption process, the system sustains a residual local tension to balance the difference between the actual mean curvature and the imposed spontaneous curvature. We also show that this change in membrane tension can have a functional impact such as altered response to pulling forces in the presence of proteins.     

Protein-Induced Membrane Curvature Alters Local Membrane Tension

Adsorption of proteins onto membranes can alter the local membrane curvature. This phenomenon has been observed in biological processes such as endocytosis, tubulation, and vesiculation. However, it is not clear how the local surface properties of the membrane, such as membrane tension, change in response to protein adsorption. In this article, we show that the partial differential equations arising from classical elastic model of lipid membranes, which account for simultaneous changes in shape and membrane tension due to protein adsorption in a local region, cannot be solved for nonaxisymmetric geometries using straightforward numerical techniques; instead, a viscous-elastic formulation is necessary to fully describe the system. Therefore, we develop a viscous-elastic model for inhomogeneous membranes of the Helfrich type. Using the newly available viscous-elastic model, we find that the lipids flow to accommodate changes in membrane curvature during protein adsorption. We show that, at the end of protein adsorption process, the system sustains a residual local tension to balance the difference between the actual mean curvature and the imposed spontaneous curvature. We also show that this change in membrane tension can have a functional impact such as altered response to pulling forces in the presence of proteins.     

Micron-scale plasma membrane curvature is recognized by the septin cytoskeleton

Cells change shape in response to diverse environmental and developmental conditions, creating topologies with micron-scale features. Although individual proteins can sense nanometer-scale membrane curvature, it is unclear if a cell could also use nanometer-scale components to sense micron-scale contours, such as the cytokinetic furrow and base of neuronal branches. Septins are filament-forming proteins that serve as signaling platforms and are frequently associated with areas of the plasma membrane where there is micron-scale curvature, including the cytokinetic furrow and the base of cell protrusions. We report here that fungal and human septins are able to distinguish between different degrees of micron-scale curvature in cells. By preparing supported lipid bilayers on beads of different curvature, we reconstitute and measure the intrinsic septin curvature preference. We conclude that micron-scale curvature recognition is a fundamental property of the septin cytoskeleton that provides the cell with a mechanism to know its local shape.     

Crenation and cupping of the red cell: A new theoretical approach. Part II. Cupping

Typical, axisymmetrical cup shaped cells have been carefully measured and the shapes analyzed mathematically. The results show that the strain energy of a cup shaped cell is always higher than that of a biconcave cell except when the two layers of the membrane involved in resistance to bending are free to slide over one another. This is true whether intrinsic curvature of the membrane is positive, negative or zero. If the two layers can slide over one another, the cup shape becomes the lower energy form. Shear resistance, if appreciable, must cause the cup cell to buckle. Photomicrographs of cup shaped cells show buckled configurations characteristic of those of a partly deflated thin-walled rubber ball, which is a similar object having a low ratio of bending/shear strength.In light of these findings, the cup shape of the red cell can no longer be considered as evidence of intrinsic membrane curvature of opposite sign to that of the crenated cell, but appears to indicate a phase change either in the hydrophobic interior of the bimolecular membrane or in some equivalent interface.     

Minimum energy analysis of membrane deformation applied to pipet aspiration and surface adhesion of red blood cells.

An experimental procedure is demonstrated which can be used to determine the interfacial free energy density for red cell membrane adhesion and membrane elastic properties. The experiment involves micropipet aspiration of a flaccid red blood cell and manipulation of the cell proximal to a surface where adhesion occurs. A minimum free energy method is developed to model the equilibrium contour of unsupported membrane regions and to evaluate the partial derivatives of the total free energy, which correspond to the micropipet suction force and the interfacial free energy density of adhesion. It is shown that the bending elasticity of the red cell membrane does not contribute significantly to the pressure required to aspirate a flaccid red cell. Based on experimental evidence, the upper bound for the bending or curvature elastic modulus of the red cell membranes is 10-12 ergs (dyn-cm). Analysis of the adhesion experiment shows that interfacial free energy densities for red cell adhesion can be measured from a lower limit of 10-4 ergs/cm2 to an upper limit established by the membrane tension for lysis of 5-10 ergs/cm2.     

Erythrocyte shape simulation by numerical optimization

In a recent paper we examined the morphology of erythrocytes in terms of the mean mean curvature (MMC) of their cell membranes. A computer simulation of these shapes based on the different geometries showed that the MMC increased from the sphero-stomatocyte to the spheroechinocyte via the discocyte. In this work we extend this analysis by using a numerical optimization method based on importance sampling and the principle of adiabatic cooling. The erythrocyte membrane is treated as a single closed fluid lamina exhibiting viscoelastic characteristics. The energy function of the lamina includes the following terms: (i) Curvature-elastic energy terms which depend on both local and global curvature. (ii) A term describing the compression elasticity of the lamina. (iii) A term which depends on the volume of the cell and which is related to the osmotic pressure across the membrane. In the simulation the cell is assumed to have axial symmetry and it can therefore be described by a finite set of conic sections. So far we have been able to obtain an energy minimum corresponding to a discocyte shape using a sphere as the initial configuration. Our results therefore imply that the well-known sequence of erythrocyte shapes could solely be governed by the above mentioned properties of an ideal fluid forming a closed singly connected lamina.     

Molecular theory of lipid-protein interaction and the Lalpha-HII transition.

We present a molecular-level theory for lipid-protein interaction and apply it to the study of lipid-mediated interactions between proteins and the protein-induced transition from the planar bilayer (Lalpha) to the inverse-hexagonal (HII) phase. The proteins are treated as rigid, membrane-spanning, hydrophobic inclusions of different size and shape, e.g., "cylinder-like," "barrel-like," or "vase-like." We assume strong hydrophobic coupling between the protein and its neighbor lipids. This means that, if necessary, the flexible lipid chains surrounding the protein will stretch, compress, and/or tilt to bridge the hydrophobic thickness mismatch between the protein and the unperturbed bilayer. The system free energy is expressed as an integral over local molecular contributions, the latter accounting for interheadgroup repulsion, hydrocarbon-water surface energy, and chain stretching-tilting effects. We show that the molecular interaction constants are intimately related to familiar elastic (continuum) characteristics of the membrane, such as the bending rigidity and spontaneous curvature, as well as to the less familiar tilt modulus. The equilibrium configuration of the membrane is determined by minimizing the free energy functional, subject to boundary conditions dictated by the size, shape, and spatial distribution of inclusions. A similar procedure is used to calculate the free energy and structure of peptide-free and peptide-rich hexagonal phases. Two degrees of freedom are involved in the variational minimization procedure: the local length and local tilt angle of the lipid chains. The inclusion of chain tilt is particularly important for studying noncylindrical (for instance, barrel-like) inclusions and analyzing the structure of the HII lipid phase; e.g., we find that chain tilt relaxation implies strong faceting of the lipid monolayers in the hexagonal phase. Consistent with experiment, we find that only short peptides (large negative mismatch) can induce the Lalpha --> HII transition. At the transition, a peptide-poor Lalpha phase coexists with a peptide-rich HII phase.     

Interaction between bending and tension forces in bilayer membranes.

A theoretical analysis is presented of the bending mechanics of a membrane consisting of two tightly-coupled leaflets, each of which shears and bends readily but strongly resists area changes. Structures of this type have been proposed to model biological membranes such as red blood cell membrane. It is shown that when such a membrane is bent, anisotropic components of resultant membrane tension (shear stresses) are induced, even when the tension in each leaflet is isotropic. The induced shear stresses increase as the square of the membrane curvature, and become significant for moderate curvatures (when the radius of curvature is much larger than the distance between the leaflets). This effect has implications for the analysis of shape and deformation of freely suspended and flowing red blood cells.     

Role of the membrane cortex in neutrophil deformation in small pipets.

The simplest model for a neutrophil in its "passive" state views the cell as consisting of a liquid-like cytoplasmic region surrounded by a membrane. The cell surface is in a state of isotropic contraction, which causes the cell to assume a spherical shape. This contraction is characterized by the cortical tension. The cortical tension shows a weak area dilation dependence, and it determines the elastic properties of the cell for small curvature deformations. At high curvature deformations in small pipets (with internal radii less than 1 micron), the measured critical suction pressure for cell flow into the pipet is larger than its estimate from the law of Laplace. A model is proposed where the region consisting of the cytoplasm membrane and the underlying cortex (having a finite thickness) is introduced at the cell surface. The mechanical properties of this region are characterized by the apparent cortical tension (defined as a free contraction energy per unit area) and the apparent bending modulus (introduced as a bending free energy per unit area) of its middle plane. The model predicts that for small curvature deformations (in pipets having radii larger than 1.2 microns) the role of the cortical thickness and the resistance for bending of the membrane-cortex complex is negligible. For high curvature deformations, they lead to elevated suction pressures above the values predicted from the law of Laplace. The existence of elevated suction pressures for pipets with radii from 1 micron down to 0.24 micron is found experimentally. The measured excess suction pressures cannot be explained only by the modified law of Laplace (for a cortex with finite thickness and negligible bending resistance), because it predicts unacceptable high cortical thicknesses (from 0.3 to 0.7 micron). It is concluded that the membrane-cortex complex has an apparent bending modulus from 1 x 10(-18) to 2 x 10(-18) J for a cortex with a thickness from 0.1 micron down to values much smaller than the radius of the smallest pipet (0.24 micron) used in this study.     

Theoretical model of reticulocyte to erythrocyte shape transformation

A theoretical model describing the kinetics of reticulocyte shape transformation was developed. The model considers the evolution of a simple cellular shape under transmembrane pressure difference, and proposes a four-parameter axisymmetric approximation of the cell surface. The mathematical analysis considers plasma membrane tension in the plane of bilayer leaflets, membrane spontaneous curvature and transmembrane transport of water. Cytoskeleton dilatational and shear rigidity, and the energetic barrier preventing the decrease of cell volume below a certain minimum are also incorporated. The set of adequate physical assumptions allowed for formulation of the equation for free energy of the investigated system. Computer simulations of cell shape changes, down to the state of free energy minimum, together with estimation of the time needed for the resulting transport of water, revealed a complex, three-phase picture of temporal alterations in cellular geometry with a wide spectrum of final results, and led to propose a standard model of reticulocyte-erythrocyte transformation. According to the model, both cell volume and surface undergo changes, and the work of the pressure, initially accumulated in the cytoskeleton, is consumed for local bending of the cell membrane. Further simulations with modified initial shape or parameters of the standard model show the trajectories of system evolution and help in better understanding the conditions for the erythro-, sphero-, ovalo-, stomato-, and leptoidal metamorphosis of maturing red blood cells. The stability of the final biconcave shape was also verified. Spherogenic modifications were discussed in the context of spherocytosis. Future development of the model was proposed.     

Mechanics of curved plasma membrane vesicles: resting shapes, membrane curvature, and in-plane shear elasticity

Highly curved cell membrane structures, such as plasmalemmal vesicles (caveolae) and clathrin-coated pits, facilitate many cell functions, including the clustering of membrane receptors and transport of specific extracellular macromolecules by endothelial cells. These structures are subject to large mechanical deformations when the plasma membrane is stretched and subject to a change of its curvature. To enhance our understanding of plasmalemmal vesicles we need to improve the understanding of the mechanics in regions of high membrane curvatures. We examine here, theoretically, the shapes of plasmalemmal vesicles assuming that they consist of three membrane domains: an inner domain with high curvature, an outer domain with moderate curvature, and an outermost flat domain, all in the unstressed state. We assume the membrane properties are the same in these domains with membrane bending elasticity as well as in-plane shear elasticity. Special emphasis is placed on the effects of membrane curvature and in-plane shear elasticity on the mechanics of vesicle during unfolding by application of membrane tension. The vesicle shapes were computed by minimization of bending and in-plane shear strain energy. Mechanically stable vesicles were identified with characteristic membrane necks. Upon stretch of the membrane, the vesicle necks disappeared relatively abruptly leading to membrane shapes that consist of curved indentations. While the resting shape of vesicles is predominantly affected by the membrane spontaneous curvatures, the membrane shear elasticity (for a range of values recorded in the red cell membrane) makes a significant contribution as the vesicle is subject to stretch and unfolding. The membrane tension required to unfold the vesicle is sensitive with respect to its shape, especially as the vesicle becomes fully unfolded and approaches a relative flat shape.     

Elastic energy of curvature-driven bump formation on red blood cell membrane.

Model calculations were performed to explore quantitative aspects of the discocyte-echinocyte shape transformation in red blood cells. The shape transformation was assumed to be driven by changes in the preferred curvature of the membrane bilayer and opposed by the elastic shear rigidity of the membrane skeleton. The energy required for echinocyte bump formation was calculated for a range of bump shapes for different preferred curvatures. Energy minima corresponding to nonzero bump heights were found when the stress-free area difference between the membrane leaflets or the spontaneous curvature of the membrane became sufficiently large, but the calculations predict that the membrane can tolerate significant differences in the resting areas of the inner and outer leaflets or significant spontaneous curvature without visible changes in shape. Thus, if the cell is near the threshold for bump formation, the calculations predict that small changes in membrane properties would produce large changes in cellular geometry. These results provide a rational framework for interpreting observations of shape transformations in red cells and for understanding the mechanism by which small changes in membrane elastic properties might lead to significant changes in geometry.     

Lateral redistribution of transmembrane proteins and liquid-ordered domains in lipid membranes with inhomogeneous curvature

A number of processes in living cells are accompanied by significant changes of the geometric curvature of lipid membranes. In turn, heterogeneity of the lateral curvature can lead to spatial redistribution of membrane components, most important of which are transmembrane proteins and liquid-ordered lipid-protein domains. These components have a so-called hydrophobic mismatch: the length of the transmembrane domain of the protein, or the thickness of the bilayer of the domain differ from the thickness of the surrounding membrane. In this work we consider redistribution of membrane components with hydrophobic mismatch in membranes with non-uniform geometric curvature. Dependence of the components’ energy on the curvature is calculated in terms of theory of elasticity of liquid crystals adapted to lipid membranes. According to the calculations, transmembrane proteins prefer regions of the membrane with zero curvature. Liquid-ordered domains having a size of a few nm distribute mainly into regions of the membrane with small negative curvature appearing in the cell plasma membrane in the process of endocytosis. The distribution of domains of a large radius is determined by a decrease of their perimeter upon bending; these domains distribute into membrane regions with relatively large curvature.     

Curling and Local Shape Changes of Red Blood Cell Membranes Driven by Cytoskeletal Reorganization

Human red blood cells (RBCs) lack the actin-myosin-microtubule cytoskeleton that is responsible for shape changes in other cells. Nevertheless, they can display highly dynamic local deformations in response to external perturbations, such as those that occur during the process of apical alignment preceding merozoite invasion in malaria. Moreover, after lysis in divalent cation-free media, the isolated membranes of ruptured ghosts show spontaneous inside-out curling motions at the free edges of the lytic hole, leading to inside-out vesiculation. The molecular mechanisms that drive these rapid shape changes are unknown. Here, we propose a molecular model in which the spectrin filaments of the RBC cortical cytoskeleton control the sign and dynamics of membrane curvature depending on whether the ends of the filaments are free or anchored to the bilayer. Computer simulations of the model reveal that curling, as experimentally observed, can be obtained either by an overall excess of weakly-bound filaments throughout the cell, or by the flux of such filaments toward the curling edges. Divalent cations have been shown to arrest the curling process, and Ca²⁺ ions have also been implicated in local membrane deformations during merozoite invasion. These effects can be replicated in our model by attributing the divalent cation effects to increased filament-membrane binding. This process converts the curl-inducing loose filaments into fully bound filaments that arrest curling. The same basic mechanism can be shown to account for Ca²⁺-induced local and dynamic membrane deformations in intact RBCs. The implications of these results in terms of RBC membrane dynamics under physiological, pathological, and experimental conditions is discussed.     

Mechanical stability of membrane nanotubular protrusions influenced by attachment of flexible rod-like proteins

It is indicated that nonhomogeneous lateral distribution of membrane attached and flexible rod-like proteins (MRPs) may stabilize nanotubular membrane protrusions. We have shown that curvature induced accumulation of MRPs in the nanotubular membrane protrusion and the corresponding reduction of the membrane free energy are possible if the decrease of the deviatoric free energy of MRPs in the nanotubular protrusions is large enough to overcome the increase of the free energy due to decrease of configurational entropy in the process of lateral sorting of MRPs. The decrease of isotropic curvature energy of MRPs in the region of membrane protrusion is usually not sufficient for substantial MRPs sorting and consequent stabilization of the nanotubular membrane protrusions.