Biophys J, November 1999, p. 2462-2469, Vol. 77, No. 5

Combined Monte Carlo and Molecular Dynamics Simulation of Fully Hydrated Dioleyl and Palmitoyl-oleyl Phosphatidylcholine Lipid Bilayers

 ^*Department of Molecular and Integrative Physiology, Department of Biochemistry, UIUC Programs in Biophysics, Neuroscience, and Bioengineering, and Beckman Institute, University of Illinois, Urbana, Illinois 61801; and  ^#Department of Physics, Oklahoma State University, Stillwater, Oklahoma 74078 USA

We have applied a new equilibration procedure for the atomic level simulation of a hydrated lipid bilayer to hydrated bilayers of dioleyl-phosphatidylcholine (DOPC) and palmitoyl-oleyl phosphatidylcholine (POPC). The procedure consists of alternating molecular dynamics trajectory calculations in a constant surface tension and temperature ensemble with configurational bias Monte Carlo moves to different regions of the configuration space of the bilayer in a constant volume and temperature ensemble. The procedure is applied to bilayers of 128 molecules of POPC with 4628 water molecules, and 128 molecules of DOPC with 4825 water molecules. Progress toward equilibration is almost three times as fast in central processing unit (CPU) time compared with a purely molecular dynamics (MD) simulation. Equilibration is complete, as judged by the lack of energy drift in 200-ps runs of continuous MD. After the equilibrium state was reached, as determined by agreement between the simulation volume per lipid molecule with experiment, continuous MD was run in an ensemble in which the lateral area was restrained to fluctuate about a mean value and a pressure of 1 atm applied normal to the bilayer surface. Three separate continuous MD runs, 200 ps in duration each, separated by 10,000 CBMC steps, were carried out for each system. Properties of the systems were calculated and averaged over the three separate runs. Results of the simulations are presented and compared with experimental data and with other recent simulations of POPC and DOPC. Analysis of the hydration environment in the headgroups supports a mechanism by which unsaturation contributes to reduced transition temperatures. In this view, the relatively horizontal orientation of the unsaturated bond increases the area per lipid, resulting in increased water penetration between the headgroups. As a result the headgroup-headgroup interactions are attenuated and shielded, and this contributes to the lowered transition temperature.

Biophys J, November 1999, p. 2462-2469, Vol. 77, No. 5
© 1999 by the Biophysical Society   0006-3495/99/11/2462/08  $2.00

This article has been cited by other articles:

				N. Kucerka, J. D. Perlmutter, J. Pan, S. Tristram-Nagle, J. Katsaras, and J. N. Sachs The Effect of Cholesterol on Short- and Long-Chain Monounsaturated Lipid Bilayers as Determined by Molecular Dynamics Simulations and X-Ray Scattering Biophys. J., September 15, 2008; 95(6): 2792 - 2805. [Abstract] [Full Text] [PDF]

				A. Dickey and R. Faller Examining the Contributions of Lipid Shape and Headgroup Charge on Bilayer Behavior Biophys. J., September 15, 2008; 95(6): 2636 - 2646. [Abstract] [Full Text] [PDF]

				R. B. M. Koehorst, R. B. Spruijt, and M. A. Hemminga Site-Directed Fluorescence Labeling of a Membrane Protein with BADAN: Probing Protein Topology and Local Environment Biophys. J., May 15, 2008; 94(10): 3945 - 3955. [Abstract] [Full Text] [PDF]

				L. C. Silva, R. F. M. de Almeida, B. M. Castro, A. Fedorov, and M. Prieto Ceramide-Domain Formation and Collapse in Lipid Rafts: Membrane Reorganization by an Apoptotic Lipid Biophys. J., January 15, 2007; 92(2): 502 - 516. [Abstract] [Full Text] [PDF]

				M. Patra, E. Salonen, E. Terama, I. Vattulainen, R. Faller, B. W. Lee, J. Holopainen, and M. Karttunen Under the Influence of Alcohol: The Effect of Ethanol and Methanol on Lipid Bilayers Biophys. J., February 15, 2006; 90(4): 1121 - 1135. [Abstract] [Full Text] [PDF]

				R. W. Benz, F. Castro-Roman, D. J. Tobias, and S. H. White Experimental Validation of Molecular Dynamics Simulations of Lipid Bilayers: A New Approach Biophys. J., February 1, 2005; 88(2): 805 - 817. [Abstract] [Full Text] [PDF]

				R. F. M. de Almeida, L. M. S. Loura, M. Prieto, A. Watts, A. Fedorov, and F. J. Barrantes Cholesterol Modulates the Organization of the {gamma}M4 Transmembrane Domain of the Muscle Nicotinic Acetylcholine Receptor Biophys. J., April 1, 2004; 86(4): 2261 - 2272. [Abstract] [Full Text] [PDF]

				T. Rog, K. Murzyn, R. Gurbiel, Y. Takaoka, A. Kusumi, and M. Pasenkiewicz-Gierula Effects of phospholipid unsaturation on the bilayer nonpolar region: a molecular simulation study J. Lipid Res., February 1, 2004; 45(2): 326 - 336. [Abstract] [Full Text] [PDF]

				A. Filippov, G. Oradd, and G. Lindblom The Effect of Cholesterol on the Lateral Diffusion of Phospholipids in Oriented Bilayers Biophys. J., May 1, 2003; 84(5): 3079 - 3086. [Abstract] [Full Text] [PDF]

				S. W. Chiu, E. Jakobsson, R. J. Mashl, and H. L. Scott Cholesterol-Induced Modifications in Lipid Bilayers: A Simulation Study Biophys. J., October 1, 2002; 83(4): 1842 - 1853. [Abstract] [Full Text] [PDF]