Molecular studies of the gel to liquid-crystalline phase transition for fully hydrated DPPC and DPPE bilayers

Molecular dynamics simulations were used for a comprehensive study of the structural properties of saturated lipid bilayers, DPPC and DPPE, near the main phase transition. Though the chemical structure of DPPC and DPPE are largely similar (they only differ in the choline and ethanolamine groups), their transformation process from a gel to a liquid-crystalline state is contrasting. For DPPC, three distinct structures can be identified relative to the melting temperature (T_m): below T_m with “mixed” domains consisting of lipids that are tilted with partial overlap of the lipid tails between leaflet; near T_m with a slight increase in the average area per lipid, resulting in a rearrangement of the lipid tails and an increase in the bilayer thickness; and above T_m with unhindered lipid tails in random motion resulting in an increase in %gauche formed and increase in the level of interdigitation between lipid leaflets. For DPPE, the structures identified were below T_m with “ordered” domains consisting of slightly tilted lipid tails and non-overlapping lipid tails between leaflets, near T_m with minimal rearrangement of the lipids as the bilayer thickness reduces slightly with increasing temperature, and above T_m with unhindered lipid tails as that for DPPC. For DPPE, most of the lipid tails do not overlap as observed to DPPC, which is due to the tight packing of the DPPE molecules. The non-overlapping behavior of DPPE above T_m is confirmed from the density profile of the terminal carbon atoms in each leaflet, which shows a narrow distribution near the center of the bilayer core. This study also demonstrates that atomistic simulations are capable of capturing the phase transition behavior of lipid bilayers, providing a rich set of molecular and structural information at and near the transition state.