Computer simulation study of intermolecular voids in unsaturated phosphatidylcholine lipid bilayers

Computer simulation of the liquid crystalline phase of five different hydrated unsaturated phosphadidylcholine (PC) lipid bilayers, i.e., membranes built up by 18:0/18:1omega9cis PC, 18:0/18:2omega6cis PC, 18:0/18:3omega3cis PC, 18:0/20:4omega6cis PC, and 18:0/22:6omega3cis PC molecules have been performed on the isothermal-isobaric ensemble at 1 atm and 303 K. (The notation n:domegapcis specifies the lipid tails: n refers to the total number of carbon atoms in the chain, d is the number of the methylene-interrupted double bonds, p denotes the number of carbons between the chain terminal CH(3) group and the nearest double bond, and cis refers to the conformation around the double bonds.) The characteristics of the free volume in these systems have been analyzed by means of a generalized version of the Voronoi-Delaunay method [M. G. Alinchenko et al., J. Phys. Chem. B 108, 19056 (2004)]. As a reference system, the hydrated bilayer of the saturated 14:014:0 PC molecules (dimyristoylphosphatidylcholine) has also been analyzed. It has been found that the profiles of the fraction of the free volume across the membrane exhibit a rather complex pattern. This fine structure of the free volume fraction profiles can be interpreted by dividing the membrane into three separate major zones (i.e., zones of the aqueous, polar, and apolar parts of the membrane) and defining five subzones within these zones according to the average position of various atomic groups in the membrane. The fraction of the free volume in the middle of the membrane is found to increase with increasing unsaturation of the sn-2 chain of the lipid molecule. This is due to the fact that with increasing number of methylene-interrupted double bonds the lipid tails become more flexible, and hence they do not extend to the middle of the membrane. It is found that there are no broad enough preformed channels in the bilayers through which small penetrants, such as water molecules, can readily go through; however, the existing channels can largely facilitate the permeation of these molecules.