Factors influencing local membrane curvature induction by N-BAR domains as revealed by molecular dynamics simulations

N-BAR domains are protein modules that bind to and induce curvature in membranes via a charged concave surface and N-terminal amphipathic helices. Recently, molecular dynamics simulations have demonstrated that the N-BAR domain can induce a strong local curvature that matches the curvature of the BAR domain surface facing the bilayer. Here we present further molecular dynamics simulations that examine in greater detail the roles of the concave surface and amphipathic helices in driving local membrane curvature. We find that the strong curvature induction observed in our previous simulations requires the stable presentation of the charged concave surface to the membrane and is not driven by the membrane-embedded amphipathic helices. Nevertheless, without these amphipathic helices embedded in the membrane, the N-BAR domain does not maintain a close association with the bilayer, and fails to drive membrane curvature. Increasing the membrane negative charge through the addition of PIP₂ facilitates closer association with the membrane in the absence of embedded helices. At sufficiently high concentrations, amphipathic helices embedded in the membrane drive membrane curvature independently of the BAR domain.     

Interplay of septin amphipathic helices in sensing membrane-curvature and filament bundling

The curvature of the membrane defines cell shape. Septins are GTP-binding proteins that assemble into heteromeric complexes and polymerize into filaments at areas of micron-scale membrane curvature. An amphipathic helix (AH) domain within the septin complex is necessary and sufficient for septins to preferentially assemble onto micron-scale curvature. Here we report that the nonessential fungal septin, Shs1, also has an AH domain capable of recognizing membrane curvature. In a septin mutant strain lacking a fully functional Cdc12 AH domain (cdc12-6), the C-terminal extension of Shs1, containing an AH domain, becomes essential. Additionally, we find that the Cdc12 AH domain is important for regulating septin filament bundling, suggesting septin AH domains have multiple, distinct functions and that bundling and membrane binding may be coordinately controlled.     

Subcellular membrane curvature mediated by the BAR domain superfamily proteins

The Bin-Amphiphysin-Rvs167 (BAR) domain superfamily consists of proteins containing the BAR domain, the extended FCH (EFC)/FCH-BAR (F-BAR) domain, or the IRSp53-MIM homology domain (IMD)/inverse BAR (I-BAR) domain. These domains bind membranes through electrostatic interactions between the negative charges of the membranes and the positive charges on the structural surface of homo-dimeric BAR domain superfamily members. Some BAR superfamily members have membrane-penetrating insertion loops, which also contribute to the membrane binding by the proteins. The membrane-binding surface of each BAR domain superfamily member has its own unique curvature that governs or senses the curvature of the membrane for BAR-domain binding. The wide range of BAR-domain surface curvatures correlates with the various invaginations and protrusions of cells. Therefore, each BAR domain superfamily member may generate and recognize the curvature of the membrane of each subcellular structure, such as clathrin-coated pits or filopodia. The BAR domain superfamily proteins may regulate their own catalytic activity or that of their binding proteins, depending on the membrane curvature of their corresponding subcellular structures.     

Mechanism of endophilin N-BAR domain-mediated membrane curvature

Endophilin-A1 is a BAR domain-containing protein enriched at synapses and is implicated in synaptic vesicle endocytosis. It binds to dynamin and synaptojanin via a C-terminal SH3 domain. We examine the mechanism by which the BAR domain and an N-terminal amphipathic helix, which folds upon membrane binding, work as a functional unit (the N-BAR domain) to promote dimerisation and membrane curvature generation. By electron paramagnetic resonance spectroscopy, we show that this amphipathic helix is peripherally bound in the plane of the membrane, with the midpoint of insertion aligned with the phosphate level of headgroups. This places the helix in an optimal position to effect membrane curvature generation. We solved the crystal structure of rat endophilin-A1 BAR domain and examined a distinctive insert protruding from the membrane interaction face. This insert is predicted to form an additional amphipathic helix and is important for curvature generation. Its presence defines an endophilin/nadrin subclass of BAR domains. We propose that N-BAR domains function as low-affinity dimers regulating binding partner recruitment to areas of high membrane curvature.     

BAR domains and membrane curvature: bringing your curves to the BAR

BAR (bin, amphiphysin and Rvs161/167) domains are a unique class of dimerization domains, whose dimerization interface is edged by a membrane-binding surface. In its dimeric form, the membrane-binding interface is concave, and this gives the ability to bind better to curved membranes, i.e. to sense membrane curvature. When present at higher concentrations, the domain can stabilize membrane curvature, generating lipid tubules. This domain is found in many contexts in a wide variety of proteins, where the dimerization and membrane-binding function of this domain is likely to have a profound effect on protein activity. If these proteins function as predicted, then there will be membrane subdomains based on curvature, and thus there is an additional layer of compartmentalization on membranes. These and other possible functions of the BAR domain are discussed.     

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.     

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.     

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.     

Roles of Amphipathic Helices and the Bin/Amphiphysin/Rvs (BAR) Domain of Endophilin in Membrane Curvature Generation

Control of membrane curvature is required in many important cellular processes, including endocytosis and vesicular trafficking. Endophilin is a bin/amphiphysin/rvs (BAR) domain protein that induces vesicle formation by promotion of membrane curvature through membrane binding as a dimer. Using site-directed spin labeling and EPR spectroscopy, we show that the overall BAR domain structure of the rat endophilin A1 dimer determined crystallographically is maintained under predominantly vesiculating conditions. Spin-labeled side chains on the concave surface of the BAR domain do not penetrate into the acyl chain interior, indicating that the BAR domain interacts only peripherally with the surface of a curved bilayer. Using a combination of EPR data and computational refinement, we determined the structure of residues 63–86, a region that is disordered in the crystal structure of rat endophilin A1. Upon membrane binding, residues 63–75 in each subunit of the endophilin dimer form a slightly tilted, amphipathic α-helix that directly interacts with the membrane. In their predominant conformation, these helices are located orthogonal to the long axis of the BAR domain. In this conformation, the amphipathic helices are positioned to act as molecular wedges that induce membrane curvature along the concave surface of the BAR domain.     

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.     

Inner but Not Outer Membrane Leaflets Control the Transition from Glycosylphosphatidylinositol-anchored Influenza Hemagglutinin-induced Hemifusion to Full Fusion

Cells that express wild-type influenza hemagglutinin (HA) fully fuse to RBCs, while cells that express the HA-ectodomain anchored to membranes by glycosylphosphatidylinositol, rather than by a transmembrane domain, only hemifuse to RBCs. Amphipaths were inserted into inner and outer membrane leaflets to determine the contribution of each leaflet in the transition from hemifusion to fusion. When inserted into outer leaflets, amphipaths did not promote the transition, independent of whether the agent induces monolayers to bend outward (conferring positive spontaneous monolayer curvature) or inward (negative curvature). In contrast, when incorporated into inner leaflets, positive curvature agents led to full fusion. This suggests that fusion is completed when a lipidic fusion pore with net positive curvature is formed by the inner leaflets that compose a hemifusion diaphragm. Suboptimal fusion conditions were established for RBCs bound to cells expressing wild-type HA so that lipid but not aqueous dye spread was observed. While this is the same pattern of dye spread as in stable hemifusion, for this “stunted” fusion, lower concentrations of amphipaths in inner leaflets were required to promote transfer of aqueous dyes. Also, these amphipaths induced larger pores for stunted fusion than they generated within a stable hemifusion diaphragm. Therefore, spontaneous curvature of inner leaflets can affect formation and enlargement of fusion pores induced by HA. We propose that after the HA-ectodomain induces hemifusion, the transmembrane domain causes pore formation by conferring positive spontaneous curvature to leaflets of the hemifusion diaphragm.     

Endophilin BAR domain drives membrane curvature by two newly identified structure-based mechanisms

The crescent-shaped BAR (Bin/Amphiphysin/Rvs-homology) domain dimer is a versatile protein module that senses and generates positive membrane curvature. The BAR domain dimer of human endophilin-A1, solved at 3.1 A, has a unique structure consisting of a pair of helix-loop appendages sprouting out from the crescent. The appendage's short helices form a hydrophobic ridge, which runs across the concave surface at its center. Examining liposome binding and tubulation in vitro using purified BAR domain and its mutants indicated that the ridge penetrates into the membrane bilayer and enhances liposome tubulation. BAR domain-expressing cells exhibited marked plasma membrane tubulation in vivo. Furthermore, a swinging-arm mutant lost liposome tubulation activity yet retaining liposome binding. These data suggested that the rigid crescent dimer shape is crucial for the tubulation. We here propose that the BAR domain drives membrane curvature by coordinate action of the crescent's scaffold mechanism and the ridge's membrane insertion in addition to membrane binding via amino-terminal amphipathic helix.     

Membrane sorting via the extracellular matrix

We consider the coupling between a membrane and the extracellular matrix. Computer simulations demonstrate that the latter coupling is able to sort lipids. It is assumed that membranes are elastic manifolds, and that this manifold is disrupted by the extracellular matrix. For a solid-supported membrane with an actin network on top, regions of positive curvature are induced below the actin fibers. A similar mechanism is conceivable by assuming that the proteins which connect the cytoskeleton to the membrane induce local membrane curvature. The regions of non-zero curvature exist irrespective of any phase transition the lipids themselves may undergo. For lipids that prefer certain curvature, the extracellular matrix thus provides a spatial template for the resulting lateral domain structure of the membrane.     

Nonlinear sorting, curvature generation, and crowding of endophilin N-BAR on tubular membranes

The curvature of biological membranes is controlled by membrane-bound proteins. For example, during endocytosis, the sorting of membrane components, vesicle budding, and fission from the plasma membrane are mediated by adaptor and accessory proteins. Endophilin is a peripherally binding membrane protein that functions as an endocytic accessory protein. Endophilin's membrane tubulation capacity is well known. However, to understand the thermodynamic and mechanical aspects of endophilin function, experimental measurements need to be compared to quantitative theoretical models. We present measurements of curvature sorting and curvature generation of the endophilin A1 N-BAR domain on tubular membranes pulled from giant unilamellar vesicles. At low concentration, endophilin functions primarily as a membrane curvature sensor; at high concentrations, it also generates curvature. We determine the spontaneous curvature induced by endophilin and observe sigmoidal curvature/composition coupling isotherms that saturate at high membrane tensions and protein solution concentrations. The observation of saturation is supported by a strong dependence of lateral diffusion coefficients on protein density on the tether membrane. We develop a nonlinear curvature/composition coupling model that captures our experimental observations. Our model predicts a curvature-induced phase transition among two states with varying protein density and membrane curvature. This transition could act as a switch during endocytosis.     

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.     

Molecular Basis of the Potent Membrane-remodeling Activity of the Epsin 1 N-terminal Homology Domain

The mechanisms by which cytosolic proteins reversibly bind the membrane and induce the curvature for membrane trafficking and remodeling remain elusive. The epsin N-terminal homology (ENTH) domain has potent vesicle tubulation activity despite a lack of intrinsic molecular curvature. EPR revealed that the N-terminal α-helix penetrates the phosphatidylinositol 4,5-bisphosphate-containing membrane at a unique oblique angle and concomitantly interacts closely with helices from neighboring molecules in an antiparallel orientation. The quantitative fluorescence microscopy showed that the formation of highly ordered ENTH domain complexes beyond a critical size is essential for its vesicle tubulation activity. The mutations that interfere with the formation of large ENTH domain complexes abrogated the vesicle tubulation activity. Furthermore, the same mutations in the intact epsin 1 abolished its endocytic activity in mammalian cells. Collectively, these results show that the ENTH domain facilitates the cellular membrane budding and fission by a novel mechanism that is distinct from that proposed for BAR domains.     

Membrane remodeling from N-BAR domain interactions: insights from multi-scale simulation

Liposome remodeling processes (e.g., vesiculation and tubulation) due to N-BAR domain interactions with the lipid bilayer are explored with a multi-scale simulation approach. Results from atomistic-level molecular dynamics simulations of membrane binding to the concave face of N-BAR domains are used along with discretized mesoscopic field-theoretic simulations to examine how the spontaneous curvature fields generated by N-BAR domains result in membrane remodeling. It is found that tubulation can be generated by anisotropic N-BAR spontaneous curvature fields, whereas vesiculation is only observed with isotropic N-BAR spontaneous curvature fields at high density. The results of the multi-scale simulations provide insight into recent experimental observations.     

Autoinhibition of Endophilin in Solution via Interdomain Interactions

Endophilin is a key protein involved in clathrin-mediated endocytosis. Previous computational and experimental work suggested that the N-terminal helix is embedded into the membrane to induce curvature; however, the role of the SH3 domain remains controversial. To address this issue, we performed computer simulations of the endophilin dimer in solution to understand the interaction between the N-BAR and SH3 domains and its effect on biological function. We predict that the helix binds to the SH3 domain through hydrophobic and salt-bridge interactions. This protects the hydrophobic residues on both domains and keeps the SH3 domain near the end of the N-BAR domain, in agreement with previous experimental results. The complex has a binding strength similar to a few hydrogen bonds (13.0 ± 0.6 kcal/mol), and the SH3 domain stabilizes the structure of the N-terminal helix in solution. Electrostatic calculations show a large region of strongly positive electrostatic potential near the N-terminal that can orient the helix toward the membrane and likely embed the helix into the membrane surface. This predicted mechanism suggests that endophilin can select for both curvature and electrostatic potential when interacting with membranes, highlighting the importance of the SH3 domain in regulating the function of endophilin.     

Protein shape and crowding drive domain formation and curvature in biological membranes

Folding, curvature, and domain formation are characteristics of many biological membranes. Yet the mechanisms that drive both curvature and the formation of specialized domains enriched in particular protein complexes are unknown. For this reason, studies in membranes whose shape and organization are known under physiological conditions are of great value. We therefore conducted atomic force microscopy and polarized spectroscopy experiments on membranes of the photosynthetic bacterium Rhodobacter sphaeroides. These membranes are densely populated with peripheral light harvesting (LH2) complexes, physically and functionally connected to dimeric reaction center-light harvesting (RC-LH1-PufX) complexes. Here, we show that even when converting the dimeric RC-LH1-PufX complex into RC-LH1 monomers by deleting the gene encoding PufX, both the appearance of protein domains and the associated membrane curvature are retained. This suggests that a general mechanism may govern membrane organization and shape. Monte Carlo simulations of a membrane model accounting for crowding and protein geometry alone confirm that these features are sufficient to induce domain formation and membrane curvature. Our results suggest that coexisting ordered and fluid domains of like proteins can arise solely from asymmetries in protein size and shape, without the need to invoke specific interactions. Functionally, coexisting domains of different fluidity are of enormous importance to allow for diffusive processes to occur in crowded conditions.     

Amphipathic motifs in BAR domains are essential for membrane curvature sensing

BAR (Bin/Amphiphysin/Rvs) domains and amphipathic α‐helices (AHs) are believed to be sensors of membrane curvature thus facilitating the assembly of protein complexes on curved membranes. Here, we used quantitative fluorescence microscopy to compare the binding of both motifs on single nanosized liposomes of different diameters and therefore membrane curvature. Characterization of members of the three BAR domain families showed surprisingly that the crescent‐shaped BAR dimer with its positively charged concave face is not able to sense membrane curvature. Mutagenesis on BAR domains showed that membrane curvature sensing critically depends on the N‐terminal AH and furthermore that BAR domains sense membrane curvature through hydrophobic insertion in lipid packing defects and not through electrostatics. Consequently, amphipathic motifs, such as AHs, that are often associated with BAR domains emerge as an important means for a protein to sense membrane curvature. Measurements on single liposomes allowed us to document heterogeneous binding behaviour within the ensemble and quantify the influence of liposome polydispersity on bulk membrane curvature sensing experiments. The latter results suggest that bulk liposome‐binding experiments should be interpreted with great caution.