We report a constant pressure and temperature molecular dynamics simulation of a fully hydrated liquid crystal (L alpha) phase bilayer of dipalmitoylphosphatidylcholine at 50 degrees C and 28 water molecules/lipid. We have shown that the bilayer is stable throughout the 1550-ps simulation and have demonstrated convergence of the system dimensions. Several important aspects of the bilayer structure have been investigated and compared favorably with experimental results. For example, the average positions of specific carbon atoms along the bilayer normal agree well with neutron diffraction data, and the electron density profile is in accord with x-ray diffraction results. The hydrocarbon chain deuterium order parameters agree reasonably well with NMR results for the middles of the chains, but the simulation predicts too much order at the chain ends. In spite of the deviations in the order parameters, the hydrocarbon chain packing density appears to be essentially correct, inasmuch as the area/lipid and bilayer thickness are in agreement with the most refined experimental estimates. The deuterium order parameters for the glycerol and choline groups, as well as the phosphorus chemical shift anisotropy, are in qualitative agreement with those extracted from NMR measurements.

We report a 1.4-ns constant-pressure molecular dynamics simulation of cholesterol at 12.5 mol% in a dipalmitoylphosphatidylcholine (DPPC) bilayer at 50 degrees C and compare the results to our previous simulation of a pure DPPC bilayer. The interlamellar spacing was increased by 2.5 A in the cholesterol-containing bilayer, consistent with x-ray diffraction results, whereas the bilayer thickness was increased by only 1 A. The bilayer/water interface was more abrupt because the lipid headgroups lie flatter to fill spaces left by the cholesterol molecules. This leads to less compensation by the lipid headgroups of the oriented water contribution to the membrane dipole potential and could explain the experimentally observed increase in the magnitude of the dipole potential by cholesterol. Our calculations suggested that 12.5 mol% cholesterol does not significantly affect the conformations and packing of the hydrocarbon chains and produces only a slight reduction in the empty free volume. However, cholesterol has a significant influence on the subnanosecond time scale lipid dynamics: the diffusion constant for the center-of-mass "rattling" motion was reduced by a factor of 3, and the reorientational motion of the methylene groups was slowed along the entire length of the hydrocarbon chains.

We report the results of a constant pressure and temperature molecular dynamics simulation of a gel-phase dipalmitoylphosphatidylcholine bilayer with nw = 11.8 water molecules/lipid at 19 degrees C. The results of the simulation were compared in detail with a variety of x-ray and neutron diffraction data. The average positions of specific carbon atoms along the bilayer normal and the interlamellar spacing and electron density profile were in very good agreement with neutron and x-ray diffraction results. The area per lipid and the details of the in-plane hydrocarbon chain structure were in excellent agreement with wide-angle x-ray diffraction results. The only significant deviation is that the chains met in a pleated arrangement at the bilayer center, although they should be parallel. Novel discoveries made in the present work include the observation of a bimodal headgroup orientational distribution. Furthermore, we found that there are a significant number of gauche conformations near the ends of the hydrocarbon chains and, in addition to verifying a previous suggestion that there is partial rotational ordering in the hydrocarbon chains, that the two chains in a given molecule are inequivalent with respect to rotations. Finally, we have investigated the lipid/water interface and found that the water penetrates beneath the headgroups, but not as far as the carbonyl groups, that the phosphates are strongly hydrated almost exclusively at the nonesterified oxygen atoms, and that the hydration of the ammonium groups is more diffuse, with some water molecules concentrated in the grooves between the methyl groups.

The conventional formula for relating CD2 average order parameters to average methylenic travel is flawed when compared to molecular dynamics simulations of dipalmitoylphosphatidylcholine. Inspired by the simulated probability distribution functions, a new formula is derived that satisfactorily relates these quantities. This formula is used to obtain the average chain length , and the result agrees with the direct simulation result for . The simulation also yields a hydrocarbon thickness 2. The result = is consistent with a model of chain packing with both early chain termination and partial interdigitation of chains from opposing monolayers. The actual simulated area per lipid is easily obtained from the order parameters. However, when this method is applied to NMR order parameter data from dimyristoylphosphatidylcholine, the resulting is 10% larger than the currently accepted value.

We report the results of constant temperature and pressure molecular dynamics calculations carried out on the liquid crystal (Lalpha) phase of dipalmitoylphosphatidylcholine with a mole fraction of 6.5% halothane (2-3 MAC). The present results are compared with previous simulations for pure dipalmitoylphosphatidylcholine under the same conditions (Tu et al., 1995. Biophys. J. 69:2558-2562) and with various experimental data. We have found subtle structural changes in the lipid bilayer in the presence of the anesthetic compared with the pure lipid bilayer: a small lateral expansion is accompanied by a modest contraction in the bilayer thickness. However, the overall increase in the system volume is found to be comparable to the molecular volume of the added anesthetic molecules. No significant change in the hydrocarbon chain conformations is apparent. The observed structural changes are in fair agreement with NMR data corresponding to low anesthetic concentrations. We have found that halothane exhibits no specific binding to the lipid headgroup or to the acyl chains. No evidence is obtained for preferential orientation of halothane molecules with respect to the lipid/water interface. The overall dynamics of the lipid-bound halothane molecules appears to be reminiscent of that of other small solutes (Bassolino-Klimas et al., 1995. J. Am. Chem. Soc. 117:4118-4129).