Molecular dynamics simulations of lipid bilayer ribbons have been performed to investigate the structures and line tensions associated with free bilayer edges. Simulations carried out for dioleoyl phosphatidylcholine with three different force-field parameter sets yielded edge line tensions of 45 ± 2 pN, over 50% greater than the most recently reported experimentally determined value for this lipid. Edge tensions obtained from simulations of a series of phosphatidylcholine lipid bilayer ribbons with saturated acyl tails of length 12-16 carbons and with monounsaturated acyl tails of length 14-18 carbons could be correlated with the excess area associated with forming the edge, through a two-parameter fit. Saturated-tail lipids underwent local thickening near the edge, producing denser packing that correlated with lower line tensions, while unsaturated-tail lipids showed little or no local thickening. In a dipalmitoyl phosphatidylcholine ribbon initiated in a tilted gel-phase structure, lipid headgroups tended to tilt toward the nearer edge producing a herringbone pattern, an accommodation that may account for the reported edge-induced stabilization of an ordered structure at temperatures near a lipid gel-fluid phase transition.

The dynamics of the gel to fluid phase transformation in 100 nm large unilamellar vesicles (LUV) of 1,2-dipalmitoyl(d62)-sn-glycero-3-phosphocholine (d62-DPPC), has been studied by laser-induced temperature-jump initiation coupled with time-resolved infrared spectroscopy and by MD simulations. The infrared transients that characterize the temperature dependent phase transformation are complex, extending from the nanosecond to the millisecond time scales. An initial fast (submicrosecond) component can be modeled by partial melting of the gel domains, initiated at pre-existing defects at the edges of the faceted structure of the gel phase. Molecular dynamics simulations support the model of fast melting from edge defects. The extent of melting during the fast phase is limited by the area expansion on melting, which leads to a surface pressure that raises the effective melting temperature. Subsequent melting is observed to follow highly stretched exponential kinetics, consistent with collective relaxation of the surface pressure through a hierarchy of surface undulations with different relaxation times. The slowest step is water diffusion through the bilayer to allow the vesicle volume to grow along with its expanded surface area. The results demonstrate that the dominant relaxation in the gel to fluid phase transformation in response to a large T-jump perturbation (compared to the transition width) is fast (submicrosecond), which has important practical and fundamental consequences.

No abstract