Membrane remodelling plays an important role in cellular tasks such as endocytosis, vesiculation and protein sorting, and in the biogenesis of organelles such as the endoplasmic reticulum or the Golgi apparatus. It is well established that the remodelling process is aided by specialized proteins that can sense as well as create membrane curvature, and trigger tubulation when added to synthetic liposomes. Because the energy needed for such large-scale changes in membrane geometry significantly exceeds the binding energy between individual proteins and between protein and membrane, cooperative action is essential. It has recently been suggested that curvature-mediated attractive interactions could aid cooperation and complement the effects of specific binding events on membrane remodelling. But it is difficult to experimentally isolate curvature-mediated interactions from direct attractions between proteins. Moreover, approximate theories predict repulsion between isotropically curving proteins. Here we use coarse-grained membrane simulations to show that curvature-inducing model proteins adsorbed on lipid bilayer membranes can experience attractive interactions that arise purely as a result of membrane curvature. We find that once a minimal local bending is realized, the effect robustly drives protein cluster formation and subsequent transformation into vesicles with radii that correlate with the local curvature imprint. Owing to its universal nature, curvature-mediated attraction can operate even between proteins lacking any specific interactions, such as newly synthesized and still immature membrane proteins in the endoplasmic reticulum.

Dissipative particle dynamics is used to extract the material parameters (bending and area stretch moduli) of a bilayer membrane patch. Some experiments indicate that the area stretch modulus of lipid vesicles varies little as the chain length of the lipids composing the bilayer increases. Here we show that making the interactions between the hydrophilic head groups of the model amphiphiles proportional to the hydrophobic tail length reproduces the above result for the area stretch modulus. We also show that the area stretch modulus of bilayers composed of amphiphiles with the same number of tail beads but with asymmetric chains is less than that of bilayers with symmetric chains. The effects on the bilayer density and lateral stress profiles of changes to the amphiphile architecture are also presented.

We investigate the interactions between lipid bilayers and amphiphilic peptides using a solvent-free coarse-grained simulation technique. In our model, each lipid is represented by one hydrophilic and three hydrophobic beads. The amphiphilic peptide is modeled as a hydrophobic-hydrophilic cylinder with hydrophilic caps. We find that with increasing peptide-lipid attraction the preferred state of the peptide changes from desorbed, to adsorbed, to inserted. A single peptide with weak attraction binds on the bilayer surface, while one with strong attraction spontaneously inserts into the bilayer. We show how several peptides, which individually bind only to the bilayer surface, cooperatively insert. Furthermore, hydrophilic strips along the peptide cylinder induce the formation of multipeptide pores, whose size and morphology depend on the peptides' overall hydrophilicity, the distribution of hydrophilic residues, and the peptide-peptide interactions. Strongly hydrophilic peptides insert less readily, but prove to be more destructive to bilayer integrity.

The material parameters (area stretch modulus and bending rigidity) of two-component amphiphilic membranes are determined from dissipative particle dynamics simulations. The preferred area per molecule for each species is varied so as to produce homogeneous mixtures or nonhomogeneous mixtures that form domains. If the latter mixtures are composed of amphiphiles with the same tail length, but different preferred areas per molecule, their material parameters increase monotonically as a function of composition. By contrast, mixtures of amphiphiles that differ in both tail length and preferred area per molecule form both homogeneous and nonhomogeneous mixtures that both exhibit smaller values of their material properties compared to the corresponding pure systems. When the same nonhomogeneous mixtures of amphiphiles are assembled into planar membrane patches and vesicles, the resulting domain shapes are different when the bending rigidities of the domains are sufficiently different. Additionally, both bilayer and monolayer domains are observed in vesicles. We conclude that the evolution of the domain shapes is influenced by the high curvature of the vesicles in the simulation, a result that may be relevant for biological vesicle membranes.