Characterization of lipid membrane dynamics by simulation: I. Torsion angle motions of the linear chains

The torsion angle motions, generated from molecular dynamics (MD) simulations, of the two aliphatic chains of 1,2-dimyristoyl-sn-glycero-3-phosphatidylcholine (DMPC) in its lipid monolayer were evaluated by comparing these motions to those of an equivalent isolated (free) n-alkane chain, and the same n-alkane chain in its crystal lattice. The time-dependent autocorrelation and (1,2)-, (1,3)-, (1,4)-, and (1,5)-cross-correlation functions were constructed to analyze the torsion angle motions. It was found that the torsion angle motions of the DMPC lipid monolayer aliphatic chains are intermediate to those of the free n-alkane chain and the same n-alkane chain in its crystal lattice, particularly for short correlation times. The torsion angle motions of the aliphatic chains of DMPC are also found to be essentially independent of the charge state on the head group. The linear aliphatic chains of a DMPC lipid monolayer behave most like the isolated n-alkane chains with respect to torsion angle flexibility, even though the pairs of aliphatic chains of each DMPC are part of an ordered monolayer assembly. The aliphatic chains of the DMPC molecules in their monolayer exhibit at least two types of wave motions. One of the wave motions is the same in form, though somewhat more diffuse, as a traveling wave found in n-alkane crystals. The other wave motion involves major torsion angle transitions, and has some characteristics of the soliton properties observed in n-alkane crystals near their respective melt transition temperatures.