Application of combined Monte Carlo and molecular dynamics method to simulation of dipalmitoyl phosphatidylcholine lipid bilayer

We describe a new equilibration procedure for the atomic level simulation of a hydrated lipid bilayer. 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 described in detail and is applied to a bilayer of 100 molecules of dipalmitoyl phosphatidylcholine (DPPC) and 3205 water molecules. We find that the hybrid simulation procedure enhances the equilibration of the bilayer as measured by the convergence of the area per molecule and the segmental order parameters, as compared with a simulation using only molecular dynamics (MD). Progress toward equilibration is almost three times as fast in CPU time, compared with a purely MD simulation. Equilibration is complete, as judged by the lack of energy drift in three separate 200-ps runs of continuous MD started from different initial states. Results of the simulation are presented and compared with experimental data and with other recent simulations of DPPC. ©1999 John Wiley & Sons, Inc. J Comput Chem 20: 1153–1164, 1999     

Membrane orientation of MSI-78 measured by sum frequency generation vibrational spectroscopy

Antimicrobial peptides (AMPs) selectively disrupt bacterial cell membranes to kill bacteria whereas they either do not or weakly interact with mammalian cells. The orientations of AMPs in lipid bilayers mimicking bacterial and mammalian cell membranes are related to their antimicrobial activity and selectivity. To understand the role of AMP-lipid interactions in the functional properties of AMPs better, we determined the membrane orientation of an AMP (MSI-78 or pexiganan) in various model membranes using sum frequency generation (SFG) vibrational spectroscopy. A solid-supported single 1,2-dipalmitoyl-an-glycero-3-[phospho-rac-(1-glycerol)] (DPPG) bilayer or 1-palmitoyl-2-oleoyl-sn-glycero-3-[phospho-rac-(1-glycerol)] (POPG) bilayer was used as a model bacterial cell membrane. A supported 1,2-dipalmitoyl-an-glycero-3-phosphocholine (DPPC) bilayer or a 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) bilayer was used as a model mammalian cell membrane. Our SFG results indicate that the helical MSI-78 molecules are associated with the bilayer surface with ～70° deviation from the bilayer normal in the negatively charged gel-phase DPPG bilayer at 400 nM peptide concentration. However, when the concentration was increased to 600 nM, MSI-78 molecules changed their orientation to make a 25° tilt from the lipid bilayer normal whereas multiple orientations were observed for an even higher peptide concentration in agreement with toroidal-type pore formation as reported in a previous solid-state NMR study. In contrary, no interaction between MSI-78 and a zwitterionic DPPC bilayer was observed even at a much higher peptide concentration (～12,000 nM). These results demonstrate that SFG can provide insights into the antibacterial activity and selectivity of MSI-78. Interestingly, the peptide exhibits a concentration-dependent membrane orientation in the lamellar-phase POPG bilayer and was also found to induce toroidal-type pore formation. The deduced lipid flip-flop from SFG signals observed from lipids also supports MSI-78-induced toroidal-type pore formation.     

Monte Carlo study of lipid membranes: Simulation of diparmitoylphosphatidylcholine bilayers in gel and liquid-crystalline phases

The statistical properties of the bilayer membranes of diparmitoylphosphatidylcholine (DPPC) in the gel and liquid-crystal phases were studied by Monte Carlo (MC) simulation using potential functions of the Lennard-Jones, the simple Coulomb, and the bond torsion. The simulation was undertaken on a two-dimensional periodic condition imposed on the bilayer model consisting of faithfully described molecules. The structure and ordering of the model bilayers accorded well with experiments, and the segment order parameters were in agreement with those of the nuclear magnetic resonance (NMR) experiments. The two kinds of lipid chains of DPPC do not equivalently behave in the bilayers, and chain 2 has lower ordering than chain 1. The order parameters of the first eight segments of chain 2 in the liquid-crystal model are particularly small and are roughly constant. From electron density analysis, it has been observed that the liquid-crystal bilayer has about one excess water molecule per one lipid molecule in comparison with the gel bilayer. The energy difference between the two bilayer models, taking account of the water contribution, is consistent with the latent heat of the phase transition. © 1995 by John Wiley & Sons, Inc.     

Coarse-grained model for phospholipid/cholesterol bilayer

We construct a coarse-grained (CG) model for dipalmitoylphosphatidylcholine (DPPC)/cholesterol bilayers and apply it to large-scale simulation studies of lipid membranes. Our CG model is a two-dimensional representation of the membrane, where the individual lipid and sterol molecules are described by pointlike particles. The effective intermolecular interactions used in the model are systematically derived from detailed atomic-scale molecular dynamics simulations using the Inverse Monte Carlo technique, which guarantees that the radial distribution properties of the CG model are consistent with those given by the corresponding atomistic system. We find that the coarse-grained model for the DPPC/cholesterol bilayer is substantially more efficient than atomistic models, providing a speedup of approximately eight orders of magnitude. The results are in favor of formation of cholesterol-rich and cholesterol-poor domains at intermediate cholesterol concentrations, in agreement with the experimental phase diagram of the system. We also explore the limits of the coarse-grained model, and discuss the general validity and applicability of the present approach.     

Phase Transitions and Structural Changes in DPPC Liposomes Induced by a 1-Carba-Alpha-Tocopherol Analogue

Steady-state emission spectroscopy of 1-anilino-8- naphthalene sulfonate (ANS) and 1,6-diphenyl-1,3,5-hexatriene (DPH), fluorescence anisotropy, and DSC methods were used to characterize the interactions of the newly synthesized 1-carba-alpha-tocopherol (CT) with a 1,2-dipalmitoyl-sn-glycero-3-phosphocholine (DPPC) membrane. The DSC results showed significant perturbations in the DPPC structure for CT concentrations as low as 2 mol%. The main phase transition peak was broadened and shifted to lower temperatures in a concentration-dependent manner, and pretransition was abolished. Increasing CT concentrations induced the formation of new phases in the DPPC structure, leading to melting at lower temperatures and, finally, disruption of the ordered DPPC structure. Hydration and structural changes of the DPPC liposomes using ANS and DPH fluorescent probes, which are selectively located at different places in the bilayer, were studied. With the increased concentration of CT molecules in the DPPC liposomes, structural changes with the simultaneous formation of different phases of such mixture were observed. Temperature studies of such mixtures revealed a decrease in the temperature of the main phase transition and fluidization at decreasing temperatures related to increasing hydration in the bilayer. Contour plots obtained from concentration–temperature data with fluorescent probes allowed for identification of different phases, such as gel, ordered liquid, disordered liquid, and liquid crystalline phases. The CT molecule with a modified chromanol ring embedded in the bilayer led to H-bonding interactions, expelling water molecules from the interphase, thus introducing disorder and structural changes to the highly ordered gel phase.     

Molecular dynamics simulations of DiI-C18(3) in a DPPC lipid bilayer

We performed a 40 ns simulation of 1,1'-dioctadecyl-3,3,3',3'-tetramethylindocarbocyanine perchlorate (DiI-C18(3)) in a 1,2-dipalmitoyl-sn-glycero-3-phosphatidyl choline (DPPC) bilayer in order to facilitate interpretation of lipid dynamics and membrane structure from fluorescence lifetime, anisotropy, and fluorescence correlations spectroscopy (FCS). Incorporation of DiI of 1.6 to 3.2 mol% induced negligible changes in area per lipid but detectable increases in bilayer thickness, each of which are indicators of membrane structural perturbation. The DiI chromophore angle was 77 +/- 17 degrees with respect to the bilayer normal, consistent with rotational diffusion inferred from polarization studies. The DiI headgroup was located 0.63 nm below the lipid head group-water interface, a novel result in contrast to some popular cartoon representations of DiI but consistent with DiI's increase in quantum yield when incorporated into lipid bilayers. Importantly, the fast component of rotational anisotropy matched published experimental results demonstrating that sufficient free volume exists at the sub-interfacial region to support fast rotations. Simulations with non-charged DiI head groups exhibited DiI flip-flop, demonstrating that the positively-charged chromophore stabilizes the orientation and location of DiI in a single monolayer. DiI induced detectable changes in interfacial properties of water ordering, electrostatic potential, and changes in P-N vector orientation of DPPC lipids. The diffusion coefficient of DiI (9.7 +/- 0.02 x 10(-8) cm2 s(-1)) was similar to the diffusion of DPPC molecules (10.7 +/- 0.04 x 10(-8) cm2 s(-1)), supporting the conclusion that DiI dynamics reflect lipid dynamics. These results provide the first atomistic level insight into DiI dynamics, results essential in elucidating lipid dynamics through single molecule fluorescence studies.     

Molecular dynamics simulation of the spontaneous formation of a small DPPC vesicle in water in atomistic detail

Molecular dynamics simulations have been used to study the spontaneous aggregation of a concentrated solution of dipalmitoylphosphatidylcholine (DPPC) molecules in water into a small vesicle. The molecules were represented in atomistic detail. Starting from a DPPC solution in water, an oblong vesicle with a long axis of 15 nm and short axes of 10 nm was formed spontaneously. After 90 ns of simulation, the vesicle contained a number of water pores. Water pores were shown to facilitate exchange of lipids between inner and outer leaflets. Lipid tails were shown to be less ordered in the inner leaflet of the vesicle, as compared to those in the outer leaflet of the vesicle and an equilibrated lamellar bilayer.     

Time-domain calculations of the infrared and polarized Raman spectra of tetraalanine in aqueous solution

The IR and polarized (isotropic and anisotropic) Raman spectra are calculated for the amide I band of tetraalanine ((Ala)4) in aqueous solution by using a time-domain computational method, which includes the effects of the diagonal frequency modulations (of individual peptide groups), the off-diagonal (interpeptide) vibrational couplings, and structural dynamics. It is shown that the calculated band profiles, especially the existence of a large negative noncoincidence effect (i.e., large frequency separations between the IR, isotropic Raman, and anisotropic Raman bands, with the isotropic Raman being higher in frequency), are in reasonable agreement with the experimental results. This negative noncoincidence effect derives from two conditions: the positive coupling between the amide I vibrations of peptide groups and the angle larger than 90 degrees between the transition dipoles of the coupled vibrations. This result means that the dynamically changing structures mainly in the polyproline II and beta-type conformations containing some repeated interconversions obtained from the molecular dynamics calculation are consistent with the existence of a large negative noncoincidence effect, as far as the structures satisfy the above two conditions. It is also shown that the electric fields from solvent water molecules induce larger frequency shifts than those of intrachain interactions, with rapid underdamped oscillatory modulations ( approximately 43 fs) due to the librational motions of water molecules that give rise to motional narrowing effect on the spectra. The reason for the difference from the behavior seen for the O-H stretching mode of liquid water is discussed. The time-domain analysis of the mode identity shows that the system proceeds halfway to complete mode mixing with a similar time scale ( approximately 60 fs), suggesting the importance of the nonadiabatic effect, which is included in a natural way in the present computational method.     

Gadolinium (III) ion in liquid water: structure, dynamics, and magnetic interactions from first principles

We applied first principles molecular dynamics (MD) technique to study structure, dynamics, and magnetic interactions of the Gd(3+) aqua ion dissolved in liquid water, a prototypical system for Gd-based complexes used as contrast agents for magnetic resonance imaging. The first coordination sphere contains eight water molecules with an average Gd-O distance of 2.37 A and an average geometric arrangement close to a square antiprism. The mean tilt angle of the electric dipole vector of these water molecules is theta=145 degrees . In our picosecond time scale simulation we observe no exchange event from the first coordination sphere but only fast "wagging" motions. The second coordination sphere is well pronounced though water molecules in this sphere are subjected to large amplitude dynamic motions. The isotropic hyperfine coupling constants for the inner sphere water molecules [A(iso)((17)O(I))=0.65+/-0.03 MHz, A(iso)((1)H(I))=0.085+/-0.005 MHz] are in good agreement with experimental data and with an earlier study using classical MD. Second sphere Fermi contact hyperfine coupling constants calculated are more than one order of magnitude smaller and of opposite sign as those of the first coordination sphere. The effect of spin polarization induced by the paramagnetic Gd(3+) ion on the dipolar hyperfine interaction was found to be sizable only for the (17)O nuclei of inner sphere water molecules and has a screening character.     

Effect of pressure on the Prodan fluorescence in bilayer membranes of phospholipids with varying acyl chain lengths

The fluorescence spectra of 6-propionyl-2-(dimethylamino)naphthalene (Prodan) were observed as a function of pressure for the bilayer membrane systems of dilauroylphosphatidylcholine (DLPC), dimyristoylphosphatidylcholine (DMPC), dipalmitoylphosphatidylcholine (DPPC), and distearoylphosphatidylcholine (DSPC). The wavelength of the emission maximum, lambdamax, was found to be 480, 430, and 500 nm for the liquid crystalline (Lalpha), ripple gel (P'beta), and pressure-induced interdigitated gel (LbetaI) phase, respectively. Since the lambdamax reflects the solvent property around the probe molecules, we could speculate on the location of the Prodan molecules in the bilayer membranes; in the Lalpha phase of the lipid bilayer, the Prodan molecules distribute around the phosphate of the lipids (i.e. the polar region). The Lalpha/P'beta phase transition caused the Prodan molecules to move into the less polar region near the glycerol backbone. The fluorescence intensity of the Prodan in the P'beta phase was dependent on the chain length of the lipids and on pressure; the shorter the chain length of the lipid, the stronger the fluorescence intensity of the Prodan. Moreover, for the DLPC bilayer membrane system, the fluorescence intensity at 430 nm increased with increasing pressure, indicating that the partition of Prodan into the DLPC bilayer membrane is promoted by applying pressure. In the case of the DPPC and DSPC bilayers, as the pressure increased further, the pressure-induced interdigitation caused the Prodan molecules to squeeze out of the glycerol backbone region and to move the hydrophilic region near the bilayer surface. The ratio of fluorescence intensity at 480 nm to that at 430 nm, F480/F430, showed a sharp change at the phase-transition pressure. In the case of the DPPC and DSPC bilayers, the values of F480/F430 showed an abrupt increase above a certain pressure higher than the Lalpha/P'beta transition pressure, which corresponds to the interdigitation from the P'beta to the LbetaI phase. The plot of F480/F430 versus pressure is available for recognition of the bilayer phase transitions, especially the bilayer interdigitation.     

Effect of different treatments of long-range interactions and sampling conditions in molecular dynamic simulations of rhodopsin embedded in a dipalmitoyl phosphatidylcholine bilayer

The present study analyzes the effect of the simulation conditions on the results of molecular dynamics simulations of G-protein coupled receptors (GPCRs) performed with an explicit lipid bilayer. Accordingly, the present work reports the analysis of different simulations of bovine rhodopsin embedded in a dipalmitoyl phosphatidylcholine (DPPC) lipid bilayer using two different sampling conditions and two different approaches for the treatment of long-range electrostatic interactions. Specifically, sampling was carried out either by using the statistical ensembles NVT or NPT (constant number of atoms, a pressure of 1 atm in all directions and fixed temperature), and the electrostatic interactions were treated either by using a twin-cutoff, or the particle mesh Ewald summation method (PME). The results of the present study suggest that the use of the NPT ensemble in combination with the PME method provide more realistic simulations. The use of NPT during the equilibration avoids the need of an a priori estimation of the box dimensions, giving the correct area per lipid. However, once the system is equilibrated, the simulations are irrespective of the sampling conditions used. The use of an electrostatic cutoff induces artifacts on both lipid thickness and the ion distribution, but has no direct effect on the protein and water molecules.     

The effect of pressure on the hydration structure around hydrophobic solute: a molecular dynamics simulation study

Molecular dynamics simulations are performed to study the effects of pressure on the hydrophobic interactions between neopentane molecules immersed in water. Simulations are carried out for five different pressure values ranging from 1 atm to 8000 atm. From potential of mean force calculations, we find that with enhancement of pressure, there is decrease in the well depth of contact minimum (CM) and the relative stability of solvent separated minimum over CM increases. Lower clustering of neopentane at high pressure is also observed in association constant and cluster-structure analysis. Selected site-site radial distribution functions suggest efficient packing of water molecules around neopentane molecules at elevated pressure. The orientational profile calculations of water molecules show that the orientation of water molecules in the vicinity of solute molecule is anisotropic and this distribution becomes flatter as we move away from the solute. Increasing pressure slightly changes the water distribution. Our hydrogen bond properties and dynamics calculations reveal pressure-induced formation of more and more number of water molecules with five and four hydrogen bond at the expense of breaking of two and three hydrogen bonded water molecules. We also find lowering of water-water continuous hydrogen bond lifetime on application of pressure. Implication of these results for relative dispersion of hydrophobic molecules at high pressure are discussed.     

Dissipative particle dynamics simulation study of the bilayer-vesicle transition

A bilayer structure is an important immediate for the vesicle formation. However,the mechanism for the bilayer-vesicle transition remains unclear. In this work,a dissipative particle dynamics(DPD) simulation method was employed to study the mechanism of the bilayer-vesicle transition. A coarse-grained model was built based on a lipid molecule termed dimyristoylphosphatidylcholine(DMPC). Simulations were performed from two different initial configurations:a random dispersed solution and a tensionless bilayer. It was found that the bilayer-vesicle transition was driven by the minimization of the water-tail hydrophobic interaction energy,and was accompanied with the increase of the position entropy due to the redistribution of water molecules. The bulk pressure was reduced during the bilayer-vesicle transition,suggesting the evolved vesicle morphology was at the relatively low free energy state. The membrane in the product vesicle was a two-dimensional fluid. It can be concluded that the membrane of a vesicle is not interdigitated and most of the bonds in lipid chains are inclined to orient along the radical axis of the vesicle.     

Effects of silver nanoparticles on the fluidity of bilayer in phospholipid liposome

This paper describes the formation and characterization of liposome entrapping the silver nanoparticles in bilayer. Silver nanoparticles were entrapped in the bilayer of dipalmitoylphosphatidylcholine (DPPC) liposome, named as silver-loaded liposome. Specifically, above the gel to liquid-crystalline phase transition temperature of this lipid (i.e., 41 degrees C), it was observed that membrane fluidities of silver-loaded liposomes were increased, and fluorescence anisotropy values were reduced from 0.114 to 0.097. This might be due to the structural modifications and interactions between DPPC molecules and silver nanoparticles within the bilayer. It was also confirmed that silver nanoparticles were entrapped in hydrophobic region of lipid bilayer with transmission electron microscopy (TEM) and electron energy loss spectroscopy (EELS) measurements.     

Molecular dynamics study on the stabilization of dehydrated lipid bilayers with glucose and trehalose

In this study, we use molecular dynamics simulations to investigate and compare the interactions of DPPC bilayers with and without saccharides (glucose or trehalose) under dehydrated conditions. Results from the simulations indicate that unilamellar bilayers lose their structural integrity under dehydrated conditions in the absence of saccharides; however, in the presence of either glucose or trehalose, the bilayers maintain their stability. Hydrogen bond analysis shows that the saccharide molecules displace a significant amount of water surrounding the lipid headgroups. At the same time, the additional hydrogen bonds formed between water and saccharide molecules help to maintain a hydration layer on the lipid bilayer interface. On the basis of the hydrogen bond distributions, trehalose forms more hydrogen bonds with the lipids than glucose, and it is less likely to interact with neighboring saccharide molecules. These results suggest that the interaction between the saccharide and lipid molecules through hydrogen bonds is an essential component of the mechanism for the stabilization of lipid bilayers.     

Quantum chemical investigation of hyperfine coupling constants on first coordination sphere water molecule of gadolinium(III) aqua complexes

Hyperfine interactions (HFI) on the nuclei of the first coordination sphere water molecules in a model [Gd(H(2)O)(8)](3+) aqua complex and in the magnetic resonance imaging contrast agent [Gd(DOTA)(H(2)O)](-) were studied theoretically. Density functional theory (DFT) calculations combined with classical molecular dynamics (MD) simulations have been used in order to take into account dynamic effects in aqueous solution. DFT relativistic calculations show a strong spin-polarization of the first coordination sphere water molecules. This spin-polarization leads to a positive (17)O isotropic hyperfine coupling constant (A(iso)((17)O) = 0.58 +/- 0.11 MHz) and to a significant increase of the effective distance (r(eff)(Gd-O) = 2.72 +/- 0.06 A) of dipolar interaction compared to the mean internuclear distance (r(Gd-O) = 2.56 +/- 0.06 A) obtained from the MD trajectory of [Gd(DOTA)(H(2)O)](-) in aqueous solution. The point-dipole model for anisotropic hyperfine interaction overestimates therefore the longitudinal relaxation rate of the (17)O nucleus by approximately 45%. The (1)H isotropic hyperfine coupling constant of the bound water molecule is predicted to be very small (A(iso)((1)H) = 0.03 +/- 0.02 MHz), and the point-dipole approximation for first coordination sphere water protons holds. The calculated hyperfine parameters are in good agreement with available experimental data.     

Membrane resistance to Triton X-100 explored by real-time atomic force microscopy

Lateral segregation of lipids and proteins in biological membranes leads to the formation of detergent-resistant domains, also called "rafts". Understanding the mechanisms governing the biomembrane's resistance to solubilization by detergents is crucial in biochemical research. Here, we used real-time atomic force microscopy (AFM) imaging to visualize the behavior of a model supported lipid bilayer in the presence of different Triton X-100 (TX-100) concentrations. Mixed dioleoylphosphatidylcholine/dipalmitoylphosphatidylcholine (DOPC/DPPC) supported bilayers were prepared by vesicle fusion. Real-time AFM imaging revealed that, at concentrations below the critical micelle concentration (CMC), TX-100 did not solubilize the bilayer, but the DPPC domains were eroded in a time-dependent manner. This effect was attributed to the DPPC molecular packing disorganization by the detergent starting from the DOPC/DPPC interface. Just above the CMC, the detergent led to a complete solubilization of the DOPC matrix, leaving the DPPC domains unaltered. At higher TX-100 concentrations, the DOPC was also immediately removed just after detergent addition, and the DPPC domains remaining on the mica surface appeared to be more swollen and were gradually solubilized. This progressive solubilization of the DPPC remaining phase did not start at the edge of the domains but from holes appearing and expanding at the center of DPPC patches. The swelling of the DPPC domains was directly correlated with TX-100 concentration above the CMC and with detergent intercalation between DPPC molecules. We are convinced that this approach will provide a key system to elucidate the physical mechanisms of membrane solubilization by nonionic detergents.     

Use of molecular dynamics simulations with ab initio SCF calculations for the determination of the deuterium quadrupole coupling constant in liquid water and bond lengths in ice

The deuterium quadrupole coupling constant and asymmetry parameter in heavy water were determined using ab initio SCF calculations. Snapshots from a molecular dynamics simulation were used to give liquid water cluster configurations and the influence of simulation parameters on the quadrupole coupling constant was investigated. The electronic potential model and the number of molecules in the molecular dynamics simulation and the pressure of the system were found to have only a small influence on the quadrupole coupling constant. The average value of the quadrupole coupling constant at room temperature, corrected for the known deficiency of the ab initio calculation in the gas phase, yields a quadrupole coupling constant of 253 kHz, in perfect agreement with the most recent experiments. The oxygen—deuterium bond lengths in ice Ih, ice II, and ice IX were determined using experimental quadrupole coupling constants and a model equation. An averaged bond length of 98.9 pm was obtained for the Ih form, which is approximately 2 pm shorter than that determined by neutron diffraction studies, whereas the bond lengths for the four deuterium sites in ice II and the three sites in ice IX are in fair agreement with experiment. © John Wiley & Sons, Inc.     

Molecular dynamics simulation of the formation, structure, and dynamics of small phospholipid vesicles

Here, we use coarse grained molecular dynamics (MD) simulations to study the spontaneous aggregation of dipalmitoylphosphatidylcholine (DPPC) lipids into small unilamellar vesicles. We show that the aggregation process occurs on a nanosecond time scale, with bicelles and cuplike vesicles formed at intermediate stages. Formation of hemifused vesicles is also observed at higher lipid concentration. With either 25% dipalmitoylphosphatidylethanolamine (DPPE) or lysoPC mixed into the system, the final stages of the aggregation process occur significantly faster. The structure of the spontaneously formed vesicles is analyzed in detail. Microsecond simulations of isolated vesicles reveal significant differences in the packing of the lipids between the inner and outer monolayers, and between PC, PE, and lysoPC. Due to the small size of the vesicles they remain almost perfectly spherical, undergoing very limited shape fluctuations or bilayer undulations. The lipid lateral diffusion rate is found to be faster in the outer than in the inner monolayer. The water permeability coefficient of the pure DPPC vesicles is of the order of 10(-)(3) cm s(-)(1), in agreement with experimental measurements.     

Structure and dynamics of a fluid phase bilayer on a solid support as observed by a molecular dynamics computer simulation

Simulations of a 1-palmitoyl-2-oleoyl- sn-glycero-3-phosphocholine lipid bilayer interacting with a solid surface of hydroxylated nanoporous amorphous silica have been carried out over a range of lipid-solid substrate distances. The porous solid surface allowed the water layer to dynamically adjust its thickness, maintaining equal pressures above and below the membrane bilayer. Qualitative estimates of the force between the surfaces leads to an estimated lipid-silicon distance in very good agreement with the results of neutron scattering experiments. Detailed analysis of the simulation at the separation suggested by experiment shows that for this type of solid support the water layer between surfaces is very narrow, consisting only of bound waters hydrating the lipid head groups and hydrophilic silica surface. The reduced hydration, however, has only minor effects on the head group hydration, the orientation of water molecules at the interface, and the membrane dipole potential. Whereas these structural properties were not sensitive to the presence of the solid substrate, the calculated diffusion coefficient for translation of the lipid molecules was altered significantly by the silica surface.