Understanding the phase behavior of coarse-grained model lipid bilayers through computational calorimetry

We study the phase behavior of saturated lipids as a function of temperature and tail length for two coarse-grained models: the soft-repulsive model typically employed with dissipative particle dynamics (DPD) and the MARTINI model. We characterize the simulated transitions through changes in structural properties, and we introduce a computational method to monitor changes in enthalpy, as is done experimentally with differential scanning calorimetry. The lipid system experimentally presents four different bilayer phases - subgel, gel, ripple, and fluid - and the DPD model describes all of these phases structurally while MARTINI describes a single order-disorder transition between the gel and the fluid phases. Given both models' varying degrees of success in displaying accurate structural and thermodynamic signatures, there is an overall satisfying extent of agreement for the coarse-grained models. We also study the lipid dynamics displayed by these models for the various phases, discussing this dynamics with relation to fidelity to experiment and computational efficiency.