Rm. Venable et al., Molecular dynamics simulations of gel (L-beta I) phase lipid bilayers in constant pressure and constant surface area ensembles, J CHEM PHYS, 112(10), 2000, pp. 4822-4832
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].