Simulation studies of pore and domain formation in a phospholipid monolayer

Despite extensive study the phase behavior of phospholipid monolayers at an air-water interface is still not fully understood. In particular recent vibrational sum-frequency generation (VSFG) spectra of DPPC monolayers as a function of area density show a sharp transition in the order of the lipid chains at 1.10 nm²/molecule. This is in a region where the lateral pressure as a function of area is effectively constant. We have investigated the nature of this transition by studying the phase behavior of DPPC monolayers as a function of area density using molecular-dynamics simulations. The changes in order within the monolayer as a function of area density correlate well with the experimental signal. At 0.58 nm²/molecule we observe the onset of lateral separation of highly ordered and disordered lipids, indicating the coexistence of a gel-like liquid condensed and a fluidlike liquid expanded phase. At 0.97 nm²/molecule the monolayer ruptures, marking the onset of the liquid-gas (G) coexistence region. This is much earlier than suggested by fluorescence microscopy results and implies that at the point of rupture, the initial pores have an equilibrium size smaller than ~500 nm in diameter. The rupture of the monolayer leads to a sharp increase in the overall lipid order that explains the sharp transition observed in the VSFG measurements. VSFG measurements thus may represent a sensitive means to determine the onset of the liquid-gas (G) coexistence region for such systems.