Self-consistent-field modeling of hydrated unsaturated lipid bilayers in the liquid-crystal phase and comparison to...

A molecular-level self-consistent-field (SCF) theory is applied to model the lipid bilayer structures composed of 1-stearoyl-2-oleoyl-sn-glycero-3-phosphatidylcholine (18:0/18:1 omega 9cis PC) and 1-stearoyl-2-docosahexaenoyl-sn-glycero-3-phosphatidylcholine (18:0/22:6 omega 3cis PC). As compared to earlier attempts to model (saturated) PC membranes several additional features are implemented: (i) A water model is used which correctly leads to low water concentration in the bilayers. (ii) Free volume is allowed for, which is important to obtain bilayers in the fluid state. (iii) A polarization term is included in the segment potentials; this new feature corrects for a minor thermodynamic inconsistency present in (all) earlier results for charged bilayers. (iv) The CH3 groups in the lipid molecules are assumed to have twice the volume of a CH2 group; this leads to stable noninterdigitated bilayers. (v) A cis double bond is simulated by forcing gauche conformations along the sn-2 acyl chain. Results of an all-atom molecular dynamics (MD) simulation, using the collision dynamics method, on the same system are presented. Both SCF and MD prove, in accordance with experimental facts, that acyl unsaturation effectively reduces the length of the chain which counteracts interdigitation. It is also found that the phosphatidylcholine head group is lying almost flat on the membrane surface and the water penetrates into the bilayer upto the glycerol backbone units. From the SCF results it further followed that the free volume is not exactly evenly distributed over the bilayer. There is a small increase in free volume in the center of the bilayer as well as in the glycerol backbone region.