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.     

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.     

Statistical thermodynamics of membrane bending-mediated protein-protein attractions

Highly wedge-shaped integral membrane proteins, or membrane-adsorbed proteins can induce long-ranged deformations. The strain in the surrounding bilayer creates relatively long-ranged forces that contribute to interactions with nearby proteins. In contrast, to direct short-ranged interactions such as van der Waal's, hydrophobic, or electrostatic interactions, both local membrane Gaussian curvature and protein ellipticity can induce forces acting at distances of up to a few times their typical radii. These forces can be attractive or repulsive, depending on the proteins' shape, height, contact angle with the bilayer, and a pre-existing local membrane curvature. Although interaction energies are not pairwise additive, for sufficiently low protein density, thermodynamic properties depend only upon pair interactions. Here, we compute pair interaction potentials and entropic contributions to the two-dimensional osmotic pressure of a collection of noncircular proteins. For flat membranes, bending rigidities of approximately 100k(B)T, moderate ellipticities, and large contact angle proteins, we find thermally averaged attractive interactions of order k(B)T. These interactions may play an important role in the intermediate stages of protein aggregation. Numerous biological processes where membrane bending-mediated interactions may be relevant are cited, and possible experiments are discussed.     

Membrane-Elasticity Model of Coatless Vesicle Budding Induced by ESCRT Complexes

The formation of vesicles is essential for many biological processes, in particular for the trafficking of membrane proteins within cells. The Endosomal Sorting Complex Required for Transport (ESCRT) directs membrane budding away from the cytosol. Unlike other vesicle formation pathways, the ESCRT-mediated budding occurs without a protein coat. Here, we propose a minimal model of ESCRT-induced vesicle budding. Our model is based on recent experimental observations from direct fluorescence microscopy imaging that show ESCRT proteins colocalized only in the neck region of membrane buds. The model, cast in the framework of membrane elasticity theory, reproduces the experimentally observed vesicle morphologies with physically meaningful parameters. In this parameter range, the minimum energy configurations of the membrane are coatless buds with ESCRTs localized in the bud neck, consistent with experiment. The minimum energy configurations agree with those seen in the fluorescence images, with respect to both bud shapes and ESCRT protein localization. On the basis of our model, we identify distinct mechanistic pathways for the ESCRT-mediated budding process. The bud size is determined by membrane material parameters, explaining the narrow yet different bud size distributions in vitro and in vivo. Our membrane elasticity model thus sheds light on the energetics and possible mechanisms of ESCRT-induced membrane budding.     

Frequency-dependent electrodeformation of giant phospholipid vesicles in AC electric field

A model of vesicle electrodeformation is described which obtains a parametrized vesicle shape by minimizing the sum of the membrane bending energy and the energy due to the electric field. Both the vesicle membrane and the aqueous media inside and outside the vesicle are treated as leaky dielectrics, and the vesicle itself is modeled as a nearly spherical shape enclosed within a thin membrane. It is demonstrated (a) that the model achieves a good quantitative agreement with the experimentally determined prolate-to-oblate transition frequencies in the kilohertz range and (b) that the model can explain a phase diagram of shapes of giant phospholipid vesicles with respect to two parameters: the frequency of the applied alternating current electric field and the ratio of the electrical conductivities of the aqueous media inside and outside the vesicle, explored in a recent paper (S. Aranda et al., Biophys J 95:L19–L21, 2008). A possible use of the frequency-dependent shape transitions of phospholipid vesicles in conductometry of microliter samples is discussed.     

The Vesicle Trafficking Protein Sar1 Lowers Lipid Membrane Rigidity

The sculpting of membranes into dynamic, curved shapes is central to intracellular cargo trafficking. Though the generation of membrane curvature during trafficking necessarily involves both lipids and membrane-associated proteins, current mechanistic views focus primarily on the formation of rigid cages and curved scaffolds by protein assemblies. Here we report on a different mechanism for the control of membrane deformation, unrelated to the imposition of predefined curvature, involving modulation of membrane material properties: Sar1, a GTPase that regulates vesicle trafficking from the endoplasmic reticulum, lowers the rigidity of the lipid bilayer membrane to which it binds. In vitro assays in which optically trapped microspheres create controlled membrane deformations revealed a monotonic decline in bending modulus as a function of Sar1 concentration, down to nearly zero rigidity, indicating a dramatic lowering of the energetic cost of curvature generation. This is the first demonstration that a vesicle trafficking protein lowers the rigidity of its target membrane, leading to a new conceptual framework for vesicle biogenesis.     

Modeling Membrane Deformations and Lipid Demixing upon Protein-Membrane Interaction: The BAR Dimer Adsorption

We use a self-consistent mean-field theory, designed to investigate membrane reshaping and lipid demixing upon interaction with proteins, to explore BAR domains interacting with large patches of lipid membranes of heterogeneous compositions. The computational model includes contributions to the system free energy from electrostatic interactions and elastic energies of the membrane, as well as salt and lipid mixing entropies. The results from our simulation of a single adsorbing Amphiphysin BAR dimer indicate that it is capable of stabilizing a significantly curved membrane. However, we predict that such deformations will occur only for membrane patches that have the inherent propensity for high curvature, reflected in the tendency to create local distortions that closely match the curvature of the BAR dimer itself. Such favorable preconditioning for BAR-membrane interaction may be the result of perturbations such as local lipid demixing induced by the interaction, or of a prior insertion of the BAR domain's amphiphatic N-helix. From our simulations it appears that local segregation of charged lipids under the influence of the BAR dimer cannot produce high enough asymmetry between bilayer leaflets to induce significant bending. In the absence of additional energy contributions that favor membrane asymmetry, the membrane will remain nearly flat upon single BAR dimer adsorption, relative to the undulation expected from thermal fluctuations. Thus, we conclude that the N-helix insertions have a critical mechanistic role in the local perturbation and curving of the membrane, which is then stabilized by the electrostatic interaction with the BAR dimer. We discuss how these results can be used to estimate the tendency of BARs to bend membranes in terms of a spatially nonisotropic spontaneous curvature.     

Membrane Curvature Induction and Tubulation Are Common Features of Synucleins and Apolipoproteins

Synucleins and apolipoproteins have been implicated in a number of membrane and lipid trafficking events. Lipid interaction for both types of proteins is mediated by 11 amino acid repeats that form amphipathic helices. This similarity suggests that synucleins and apolipoproteins might have comparable effects on lipid membranes, but this has not been shown directly. Here, we find that α-synuclein, β-synuclein, and apolipoprotein A-1 have the conserved functional ability to induce membrane curvature and to convert large vesicles into highly curved membrane tubules and vesicles. The resulting structures are morphologically similar to those generated by amphiphysin, a curvature-inducing protein involved in endocytosis. Unlike amphiphysin, however, synucleins and apolipoproteins do not require any scaffolding domains and curvature induction is mediated by the membrane insertion and wedging of amphipathic helices alone. Moreover, we frequently observed that α-synuclein caused membrane structures that had the appearance of nascent budding vesicles. The ability to function as a minimal machinery for vesicle budding agrees well with recent findings that α-synuclein plays a role in vesicle trafficking and enhances endocytosis. Induction of membrane curvature must be under strict regulation in vivo; however, as we find it can also cause disruption of membrane integrity. Because the degree of membrane curvature induction depends on the concerted action of multiple proteins, controlling the local protein density of tubulating proteins may be important. How cellular safeguarding mechanisms prevent such potentially toxic events and whether they go awry in disease remains to be determined.     

Theory of cation-phospholipid-induced shape changes in a lipid bilayer couple

A quantitative model of ion binding and molecular interactions in the lipid bilayer membrane is proposed and found to be useful in examining the factors underlying such membrane characteristics as shape, sidedness, stability and vesicle size at various cation concentrations. The lipid membrane behaves as a bilayer couple whose preferential radius of curvature depends on the expansion or contraction of one monolayer relative to the other. It is proposed that molecular packing may be altered by electrostatic repulsion of adjacent like-charged phospholipid headgroups, or by bringing two headgroups closer together by divalent cation crossbridging. The surface concentrations of each type of cation-phospholipid complex can be described by simple binding equilibria and the Gouy-Chapman-Stern formulation for the surface potential in a diffuse double layer. The asymmetric distribution of acidic phospholipids in most biological membranes can account for the differential effects of identical ionic environments on either side of the bilayer. The fraction of vesicle material which tends to have a right-side-out orientation may be approximated by a normal distribution about the mean curvature. The theory generates vesicle sidedness distributions that, when fitted to experimental results from human erythrocyte membranes, provide an alternative method of estimating intrinsic cationphospholipid dissociation constants and other molecular parameters of the bilayer. The results also corroborate earlier suggestions that the Gouy-Chapman theory tends to overestimate free counter-ion concentrations at the surface under large surface potentials.     

Membrane Curvature Sensing by Amphipathic Helices Is Modulated by the Surrounding Protein Backbone

Membrane curvature is involved in numerous biological pathways like vesicle trafficking, endocytosis or nuclear pore complex assembly. In addition to its topological role, membrane curvature is sensed by specific proteins, enabling the coordination of biological processes in space and time. Amongst membrane curvature sensors are the ALPS (Amphipathic Lipid Packing Sensors). ALPS motifs are short peptides with peculiar amphipathic properties. They are found in proteins targeted to distinct curved membranes, mostly in the early secretory pathway. For instance, the ALPS motif of the golgin GMAP210 binds trafficking vesicles, while the ALPS motif of Nup133 targets nuclear pores. It is not clear if, besides curvature sensitivity, ALPS motifs also provide target specificity, or if other domains in the surrounding protein backbone are involved. To elucidate this aspect, we studied the subcellular localization of ALPS motifs outside their natural protein context. The ALPS motifs of GMAP210 or Nup133 were grafted on artificial fluorescent probes. Importantly, ALPS motifs are held in different positions and these contrasting architectures were mimicked by the fluorescent probes. The resulting chimeras recapitulated the original proteins localization, indicating that ALPS motifs are sufficient to specifically localize proteins. Modulating the electrostatic or hydrophobic content of Nup133 ALPS motif modified its avidity for cellular membranes but did not change its organelle targeting properties. In contrast, the structure of the backbone surrounding the helix strongly influenced targeting. In particular, introducing an artificial coiled-coil between ALPS and the fluorescent protein increased membrane curvature sensitivity. This coiled-coil domain also provided membrane curvature sensitivity to the amphipathic helix of Sar1. The degree of curvature sensitivity within the coiled-coil context remains correlated to the natural curvature sensitivity of the helices. This suggests that the chemistry of ALPS motifs is a key parameter for membrane curvature sensitivity, which can be further modulated by the surrounding protein backbone.     

An Electrostatically Preferred Lateral Orientation of SNARE Complex Suggests Novel Mechanisms for Driving Membrane Fusion

Biological membrane fusion is a basic cellular process catalyzed by SNARE proteins and additional auxiliary factors. Yet, the critical mechanistic details of SNARE-catalyzed membrane fusion are poorly understood, especially during rapid synaptic transmission. Here, we systematically assessed the electrostatic forces between SNARE complex, auxiliary proteins and fusing membranes by the nonlinear Poisson-Boltzmann equation using explicit models of membranes and proteins. We found that a previously unrecognized, structurally preferred and energetically highly favorable lateral orientation exists for the SNARE complex between fusing membranes. This preferred orientation immediately suggests a novel and simple synaptotagmin-dependent mechanistic trigger of membrane fusion. Moreover, electrostatic interactions between membranes, SNARE complex, and auxiliary proteins appear to orchestrate a series of membrane curvature events that set the stage for rapid synaptic vesicle fusion. Together, our electrostatic analyses of SNAREs and their regulatory factors suggest unexpected and potentially novel mechanisms for eukaryotic membrane fusion proteins.     

Water under the BAR

Many cellular processes require the generation of highly curved regions of cell membranes by interfacial membrane proteins. A number of such proteins are now known, and several mechanisms of curvature generation have been suggested, but so far a quantitative understanding of the importance of the various potential mechanisms remains elusive. Following previous theoretical work, we consider the electrostatic attraction that underlies the scaffold mechanism of membrane bending in the context of the N-BAR domain of amphiphysin. Analysis of atomistic molecular dynamics simulations reveals considerable water between the membrane and the positively charged concave face of the BAR, even when it is tightly bound to highly curved membranes. This results in significant screening of electrostatic interactions, suggesting that electrostatic attraction is not the main driving force behind curvature sensing, supporting recent experimental work. These results also emphasize the need for care when building coarse-grained models of protein-membrane interactions. These results are emphasized by simulations of oligomerized amphiphysin N-BARs at the atomistic and coarse-grained level. In the coarse-grained simulations, we find a strong dependence of the induced curvature on the dielectric screening.     

Regulation of Membrane-Shape Transitions Induced by I-BAR Domains

I-BAR proteins are well-known actin-cytoskeleton adaptors and have been observed to be involved in the formation of plasma membrane protrusions (filopodia). I-BAR proteins contain an all-helical, crescent-shaped IRSp53-MIM domain (IMD) dimer that is believed to be able to couple with a membrane shape. This coupling could involve the sensing and even the generation of negative plasma membrane curvature. Indeed, the in vitro studies have shown that IMDs can induce inward tubulation of liposomes. While N-BAR domains, which generate positive membrane curvature, have received a considerable amount of attention from both theory and experiments, the mechanisms of curvature coupling through IMDs are comparatively less studied and understood. Here we used a membrane-shape stability assay developed recently in our lab to quantitatively characterize IMD-induced membrane-shape transitions. We determined a membrane-shape stability diagram for IMDs that reveals how membrane tension and protein density can comodulate the generation of IMD-induced membrane protrusions. From comparison to analytical theory, we determine three key parameters that characterize the curvature coupling of IMD. We find that the curvature generation capacity of IMDs is significantly stronger compared to that of endophilin, an N-BAR protein known to be involved in plasma membrane shape transitions. Contrary to N-BAR domains, where amphipathic helix insertion is known to promote its membrane curvature generation, for IMDs we find that amphipathic helices inhibit membrane shape transitions, consistent with the inverse curvature that IMDs generate. Importantly, in both of these types of BAR domains, electrostatic interactions affect membrane-binding capacity, but do not appear to affect the curvature generation capacity of the protein. These two types of BAR domain proteins show qualitatively similar membrane shape stability diagrams, suggesting an underlying ubiquitous mechanism by which peripheral proteins regulate membrane curvature.     

Determination of bilayer membrane bending stiffness by tether formation from giant, thin-walled vesicles.

The curvature elastic modulus (bending stiffness) of stearoyloleoyl phosphatidylcholine (SOPC) bilayer membrane is determined from membrane tether formation experiments. R. E. Waugh and R. M. Hochmuth 1987. Biophys. J. 52:391-400) have shown that the radius of a bilayer cylinder (tether) is inversely related to the force supported along its axis. The coefficient that relates the axial force on the tether to the tether radius is the membrane bending stiffness. Thus, the bending stiffness can be calculated directly from measurements of the tether radius as a function of force. Giant (10-50-microns diam) thin-walled vesicles were aspirated into a micropipette and a tether was pulled out of the surface by gravitational forces on small glass beads that had adhered to the vesicle surface. Because the vesicle keeps constant surface area and volume, formation of the tether requires displacement of material from the projection of the vesicle in the pipette. Tethers can be made to grow longer or shorter or to maintain equilibrium by adjusting the aspiration pressure in the micropipette at constant tether force. The ratio of the change in the length of the tether to the change in the projection length is proportional to the ratio of the pipette radius to the tether radius. Thus, knowing the density and diameter of the glass beads and measuring the displacement of the projection as a function of tether length, independent determinations of the force on the tether and the tether radius were obtained. The bending stiffness for an SOPC bilayer obtained from these data is approximately 2.0 x 10(-12) dyn cm, for tether radii in the range of 20-100 nm.(ABSTRACT TRUNCATED AT 250 WORDS)     

Membrane bending by protein-protein crowding

Curved membranes are an essential feature of dynamic cellular structures, including endocytic pits, filopodia protrusions and most organelles. It has been proposed that specialized proteins induce curvature by binding to membranes through two primary mechanisms: membrane scaffolding by curved proteins or complexes; and insertion of wedge-like amphipathic helices into the membrane. Recent computational studies have raised questions about the efficiency of the helix-insertion mechanism, predicting that proteins must cover nearly 100% of the membrane surface to generate high curvature, an improbable physiological situation. Thus, at present, we lack a sufficient physical explanation of how protein attachment bends membranes efficiently. On the basis of studies of epsin1 and AP180, proteins involved in clathrin-mediated endocytosis, we propose a third general mechanism for bending fluid cellular membranes: protein-protein crowding. By correlating membrane tubulation with measurements of protein densities on membrane surfaces, we demonstrate that lateral pressure generated by collisions between bound proteins drives bending. Whether proteins attach by inserting a helix or by binding lipid heads with an engineered tag, protein coverage above ~20% is sufficient to bend membranes. Consistent with this crowding mechanism, we find that even proteins unrelated to membrane curvature, such as green fluorescent protein (GFP), can bend membranes when sufficiently concentrated. These findings demonstrate a highly efficient mechanism by which the crowded protein environment on the surface of cellular membranes can contribute to membrane shape change.     

Spontaneous curvature of fluid vesicles induced by trans-bilayer sugar asymmetry

We present measurements of the effective spontaneous curvature of fluid lipid bilayers as a function of trans-bilayer asymmetry. Experiments are performed on micrometer-scale vesicles in sugar solutions with varying species across the membrane. There are two effects leading to a preferred curvature of such a vesicle. The spontaneous curvatures of the two monolayers as well as their area difference combine into an effective spontaneous curvature of the membrane. Our technique for measuring this parameter allows us to use vesicle morphology as a probe for general membrane-solute interactions affecting elasticity. Received: 3 June 1998 / Revised version: 18 August 1998 / Accepted: 21 August 1998     

The influence of curvature on membrane domains

An interdependence between local curvature and domain formation has been observed in both cell and model membranes. An implication of this observation is that domain formation in model membranes may be modulated by membrane curvature. In this paper, small-angle neutron scattering (SANS) is used to examine the influence of membrane curvature (i.e., vesicle size) on the formation of membrane domains. It is found that, although vesicle size and polydispersity are not significantly altered by the formation of membrane domains, the area fraction occupied by domains depends on the overall vesicle size. In particular, increasing membrane curvature (i.e., decreasing vesicle size) results in increased area fractions of membrane domains.     

Membrane curvature affects the stability and folding kinetics of bacteriorhodopsin

Biological membrane is crucial for the function, stability and folding of membrane proteins. By studying the stability and folding kinetics of bacteriorhodopsin (bR) in lipid vesicles with different sizes, here we report the influence of membrane curvature (vesicle size) on the stability and folding kinetics of bR. The results show that both the stability and folding kinetics of bR can be significantly changed when reconstituted into mimic membranes with different curvatures. The stability of bR decreases and unfolding rate of bR increases with the growth of vesicle size, i.e. decrease of membrane curvature. Our results suggest that it is possible to regulate the properties of membrane proteins by changing the curvature of membranes.     

Membrane Bending Energy and Fusion Pore Kinetics in Ca-Triggered Exocytosis

A fusion pore composed of lipid is an obligatory kinetic intermediate of membrane fusion, and its formation requires energy to bend membranes into highly curved shapes. The energetics of such deformations in viral fusion is well established, but the role of membrane bending in Ca²⁺-triggered exocytosis remains largely untested. Amperometry recording showed that during exocytosis in chromaffin and PC12 cells, fusion pores formed by smaller vesicles dilated more rapidly than fusion pores formed by larger vesicles. The logarithm of 1/(fusion pore lifetime) varied linearly with vesicle curvature. The vesicle size dependence of fusion pore lifetime quantitatively accounted for the nonexponential fusion pore lifetime distribution. Experimentally manipulating vesicle size failed to alter the size dependence of fusion pore lifetime. Manipulations of membrane spontaneous curvature altered this dependence, and applying the curvature perturbants to the opposite side of the membrane reversed their effects. These effects of curvature perturbants were opposite to those seen in viral fusion. These results indicate that during Ca²⁺-triggered exocytosis membrane bending opposes fusion pore dilation rather than fusion pore formation. Ca²⁺-triggered exocytosis begins with a proteinaceous fusion pore with less stressed membrane, and becomes lipidic as it dilates, bending membrane into a highly curved shape.     

Bending membranes

It is widely assumed that peripheral membrane proteins induce intracellular membrane curvature by the asymmetric insertion of a protein segment into the lipid bilayer, or by imposing shape by adhesion of a curved protein domain to the membrane surface. Two papers now provide convincing evidence challenging these views. The first shows that specific assembly of a clathrin protein scaffold, coupled to the membrane, seems to be the most prevalent mechanism for bending a lipid bilayer in a cell. The second reports that membrane crowding, driven by protein-protein interactions, can also drive membrane bending, even in the absence of any protein insertion into the bilayer.