In this paper, we present a coarse-grained model of a hydrated saturated phospholipid bilayer (dimyristoylphosphatidylcholine, DMPC) containing cholesterol that we study using a hybrid dissipative particle dynamics-Monte Carlo method. This approach allows us to reach the time and length scales necessary to study structural and mechanical properties of the bilayer at various temperatures and cholesterol concentrations. The properties studied are the area per lipid, condensation, bilayer thickness, tail order parameters, bending modulus, and area compressibility. Our model quantitatively reproduces most of the experimental effects of cholesterol on these properties and reproduces the main features of the experimental phase and structure diagrams. We also present all-atom simulation results of the system and use these results to further validate the structure of our coarse-grained bilayer. On the basis of the changes in structural properties, we propose a temperature-composition structure diagram, which we compare with the experimental phase and structure diagrams. Attention is paid to the reliability and interpretation of the model and simulation method and of the different experimental techniques. The lateral organization of cholesterol in the bilayer is discussed.

First principles molecular dynamics simulation protocol is established using revised functional of Perdew-Burke-Ernzerhof (revPBE) in conjunction with Grimme's third generation of dispersion (D3) correction to describe properties of water at ambient conditions. This study also demonstrates the consistency of the structure of water across both isobaric (NpT) and isothermal (NVT) ensembles. Going beyond the standard structural benchmarks for liquid water, we compute properties that are connected to both local structure and mass density fluctuations that are related to concepts of solvation and hydrophobicity. We directly compare our revPBE results to the Becke-Lee-Yang-Parr (BLYP) plus Grimme dispersion corrections (D2) and both the empirical fixed charged model (SPC/E) and many body interaction potential model (MB-pol) to further our understanding of how the computed properties herein depend on the form of the interaction potential.

No abstract