In this paper we report on the molecular dynamics simulation of a fluid phase hydrated dimyristoylphosphatidylcholine bilayer. The initial configuration of the lipid was the x-ray crystal structure. A distinctive feature of this simulation is that, upon heating the system, the fluid phase emerged from parameters, initial conditions, and boundary conditions determined independently of the collective properties of the fluid phase. The initial conditions did not include chain disorder characteristic of the fluid phase. The partial charges on the lipids were determined by ab initio self-consistent field calculations and required no adjustment to produce a fluid phase. The boundary conditions were constant pressure and temperature. Thus the membrane was not explicitly required to assume an area/phospholipid molecule thought to be characteristic of the fluid phase, as is the case in constant volume simulations. Normal to the membrane plane, the pressure was 1 atmosphere, corresponding to the normal laboratory situation. Parallel to the membrane plane a negative pressure of -100 atmospheres was applied, derived from the measured surface tension of a monolayer at an air-water interface. The measured features of the computed membrane are generally in close agreement with experiment. Our results confirm the concept that, for appropriately matched temperature and surface pressure, a monolayer is a close approximation to one-half of a bilayer. Our results suggest that the surface area per phospholipid molecule for fluid phosphatidylcholine bilayer membranes is smaller than has generally been assumed in computational studies at constant volume. Our results confirm that the basis of the measured dipole potential is primarily water orientations and also suggest the presence of potential barriers for the movement of positive charges across the water-headgroup interfacial region of the phospholipid.

We introduce a new force field (43A1-S3) for simulation of membranes by the Gromacs simulation package. Construction of the force fields is by standard methods of electronic structure computations for bond parameters and charge distribution and specific volumes and heats of vaporization for small-molecule components of the larger lipid molecules for van der Waals parameters. Some parameters from the earlier 43A1 force field are found to be correct in the context of these calculations, while others are modified. The validity of the force fields is demonstrated by correct replication of X-ray form factors and NMR order parameters over a wide range of membrane compositions in semi-isotropic NTP 1 atm simulations. 43-A1-S3 compares favorably with other force fields used in conjunction with the Gromacs simulation package with respect to the breadth of phenomena that it accurately reproduces.

We present analysis of new configurational bias Monte Carlo and molecular dynamics simulation data for bilayers of dipalmitoyl phosphatidyl choline and cholesterol for dipalmitoyl phosphatidyl choline:cholesterol ratios of 24:1, 47:3, 11.5:1, 8:1, 7:1, 4:1, 3:1, 2:1, and 1:1, using long molecular dynamics runs and interspersed configurational bias Monte Carlo to ensure equilibration and enhance sampling. In all cases with cholesterol concentrations above 12.5% the area per molecule of the heterogeneous membrane varied linearly with cholesterol fraction. By extrapolation to pure cholesterol, we find the cross-sectional area of cholesterol in these mixtures is approximately 22.3 A(2). From the slope of the area/molecule relationship, we also find that the phospholipid in these mixtures is in a liquid ordered state with an average cross-sectional area per lipid of 50.7 A(2), slightly above the molecular area of a pure phospholipid gel. For lower concentrations of cholesterol, the molecular area rises above the straight line, indicating the "melting" of at least some of the phospholipid into a fluid state. Analysis of the lateral distribution of cholesterol molecules in the leaflets reveals peaks in radial distributions of cholesterols at multiples of approximately 5 A. These peaks grow in size as the simulation progresses, suggesting a tendency for small subunits of one lipid plus one cholesterol, hydrogen bonded together, to act as one composite particle, and perhaps to aggregate with other composites. Our results are consistent with experimentally observed effects of cholesterol, including the condensation effect of cholesterol in phospholipid monolayers and the tendency of cholesterol-rich domains to form in cholesterol-lipid bilayers. We are continuing to analyze this tendency on longer timescales and for larger bilayer patches.

We have carried out an atomic-level molecular dynamics simulation of a system of nanoscopic size containing a domain of 18:0 sphingomyelin and cholesterol embedded in a fully hydrated dioleylposphatidylcholine (DOPC) bilayer. To analyze the interaction between the domain and the surrounding phospholipid, we calculate order parameters and area per molecule as a function of molecule type and proximity to the domain. We propose an algorithm based on Voronoi tessellation for the calculation of the area per molecule of various constituents in this ternary mixture. The calculated areas per sphingomyelin and cholesterol are in agreement with previous simulations. The simulation reveals that the presence of the liquid-ordered domain changes the packing properties of DOPC bilayer at a distance as large as approximately 8 nm. We calculate electron density profiles and also calculate the difference in the thickness between the domain and the surrounding DOPC bilayer. The calculated difference in thickness is consistent with data obtained in atomic force microscopy experiments.

We have carried out a molecular dynamics simulation of a hydrated 18:0 sphingomyelin lipid bilayer. The bilayer contained 1600 sphingomyelin (SM) molecules, and 50,592 water molecules. After construction and initial equilibration, the simulation was run for 3.8 ns at a constant temperature of 50 degrees C and a constant pressure of 1 atm. We present properties of the bilayer calculated from the simulation, and compare with experimental data and with properties of dipalmitoyl phosphatidylcholine (DPPC) bilayers. The SM bilayers are significantly more ordered and compact than DPPC bilayers at the same temperature. SM bilayers also exhibit significant intramolecular hydrogen bonding between phosphate ester oxygen and hydroxyl hydrogen atoms. This results in a decreased hydration in the polar region of the SM bilayer compared with DPPC. Since our simulation system is very large we have calculated the power spectrum of bilayer undulation and peristaltic modes, and we compare these data with similar calculations for DPPC bilayers. We find that the SM bilayer has significantly larger bending modulus and area compressibility compared to DPPC.

The fundamental role of cholesterol in the regulation of eukaryotic membrane structure is wellestablished. However the manner in which atomic level interactions between cholesterol and lipids, with varying degrees of chain unsaturation and polar groups, affect the overall structure and organization of the bilayer is only beginning to be understood. In this paper we describe a series of Molecular Dynamics simulations designed to provide new insights into lipid-cholesterol interactions as a function of chain unsaturation. We have run simulations of varying concentrations of cholesterol in dipalmitoyl phosphatidylcholine (DPPC), palmitoyl-oleyol phosphatidylcholine (POPC), and dioleyol phosphatidylcholine (DOPC) bilayers. Structural analysis of the simulations reveals both atomistic and systemic details of the interactions and are presented here. In particular, we find that the minimum partial molecular area of cholesterol occurs in POPC-Chol mixtures implying the most favorable packing. Physically, this appears to be related to the fact that the two faces of the cholesterol molecule are different from each other and that the steric cross section of cholesterol molecules drops sharply near the small chain tails.

Aβ oligomers are potential targets for the diagnosis and therapy of Alzheimer's disease (AD). On the other hand, the molecule curcumin has been shown to possess significant therapeutic potential in many areas. In this paper, we use all-atom explicit solvent molecular dynamics simulations to study the effect of curcumin on the stability of Aβ amyloid protein oligomers. We observed that curcumin decreases the β-sheet secondary structural content within the Aβ oligomers without reducing the contacts between the monomers. The breaking of the β-sheet is found to be preceded by a deformation of the β-sheet structure due to hydrophobic interaction from the nearby curcumin. Furthermore, the π-stacking interaction between curcumin (keto ring and enol ring) and the aromatic residues of Aβ, which exists throughout the simulations, has also contributed to the diminishing of the β-sheet structure. Our analysis of the underwrapped amide-carbonyl hydrogen bonds reveals several stable dehydrons of the oligomer, especially the dehydron pair 34L and 41I, which curcumin tends to hover over. We have examined the paths of curcumin on the Aβ proteins and determined the common routes where curcumin lingers as it traverses around the Aβ. In consequence, our study has provided a detailed interaction picture between curcumin and the Aβ oligomers.

This paper reports on a simulation of a gramicidin channel inserted into a fluid phase DMPC bilayer with 100 lipid molecules. Two lipid molecules per leaflet were removed to insert the gramicidin, so the resulting preparation had 96 lipid molecules and 3209 water molecules. Constant surface tension boundary conditions were employed. Like previous simulations with a lower lipid/gramicidin ratio (Woolf, T. B., and B. Roux. 1996. Proteins: Struct., Funct., Genet. 24:92-114), it is found that tryptophan-water hydrogen bonds are more common than tryptophan-phospholipid hydrogen bonds. However, one of the tryptophan NH groups entered into an unusually long-lived hydrogen bonding pattern with two glycerol oxygens of one of the phospholipid molecules. Comparisons were made between the behavior of the lipids adjacent to the channel with those farther away. It was found that hydrocarbon chains of lipids adjacent to the channel had higher-order parameters than those farther away. The thickness of the lipid bilayer immediately adjacent to the channel was greater than it was farther away. In general, the lipids adjacent to the membrane had similar orientations to those seen by Woolf and Roux, while those farther away had similar orientations to those pertaining before the insertion of the gramicidin. A corollary to this observation is that the thickness of the hydrocarbon region adjacent to the gramicidin was much thicker than what other studies have identified as the "hydrophobic length" of the gramicidin channel.