The window of a few tens to a few hundred nanometers in length scale is a booming field in lipid membrane research, owing largely to two reasons. First, many exciting biophysical and cell biological ...processes take place within it. Second, experimental techniques manage to zoom in on this sub‐optical scale, while computer simulations zoom out to system sizes previously unattainable, and both will be meeting soon. This paper reviews a selection of questions and concepts in this field and demonstrates that they can often be favorably addressed with highly simplified simulation models. Among the topics discussed are membrane adhesion to substrates, mixed lipid bilayers, lipid curvature coupling, pore formation by antimicrobial peptides, composition‐driven protein aggregation, and curvature driven vesiculation.
Biomembranes exhibit fascinating physics covering orders of magnitude in length and time scales. Computer simulations within the mesoscopic window require systematically reduced descriptions of lipids and proteins to overcome otherwise insurmountable barriers in computational feasibility. This article illustrates recent advances gained by remarkably simple coarse grained models.
Lipid tilt affects the energy of membrane deformations on scales comparable to the membrane’s thickness. The surface divergence of the tilt field creates a local spontaneous curvature, while tilt ...itself is quadratically penalized with a strength given by the tilt modulus. Traditionally, this modulus is determined by measuring the power spectrum of lipid orientation fluctuations. Here we present a novel approach which does not rely on fluctuations but instead exploits the fact that curvature gradients induce a tilt field. Its implementation extends a technique previously developed by us for localizing the position of the pivotal plane in buckling simulations, which quantifies the lipid imbalance across segments cut out from a complete buckle. Lipid tilt affects this count in a predictable way, and the signal can be quantified well enough to back out the tilt modulusat no additional cost and with about 5% precision for not too coarse models. We apply our technique to three lipid models of very different resolution: the highly coarse grained Cooke model, and two versions of DMPC, using both the (less highly coarse grained) MARTINI and the (united atom) Berger force field. For Cooke, we find an effective bilayer tilt modulus of 29 ± 9 pN/nm, and for the less generic DMPC lipid, we find 115 ± 6 pN/nm for MARTINI and 39 ± 2 pN/nm for Berger, both in reasonable agreement with existing studies for these models. We also show that the position of the pivotal plane for Berger DMPC lies just below the glycerol backboneunlike for MARTINI DMPC, where this plane lies closer to the middle of the lipid.
Display omitted
•Lipid membranes can be described by a curvature elastic Hamiltonian.•Its covariant differential geometric description is extensively reviewed.•A variation under geometric constraints ...reveals conserved stresses and torques.•Examples are provided for how to efficiently calculate with these objects.
A fluid lipid membrane transmits stresses and torques that are fully determined by its geometry. They can be described by a stress- and torque-tensor, respectively, which yield the force or torque per length through any curve drawn on the membrane's surface. In the absence of external forces or torques the surface divergence of these tensors vanishes, revealing them as conserved quantities of the underlying Euler–Lagrange equation for the membrane's shape. This review provides a comprehensive introduction into these concepts without assuming the reader's familiarity with differential geometry, which instead will be developed as needed, relying on little more than vector calculus. The Helfrich Hamiltonian is then introduced and discussed in some depth. By expressing the quest for the energy-minimizing shape as a functional variation problem subject to geometric constraints, as proposed by Guven (2004), stress- and torque-tensors naturally emerge, and their connection to the shape equation becomes evident. How to reason with both tensors is then illustrated with a number of simple examples, after which this review concludes with four more sophisticated applications: boundary conditions for adhering membranes, corrections to the classical micropipette aspiration equation, membrane buckling, and membrane mediated interactions.
We present an implicit solvent coarse-grained (CG) model for quantitative simulations of 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) bilayers. The absence of explicit solvent enables ...membrane simulations on large length and time scales at moderate computational expense. Despite improved computational efficiency, the model preserves chemical specificity and quantitative accuracy. The bonded and nonbonded interactions together with the effective cohesion mimicking the hydrophobic effect were systematically tuned by matching structural and mechanical properties from experiments and all-atom bilayer simulations, such as saturated area per lipid, radial distribution functions, density and pressure profiles across the bilayer, P 2 order, etc. The CG lipid model is shown to self-assemble into a bilayer starting from a random dispersion. Its line tension and elastic properties, such as bending and stretching modulus, are semiquantitatively consistent with experiments. The effects of (i) reduced molecular friction and (ii) more efficient integration combine to an overall speed-up of 3−4 orders of magnitude compared to all-atom bilayer simulations. Our CG lipid model is especially useful for studies of large-scale phenomena in membranes that nevertheless require a fair description of chemical specificity, e.g., membrane patches interacting with movable and transformable membrane proteins and peptides.
An approach based on effective field theory (EFT) is discussed and applied to the problem of surface-mediated interactions between rigid inclusions of circular footprint on a membrane. Instead of ...explicitly constraining the surface fluctuations in accord with the boundary conditions around the inclusions, the EFT formalism rewrites the theory; the Hamiltonian of a freely fluctuating surface is augmented by pointwise localized terms that capture the same constraints. This allows one to compute the interaction free energy as an asymptotic expansion in inverse separations in a systematic, efficient, and transparent way. Both entropic (fluctuation-induced, Casimir-like) and curvature-elastic (ground-state) forces are considered. Our findings include higher-order corrections to known asymptotic results, on both the pair and the multibody levels. We also show that the few previous attempts in the literature at predicting subleading orders missed some terms due to an uncontrolled point-particle approximation.
The importance of curvature as a structural feature of biological membranes has been recognized for many years and has fascinated scientists from a wide range of different backgrounds. On the one ...hand, changes in membrane morphology are involved in a plethora of phenomena involving the plasma membrane of eukaryotic cells, including endo- and exocytosis, phagocytosis and filopodia formation. On the other hand, a multitude of intracellular processes at the level of organelles rely on generation, modulation, and maintenance of membrane curvature to maintain the organelle shape and functionality. The contribution of biophysicists and biologists is essential for shedding light on the mechanistic understanding and quantification of these processes. Given the vast complexity of phenomena and mechanisms involved in the coupling between membrane shape and function, it is not always clear in what direction to advance to eventually arrive at an exhaustive understanding of this important research area. The 2018 Biomembrane Curvature and Remodeling Roadmap of Journal of Physics D: Applied Physics addresses this need for clarity and is intended to provide guidance both for students who have just entered the field as well as established scientists who would like to improve their orientation within this fascinating area.
The force needed to buckle a thin elastic surface is proportional to its bending rigidity. This fact suggests using a buckling setup to measure the bending modulus of lipid membranes. Extending the ...work of Noguchi Phys. Rev. E 83, 061919 (2011), we systematically derive highly accurate analytical expressions for the forces along and perpendicular to the buckle, and we elucidate some of their counterintuitive properties using the framework of a surface stress tensor. Furthermore, we estimate the corrections to buckling forces due to thermal fluctuations and find them significant only for stresses along the ridges. We then apply this buckling protocol to four different lipid membrane models, which widely differ in their level of resolution and the treatment of solvent, and show that in all cases buckling is a reliable and accurate means for measuring their rigidity. Finally, we show that monitoring both stresses and energies during a simulation offers additional insights into the thermodynamics of curvature elasticity and permits one to predict the bending rigidity for a range of temperatures around the actual simulation temperature.
Lipid bilayers can exhibit asymmetric states, in which the physical characteristics of one leaflet differ from those of the other. This most visibly manifests in a different lipid composition, but it ...can also involve opposing lateral stresses in each leaflet that combine to an overall vanishing membrane tension. Here, we use theoretical modeling and coarse-grained simulation to explore the interplay between a compositional asymmetry and a nonvanishing differential stress. Minimizing the total elastic energy leads to a preferred spontaneous curvature that balances torques due to both bending moments and differential stress, with sometimes unexpected consequences. For instance, asymmetric flat bilayers, whose specific areas in each leaflet are matched to those of corresponding tensionless symmetric flat membranes, still exhibit a residual differential stress because the conditions of vanishing area strain and vanishing bending moment differ. We also measure the curvature rigidity of asymmetric bilayers and find that a sufficiently strong differential stress, but not compositional asymmetry alone, can increase the bending modulus. The likely cause is a stiffening of the compressed leaflet, which appears to be related to its gel transition but not identical with it. We finally show that the impact of cholesterol on differential stress depends on the relative strength of elastic and thermodynamic driving forces: if cholesterol solvates equally well in both leaflets, it will redistribute to cancel both leaflet tensions almost completely, but if its partitioning free energy prefers one leaflet over the other, the resulting distribution bias may even create differential stress. Because cells keep most of their lipid bilayers in an asymmetric nonequilibrium steady state, our findings suggest that biomembranes are elastically more complex than previously thought: besides a spontaneous curvature, they might also exhibit significant differential stress, which could strongly affect their curvature energetics.