N-BAR Domain Research Articles

Membrane Morphologies Induced by Arc-Shaped Scaffolds Are Determined by Arc Angle and Coverage

The intricate shapes of biological membranes such as tubules and membrane stacks are induced by proteins. In this article, we systematically investigate the membrane shapes induced by arc-shaped scaffolds such as proteins and protein complexes with coarse-grained modeling and simulations. We find that arc-shaped scaffolds induce membrane tubules at membrane coverages larger than a threshold of ∼40%, irrespective of their arc angle. The membrane morphologies at intermediate coverages below this tubulation threshold, in contrast, strongly depend on the arc angle. Scaffolds with arc angles of about 60°, akin to N-BAR domains, do not change the membrane shape at coverages below the tubulation threshold, whereas scaffolds with arc angles larger than about 120° induce double-membrane stacks at intermediate coverages. The scaffolds stabilize the curved membrane edges that connect the membrane stacks, as suggested for complexes of reticulon proteins. Our results provide general insights on the determinants of membrane shaping by arc-shaped scaffolds.

Regulation of neuromuscular junction organization by Rab2 and its effector ICA69 in Drosophila

The mechanisms underlying synaptic differentiation, which involves neuronal membrane and cytoskeletal remodeling, are not completely understood. We performed a targeted RNAi-mediated screen of Drosophila BAR-domain proteins and identified islet cell autoantigen 69 kDa (ICA69) as one of the key regulators of morphological differentiation of the larval neuromuscular junction (NMJ). We show that Drosophila ICA69 colocalizes with α-Spectrin at the NMJ. The conserved N-BAR domain of ICA69 deforms liposomes in vitro Full-length ICA69 and the ICAC but not the N-BAR domain of ICA69 induce filopodia in cultured cells. Consistent with its cytoskeleton regulatory role, ICA69 mutants show reduced α-Spectrin immunoreactivity at the larval NMJ. Manipulating levels of ICA69 or its interactor PICK1 alters the synaptic level of ionotropic glutamate receptors (iGluRs). Moreover, reducing PICK1 or Rab2 levels phenocopies ICA69 mutation. Interestingly, Rab2 regulates not only synaptic iGluR but also ICA69 levels. Thus, our data suggest that: (1) ICA69 regulates NMJ organization through a pathway that involves PICK1 and Rab2, and (2) Rab2 functions genetically upstream of ICA69 and regulates NMJ organization and targeting/retention of iGluRs by regulating ICA69 levels.

The N-Terminal Amphipathic Helix of Endophilin Does Not Contribute to Its Molecular Curvature Generation Capacity.

N-BAR proteins such as endophilin are thought to bend lipid membranes via scaffolding (the molding of membranes through the crescent protein shape) and membrane insertion (also called wedging) of amphipathic helices. However, the contributions from these distinct mechanisms to membrane curvature generation and sensing have remained controversial. Here we quantitatively demonstrate that the amphipathic N-terminal H0 helix of endophilin is important for recruiting this protein to the membrane, but does not contribute significantly to its intrinsic membrane curvature generation capacity. These observations elevate the importance of the scaffolding mechanism, rather than H0 insertion, for the membrane curvature generation by N-BAR domains. Furthermore, consistent with the thermodynamically required coupling between curvature generation and sensing, we observed that the H0-truncated N-BAR domain is capable of sensing membrane curvature. Overall, our contribution clarifies an important mechanistic controversy in the function of N-BAR domain proteins.

Calmodulin Promotes N-BAR Domain-Mediated Membrane Constriction and Endocytosis

Membrane remodeling by BAR (Bin, Amphiphysin, RVS) domain-containing proteins, such as endophilins and amphiphysins, is integral to the process of endocytosis. However, little is known about the regulation of endocytic BAR domain activity. We have identified an interaction between the yeast Rvs167 N-BAR domain and calmodulin. Calmodulin-binding mutants of Rvs167 exhibited defects in endocytic vesicle release. In vitro, calmodulin enhanced membrane tubulation and constriction by wild-type Rvs167 but not calmodulin-binding-defective mutants. A subset of mammalian N-BAR domains bound calmodulin, and co-expression of calmodulin with endophilin A2 potentiated tubulation in vivo. These studies reveal a conserved role for calmodulin in regulating the intrinsic membrane-sculpting activity of endocytic N-BAR domains.

Supramolecular organization of the human N-BAR domain in shaping the sarcolemma membrane.

The 30kDa N-BAR domain of the human Bin1 protein is essential for the generation of skeletal muscle T-tubules. By electron cryo-microscopy and electron cryo-tomography with a direct electron detector, we found that Bin1-N-BAR domains assemble into scaffolds of low long-range order that form flexible membrane tubules. The diameter of the tubules closely matches the curved shape of the N-BAR domain, which depends on the composition of the target membrane. These insights are fundamental to our understanding of T-tubule formation and function in human skeletal muscle.

Bar Domain Proteins can Differ Substantially in their Capacity to Generate Membrane Curvature

BAR (Bin/amphiphysin/Rvs) domains bind to and reshape cell membranes. Numerous BAR proteins are known to participate in endocytosis pathways. Whether their function is redundant or specifically different has remained unclear. Here we first compared the membrane curvature generation abilities of three N-BAR domains involved in endocytosis: endophilin, amphiphysin and SNX9. For this purpose, we used a membrane shape stability assay recently developed in our lab that allows determination of shape stability diagrams obtained via pipette aspirated giant unilamellar vesicles (GUVs). Significantly different membrane stability phase diagrams were observed for these three proteins, suggesting their functional non-redundant. We further compared the membrane binding and curvature generating of endophilin with an F-BAR protein: FCHo2, which arrives at clathrin coated pits (CCPs) earlier than endophilin. We observed interesting competitive membrane binding between FCHo2 and endophilin, as well as stronger curvature induction abilities of FCHo2 when compared with endophilin. These observations suggest a membrane binding regulation mechanism by the BAR domains, and promising to explain why FCHo2 act as the nucleator of the CCPs.

Osmotically-induced tension and the binding of N-BAR protein to lipid vesicles.

The binding affinity of a curvature-sensing protein domain (N-BAR) is measured as a function of applied osmotic stress while the membrane curvature is nearly constant. Varying the osmotic stress allows us to control membrane tension, which provides a probe of the mechanism of binding. We study the N-BAR domain of the Drosophila amphiphysin and monitor its binding on 50 nm-radius vesicles composed of 90 mol% DOPC and 10 mol% PIP. We find that the bound fraction of N-BAR is enhanced by a factor of approximately 6.5 when the tension increases from zero to 2.6 mN m(-1). This tension-induced response can be explained by the hydrophobic insertion mechanism. From the data we extract a hydrophobic domain area that is consistent with known structure. These results indicate that membrane stress and strain could play a major role in the previously reported curvature-affinity of N-BAR.

Structural insights into the cooperative remodeling of membranes by amphiphysin/BIN1

Amphiphysin2/BIN1 is a crescent-shaped N-BAR protein playing a key role in forming deeply invaginated tubes in muscle T-tubules. Amphiphysin2/BIN1 structurally stabilizes tubular formations in contrast to other N-BAR proteins involved in dynamic membrane scission processes; however, the molecular mechanism of the stabilizing effect is poorly understood. Using cryo-EM, we investigated the assembly of the amphiphysin/BIN1 on a membrane tube. We found that the N-BAR domains self-assemble on the membrane surface in a highly cooperative manner. Our biochemical assays and 3D reconstructions indicate that the N-terminal amphipathic helix H0 plays an important role in the initiation of the tube assembly and further in organizing BAR-mediated polymerization by locking adjacent N-BAR domains. Mutants that lack H0 or the tip portion, which is also involved in interactions of the neighboring BAR unit, lead to a disruption of the polymer organization, even though tubulation can still be observed. The regulatory region of amphiphysin/BIN1 including an SH3 domain does not have any apparent involvement in the polymer lattice. Our study indicates that the H0 helix and the BAR tip are necessary for efficient and organized self-assembly of amphiphysin/N-BAR.

Protein Folding Mechanism of the Dimeric AmphiphysinII/Bin1 N-BAR Domain.

The human AmphyphisinII/Bin1 N-BAR domain belongs to the BAR domain superfamily, whose members sense and generate membrane curvatures. The N-BAR domain is a 57 kDa homodimeric protein comprising a six helix bundle. Here we report the protein folding mechanism of this protein as a representative of this protein superfamily. The concentration dependent thermodynamic stability was studied by urea equilibrium transition curves followed by fluorescence and far-UV CD spectroscopy. Kinetic unfolding and refolding experiments, including rapid double and triple mixing techniques, allowed to unravel the complex folding behavior of N-BAR. The equilibrium unfolding transition curve can be described by a two-state process, while the folding kinetics show four refolding phases, an additional burst reaction and two unfolding phases. All fast refolding phases show a rollover in the chevron plot but only one of these phases depends on the protein concentration reporting the dimerization step. Secondary structure formation occurs during the three fast refolding phases. The slowest phase can be assigned to a proline isomerization. All kinetic experiments were also followed by fluorescence anisotropy detection to verify the assignment of the dimerization step to the respective folding phase. Based on these experiments we propose for N-BAR two parallel folding pathways towards the homodimeric native state depending on the proline conformation in the unfolded state.

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.

How Synaptotagmin I, N-BAR and F-BAR Domains Generate Membrane Curvature

Protein-induced membrane curving governs many cellular processes, including cell division, growth and cell-cell communication. In order to unravel the mechanism for how membrane curvature occurs driven by the proteins, one needs to have detailed molecular pictures on the key membrane-protein interactions and protein-protein interactions. Computer simulation serves as a great tool in providing such molecular pictures and thus in rationalizing how proteins curve membrane. For example, it is known that neurotransmitter release involves a Ca++ ion-regulated fusion process between synaptic vesicles and the presynaptic membrane, but little is known for how such small ions function. With the help of computer simulation, we realized the Ca++-sensor protein, synaptotagmin I, undergoes a conformational transition after binding to both the Ca++ and membrane that differs from the membrane-free crystal/NMR structure. The new synaptotagmin conformation has its C-terminal helix to interact with the presynaptic membrane, and thereby, causes the presynaptic membrane to curve and facilitates membrane-vesicle fusion. Another example in showing the power of computer simulation lies in resolving how N-BAR and F-BAR domains sculpt a flat membrane. As these domains form a lattice in the sculpting process, we identified the optimal lattice type and key protein-protein interactions within the lattice. The whole process of membrane tubulation from a flat membrane was resolved from the simulations, and agreements were found between the lattice structures observed via cryo-electron microscopy and the simulations.

Finely-Fabricated Nanomaterials to Reveal Curvature Roles in Neuronal Network Formation

In the vertebrate central nervous system, exploratory filopodia transiently form on dendritic branches to sample the neuronal environment and initiate new trans-neuronal contacts. While much is known about the molecules that control filopodia extension and subsequent maturation into functional synapses, the mechanisms that regulate initiation of these dynamic, actin-rich structures have remained elusive.Using nanomaterials that deform the plasma membrane of neurons, correlative light-electron microscopy (CLEM) and quantitative live cell microscopy, we find that filopodia initiation is suppressed by recruitment of ArhGAP44 to actin-patches that seed filopodia in dendritic branches. Recruitment is mediated by binding of a membrane curvature-sensing ArhGAP44 N-BAR domain to convoluted plasma membrane sections that were deformed inward by acto-myosin mediated contractile forces. A GAP domain in the enriched ArhGAP44 triggers local Rac-GTP hydrolysis, thus reducing local actin polymerization required for filopodia formation. ArhGAP44 expression increases as the neuronal network is established and the frequency of exploratory filopodia formation is diminished, suggesting that ArhGAP44 may facilitate the transition of neurons from a dynamic exploratory mode to a mature more static state, a hallmark of nervous system development.Together, our data reveals a local and receptor-independent auto-regulatory mechanism that limits initiation of exploratory filopodia in neurons via protein recruitment to nanoscale membrane deformations.

Dynamic recruitment of the curvature-sensitive protein ArhGAP44 to nanoscale membrane deformations limits exploratory filopodia initiation in neurons.

In the vertebrate central nervous system, exploratory filopodia transiently form on dendritic branches to sample the neuronal environment and initiate new trans-neuronal contacts. While much is known about the molecules that control filopodia extension and subsequent maturation into functional synapses, the mechanisms that regulate initiation of these dynamic, actin-rich structures have remained elusive. Here, we find that filopodia initiation is suppressed by recruitment of ArhGAP44 to actin-patches that seed filopodia. Recruitment is mediated by binding of a membrane curvature-sensing ArhGAP44 N-BAR domain to plasma membrane sections that were deformed inward by acto-myosin mediated contractile forces. A GAP domain in ArhGAP44 triggers local Rac-GTP hydrolysis, thus reducing actin polymerization required for filopodia formation. Additionally, ArhGAP44 expression increases during neuronal development, concurrent with a decrease in the rate of filopodia formation. Together, our data reveals a local auto-regulatory mechanism that limits initiation of filopodia via protein recruitment to nanoscale membrane deformations.

BIN1 membrane curvature sensing and generation show autoinhibition regulated by downstream ligands and PI(4,5)P2.

In striated muscles, invaginations from the plasma membrane, termed transverse tubules (T-tubule), function in the excitation–contraction coupling machinery. BIN1 (isoform8) plays a critical role in the biogenesis of T-tubules. BIN1 contains an N-terminal BAR domain to sense and induce membrane curvature, an isoform8-specific polybasic motif (exon10) as the phosphoinositide binding module and a C-terminal Src homology 3 (SH3) domain for the recruitment of downstream proteins such as dynamin 2. Previous studies of N-BAR domains focused on elucidating mechanisms of membrane curvature sensing and generation (MC-S&G). Less is known about how MC-S&G is regulated. We found that the SH3 domain binds to the exon10 motif more strongly compared to the proline-rich domain (PRD) of dynamin 2. Furthermore, we found that the MC-S&G ability of full-length BIN1 is inhibited on membranes lacking PI(4,5)P2. Addition of PI(4,5)P2 in the membrane activates BIN1 to sense and induce membrane curvature. Co-presence of the SH3 domain and exon10 motif leads to the strongest phosphoinositide-mediated control of BIN1 function. Addition of SH3 domain ligand (such as PRD peptides), as well as addition of the water-soluble PI(4,5)P2 analogue, can both enhance the MC-S&G ability of BIN1 on membranes without PI(4,5)P2, indicating that the key to activate BIN1 is to disrupt the exon10–SH3 interaction. The nonsense mutation K436X, found in centronuclear myopathy (CNM) patients, abolishes SH3 domain binding with either exon10 or the PRD motif, resulting in increased membrane deformation capacity. Our results suggest an autoinhibition model for BIN1 that involves a synergistic regulation by membrane composition and protein–protein interactions.

Arf6 exchange factor EFA6 and endophilin directly interact at the plasma membrane to control clathrin-mediated endocytosis.

Members of the Arf family of small G proteins are involved in membrane traffic and organelle structure. They control the recruitment of coat proteins, and modulate the structure of actin filaments and the lipid composition of membranes. The ADP-ribosylation factor 6 (Arf6) isoform and the exchange factor for Arf6 (EFA6) are known to regulate the endocytic pathway of many different receptors. To determine the molecular mechanism of the EFA6/Arf6 function in vesicular transport, we searched for new EFA6 partners. In a two-hybrid screening using the catalytic Sec7 domain as a bait, we identified endophilin as a new partner of EFA6. Endophilin contains a Bin/Amphiphysin/Rvs (BAR) domain responsible for membrane bending, and an SH3 domain responsible for the recruitment of dynamin and synaptojanin, two proteins involved, respectively, in the fission and uncoating of clathrin-coated vesicles. By using purified proteins, we confirmed the direct interaction, and identified the N-BAR domain as the binding motif to EFA6A. We showed that endophilin stimulates the catalytic activity of EFA6A on Arf6. In addition, we observed that the Sec7 domain competes with flat but not with highly curved lipid membranes to bind the N-BAR. In cells, expression of EFA6A recruits endophilin to EFA6A-positive plasma membrane ruffles, whereas expression of endophilin rescues the EFA6A-mediated inhibition of transferrin internalization. Overall, our results support a model whereby EFA6 recruits endophilin on flat areas of the plasma membrane to control Arf6 activation and clathrin-mediated endocytosis.

Toward the assembly of a minimal divisome.

The construction of an irreducible minimal cell having all essential attributes of a living system is one of the biggest challenges facing synthetic biology. One ubiquitous task accomplished by any living systems is the division of the cell envelope. Hence, the assembly of an elementary, albeit sufficient, molecular machinery that supports compartment division, is a crucial step towards the realization of self-reproducing artificial cells. Looking backward to the molecular nature of possible ancestral, supposedly more rudimentary, cell division systems may help to identify a minimal divisome. In light of a possible evolutionary pathway of division mechanisms from simple lipid vesicles toward modern life, we define two approaches for recapitulating division in primitive cells: the membrane deforming protein route and the lipid biosynthesis route. Having identified possible proteins and working mechanisms participating in membrane shape alteration, we then discuss how they could be integrated into the construction framework of a programmable minimal cell relying on gene expression inside liposomes. The protein synthesis using recombinant elements (PURE)system, a reconstituted minimal gene expression system, is conceivably the most versatile synthesis platform. As a first step towards the de novo synthesis of a divisome, we showed that the N-BAR domain protein produced from its gene could assemble onto the outer surface of liposomes and sculpt the membrane into tubular structures. We finally discuss the remaining challenges for building up a self-reproducing minimal cell, in particular the coupling of the division machinery with volume expansion and genome replication.

Mutations in BIN1 associated with centronuclear myopathy disrupt membrane remodeling by affecting protein density and oligomerization.

The regulation of membrane shapes is central to many cellular phenomena. Bin/Amphiphysin/Rvs (BAR) domain-containing proteins are key players for membrane remodeling during endocytosis, cell migration, and endosomal sorting. BIN1, which contains an N-BAR domain, is assumed to be essential for biogenesis of plasma membrane invaginations (T-tubules) in muscle tissues. Three mutations, K35N, D151N and R154Q, have been discovered so far in the BAR domain of BIN1 in patients with centronuclear myopathy (CNM), where impaired organization of T-tubules has been reported. However, molecular mechanisms behind this malfunction have remained elusive. None of the BIN1 disease mutants displayed a significantly compromised curvature sensing ability. However, two mutants showed impaired membrane tubulation both in vivo and in vitro, and displayed characteristically different behaviors. R154Q generated smaller membrane curvature compared to WT N-BAR. Quantification of protein density on membranes revealed a lower membrane-bound density for R154Q compared to WT and the other mutants, which appeared to be the primary reason for the observation of impaired deformation capacity. The D151N mutant was unable to tubulate liposomes under certain experimental conditions. At medium protein concentrations we found ‘budding’ structures on liposomes that we hypothesized to be intermediates during the tubulation process except for the D151N mutant. Chemical crosslinking assays suggested that the D151N mutation impaired protein oligomerization upon membrane binding. Although we found an insignificant difference between WT and K35N N-BAR in in vitro assays, depolymerizing actin in live cells allowed tubulation of plasma membranes through the K35N mutant. Our results provide insights into the membrane-involved pathophysiological mechanisms leading to human disease.

Intradimer/Intermolecular Interactions Suggest Autoinhibition Mechanism in Endophilin A1

Endophilin A1 is a homodimeric membrane-binding endocytic accessory protein with a high dimerization affinity. Its function has been hypothesized to involve autoinhibition. However, the autoinhibition mechanism, as well as the physicochemical basis for the high dimerization affinity of endophilin in solution, have remained unclear. In this contribution, we use a Förster resonance energy transfer (FRET) method to investigate the homodimerization mechanism and intradimer molecular interactions in endophilin. For the endophilin N-BAR domain (which lacks the SH3 domain including a linker region of the full length protein), we observe a large temperature dependence of the dimerization affinity and dimer dissociation kinetics, implying large dimerization enthalpy and dissociation activation enthalpy, respectively. Our evaluation of the protein concentration dependence of dimer dissociation kinetics implies that endophilin reversibly forms monomers via a dissociation/reassociation mechanism. Furthermore, we use a kinetic method that allows us to compare the dissociation kinetics of full-length endophilin to that of truncated mutants. We find that mutants that lack either H0 helix or SH3 domain show significantly faster dissociation kinetics relative to full-length endophilin. This observation supports the presence of an intradimer, intermonomer cross-interaction between H0 helix and SH3 domain from different subunits within a homodimer. Because the H0 helix is known to play a significant role in endophilin’s membrane interactions, our measurements support a syngergistic model where these interactions are inhibited in the absence of SH3 domain binding ligands such as dynamin’s prolin rich domains, and where the binding of these ligands may be suppressed for non-membrane-bound endophilin.

Endophilin B is required for the Drosophila oocyte to endocytose yolk downstream of Oskar

The nutritional environment is crucial for Drosophila oogenesis in terms of controlling hormonal conditions that regulate yolk production and the progress of vitellogenesis. Here, we discovered that Drosophila Endophilin B (D-EndoB), a member of the endophilin family, is required for yolk endocytosis as it regulates membrane dynamics in developing egg chambers. Loss of D-EndoB leads to yolk content reduction, similar to that seen in yolkless mutants, and also causes poor fecundity. In addition, mutant egg chambers exhibit an arrest at the previtellogenic stage. D-EndoB displayed a crescent localization at the oocyte posterior pole in an Oskar-dependent manner; however, it did not contribute to pole plasm assembly. D-EndoB was found to partially colocalize with Long Oskar and Yolkless at the endocytic membranes in ultrastructure analysis. Using an FM4-64 dye incorporation assay, D-EndoB was also found to promote endocytosis in the oocyte. When expressing the full-length D-endoB(FL) or D-endoB(ΔSH3) mutant transgenes in oocytes, the blockage of vitellogenesis and the defect in fecundity in D-endoB mutants was restored. By contrast, a truncated N-BAR domain of the D-EndoB only partially rescued these defects. Taken together, these results allow us to conclude that D-EndoB contributes to the endocytic activity downstream of Oskar by facilitating membrane dynamics through its N-BAR domain in the yolk uptake process, thereby leading to normal progression of vitellogenesis.

Cellular Scale Biophysical Models of Membrane Sculpting by the Proteins During Endocytosis and Exocytosis

We employ coarse-grained molecular dynamics (CGMD) simulations and continuum mechanics methods to quantify the membrane-bending action of exo70. CGMD simulations show that exo70 produces the strongest curvature when bound to a key lipid PIP2, and when dimerized. To relate these simulations to experiments, we first quantify the deformation field generated by exo70 using CGMD simulation, and then investigate the resulting cellular morphologies using a continuum mechanics model for membrane bending. In the first step, we characterize the deformation field by directly analyzing curved membrane conformations generated by CGMD simulations of mutant and wild-type exo70 on mixed lipid bilayers. Monte Carlo simulations yield membrane morphologies which we relate to cryo-transmission electron microscopy images of exo70 remodeling a GUV. This model agrees with experiments which show that PIP2 binding and dimerization are necessary for generating membrane tubulues and protrusions. We show that this multi-scale model allows us to relate the structural as well as thermodynamic details of exo70 membrane-binding and self-association to the induced deformation field and resulting cellular morphologies. We also extend our studies to ENTH, BAR, and ESCRT family proteins. We show that thermal undulations in the membrane and cooperativity in the curvature fields, collectively drive the membrane into different morphological states (buds, tubules, etc.) that resemble those in cellular experiments in vivo and vesicle experiments in vitro. We determine the relative stability of the above mentioned shapes based on the free energy of these membrane configurations, determined using the methods of Thermodynamic Integration (TI) and Bennett Acceptance (BA), and Widom Insertion techniques. Results are shown for the case of Exo70 protein and ENTH, and N-BAR domains and also compared against measurements determined from experiments.