Interbilayer repulsion forces between tension-free lipid bilayers from simulation

Abstract

Here we report studies on biologically important intermembrane repulsion forces using molecular dynamics (MD) simulations and experimental (osmotic stress) investigations of repulsion forces between 1-palmitoyl-2-oleyl-sn-glycero-3-phosphocholine bilayers. We show that the repulsion between tensionfree membranes can be determined from MD simulations by either (i) simulating membrane stacks under different hydration conditions (unrestrained setup) and monitoring the change in the area per lipid upon dehydration or (ii) simulating two single punctured membranes immersed in a water reservoir and controlling the center-of-mass distance between the bilayers using an external potential (umbrella sampling setup). Despite the coarse-grained nature of the (MARTINI) model employed, the disjoining pressure profiles obtained from the simulations are in good agreement with our experiments.Remarkably, the two setups behave very differently in terms of membrane structure, as explained by considerations using elasticity theory, and the balance of interactions. In the unrestrained setup, dehydration decreases the area per lipid and lipid entropy. Dehydration in the umbrella sampling setup, in contrast, leads to an increase in area per lipid and lipid entropy. Hence, in the latter case, entropic effects from protrusion and zippering forces appear to be overcompensated by the entropy gain due to the disorder emerging from the expansion of the bilayers. The balance of interactions involves near cancellation of large opposing terms, for which also intramembrane and water-water interactions are important, and which appears to be largely a consequence, rather than the cause, of the intermembrane repulsion. Hence, care must be taken when drawing conclusions on the origin of intermembrane repulsion from thermodynamic analyses.

SNARE ‐mediated membrane fusion arrests at pore expansion to regulate the volume of an organelle

“…Fusion pore closure is opposed by an increase in bending stress (Ryham et al , ) and the unfavorable dehydration of lipid head groups (hydration repulsion) (Smirnova et al , ). The resulting, here‐predicted minimal pore size (~1.5 nm) should in principle allow free passage of small molecules, such as the fluorescent lumenal markers that we had used.…”

Effect of Sodium and Chloride Binding on a Lecithin Bilayer. A Molecular Dynamics Study

“…The area per lipid at the infinite bilayer separation may be expressed in terms of the disjoining pressure [38]: $P_{\mathrm{dis}}={\kappa }_{\mathrm{A}}{d}_{\mathrm{normalw}}^{-1}(1-a_{\mathrm{L},\mathrm{o}}^{-1}{a}_{\mathrm{L},\mathrm{o},\infty }),$ where $d_{\mathrm{normalw}}$ is the water layer thickness. Setting $P_{\mathrm{dis}}$ equal to the disjoining pressure present in an MD simulation, i.e., [38]: $P_{\mathrm{MD}}={\kappa }_{\mathrm{A}}{d}_{\mathrm{normalw}}^{-1}(1-a_{\mathrm{L},\mathrm{o}}^{-1}{a}_{\mathrm{normalL}}),$ allows solving Equation (5) for $a_{\mathrm{L},\mathrm{o},\infty }$. Comparison of $a_{\mathrm{L},\mathrm{o},\infty }$ with the actual area per lipid observed during the simulation gives an estimate of the magnitude of finite-size effects on the simulated membrane structure.…”

The Interaction between Soft Interfaces: Forces and Structural Aspects

“…[101] The monomer profiles revealed considerable brush interpenetration, as shown in Figure 7B. [103] Specific binding between opposing membranes was considered by Hu et al using solvent-implicit CG simulations. [102] Also the interactions between lipid bilayers have been addressed using CG simulations.…”