Molecular dynamics simulation of dipalmitoylphosphatidylcholine membrane with cholesterol sulfate

Using the molecular dynamics simulation technique, we studied the changes occurring in a dipalmitoylphosphatidylcholine (DPPC):cholesterol (CH) membrane at 50 mol% sterol when cholesterol is replaced with cholesterol sulfate (CS). Our simulations were performed at constant pressure and temperature on a nanosecond time scale. We found that 1) the area per DPPC:CS heterodimer is greater than the area of the DPPC:CH heterodimer; 2) CS increases ordering of DPPC acyl chains, but to a lesser extent than CH; 3) the number of hydrogen bonds between DPPC and water is decreased in a CS-containing membrane, but CS forms more water hydrogen bonds than CH; and 4) the membrane dipole potential reverses its sign for a DPPC-CS membrane compared to a DPPC-CH bilayer. We also studied the changes occurring in lipid headgroup conformations and determined the location of CS molecules in the membrane. Our results are in good agreement with the data available from experiments.     

Influence of stigmastanol and stigmastanyl-phosphorylcholine, two plasma cholesterol lowering substances, on synthetic phospholipid membranes. A 2H- and 31P-NMR study.

Cholesterol, stigmastanol, and stigmastanyl-phosphorylcholine (ST-PC) were incorporated into model membranes composed of 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) or 1,2-dioleoyl-sn-glycero-3-phosphocholine (DOPC). POPC and ST-PC were deuterated at the lipid headgroup, DOPC at the cis-double bonds. The influence of the three sterols on the motion and conformation of the lipid headgroups and the hydrocarbon chains was monitored with 2H- and 31P-NMR. All three sterols were freely miscible with the lipid matrix in concentrations of up to 50 mol% without inducing phase separations or nonbilayer structures. However, the molecules exert quite different effects on the phospholipid bilayer. Cholesterol and stigmastanol are largely buried in the hydrocarbon part of the membrane, distinctly restricting the flexing motions of the fatty acyl chains whereas the conformation of the phospholipid headgroups is little affected. In contrast, ST-PC is anchored with its headgroup in the layer of phospholipid dipoles, preventing an extensive penetration of the sterol ring into the hydrocarbon layer. Hence ST-PC has almost no effect on the hydrocarbon chains but induces a characteristic conformational change of the phospholipid headgroups. The 2H- and 31P-NMR spectra of mixed phospholipid/ST-PC membranes further demonstrate that the PC headgroup of ST-PC has a similar orientation as the surrounding phosphatidylcholine headgroups. For both types of molecules the -P-N+ dipole is essentially parallel to the membrane surface. Addition of ST-PC induces a small rotation of the POPC headgroup towards the water phase.     

The influence of cholesterol on interactions and dynamics of ibuprofen in a lipid bilayer

In this work, molecular dynamics (MD) simulations with atomistic details were performed to examine the influence of the cholesterol on the interactions and the partitioning of the hydrophobic drug ibuprofen in a fully hydrated 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC) bilayer. Analysis of MD simulations indicated that ibuprofen molecules prefer to be located in the hydrophobic acyl chain region of DMPC/cholesterol bilayers. This distribution decreases the lateral motion of lipid molecules. The presence of ibuprofen molecules in the bilayers with 0 and 25 mol% cholesterol increases the ordering of hydrocarbon tails of lipids whereas for the bilayers with 50 mol% cholesterol, ibuprofen molecules perturb the flexible chains of DMPC lipids which leads to the reduction of the acyl chain order parameter. The potential of the mean force (PMF) method was used to calculate the free energy profile for the transferring of an ibuprofen molecule from the bulk water into the DMPC/cholesterol membranes. The PMF studies indicated that the presence of 50 mol% cholesterol in the bilayers increases the free energy barrier and slows down the permeation of the ibuprofen drug across the DMPC bilayer. This can be due to the condensing and ordering effects of the cholesterol on the bilayer.     

Molecular-dynamic simulation of DPPC bilayer in different phase state: Hydration and electric field distribution in the presence of Be<Superscript>2+</Superscript> cations

New molecular-dynamic topology of phosphatidylcholine bilayer (DPPC) in total atomic OPLS force field was developed and used to study the structural characteristics of liquid-crystalline and gel state of lipid bilayer in the absence and in the presence of Na⁺ and Be²⁺ cations adsorbed at the interface and different in their affinity. The parameters of bilayer geometry, the amount of surface water, and the electrostatic potential distribution were estimated quantitatively from the simulation in both phase states. The azimuthal angle of hydrocarbon chains was found nearly the same in the region of each monolayer in gel state. The amount of surface water decreases upon bilayer “freezing” mainly by loss of water molecules not participating in H-bonds between lipid headgroups. The cation adsorption was shown to have a small effect on these H-bonded water molecules, whereas Be²⁺ appeared to retain surface water in the bilayer upon its freezing. The electric potential distribution in the normal direction to the membrane-water interface had a similar shape in any bilayer phase state regardless of the presence of the adsorbed cations. Analysis of the microscopic nature of the electric potential revealed a mutual compensation of the contributions of the main structural components of the system (lipids, water, and ions). The boundary potential increased by 116 mV for pure DPPC, by 212 mV in the presence of Na⁺, and by 133 mV in the presence of Be²⁺ upon the phase transition of bilayer to the gel state. The boundary potential difference in the presence of Na⁺ and Be²⁺ and its change at the bilayer phase transition are in a good agreement with the experimental data published earlier [Ermakov Yu.A., 1993].     

Structure of dipalmitoylphosphatidylcholine/cholesterol bilayer at low and high cholesterol concentrations: molecular dynamics simulation.

By using molecular dynamics simulation technique we studied the changes occurring in membranes constructed of dipalmitoylphosphatidylcholine (DPPC) and cholesterol at 8:1 and 1:1 ratios. We tested two different initial arrangements of cholesterol molecules for a 1:1 ratio. The main difference between two initial structures is the average number of nearest-neighbor DPPC molecules around the cholesterol molecule. Our simulations were performed at constant temperature (T = 50 degrees C) and pressure (P = 0 atm). Durations of the runs were 2 ns. The structure of the DPPC/cholesterol membrane was characterized by calculating the order parameter profiles for the hydrocarbon chains, atom distributions, average number of gauche defects, and membrane dipole potentials. We found that adding cholesterol to membranes results in a condensing effect: the average area of membrane becomes smaller, hydrocarbon chains of DPPC have higher order, and the probability of gauche defects in DPPC tails is lower. Our results are in agreement with the data available from experiments.     

The interfacial structure of phospholipid bilayers: differential scanning calorimetry and Fourier transform infrared spectroscopic studies of 1,2-dipalmitoyl-sn-glycero-3-phosphorylcholine and its dialkyl and acyl-alkyl analogs.

The thermotropic phase behavior of aqueous dispersions of dipalmitoylphosphatidylcholine (DPPC) and its 1,2-dialkyl, 1-acyl 2-alkyl and 1-alkyl 2-acyl analogs was examined by differential scanning calorimetry, and the organization of these molecules in those hydrated bilayers was studied by Fourier transform infrared spectroscopy. The calorimetric data indicate that substitution of either or both of the acyl chains of DPPC with the corresponding ether-linked hydrocarbon chain results in relatively small increases in the temperature (< 4 degrees C) and enthalpy (< 1 kcal/mol) of the lipid chain-melting phase transition. The spectroscopic data reveal that replacement of one or both of the ester-linked hydrocarbon chains of DPPC with its ether-linked analog causes structural changes in the bilayer assembly, which result in an increase in the polarity of the local environments of the phosphate headgroups and of the ester carbonyl groups at the bilayer polar/apolar interface. The latter observation is unexpected, given that ester linkages are considered to be intrinsically more polar that ether linkages. This finding cannot be satisfactorily rationalized unless the conformation of the glycerol backbones of the analogs containing ether-linked hydrocarbon chains differs significantly from that of diacyl glycerolipids such as DPPC. A comparison of the alpha-methylene scissoring bands and the methylene wagging band progressions of these lipids with the corresponding absorption bands of specifically chain-perdeuterated analogs of DPPC also supports the conclusion that replacement of the ester-linked hydrocarbon chains of DPPC with the corresponding ether-linked analog induces conformational changes in the lipid glycerol backbone. The suggestion that the conformation of glycerol backbones in the alkyl-acyl and dialkyl derivatives of DPPC differs from that of the naturally occurring 1,2-diacyl glycerolipid suggests that mono- and di-alkyl glycerolipids may not be good models of their diacyl analogs. These results, and previously published evidence that DPPC analogs with ether-linked hydrocarbon chains spontaneously form chain-interdigitated gel phases at low temperatures, clearly indicate that the properties of lipid bilayers can be substantially altered by small changes in the chemical structures of their polar/polar interfaces, and highlight the critical role of the interfacial region as a determinant of the structure and organization of lipid assemblies.     

The effect of cholesterol on measured interaction and compressibility of dipalmitoylphosphatidylcholine bilayers

We have examined the phase diagram of dipalmitoylphosphatidylcholine (DPPC)--cholesterol-water mixtures at low cholesterol content, and report phase separation between 3 and 10 mol% cholesterol. The two lamellar phases at equilibrium in this region appear to be pure DPPC and 11 mol% cholesterol in DPPC. For these two lamellar phases, which are made up of alternating layers of water and bimolecular lipid leaflets, we have measured the forces of interaction between leaflets and the lateral pressure and compressibility of the leaflets. Both bilayers experience a strong repulsive force when forced together only a few ?ngstr?ms (1 A = 0.1 nm) closer than their maximum separation in excess water. However, the presence of 11 mol% cholesterol causes the bilayers to move apart of 35-A separation from the 19-A characteristic of pure DPPC in excess water. This swelling may result from a decrease in van der Waals attraction between bilayers or from an increase in bilayer repulsion. Differences in bilayer interaction can be a cause for phase separation. More importantly these differences can cause changes in the composition of regions of membranes approaching contact. At 11 mol%, cholesterol substantially increases the lateral compressibility of DPPC bilayers leading to higher lateral density fluctuations and potentially higher bilayer permeability.     

Behavior of cholesterol and its effect on head group and chain conformations in lipid bilayers: a molecular dynamics study.

Cholesterol molecules were put into a computer-modeled hydrated bilayer of dimyristoyl phosphatidyl choline molecules, and molecular dynamics simulations were run to characterize the effect of this important molecule on membrane structure and dynamics. The effect was judged by observing differences in order parameters, tilt angles, and the fraction of gauche bonds along the hydrocarbon chains between lipids adjacent to cholesterol molecules and comparing them with those further away. It was observed that cholesterol causes an increase in the fraction of trans dihedrals and motional ordering of chains close to the rigid steroid ring system with a decrease in the kink population. The hydrogen-bonding interactions between cholesterol and lipid molecules were determined from radial distribution calculations and showed the cholesterol hydroxyl groups either solvated by water, or forming hydrogen bond contacts with the oxygens of lipid carbonyl and phosphate groups. The dynamics and conformation of the cholesterol molecules were investigated and it was seen that they had a smaller tilt with respect to the bilayer normal than the lipid chains and furthermore that the hydrocarbon tail of the cholesterol was conformationally flexible.     

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.     

Molecular dynamics simulations predict a tilted orientation for the helical region of dynorphin A(1-17) in dimyristoylphosphatidylcholine bilayers

The structural properties of the endogenous opioid peptide dynorphin A(1-17) (DynA), a potential analgesic, were studied with molecular dynamics simulations in dimyristoylphosphatidylcholine bilayers. Starting with the known NMR structure of the peptide in dodecylphosphocholine micelles, the N-terminal helical segment of DynA (encompassing residues 1-10) was initially inserted in the bilayer in a perpendicular orientation with respect to the membrane plane. Parallel simulations were carried out from two starting structures, systems A and B, that differ by 4 A in the vertical positioning of the peptide helix. The complex consisted of approximately 26,400 atoms (dynorphin + 86 lipids + approximately 5300 waters). After >2 ns of simulation, which included >1 ns of equilibration, the orientation of the helical segment of DynA had undergone a transition from parallel to tilted with respect to the bilayer normal in both the A and B systems. When the helix axis achieved a approximately 50 degrees angle with the bilayer normal, it remained stable for the next 1 ns of simulation. The two simulations with different starting points converged to the same final structure, with the helix inserted in the bilayer throughout the simulations. Analysis shows that the tilted orientation adopted by the N-terminal helix is due to specific interactions of residues in the DynA sequence with phospholipid headgroups, water, and the hydrocarbon chains. Key elements are the "snorkel model"-type interactions of arginine side chains, the stabilization of the N-terminal hydrophobic sequence in the lipid environment, and the specific interactions of the first residue, Tyr. Water penetration within the bilayer is facilitated by the immersed DynA, but it is not uniform around the surface of the helix. Many water molecules surround the arginine side chains, while water penetration near the helical surface formed by hydrophobic residues is negligible. A mechanism of receptor interaction is proposed for DynA, involving the tilted orientation observed from these simulations of the peptide in the lipid bilayer.     

Constant normal pressure, constant surface tension, and constant temperature molecular dynamics simulation of hydrated 1,2-dilignoceroylphosphatidylcholine monolayer

A constant normal pressure, constant surface tension, and constant temperature (NP(N)gammaT) molecular dynamics (MD) simulation of the liquid condensed phase of a 1,2-dilignoceroylphosphatidylcholine (DLGPC) monolayer has been performed at 293.15 K. A DLGPC molecule has two saturated 24-carbon acyl chains, giving the hydrocarbon core thickness of the monolayer approximately 28 A, which is close to the hydrocarbon core thickness of a membrane of a living system. NP(N)gammaT ensemble was used to reproduce the experimental observations, such as area/lipid, because surface tension is an essential factor in determining the monolayer structure. Data analysis on DLGPC/water monolayer shows that various liquid condensed-phase properties of the monolayer have been well reproduced from the simulation, indicating that surface tension 22.9 mN/M used in the simulation is an appropriate condition for the condensed-phase NP(N)gammaT simulation. The simulation results suggest that this long-chain phospholipid monolayer shares many structural characteristics with typical short-chain 1,2-diacylphosphatidylcholine systems, such as DPPC/water monolayer in the condensed phase and DPPC/water bilayer in the gel phase. Furthermore, it was found that DLGPC/water monolayer has almost completely rotationally disordered acyl chains, which have not been observed so far in short-chain 1,2-diacylphosphatidylcholine/water bilayers. This study indicates the good biological relevance of the DLGPC/water monolayer which might be useful in protein/lipid studies to reveal protein structure and protein/lipid interactions in a membrane environment.     

Differential dynamic and structural behavior of lipid-cholesterol domains in model membranes

Changes in the cholesterol (Chol) content of biological membranes are known to alter the physicochemical properties of the lipid lamella and consequently the function of membrane-associated enzymes. To characterize these changes, we used steady-state and time resolved fluorescence spectroscopy and two photon-excitation microscopy techniques. The membrane systems were chosen according to the techniques that were used: large unilamellar vesicles (LUVs) for cuvette and giant unilamellar vesicles (GUVs) for microscopy measurements; they were prepared from dipalmitoyl phosphatidylcholine (DPPC) and dioctadecyl phosphatidylcholine (DOPC) in mixtures that are well known to form lipid domains. Two fluorescent probes, which insert into different regions of the bilayer, were selected: 1,6-diphenyl-1,3,5-hexatriene (DPH) was located at the deep hydrophobic core of the acyl chain regions and 2-dimethylamino-6-lauroylnaphthalene (Laurdan) at the hydrophilic-hydrophobic membrane interface. Our spectroscopy results show that (i) the changes induced by cholesterol in the deep hydrophobic phospholipid acyl chain domain are different from the ones observed in the superficial region of the hydrophilic-hydrophobic interface, and these changes depend on the state of the lamella and (ii) the incorporation of cholesterol into the lamella induces an increase in the orientation dynamics in the deep region of the phospholipid acyl chains with a corresponding decrease in the orientation at the region close to the polar lipid headgroups. The microscopy data from DOPC/DPPC/Chol GUVs using Laurdan generalized polarization (Laurdan GP) suggest that a high cholesterol content in the bilayer weakens the stability of the water hydrogen bond network and hence the stability of the liquid-ordered phase (Lo).     

Effect of cholesterol on the structure and dynamic properties of unsaturated phospholipid bilayers

Molecular dynamics (MD) computer simulations of five different hydrated unsaturated phosphatidylcholine lipid bilayers built up by 18:0/18:1(n-9)cis PC, 18:0/18:2(n-6)cis PC, 18:0/18:3(n-3)cis PC, 18:0/20:4(n-6)cis PC, and 18:0/22:6(n-3)cis PC molecules with 40 mol% cholesterol, and the same five pure phosphatidylcholine bilayers have been performed at 303 K. The simulation box of a lipid bilayer contained 96 phosphatidylcholines, 64 cholesterols, and 3840 water molecules (48 phosphatidylcholine molecules and 32 cholesterols per layer and 24 water molecules per phospholipid or cholesterol in each case). The lateral self-diffusion coefficients of the lipids in these systems and mass density profiles with respect to the bilayer normal have been analyzed. It has been found that the lateral diffusion coefficients of phosphatidylcholine molecules increase with increasing number of double bonds in one of the lipid chains, both in pure bilayers and in bilayers with cholesterol. It has been found as well that the lateral diffusion coefficient of phosphatidylcholine molecules of a lipid bilayer with 40 mol% cholesterol is smaller than that for the corresponding pure phosphatidylcholine bilayer.     

Effect of hydrostatic pressure on water penetration and rotational dynamics in phospholipid-cholesterol bilayers.

The effect of high hydrostatic pressure on the lipid bilayer hydration, the mean order parameter, and rotational dynamics of 1,2-dipalmitoyl-sn-glycero-3-phosphocholine (DPPC) and 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) cholesterol vesicles has been studied by time-resolved fluorescence spectroscopy up to 1500 bar. Whereas the degree of hydration in the lipid headgroup and interfacial region was assessed from fluorescence lifetime data using the probe 1-(4-trimethylammonium-phenyl)-6-phenyl-1,3,5-hexatriene (TMA-DPH), the corresponding information in the upper acyl chain region was estimated from its effect on the fluorescence lifetime of and 3-(diphenylhexatrienyl)propyl-trimethylammonium (TMAP-DPH). The lifetime data indicate a greater level of interfacial hydration for DPPC bilayers than for POPC bilayers, but there is no marked difference in interchain hydration of the two bilayer systems. The addition of cholesterol at levels from 30 to 50 mol% to DPPC has a greater effect on the increase of hydrophobicity in the interfacial region of the bilayer than the application of hydrostatic pressure of several hundred to 1000 bar. Although the same trend is observed in the corresponding system, POPC/30 mol% cholesterol, the observed effects are markedly less pronounced. Whereas the rotational correlation times of the fluorophores decrease in passing the pressure-induced liquid-crystalline to gel phase transition of DPPC, the wobbling diffusion coefficient remains essentially unchanged. The wobbling diffusion constant of the two fluorophores changes markedly upon incorporation of 30 mol% cholesterol, and increases at higher pressures, also in the case of POPC/30 mol% cholesterol. The observed effects are discussed in terms of changes in the rotational characteristics of the fluorophores and the phase-state of the lipid mixture. The results demonstrate the ability of cholesterol to adjust the structural and dynamic properties of membranes composed of different phospholipid components, and to efficiently regulate the motional freedom and hydrophobicity of membranes, so that they can withstand even drastic changes in environmental conditions, such as high external hydrostatic pressure.     

The effect of cholesterol in a liposomal Muc1 vaccine

A liposomal Muc1 mucin vaccine for treatment of adenocarcinomas was formulated by incorporating a synthetic Muc1 mucin-based lipopeptide and Lipid A into a DPPC/cholesterol bilayer. Vaccination of mice with the liposomal formulation produced a peptide-specific immune response dependent on the cholesterol content. The response occurred at a threshold of 20-23 mol% cholesterol, and was optimal at cholesterol levels of > or =30 mol%. To understand this cholesterol dependency, we studied the effect of cholesterol on the liposomal bilayer and surface properties. Freeze-fracture electron microscopy showed a unique surface texture that was codependent upon cholesterol (> or =20 mol%) and lipopeptide content. Fluorescence anisotropy measurements exhibited a significant decrease in the rotational motion of 1,6-diphenyl-1,3,5-hexatriene in formulations containing >20 mol% cholesterol and only in the presence of the lipopeptide. At 20 mol% cholesterol and with lipopeptide, DSC showed a significant increase in the main phase transition of the DPPC bilayers, while Raman spectroscopy indicated a more ordered arrangement of DPPC molecules compared to control liposomes containing DPPC/cholesterol alone. Taken together, the data suggest the presence of lipopeptide-rich microdomains at and above a threshold of 20 mol% cholesterol that may play a role in the induction of a peptide-specific immunological response.     

ESR and monolayer study of the localization of coenzyme Q10 in artificial membranes

The data obtained from the ESR experiments show a complex, depth dependent effect of CoQ10 on the lipid molecules mobility in the bilayer. These effects depend both on its concentration and the temperature. CoQ10 disturbs not only the hydrophobic core of the membrane but also the region close to the hydrophilic headgroups of phospholipids. Both these effects could be explained by the fact that the high hydrophobicity of CoQ10 causes the molecules to position itself in the interior of the bilayer, but at the same time its water seeking headgroup is located close to the region of the polar headgrops of membrane lipids. The presence of CoQ10 in the hydrophobic core has further implications on the properties of membrane intrinsic domain. Results of monolayer experiments indicate that CoQ10 may form aggregates when mixed with PC molecules in the lipid hydrocarbon chain-length dependent manner. CoQ10 is not fully miscible with DMPC or DPPC but it is well miscible with the long-chain DSPC molecules. Our suggestion is that CoQ10 when present in long-chain phospholipid bilayer, interacts with saturated fatty acyl-chains and adapt the structure which allows such interactions: either parallel to the saturated acyl chains or "pseudo-ring" conformation resembling sterol structure.     

Effects of oxygenated sterol on phospholipid bilayer properties: a molecular dynamics simulation

The properties of dipalmitoylphosphatidylcholine (DPPC):6-ketocholestanol bilayer at 50 mol% sterol were studied using the molecular dynamics simulation technique. Our simulations were performed at constant pressure and temperature on a nanosecond time scale. Data from this simulation were compared to the results of our previous simulations on DPPC and DPPC-cholesterol bilayers. We conclude that the differences in the properties of membranes with cholesterol and ketocholestanol are due to the difference in 6-ketocholestanol and cholesterol location in the bilayer. The presence of the keto group in ketocholestanol moves the sterol towards the polar region closer to interface with water. We predict that similar mechanisms would govern the properties of membranes with other oxygenated sterols, such as for example 7-ketocholesterol. Results of our simulations are in a good agreement with the data available from the experiment.     

Transmembrane helix structure, dynamics, and interactions: multi-nanosecond molecular dynamics simulations.

To probe the fundamentals of membrane/protein interactions, all-atom multi-nanosecond molecular dynamics simulations were conducted on a single transmembrane poly(32)alanine helix in a fully solvated dimyristoyphosphatidylcholine (DMPC) bilayer. The central 12 residues, which interact only with the lipid hydrocarbon chains, maintained a very stable helical structure. Helical regions extended beyond these central 12 residues, but interactions with the lipid fatty-acyl ester linkages, the lipid headgroups, and water molecules made the helix less stable in this region. The C and N termini, exposed largely to water, existed as random coils. As a whole, the helix tilted substantially, from perpendicular to the bilayer plane (0 degree) to a 30 degrees tilt. The helix experienced a bend at its middle, and the two halves of the helix at times assumed substantially different tilts. Frequent hydrogen bonding, of up to 0.7 ns in duration, occurred between peptide and lipid molecules. This resulted in correlated translational diffusion between the helix and a few lipid molecules. Because of the large variation in lipid conformation, the lipid environment of the peptide was not well defined in terms of "annular" lipids and on average consisted of 18 lipid molecules. When compared with a "neat" bilayer without peptide, no significant difference was seen in the bilayer thickness, lipid conformations or diffusion, or headgroup orientation. However, the lipid hydrocarbon chain order parameters showed a significant decrease in order, especially in those methylene groups closest to the headgroup.     

The partition and transport behavior of cytotoxic ionic liquids (ILs) through the DPPC bilayer: Insights from molecular dynamics simulation

A molecular dynamics (MD) simulation with atomistic details was performed to examine the partitioning and transport behavior of moderately cytotoxic ionic liquids (ILs), namely choline bis(2-ethylhexyl) phosphate (CBEH), choline bis(2,4,4-trimethylpentyl) phosphinate (CTMP) and choline O,O-diethyl dithiophosphate (CDEP) in a fully hydrated dipalmitoylphosphatidylcholine (DPPC) bilayer in the fluid phase at 323?K. The structure of ILs was so selected to understand if the role of dipole and dispersion forces in the ILs distribution in the membrane can be possible. Several analyses including mass density, electrostatic potential, order parameter, diffusion coefficients and hydrogen bond formation, was carried out to determine the precise location of the anionic species inside the membrane. Moreover, the potential of the mean force (PMF) method was used to calculate free energy profile for transferring anionic species from the DPPC membrane into the bulk water. While less cytotoxic DEP is located within the bulk water, more cytotoxic TMP and BEH ILs were found to remain in the membrane and the energy barrier for crossing through the bilayer center of BEH was higher. Various ILs have no significant effect on P–N vector. The thickness of lipid bilayer decreased in all systems comprising ILs, while area per lipid increased.