Molecular dynamics simulations of gel (L-beta I) phase lipid bilayers in constant pressure and constant surface area ensembles

The results of a series of molecular dynamics simulations of the gel state of a dipalmitoylphosphatidylcholine bilayer at 293 K are described. The sim ulations, ranging from 40 ps to 2.5 ns, show clearly that: a flexible cell geometry is essential during equilibration; Ewald summation of electrostati cs is superior to spherical cutoff methods; water exchange with the carbony l group of chain 2 takes place on the ns time scale, while there is almost no hydration of chain 1. There is overall good agreement (D-spacing, chain tilt, fraction gauche, and area compressibility modulus) with experiment, t hough the surface area per lipid is slightly underestimated. The randomizat ion of torsion 1 of chain 2 from exclusively gauche minus (as specified in the initial condition modeled from the crystal structure of a related lipid ) to a mixture of g+/g- over the course of approximately 2 ns is a critical feature of the study. The torsional equilibration proceeded steadily when simulating at constant surface tension, but was effectively quenched by sim ulation at constant area. The associated presence of conformational degener acy of this torsion, and conformational disorder in the upper region of cha in 2, is most likely associated with the seemingly anomalous infrared (IR) results for gauche bonds in the upper region of the chains. It may also be a characteristic of the gel phase, and be related to the long time required for the gel to subgel transition. [S0021-9606(00)50210-4].