Abstract

Glycolipids such as GM1 have bulky head groups consisting of several monosaccharides. When these lipids are added to phospholipid bilayers, they generate large membrane curvatures even for small compositional asymmetries between the two leaflets of the bilayers. On the micrometer scale, these bilayer asymmetries lead to the spontaneous tubulation of giant vesicles as recently observed by optical microscopy. Here, we study these mixed membranes on the nanometer scale using coarse-grained molecular simulations. The membrane composition is defined by the mole fractions ϕ1 and ϕ2 of the large-head lipid in the two leaflets of the bilayer. Symmetric membranes are obtained for ϕ1 = ϕ2 ≡ ϕle, and asymmetric ones for ϕ1 ≠ ϕ2. In both cases, we compute the density and stress profiles across the membranes. The stress profiles are used to identify the tensionless states of the membranes. Symmetric and tensionless bilayers are found to be stable within the whole composition range 0 ≤ ϕle ≤ 1. For these symmetric bilayers, both the area compressibility modulus and the bending rigidity are found to vary non-monotonically with the leaflet mole fraction ϕle. For asymmetric bilayers, we compute the product of bending rigidity and spontaneous curvature from the first moment of the stress profile and determine the bending rigidities of the asymmetric membranes using the ϕle-dependent rigidities of the single leaflets. When we combine these results, the compositional asymmetry ϕ1 - ϕ2 is found to generate the spontaneous curvature (ϕ1 - ϕ2)/(0.63 ℓme) with the membrane thickness ℓme ≃ 4 nm. Therefore, the spontaneous curvature increases linearly with the compositional asymmetry. Furthermore, the small compositional asymmetry ϕ1 - ϕ2 = 0.04 leads to the large spontaneous curvature 1/(63 nm) and the increased asymmetry ϕ1 - ϕ2 = 0.2 generates the huge spontaneous curvature 1/(13 nm). These large values of the spontaneous curvature will facilitate future simulation studies of various membrane processes such as bud formation and nanoparticle engulfment.

Highlights

  • In spite of their diversity and molecular complexity, biomembranes have a universal architecture which is based on fluid bilayers of lipids and membrane proteins.1 Within these bilayers, the hydrophilic head groups of the lipids are positioned between their hydrophobic chains and the aqueous solutions

  • We compute the product of bending rigidity and spontaneous curvature from the first moment of the stress profile and determine the bending rigidities of the asymmetric membranes using the φle-dependent rigidities of the single leaflets

  • We study these bilayer membranes by molecular simulations based on Dissipative Particle Dynamics (DPD)

Read more

Summary

INTRODUCTION

In spite of their diversity and molecular complexity, biomembranes have a universal architecture which is based on fluid bilayers of lipids and membrane proteins. Within these bilayers, the hydrophilic head groups of the lipids are positioned between their hydrophobic chains and the aqueous solutions. Generate large spontaneous curvatures for small compositional asymmetries and are useful to elucidate the influence of this curvature on various membrane processes such as bud formation and nanoparticle engulfment We study these bilayer membranes by molecular simulations based on Dissipative Particle Dynamics (DPD).. In our DPD study, the membrane composition of the binary mixtures is described by the mole fractions φ1 and φ2 of the large-head lipids in the upper and lower leaflets of the bilayer. Various quantities such as area compressibility, bilayer thickness, and bending rigidity are related to nanoscopic length and time scales We study these quantities using the coarse-grained DPD method to simulate bilayer membranes.. After combining all forces, we integrate Newton’s equation of motion for each bead using a modified version of the velocity-Verlet algorithm and study the time evolution of the system

Leaflet compositions and mole fractions
Simulation parameters
Stress profiles and tensionless membranes
Density profiles of symmetric bilayers
Stress profiles of symmetric bilayers
Density profiles of asymmetric bilayers
Stress profiles of asymmetric bilayers
Molecular areas within individual leaflets
Symmetric and tensionless bilayers
Asymmetric and tensionless bilayers
Area compressibility modulus
Thickness of bilayer membranes
Bending rigidity of symmetric bilayers
Bending rigidities of bilayer leaflets
Composition dependence of bending rigidity
Negative first moment of stress profile
Composition dependence of spontaneous curvature
Comparison with previous simulation studies
Outlook on membrane buds and necks
Findings
SUMMARY AND CONCLUSION
Full Text
Published version (Free)

Talk to us

Join us for a 30 min session where you can share your feedback and ask us any queries you have

Schedule a call