Curvature sorting of proteins on a cylindrical lipid membrane tether connected to a reservoir

Biophysical modeling of membrane curvature generation and curvature sensing by the glycocalyx

Generation of membrane curvature is fundamental to cellular function. Recent studies have established that the glycocalyx, a sugar-rich polymer layer at the cell surface, can generate membrane curvature. While there have been some theoretical efforts to understand the interplay between the glycocalyx and membrane bending, there remain open questions about how the properties of the glycocalyx affect membrane bending. For example, the relationship between membrane curvature and the density of glycosylated proteins on its surface remains unclear. In this work, we use polymer brush theory to develop a detailed biophysical model of the energetic interactions of the glycocalyx with the membrane. Using this model, we identify the conditions under which the glycocalyx can both generate and sense curvature. Our model predicts that the extent of membrane curvature generated depends on the grafting density of the glycocalyx and the backbone length of the polymers constituting the glycocalyx. Furthermore, when coupled with the intrinsic membrane properties such as spontaneous curvature and a line tension along the membrane, the curvature generation properties of the glycocalyx are enhanced. These predictions were tested experimentally by examining the propensity of glycosylated transmembrane proteins to drive the assembly of highly curved filopodial protrusions at the plasma membrane of adherent mammalian cells. Our model also predicts that the glycocalyx has curvature-sensing capabilities, in agreement with the results of our experiments. Thus, our study develops a quantitative framework for mapping the properties of the glycocalyx to the curvature generation capability of the membrane.

Molecular mechanisms of steric pressure generation and membrane remodeling by disordered proteins

Cellular membranes, which are densely crowded by proteins, take on an elaborate array of highly curved shapes. Steric pressure generated by protein crowding plays a significant role in shaping membrane surfaces. It is increasingly clear that many proteins involved in membrane remodeling contain substantial regions of intrinsic disorder. These domains have large hydrodynamic radii, suggesting that they may contribute significantly to steric congestion on membrane surfaces. However, it has been unclear to what extent they are capable of generating steric pressure, owing to their conformational flexibility. To address this gap, we use a recently developed sensor based on Förster resonance energy transfer to measure steric pressure generated at membrane surfaces by the intrinsically disordered domain of the endocytic protein, AP180. We find that disordered domains generate substantial steric pressure that arises from both entropic and electrostatic components. Interestingly, this steric pressure is largely invariant with the molecular weight of the disordered domain, provided that coverage of the membrane surface is held constant. Moreover, equivalent levels of steric pressure result in equivalent degrees of membrane remodeling, regardless of protein molecular weight. This result, which is consistent with classical polymer scaling relationships for semi-dilute solutions, helps to explain the molecular and physical origins of steric pressure generation by intrinsically disordered domains. From a physiological perspective, these findings suggest that a broad range of membrane-associated disordered domains are likely to play a significant and previously unknown role in controlling membrane shape.

Sorting of proteins with shape and curvature anisotropy on a lipid bilayer tube

Curvature induced sorting of lipid membrane bound proteins has been widely studied through experiments that induce curvature variation in a giant unilamellar lipid-bilayer vesicle with adsorbed proteins by pulling thin cylindrical tethers. In the theoretical space, this has been supplemented with models that capture curvature dependent interaction between membrane and idealized protein particles, through free energy contributions. Many membrane proteins such as the BAR domain proteins are known to have extremely anisotropic shapes and soft interacting potentials, whereas the idealizations of protein particles explored in models have only assumed them as hard disk-like particles with curvature anisotropy. Here, we present a model of sorting of the proteins while including the effects of softness in their interaction potentials, shape anisotropy in the protein structure, and curvature anisotropy in the interactions with the membrane. This is based on a clean separation of free energy contributions from non-ideal fluid behavior of soft anisotropic particles and curvature interactions between proteins and membranes. We probe the behavior of the sorting function under limiting conditions and show that it converges to the previously derived models. In addition to this, we present a comparison of the variation in sorting ratio due to the observed variation in the shape parameter values in known membrane proteins. Finally, using published experimental data for membrane proteins, we perform fitting and derive model parameters. We observe that shape anisotropy adversely affects the sorting of proteins to a high curvature region, whereas curvature anisotropy and softer interaction between proteins favor sorting.

Recent developments in membrane curvature sensing and induction by proteins

BACKGROUND

Membrane-bound intracellular organelles have characteristic shapes attributed to different local membrane curvatures, and these attributes are conserved across species. Over the past decade, it has been confirmed that specific proteins control the large curvatures of the membrane, whereas many others due to their specific structural features can sense the curvatures and bind to the specific geometrical cues. Elucidating the interplay between sensing and induction is indispensable to understand the mechanisms behind various biological processes such as vesicular trafficking and budding.

SCOPE OF REVIEW

We provide an overview of major classes of membrane proteins and the mechanisms of curvature sensing and induction. We then discuss the importance of membrane elastic characteristics to induce the membrane shapes similar to intracellular organelles. Finally, we survey recently available assays developed for studying the curvature sensing and induction by many proteins.

MAJOR CONCLUSIONS

Recent theoretical/computational modeling along with experimental studies have uncovered fascinating connections between lipid membrane and protein interactions. However, the phenomena of protein localization and synchronization to generate spatiotemporal dynamics in membrane morphology are yet to be fully understood.

GENERAL SIGNIFICANCE

The understanding of protein-membrane interactions is essential to shed light on various biological processes. This further enables the technological applications of many natural proteins/peptides in therapeutic treatments. The studies of membrane dynamic shapes help to understand the fundamental functions of membranes, while the medicinal roles of various macromolecules (such as proteins, peptides, etc.) are being increasingly investigated.

Binding of thermalized and active membrane curvature-inducing proteins

The phase behavior of a membrane induced by the binding of curvature-inducing proteins is studied by a combination of analytical and numerical approaches. In thermal equilibrium under the detailed balance between binding and unbinding, the membrane exhibits three phases: an unbound uniform flat phase (U), a bound uniform flat phase (B), and a separated/corrugated phase (SC). In the SC phase, the bound proteins form hexagonally-ordered bowl-shaped domains. The transitions between the U and SC phases and between the B and SC phases are second order and first order, respectively. At a small spontaneous curvature of the protein or high surface tension, the transition between B and SC phases becomes continuous. Moreover, a first-order transition between the U and B phases is found at zero spontaneous curvature driven by the Casimir-like interactions between rigid proteins. Furthermore, nonequilibrium dynamics is investigated by the addition of active binding and unbinding at a constant rate. The active binding and unbinding processes alter the stability of the SC phase.

Microfluidic chip with pillar arrays for controlled production and observation of lipid membrane nanotubes

Lipid membrane nanotubes (NTs) are a widespread template for in vitro studies of cellular processes happening at high membrane curvature. Traditionally NTs are manufactured one by one, using sophisticated membrane micromanipulations, while simplified methods for controlled batch production of NTs are in growing demand. Here we propose a lab-on-a-chip (LOC) approach to the simultaneous formation of multiple NTs with length and radius controlled by the chip design. The NTs form upon rolling silica microbeads covered by lipid lamellas over the pillars of a polymer micropillar array. The array's design and surface chemistry set the geometry of the resulting free-standing NTs. The integration of the array inside a microfluidic chamber further enables fast and turbulence-free addition of components, such as proteins, to multiple preformed NTs. This LOC approach to NT production is compatible with the use of high power objectives of a fluorescence microscope, making real-time quantification of the different modes of the protein activity in a single experiment possible.

Noncooperative thermodynamics and kinetic models of ligand binding to polymers: Connecting McGhee-von Hippel model with the Tonks gas model

Ligand binding to polymers modifies the physical and chemical properties of the polymers, leading to physical, chemical, and biological implications. McGhee and von Hippel obtained the equilibrium coverage as a function of the ligand affinity, through the computation of the possible binding sites for the ligand. Here, we complete this theory deriving the kinetic model for the ligand-binding dynamics and the associated equilibrium chemical potential, which turns out to be of the Tonks gas model type. At low coverage, the Tonks chemical potential becomes the Fermi chemical potential and even the ideal gas chemical potential. We also discuss kinetic models associated with these chemical potentials. These results clarify the kinetic models of ligand binding, their relations with the chemical potentials, and their range of validity. Our results highlight the inaccuracy of ideal and simplified kinetic approaches for medium and high coverages.

Membrane Thinning Induces Sorting of Lipids and the Amphipathic Lipid Packing Sensor (ALPS) Protein Motif

Heterogeneities (e.g., membrane proteins and lipid domains) and deformations (e.g., highly curved membrane regions) in biological lipid membranes cause lipid packing defects that may trigger functional sorting of lipids and membrane-associated proteins. To study these phenomena in a controlled and efficient way within molecular simulations, we developed an external field protocol that artificially enhances packing defects in lipid membranes by enforcing local thinning of a flat membrane region. For varying lipid compositions, we observed strong thinning-induced depletion or enrichment, depending on the lipid's intrinsic shape and its effect on a membrane's elastic modulus. In particular, polyunsaturated and lysolipids are strongly attracted to regions high in packing defects, whereas phosphatidylethanolamine (PE) lipids and cholesterol are strongly repelled from it. Our results indicate that externally imposed changes in membrane thickness, area, and curvature are underpinned by shared membrane elastic principles. The observed sorting toward the thinner membrane region is in line with the sorting expected for a positively curved membrane region. Furthermore, we have demonstrated that the amphipathic lipid packing sensor (ALPS) protein motif, a known curvature and packing defect sensor, is effectively attracted to thinner membrane regions. By extracting the force that drives amphipathic molecules toward the thinner region, our thinning protocol can directly quantify and score the lipid packing sensing of different amphipathic molecules. In this way, our protocol paves the way toward high-throughput exploration of potential defect- and curvature-sensing motifs, making it a valuable addition to the molecular simulation toolbox.

Thermodynamic equilibrium of binary mixtures on curved surfaces

We study the global influence of curvature on the free energy landscape of two-dimensional binary mixtures confined on closed surfaces. Starting from a generic effective free energy, constructed on the basis of symmetry considerations and conservation laws, we identify several model-independent phenomena, such as a curvature-dependent line tension and local shifts in the binodal concentrations. To shed light on the origin of the phenomenological parameters appearing in the effective free energy, we further construct a lattice-gas model of binary mixtures on nontrivial substrates, based on the curved-space generalization of the two-dimensional Ising model. This allows us to decompose the interaction between the local concentration of the mixture and the substrate curvature into four distinct contributions, as a result of which the phase diagram splits into critical subdiagrams. The resulting free energy landscape can admit, as stable equilibria, strongly inhomogeneous mixed phases, which we refer to as "antimixed" states below the critical temperature. We corroborate our semianalytical findings with phase-field numerical simulations on realistic curved lattices. Despite this work being primarily motivated by recent experimental observations of multicomponent lipid vesicles supported by colloidal scaffolds, our results are applicable to any binary mixture confined on closed surfaces of arbitrary geometry.

Transition from curvature sensing to generation in a vesicle driven by protein binding strength and membrane tension

The ability of proteins to sense and/or generate membrane curvature is crucial for many biological processes inside the cell. We introduce a model for the binding and unbinding of curvature inducing proteins on vesicles using Dynamic Triangulation Monte Carlo (DTMC) simulations. In our study, the interaction between membrane curvature and protein binding is characterised by the binding affinity parameter μ, which indicates the interaction strength. We demonstrate that both sensing and generation of curvature can be observed in the same system as a function of the protein binding affinity on the membrane. Our results show that at low μ values, proteins only sense membrane curvature, whereas at high μ values, they induce curvature. The transition between sensing and generation regimes is marked by a sharp change in the μ-dependence of the protein bound fraction. We present ways to quantitatively characterise these two regimes. We also observe that imposing tension on the membrane (through internal excess pressure for liposomes) extends the region of curvature sensing in the parameter space.

Scaling relationships for translational diffusion constants applied to membrane domain dissolution and growth

We compare the way that relationships for diffusion constants scale with the size of diffusing membrane domains and the geometry of their environments. Then, we review our experimental work on the dynamics of dissolution/growth of membrane domains in crowding induced mixing, phase separation, and Ostwald ripening in a highly confined environment. Overall, the scaling relationships applied to diffusion constants obtained by fits to our dynamic data indicate that dissolution and growth is influenced by the diffusion of clusters or small domains of lipids, in addition to kinetic processes and geometrical constraints.

Hydrodynamics govern the pre-fusion docking time of synaptic vesicles

Synaptic vesicle fusion is a crucial step in the neurotransmission process. Neurotransmitter-filled vesicles are pre-docked at the synapse by the mediation of ribbon structures and SNARE proteins at the ribbon synapses. An electrical impulse triggers the fusion process of pre-docked vesicles, leading to the formation of a fusion pore and subsequently resulting in the release of neurotransmitter into the synaptic cleft. In this study, a continuum model of lipid membrane along with lubrication theory is used to determine the traverse time of the synaptic vesicle under the influence of hydrodynamic forces. We find that the traverse time is strongly dependent on how fast the driving force decays or grows with closure of the gap between the vesicle and the plasma membrane. If the correct behaviour is chosen, the traverse time obtained is of the order of a few hundred milliseconds and lies within the experimentally obtained value of approximately 250 ms (Zenisek D, Steyer JA, Almers W. 2000 Nature406, 849-854 (doi:10.1038/35022500)). We hypothesize that there are two different force behaviours, which complies with the experimental findings of pre-fusion docking of synaptic vesicles at the ribbon synapses. The common theme in the proposed force models is that the driving force has to very rapidly increase or decrease with the amount of clamping.

Excess area dependent scaling behavior of nano-sized membrane tethers

Thermal fluctuations in cell membranes manifest as an excess area ([Formula: see text]) which governs a multitude of physical process at the sub-micron scale. We present a theoretical framework, based on an in silico tether pulling method, which may be used to reliably estimate [Formula: see text] in live cells. We perform our simulations in two different thermodynamic ensembles: (i) the constant projected area and (ii) the constant frame tension ensembles and show the equivalence of our results in the two. The tether forces estimated from our simulations compare well with our experimental measurements for tethers extracted from ruptured GUVs and HeLa cells. We demonstrate the significance and validity of our method by showing that all our calculations performed in the initial tether formation regime (i.e. when the length of the tether is comparable to its radius) along with experiments of tether extraction in 15 different cell types collapse onto two unified scaling relationships mapping tether force, tether radius, bending stiffness κ, and membrane tension σ. We show that [Formula: see text] is an important determinant of the radius of the extracted tether, which is equal to the characteristic length [Formula: see text] for [Formula: see text], and is equal to [Formula: see text] for [Formula: see text]. We also find that the estimated excess area follows a linear scaling behavior that only depends on the true value of [Formula: see text] for the membrane, based on which we propose a self-consistent technique to estimate the range of excess membrane areas in a cell.

Analysis of diffusion in curved surfaces and its application to tubular membranes

Diffusion of particles in curved surfaces is inherently complex compared with diffusion in a flat membrane, owing to the nonplanarity of the surface. The consequence of such nonplanar geometry on diffusion is poorly understood but is highly relevant in the case of cell membranes, which often adopt complex geometries. To address this question, we developed a new finite element approach to model diffusion on curved membrane surfaces based on solutions to Fick's law of diffusion and used this to study the effects of geometry on the entry of surface-bound particles into tubules by diffusion. We show that variations in tubule radius and length can distinctly alter diffusion gradients in tubules over biologically relevant timescales. In addition, we show that tubular structures tend to retain concentration gradients for a longer time compared with a comparable flat surface. These findings indicate that sorting of particles along the surfaces of tubules can arise simply as a geometric consequence of the curvature without any specific contribution from the membrane environment. Our studies provide a framework for modeling diffusion in curved surfaces and suggest that biological regulation can emerge purely from membrane geometry.

Physical basis of some membrane shaping mechanisms

In vesicular transport pathways, membrane proteins and lipids are internalized, externalized or transported within cells, not by bulk diffusion of single molecules, but embedded in the membrane of small vesicles or thin tubules. The formation of these 'transport carriers' follows sequential events: membrane bending, fission from the donor compartment, transport and eventually fusion with the acceptor membrane. A similar sequence is involved during the internalization of drug or gene carriers inside cells. These membrane-shaping events are generally mediated by proteins binding to membranes. The mechanisms behind these biological processes are actively studied both in the context of cell biology and biophysics. Bin/amphiphysin/Rvs (BAR) domain proteins are ideally suited for illustrating how simple soft matter principles can account for membrane deformation by proteins. We review here some experimental methods and corresponding theoretical models to measure how these proteins affect the mechanics and the shape of membranes. In more detail, we show how an experimental method employing optical tweezers to pull a tube from a giant vesicle may give important quantitative insights into the mechanism by which proteins sense and generate membrane curvature and the mechanism of membrane scission.This article is part of the themed issue 'Soft interfacial materials: from fundamentals to formulation'.

When Physics Takes Over: BAR Proteins and Membrane Curvature

Cell membranes become highly curved during membrane trafficking, cytokinesis, infection, immune response, or cell motion. Bin/amphiphysin/Rvs (BAR) domain proteins with their intrinsically curved and anisotropic shape are involved in many of these processes, but with a large spectrum of modes of action. In vitro experiments and multiscale computer simulations have contributed in identifying a minimal set of physical parameters, namely protein density on the membrane, membrane tension, and membrane shape, that control how bound BAR domain proteins behave on the membrane. In this review, we summarize the multifaceted coupling of BAR proteins to membrane mechanics and propose a simple phase diagram that recapitulates the effects of these parameters.

IRSp53 senses negative membrane curvature and phase separates along membrane tubules

BAR domain proteins contribute to membrane deformation in diverse cellular processes. The inverted-BAR (I-BAR) protein IRSp53, for instance, is found on the inner leaflet of the tubular membrane of filopodia; however its role in the formation of these structures is incompletely understood. Here we develop an original assay in which proteins are encapsulated in giant unilamellar vesicles connected to membrane nanotubes. Our results demonstrate that I-BAR dimers sense negative membrane curvature. Experiment and theory reveal that the I-BAR displays a non-monotonic sorting with curvature, and expands the tube at high imposed tension while constricting it at low tension. Strikingly, at low protein density and tension, protein-rich domains appear along the tube. This peculiar behaviour is due to the shallow intrinsic curvature of I-BAR dimers. It allows constriction of weakly curved membranes coupled to local protein enrichment at biologically relevant conditions. This might explain how IRSp53 contributes in vivo to the initiation of filopodia.

The inverted-BAR domain protein IRSp53 associates with the inner leaflet of tubular membranes such as filopodia. Here, Prévost et al. demonstrate that the I-BAR domain of IRSp53 senses negative membrane curvature, and undergoes phase separation which may aid its clustering upon filopodia generation.

Application of a free-energy-landscape approach to study tension-dependent bilayer tubulation mediated by curvature-inducing proteins

We investigate the phenomenon of protein-induced tubulation of lipid bilayer membranes within a continuum framework using Monte Carlo simulations coupled with the Widom insertion technique to compute excess chemical potentials. Tubular morphologies are spontaneously formed when the density and the curvature-field strength of the membrane-bound proteins exceed their respective thresholds and this transition is marked by a sharp drop in the excess chemical potential. We find that the planar to tubular transition can be described by a micellar model and that the corresponding free-energy barrier increases with an increase in the curvature-field strength (i.e., of protein-membrane interactions) and also with an increase in membrane tension.

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.

Sorting of integral membrane proteins mediated by curvature-dependent protein-lipid bilayer interaction

Cell membrane proteins, both bound and integral, are known to preferentially accumulate at membrane locations with curvatures favorable to their shape. This is mainly due to the curvature dependent interaction between membrane proteins and their lipid environment. Here, we analyze the effects of the protein-lipid bilayer interaction energy due to mismatch between the protein shape and the principal curvatures of the surrounding bilayer. The role of different macroscopic parameters that define the interaction energy term is elucidated in relation to recent experiment in which the lateral distribution of a membrane embedded protein potassium channel KvAP is measured on a giant unilamellar lipid vesicle (reservoir) and a narrow tubular extension - a tether - kept at constant length. The dependence of the sorting ratio, defined as the ratio between the areal density of the protein on the tether and on the vesicle, on the inverse tether radius is influenced by the strength of the interaction, the intrinsic shape of the membrane embedded protein, and its abundance in the reservoir. It is described how the values of these constants can be extracted from experiments. The intrinsic principal curvatures of a protein are related to the tether radius at which the sorting ratio attains its maximum value. The estimate of the principal intrinsic curvature of the protein KvAP, obtained by comparing the experimental and theoretical sorting behavior, is consistent with the available information on its structure.

Modeling curvature-dependent subcellular localization of the small sporulation protein SpoVM in Bacillus subtilis

Recent in vivo experiments suggest that in the bacterium, Bacillus subtilis, the cue for the localization of the small sporulation protein, SpoVM, an essential factor in spore coat formation, is curvature of the bacterial plasma membrane. In vitro measurements of SpoVM adsorption to vesicles of varying sizes also find high sensitivity of adsorption to vesicle radius. This curvature-dependent adsorption is puzzling given the orders of magnitude difference in length scale between an individual protein and the radius of curvature of the cell or vesicle, suggesting protein clustering on the membrane. Here we develop a minimal model to study the relationship between curvature-dependent membrane adsorption and clustering of SpoVM. Based on our analysis, we hypothesize that the radius dependence of SpoVM adsorption observed in vitro is governed primarily by membrane tension, while for in-vivo localization of SpoVM, we propose a highly sensitive mechanism for curvature sensing based on the formation of macroscopic protein clusters on the membrane.

Two-dimensional convex-molecule fluid model for surface adsorption of proteins: effect of soft interaction on adsorption equilibria

Adsorption of proteins on membrane surfaces plays an important role in cell biological processes. In this work, we develop a two-dimensional fluid model for proteins. The protein molecules have been modeled as two-dimensional convex and soft particles. The Lennard-Jones potential for circular particles and Kihara (12,6) potential for elliptical particles with hard core have been used to model pairwise intermolecular interactions. The equation of state of the fluid model has been derived using Weeks-Chandler-Andersen decomposition and it involves three parameters, an attraction, a repulsion, and a size parameter, which depend on the shape and core size of the molecules. For validation of the model, a two-dimensional molecular dynamics simulation has been performed. Finally, the model has been applied to study the adsorption of proteins on a flat membrane. In comparison with the existing model of hard and convex particles for protein adsorption, our model predicts a higher packing fraction for the adsorption equilibria. Although the present work is based on Lennard-Jones-type interaction, it can be extended for other specific soft interactions between convex molecules. Thus the model has general applicability for any other two-dimensional adsorption systems of molecules with soft interaction.

Mesoscale computational studies of membrane bilayer remodeling by curvature-inducing proteins

Biological membranes constitute boundaries of cells and cell organelles. These membranes are soft fluid interfaces whose thermodynamic states are dictated by bending moduli, induced curvature fields, and thermal fluctuations. Recently, there has been a flood of experimental evidence highlighting active roles for these structures in many cellular processes ranging from trafficking of cargo to cell motility. It is believed that the local membrane curvature, which is continuously altered due to its interactions with myriad proteins and other macromolecules attached to its surface, holds the key to the emergent functionality in these cellular processes. Mechanisms at the atomic scale are dictated by protein-lipid interaction strength, lipid composition, lipid distribution in the vicinity of the protein, shape and amino acid composition of the protein, and its amino acid contents. The specificity of molecular interactions together with the cooperativity of multiple proteins induce and stabilize complex membrane shapes at the mesoscale. These shapes span a wide spectrum ranging from the spherical plasma membrane to the complex cisternae of the Golgi apparatus. Mapping the relation between the protein-induced deformations at the molecular scale and the resulting mesoscale morphologies is key to bridging cellular experiments across the various length scales. In this review, we focus on the theoretical and computational methods used to understand the phenomenology underlying protein-driven membrane remodeling. Interactions at the molecular scale can be computationally probed by all atom and coarse grained molecular dynamics (MD, CGMD), as well as dissipative particle dynamics (DPD) simulations, which we only describe in passing. We choose to focus on several continuum approaches extending the Canham - Helfrich elastic energy model for membranes to include the effect of curvature-inducing proteins and explore the conformational phase space of such systems. In this description, the protein is expressed in the form of a spontaneous curvature field. The approaches include field theoretical methods limited to the small deformation regime, triangulated surfaces and particle-based computational models to investigate the large-deformation regimes observed in the natural state of many biological membranes. Applications of these methods to understand the properties of biological membranes in homogeneous and inhomogeneous environments of proteins, whose underlying curvature fields are either isotropic or anisotropic, are discussed. The diversity in the curvature fields elicits a rich variety of morphological states, including tubes, discs, branched tubes, and caveola. Mapping the thermodynamic stability of these states as a function of tuning parameters such as concentration and strength of curvature induction of the proteins is discussed. The relative stabilities of these self-organized shapes are examined through free-energy calculations. The suite of methods discussed here can be tailored to applications in specific cellular settings such as endocytosis during cargo trafficking and tubulation of filopodial structures in migrating cells, which makes these methods a powerful complement to experimental studies.

Defining the free-energy landscape of curvature-inducing proteins on membrane bilayers

Curvature-sensing and curvature-remodeling proteins, such as Amphiphysin, Epsin, and Exo70, are known to reshape cell membranes, and this remodeling event is essential for key biophysical processes such as tubulation, exocytosis, and endocytosis. Curvature-inducing proteins can act as curvature sensors; they aggregate to membrane regions matching their intrinsic curvature; as well as induce curvature in cell membranes to stabilize emergent high curvature, nonspherical, structures such as tubules, discs, and caveolae. A definitive understanding of the interplay between protein recruitment and migration, the evolution of membrane curvature, and membrane morphological transitions is emerging but remains incomplete. Here, within a continuum framework and using the machinery of Monte Carlo simulations, we introduce and compare three free-energy methods to delineate the free-energy landscape of curvature-inducing proteins on bilayer membranes. We demonstrate the utility of the Widom test particle (or field) insertion methodology in computing the excess chemical potentials associated with curvature-inducing proteins on the membrane-in particular, we use this method to track the onset of morphological transitions in the membrane at elevated protein densities. We validate this approach by comparing the results from the Widom method with those of thermodynamic integration and Bennett acceptance ratio methods. Furthermore, the predictions from the Widom method have been tested against analytical calculations of the excess chemical potential at infinite dilution. Our results are useful in precisely quantifying the free-energy landscape, and also in determining the phase boundaries associated with curvature-induction, curvature-sensing, and morphological transitions. This approach can be extended to studies exploring the role of thermal fluctuations and other external (control) variables, such as membrane excess area, in shaping curvature-mediated interactions on bilayer membranes.

Bending lipid membranes: experiments after W. Helfrich's model

Current description of biomembrane mechanics originates for a large part from W. Helfrich's model. Based on his continuum theory, many experiments have been performed in the past four decades on simplified membranes in order to characterize the mechanical properties of lipid membranes and the contribution of polymers or proteins. The long-term goal was to develop a better understanding of the mechanical properties of cell membranes. In this paper, we will review representative experimental approaches that were developed during this period and the main results that were obtained.

Membrane shape modulates transmembrane protein distribution

Although membrane shape varies greatly throughout the cell, the contribution of membrane curvature to transmembrane protein targeting is unknown because of the numerous sorting mechanisms that take place concurrently in cells. To isolate the effect of membrane shape, we used cell-sized giant unilamellar vesicles (GUVs) containing either the potassium channel KvAP or the water channel AQP0 to form membrane nanotubes with controlled radii. Whereas the AQP0 concentrations in flat and curved membranes were indistinguishable, KvAP was enriched in the tubes, with greater enrichment in more highly curved membranes. Fluorescence recovery after photobleaching measurements showed that both proteins could freely diffuse through the neck between the tube and GUV, and the effect of each protein on membrane shape and stiffness was characterized using a thermodynamic sorting model. This study establishes the importance of membrane shape for targeting transmembrane proteins and provides a method for determining the effective shape and flexibility of membrane proteins.

Bending "on the rocks"--a cocktail of biophysical modules to build endocytic pathways

Numerous biological processes rely on endocytosis. The construction of endocytic pits is achieved by a bewildering complexity of biochemical factors that function in clathrin-dependent and -independent pathways. In this review, we argue that this complexity can be conceptualized by a deceptively small number of physical principles that fall into two broad categories: passive mechanisms, such as asymmetric transbilayer stress, scaffolding, line tension, and crowding, and active mechanisms driven by mechanochemical enzymes and/or cytoskeleton. We illustrate how the functional identity of biochemical modules depends on system parameters such as local protein density on membranes, thus explaining some of the controversy in the field. Different modules frequently operate in parallel in the same step and often are shared by apparently divergent uptake processes. The emergence of a novel endocytic classification system may thus be envisioned in which functional modules are the elementary bricks.

FBAR syndapin 1 recognizes and stabilizes highly curved tubular membranes in a concentration dependent manner

Syndapin 1 FBAR, a member of the Bin-amphiphysin-Rvs (BAR) domain protein family, is known to induce membrane curvature and is an essential component in biological processes like endocytosis and formation and growth of neurites. We quantify the curvature sensing of FBAR on reconstituted porcine brain lipid vesicles and show that it senses membrane curvature at low density whereas it induces and reinforces tube stiffness at higher density. FBAR strongly up-concentrates on the high curvature tubes pulled out of Giant Unilamellar lipid Vesicles (GUVs), this sorting behavior is strongly amplified at low protein densities. Interestingly, FBAR from syndapin 1 has a large affinity for tubular membranes with curvatures larger than its own intrinsic concave curvature. Finally, we studied the effect of FBAR on membrane relaxation kinetics with high temporal resolution and found that the protein increases relaxation time of the tube holding force in a density-dependent fashion.

Steric pressure between membrane-bound proteins opposes lipid phase separation

Cellular membranes are densely crowded with a diverse population of integral and membrane-associated proteins. In this complex environment, lipid rafts, which are phase-separated membrane domains enriched in cholesterol and saturated lipids, are thought to organize the membrane surface. Specifically, rafts may help to concentrate proteins and lipids locally, enabling cellular processes such as assembly of caveolae, budding of enveloped viruses, and sorting of lipids and proteins in the Golgi. However, the ability of rafts to concentrate protein species has not been quantified experimentally. Here we show that when membrane-bound proteins become densely crowded within liquid-ordered membrane regions, steric pressure arising from collisions between proteins can destabilize lipid phase separations, resulting in a homogeneous distribution of proteins and lipids over the membrane surface. Using a reconstituted system of lipid vesicles and recombinant proteins, we demonstrate that protein-protein steric pressure creates an energetic barrier to the stability of phase-separated membrane domains that increases in significance as the molecular weight of the proteins increases. Comparison with a simple analytical model reveals that domains are destabilized when the steric pressure exceeds the approximate enthalpy of membrane mixing. These results suggest that a subtle balance of free energies governs the stability of phase-separated cellular membranes, providing a new perspective on the role of lipid rafts as concentrators of membrane proteins.