Molecular dynamics investigation of bond ordering of unsaturated lipids in monolayers

Molecular dynamics simulations of three model lipid monolayers of 2,3-diacyl-D-glycerolipids, that contained stearoyl (18:0) in the position 3 and oleoyl (18:ω9cis), linoleoyl (18:2ω6cis), or linolenoyl (18:3ω3cis) in the position 2, have been carried out. The simulation systems consisted of 24 lipid molecules arranged in a rectangular simulation cell, with periodic boundary conditions in the surface plane. 1 nanosecond simulations were performed at T = 295 K. C-C and C-H bond order parameter profiles and the bond orientation distributions about the monolayer normal have been calculated. The relation of the distributions to the order parameters was analyzed in terms of maxima and widths of the distributions. The cis double bond order parameter is found to be higher than those of adjacent single C-C bonds. The widths of the two distributions of C-H bonds of the cis double bond segment in di- and triunsaturated molecules are much smaller than that obtained for methylene group located between the double bonds. The bond orientation distribution function widths depend on both the segment location in the chain and the segment chemical structure.