Self-consistent mean-field model for palmitoyloleoylphosphatidylcholine-palmitoyl sphingomyelin-cholesterol lipid bilayers

Cholesterol Distribution in Small Unilamellar Vesicles

The plasma membrane has heterogeneity, where the constituents are not only asymmetrically distributed between the inner leaflet and the outer leaflet but also laterally organized within each leaflet. There is still an ongoing controversy over the cholesterol distribution between the two leaflets, and it is also significant to explore the lateral flow and localization of cholesterol. The unilamellar vesicle is extensively employed as a simplified model of the plasma membrane in research studies. In this work, we study the spontaneous spatial distribution of cholesterol in the small unilamellar vesicles with a single type of phospholipid as the bilayer backbone through coarse-grained molecular dynamics simulations. The results show that in a spherical vesicle under ambient pressure, cholesterol is more abundant in the inner leaflet than in the outer leaflet. As the vesicle is deformed under pressurization, the net lateral flow of cholesterol in the two leaflets is in the exactly opposite directions, finally leading to a distribution strongly associated with the monolayer curvature. One of the possible explanations for our results is from the point of view of curvature elastic energy. Another possibility is from the point of view of stress: according to the correlation between the cholesterol distribution and the tail angles of lipids, it is suggested that the possible governing mechanism for the distribution of cholesterol in a membrane is to alleviate the mismatch of stress between the two leaflets. Additionally, we obtain the effect of cholesterol infiltration on the bending modulus and the spontaneous curvature of the vesicles.

Critical size dependence of domain formation observed in coarse-grained simulations of bilayers composed of ternary lipid mixtures

Model cellular membranes are known to form micro- and macroscale lipid domains dependent on molecular composition. The formation of macroscopic lipid domains by lipid mixtures has been the subject of many simulation investigations. We present a critical study of system size impact on lipid domain phase separation into liquid-ordered and liquid-disordered macroscale domains in ternary lipid mixtures. In the popular di-C16:0 PC:di-C18:2 PC:cholesterol at 35:35:30 ratio mixture, we find systems with a minimum of 1480 lipids to be necessary for the formation of macroscopic phase separated domains and systems of 10 000 lipids to achieve structurally converged conformations similar to the thermodynamic limit. To understand these results and predict the behavior of any mixture forming two phases, we develop and investigate an analytical Flory-Huggins model which is recursively validated using simulation and experimental data. We find that micro- and macroscale domains can coexist in ternary mixtures. Additionally, we analyze the distributions of specific lipid-lipid interactions in each phase, characterizing domain structures proposed based on past experimental studies. These findings offer guidance in selecting appropriate system sizes for the study of phase separations and provide new insights into the nature of domain structure for a popular ternary lipid mixture.

Cholesterol, sphingolipids, and glycolipids: what do we know about their role in raft-like membranes?

Lipids rafts are considered to be functional nanoscale membrane domains enriched in cholesterol and sphingolipids, characteristic in particular of the external leaflet of cell membranes. Lipids, together with membrane-associated proteins, are therefore considered to form nanoscale units with potential specific functions. Although the understanding of the structure of rafts in living cells is quite limited, the possible functions of rafts are widely discussed in the literature, highlighting their importance in cellular functions. In this review, we discuss the understanding of rafts that has emerged based on recent atomistic and coarse-grained molecular dynamics simulation studies on the key lipid raft components, which include cholesterol, sphingolipids, glycolipids, and the proteins interacting with these classes of lipids. The simulation results are compared to experiments when possible.

Can adding oil control domain formation in binary amphiphile bilayers?

Bilayers formed of two species of amphiphile of different chain lengths may segregate into thinner and thicker domains composed predominantly of the respective species. Using a coarse-grained mean-field model, we investigate how mixing oil with the amphiphiles affects the structure and thickness of the bilayer at and on either side of the boundary between two neighbouring domains. In particular, we find that oil molecules whose chain length is close to that of the shorter amphiphiles segregate to the thicker domain. This smooths the surface of the hydrophobic bilayer core on this side of the boundary, reducing its area and curvature and their associated free-energy penalties. The smoothing effect is weaker for oil molecules that are shorter or longer than this optimum value: short molecules spread evenly through the bilayer, while long molecules swell the thicker domain, increasing the surface area and curvature of the bilayer core in the interfacial region. Our results show that adding an appropriate oil could make the formation of domain boundaries more or less favourable, raising the possibility of controlling the domain size distribution.