Multibilayer structure of dimyristoylphosphatidylcholine dihydrate as determined by energy minimization

Complete energy minimization was carried out on the multibilayer crystal structure of 1,2-dimyristoyl-sn-glycero-3-phosphocholine dihydrate (DMPC.2H2O), starting from the X-ray structure determination reported by Pearson and Pascher (1979) Nature 281, 499-501. The asymmetric unit contains two nonidentical DMPC molecules and four water molecules. Minimization removed the acyl chain disorder present in the X-ray structure and caused the carbon planes of the acyl chains to become mutually parallel. Two energy-minimized structures (structures I and II) were found which mainly differed in the hydrogen-bonding arrangement of the waters of hydration. In structure I as in the X-ray structure, one of the water molecules forms a hydrogen-bonded bridge between successive bilayers; but in structure II, all hydrogen bonds are satisfied on the same bilayer. Structure II corresponds to the global energy minimum and is also a suitable structure for single bilayers. The lattice constants and cell volume of the minimized structures are close to the experimental values. The electrostatic force between DMPC bilayers is attractive. The mean hydration energy of the water is -14.2 kcal/mol, which is 2.5 kcal/mol lower than the binding energy of ice.