Kechuan Tu

Constant pressure and temperature molecular dynamics simulation of a fully hydrated liquid crystal phase dipalmitoylphosphatidylcholine bilayer.

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.

Atomic-scale molecular dynamics simulations of lipid membranes

In the past few years there have been dozens of reports of atomic-scale molecular dynamics simulations of lipid membranes. Regarding recent methodological developments, it is now generally accepted that there are advantages to doing membrane simulations at constant pressure, but there is also a debate on whether the external pressure should be isotropic or whether tension should be applied in the plane of the bilayer. Also, it has become clear that serious artifacts can be introduced into computer simulations if all electrostatic interactions are not taken into account. Simulations can now faithfully reproduce many known aspects of the structures of gel and liquid crystal phase bilayers, and simulators have begun to include cholesterol and proteins in membrane simulations. There is some disagreement between simulations on the details of the bilayer/water interface, and there has not been enough effort to analyze the dynamics afforded by simulations.

Molecular dynamics investigation of the structure of a fully hydrated gel-phase dipalmitoylphosphatidylcholine bilayer

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.

Analysis of Simulated NMR Order Parameters for Lipid Bilayer Structure Determination

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.

Effects of Anesthetics on the Structure of a Phospholipid Bilayer: Molecular Dynamics Investigation of Halothane in the Hydrated Liquid Crystal Phase of Dipalmitoylphosphatidylcholine

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).

Molecular Dynamics Simulations of Supported Phospholipid/Alkanethiol Bilayers on a Gold(111) Surface

Molecular dynamics simulations have been used to investigate the structure of hybrid bilayers (HB) formed by dipalmitoylphosphatidylcholine (DPPC) lipid monolayers adsorbed on a hydrophobic alkanethiol self-assembled monolayer (SAM). The HB system was studied at 20 degrees C and 60 degrees C, and the results were compared with recent neutron reflectivity measurements (Meuse, C. W., S. Krueger, C. F. Majkrzak, J. A. Dura, J. Fu, J. T. Connor, and A. L. Plant. 1998. Biophys. J. 74:1388) and previous simulations of hydrated multilamellar bilayers (MLB) of DPPC (Tu, K., D. J. Tobias, and M. L. Klein. 1995. Biophys. J. 69:2558; and 1996. 70:595). The overall structures of the HBs are in very good agreement with experiment. The structure of the SAM monolayer is hardly perturbed by the presence of the DPPC overlayer. The DPPC layer presents characteristics very similar to the MLB gel phase at low temperature and to the liquid crystal phase at high temperature. Subtle changes have been found for the lipid/water interface of the HBs compared to the MLBs. The average phosphatidylcholine headgroup orientation is less disordered, and this produces changes in the electric properties of the HB lipid/water interface. These changes are attributed to the fact that the aqueous environment of the lipids in these unilamellar films is different from that of MLB stacks. Finally, examination of the intramolecular and whole-molecule dynamics of the DPPC molecules in the fluid phase HB and MLB membranes revealed that the reorientations of the upper part of the acyl chains (near the acyl ester linkage) are slower, the single molecule protrusions are slightly damped, and the lateral rattling motions are significantly reduced in the HB compared with the MLB.