Molecular dynamics simulations of gel (LβI) phase lipid bilayers in constant pressure and constant surface area ensembles
Molecular dynamics simulations of gel (LβI) phase lipid bilayers in constant pressure and constant surface area ensembles
复制标题
DOI:
10.1063/1.481085
复制
发表时间:
2000-02
影响因子:
4.4
通讯作者:
R. Venable;B. Brooks;R. Pastor
中科院分区:
文献类型:
--
作者:
R. Venable;B. Brooks;R. Pastor
The results of a series of molecular dynamics simulations of the gel state of a dipalmitoylphosphatidylcholine bilayer at 293 K are described. The simulations, ranging from 40 ps to 2.5 ns, show clearly that: a flexible cell geometry is essential during equilibration; Ewald summation of electrostatics is superior to spherical cutoff methods; water exchange with the carbonyl 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, though the surface area per lipid is slightly underestimated. The randomization 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 w...