Examining the effects of cholesterol on model membranes at high temperatures: Laurdan and Patman see it differently

At high temperature, the presence of cholesterol in phospholipid membranes alters the influence of membrane dipoles, including water molecules, on naphthalene-based fluorescent probes such as Laurdan and Patman (solvatochromism). Although both of these probes report identical changes to their emission spectra as a function of temperature in pure phosphatidylcholine bilayers, they differ in their response to cholesterol. Computer simulations of the spectra based on a simple model of solvatochromism indicated that the presence of cholesterol reduces the probability of bilayer dipole relaxation and also blunts the tendency of heat to enhance that probability. While the overall effect of cholesterol on membrane dipoles was detected identically by the two probes, Laurdan was influenced much more by the additional effect on temperature sensitivity than was Patman. A comparison of the fluorescence data with simulations using a coarse-grained bilayer model (de Meyer et al., 2010) suggested that these probes may be differentially sensitive to two closely related properties distinguishable in the presence of cholesterol. Specifically, Patman fluorescence correlated best with the average phospholipid acyl chain order. On the other hand, Laurdan fluorescence tracked more closely with the area per lipid molecule which, although affected generally by chain order, is also impacted by additional membrane-condensing effects of cholesterol. We postulate that this difference between Laurdan and Patman may be attributed to the bulkier charged headgroup of Patman which may cause the probe to preferentially locate in juxtaposition to the diminutive headgroup of cholesterol as the membrane condenses.     

Membrane packing geometry of diphytanoylphosphatidylcholine is highly sensitive to hydration: phospholipid polymorphism induced by molecular rearrangement in the headgroup region.

Diphytanoylphosphatidylcholine (DPhPC) has often been used in the study of protein-lipid interaction and membrane channel activity, because of the general belief that it has high bilayer stability, low ion leakage, and fatty acyl packing comparable to that of phospholipid bilayers in the liquid-crystalline state. In this solid-state 31P and 2H NMR study, we find that the membrane packing geometry and headgroup orientation of DPhPC are highly sensitive to the temperature studied and its water content. The phosphocholine headgroup of DPhPC starts to change its orientation at a water content as high as approximately 16 water molecules per lipid, as evidenced by hydration-dependent 2H NMR study at room temperature. In addition, a temperature-induced structural transition in the headgroup orientation is detected in the temperature range of approximately 20-60 degrees C for lipids with approximately 8-11 water molecules per DPhPC. Dehydration of the lipid by one more water molecule leads to a nonlamellar, presumably cubic, phase formation. The lipid packing becomes a hexagonal phase at approximately 6 water molecules per lipid. A phase diagram of DPhPC in the temperature range of -40 degrees C to 80 degrees C is thus constructed on the basis of NMR results. The newly observed hydration-dependent DPhPC lipid polymorphism emphasizes the importance of molecular packing in the headgroup region in modulating membrane structure and protein-induced pore formation of the DPhPC bilayer.     

A molecular dynamics study of the response of lipid bilayers and monolayers to trehalose

Surface tensions evaluated from molecular dynamics simulations of fully hydrated dipalmitoylphosphatidylcholine bilayers and monolayers at surface areas/lipid of 54, 64, and 80 A2 are uniformly lowered 4-8 dyn/cm upon addition of trehalose in a 1:2 trehalose/lipid ratio. Constant surface tension simulations of bilayers yield the complementary result: an increase in surface area consistent with the surface pressure-surface area (pi-A) isotherms. Hydrogen bonding by trehalose, replacement of waters in the headgroup region, and modulation of the dipole potential are all similar in bilayers and monolayers at the same surface area. These results strongly support the assumption that experimental measurements on the interactions of surface active components such as trehalose with monolayers can yield quantitative insight to their effects on bilayers. The simulations also indicate that the 20-30 dyn/cm difference in surface tension of the bilayer leaflet and monolayer arises from differences in the chain regions, not the headgroup/water interfaces.     

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.     

A Molecular Dynamics Study of DMPC Lipid Bilayers Interacting with Dimethylsulfoxide–Water Mixtures

The diffusion process of dimethylsulfoxide (DMSO) through zwitterionic dimyristoylphosphatidylcholine (DMPC) lipid bilayer was studied by means of molecular dynamics (MD) simulations. To account for the cryoprotectant concentration difference between the inside and the outside of the cell, dual DMPC lipid bilayers which separate two aqueous reservoirs with and without DMSO were modeled. The initial configuration of the simulation model had DMSO molecules present in one of the aqueous phases (outside the cell) at two different concentrations of ~3 and ~6?mol%. MD simulations were performed on the systems for 50?ns at 323?K and 1?bar. Although the simulation time considered in the study was insufficient for the DMSO molecules to reach the other aqueous phase and equilibrium, early stages of the diffusion process indicated that DMSO molecules had a tendency to diffuse towards the other aqueous phase. The effects of DMSO on bilayer structural characteristics during the diffusion process were investigated. Simulations were analyzed to correlate the following properties of lipid bilayers in the presence of two different aqueous phases: area per lipid, lipid thickness, mass density profiles, lipid tail order parameter and water dipole orientation. Area per lipid calculated for the leaflet facing the aqueous DMSO?Cwater mixture did not show any significant difference compared to area per lipid for the DMSO-free pure DMPC bilayer. Mass density profiles revealed that DMSO molecules had a strong tendency to diffuse toward the aqueous phase with pure water. The lipid tail order parameter calculated for the sn-1 tail of the leaflet facing the aqueous DMSO?Cwater mixture showed that the ordering of lipid tails decreased compared to the leaflet exposed to pure water. However, the ordering of lipid tails in a system where a single bilayer is hydrated by an aqueous DMSO?Cwater mixture is far lower.     

Molecular dynamics simulations and 2H NMR study of the GalCer/DPPG lipid bilayer

Molecular dynamics simulations were performed on a two-component lipid bilayer system in the liquid crystalline phase at constant pressure and constant temperature. The lipid bilayers were composed of a mixture of neutral galactosylceramide (GalCer) and charged dipalmitoylphosphatidylglycerol (DPPG) lipid molecules. Two lipid bilayer systems were prepared with GalCer:DPPG ratio 9:1 (10%-DPPG system) and 3:1 (25%-DPPG system). The 10%-DPPG system represents a collapsed state lipid bilayer, with a narrow water space between the bilayers, and the 25%-DPPG system represents an expanded state with a fluid space of approximately 10 nm. The number of lipid molecules used in each simulation was 1024, and the length of the production run simulation was 10 ns. The simulations were validated by comparing the results with experimental data for several important aspects of the bilayer structure and dynamics. Deuterium order parameters obtained from (2)H NMR experiments for DPPG chains are in a very good agreement with those obtained from molecular dynamics simulations. The surface area per GalCer lipid molecule was estimated to be 0.608 +/- 0.011 nm(2). From the simulated electron density profiles, the bilayer thickness defined as the distance between the phosphorus peaks across the bilayer was calculated to be 4.21 nm. Both simulation systems revealed a tendency for cooperative bilayer undulations, as expected in the liquid crystalline phase. The interaction of water with the GalCer and DPPG oxygen atoms results in a strong water ordering in a spherical hydration shell and the formation of hydrogen bonds (H-bonds). Each GalCer lipid molecule makes 8.6 +/- 0.1 H-bonds with the surrounding water, whereas each DPPG lipid molecule makes 8.3 +/- 0.1 H-bonds. The number of water molecules per GalCer or DPPG in the hydration shell was estimated to be 10-11 from an analysis of the radial distribution functions. The formation of the intermolecular hydrogen bonds was observed between hydroxyl groups from the opposing GalCer sugar headgroups, giving an energy of adhesion in the range between -1.0 and -3.4 erg/cm(2). We suggest that this value is the contribution of the hydrogen-bond component to the net adhesion energy between GalCer bilayers in the liquid crystalline phase.     

Giant Vesicles under Oxidative Stress Induced by a Membrane-Anchored Photosensitizer

We have synthesized the amphiphile photosensitizer PE-porph consisting of a porphyrin bound to a lipid headgroup. We studied by optical microscopy the response to light irradiation of giant unilamellar vesicles of mixtures of unsaturated phosphatidylcholine lipids and PE-porph. In this configuration, singlet oxygen is produced at the bilayer surface by the anchored porphyrin. Under irradiation, the PE-porph decorated giant unilamellar vesicles exhibit a rapid increase in surface area with concomitant morphological changes. We quantify the surface area increase of the bilayers as a function of time and photosensitizer molar fraction. We attribute this expansion to hydroperoxide formation by the reaction of the singlet oxygen with the unsaturated bonds. Considering data from numeric simulations of relative area increase per phospholipid oxidized (15%), we measure the efficiency of the oxidative reactions. We conclude that for every 270 singlet oxygen molecules produced by the layer of anchored porphyrins, one eventually reacts to generate a hydroperoxide species. Remarkably, the integrity of the membrane is preserved in the full experimental range explored here, up to a hydroperoxide content of 60%, inducing an 8% relative area expansion.     

Alcohols reduce lateral membrane pressures: predictions from molecular theory

We explore the effects of alcohols on fluid lipid bilayers using a molecular theory with a coarse-grained model. We show that the trends predicted from the theory in the changes in area per lipid, alcohol concentration in the bilayer, and area compressibility modulus, as a function of alcohol chain length and of the alcohol concentration in the solvent far from the bilayer, follow those found experimentally. We then use the theory to study the effect of added alcohol on the lateral pressure profile across the membrane, and find that added alcohol reduces the surface tensions at both the headgroup/solvent and headgroup/tailgroup interfaces, as well as the lateral pressures in the headgroup and tailgroup regions. These changes in lateral pressures could affect the conformations of membrane proteins, providing a nonspecific mechanism for the biological effects of alcohols on cells.     

Effects of cholesterol concentration on the interaction of cytarabine with lipid membranes: a molecular dynamics simulation study

Liposomal cytarabine, DepoCyt, is a chemotherapy agent which is used in cancer treatment. This form of cytarabine has more efficacy and fewer side effects relative to the other forms. Since DepoCyt contains the cytarabine encapsulated within phosphatidylcholine and the sterol molecules, we modeled dioleoylphosphatidylcholine (DOPC)/cholesterol bilayer membrane as a carrier for cytarabine to study drug–bilayer interactions. For this purpose, we performed a series of united-atom molecular dynamics (MD) simulations for 25?ns to investigate the interactions between cytarabine and cholesterol-containing DOPC lipid bilayers. Only the uncharged form of cytarabine molecule was investigated. In this study, different levels of the cholesterol content (0, 20, and 40%) were used. MD simulations allowed us to determine dynamical and structural properties of the bilayer membrane and to estimate the preferred location and orientation of the cytarabine molecule inside the bilayer membrane. Properties such as membrane thickness, area per lipid, diffusion coefficient, mass density, bilayer packing, order parameters, and intermolecular interactions were examined. The results show that by increasing the cholesterol concentration in the lipid bilayers, the bilayer thickness increases and area per lipid decreases. Moreover, in accordance with the experiments, our calculations show that cholesterol molecules have ordering effect on the hydrocarbon acyl chains. Furthermore, the cytarabine molecule preferentially occupies the polar region of the lipid head groups to form specific interactions (hydrogen bonds). Our results fully support the experimental data. Our finding about drug–bilayer interaction is crucial for the liposomal drug design.     

Combined monte carlo and molecular dynamics simulation of fully hydrated dioleyl and palmitoyl-oleyl phosphatidylcholine lipid bilayers

We have applied a new equilibration procedure for the atomic level simulation of a hydrated lipid bilayer to hydrated bilayers of dioleyl-phosphatidylcholine (DOPC) and palmitoyl-oleyl phosphatidylcholine (POPC). 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 applied to bilayers of 128 molecules of POPC with 4628 water molecules, and 128 molecules of DOPC with 4825 water molecules. Progress toward equilibration is almost three times as fast in central processing unit (CPU) time compared with a purely molecular dynamics (MD) simulation. Equilibration is complete, as judged by the lack of energy drift in 200-ps runs of continuous MD. After the equilibrium state was reached, as determined by agreement between the simulation volume per lipid molecule with experiment, continuous MD was run in an ensemble in which the lateral area was restrained to fluctuate about a mean value and a pressure of 1 atm applied normal to the bilayer surface. Three separate continuous MD runs, 200 ps in duration each, separated by 10,000 CBMC steps, were carried out for each system. Properties of the systems were calculated and averaged over the three separate runs. Results of the simulations are presented and compared with experimental data and with other recent simulations of POPC and DOPC. Analysis of the hydration environment in the headgroups supports a mechanism by which unsaturation contributes to reduced transition temperatures. In this view, the relatively horizontal orientation of the unsaturated bond increases the area per lipid, resulting in increased water penetration between the headgroups. As a result the headgroup-headgroup interactions are attenuated and shielded, and this contributes to the lowered transition temperature.     

Structural determinants of water permeability through the lipid membrane

Despite intense study over many years, the mechanisms by which water and small nonelectrolytes cross lipid bilayers remain unclear. While prior studies of permeability through membranes have focused on solute characteristics, such as size, polarity, and partition coefficient in hydrophobic solvent, we focus here on water permeability in seven single component bilayers composed of different lipids, five with phosphatidylcholine headgroups and different chain lengths and unsaturation, one with a phosphatidylserine headgroup, and one with a phosphatidylethanolamine headgroup. We find that water permeability correlates most strongly with the area/lipid and is poorly correlated with bilayer thickness and other previously determined structural and mechanical properties of these single component bilayers. These results suggest a new model for permeability that is developed in the accompanying theoretical paper in which the area occupied by the lipid is the major determinant and the hydrocarbon thickness is a secondary determinant. Cholesterol was also incorporated into DOPC bilayers and X-ray diffuse scattering was used to determine quantitative structure with the result that the area occupied by DOPC in the membrane decreases while bilayer thickness increases in a correlated way because lipid volume does not change. The water permeability decreases with added cholesterol and it correlates in a different way from pure lipids with area per lipid, bilayer thickness, and also with area compressibility.     

How type II diabetes-related islet amyloid polypeptide damages lipid bilayers

A leading hypothesis for the decimation of insulin-producing β-cells in type 2 diabetes attributes the cause to islet amyloid polypeptide (IAPP) for its deleterious effects on the cell membranes. This idea has produced extensive investigations on human IAPP (hIAPP) and its interactions with lipid bilayers. However, it is still difficult to correlate the peptide-lipid interactions with its effects on islet cells in culture. The hIAPP fibrils have been shown to interact with lipids and damage lipid bilayers, but appear to have no effect on islet cells in culture. Thus, a modified amyloid hypothesis assumes that the toxicity is caused by hIAPP oligomers, which are not preamyloid fibrils or protofibrils. However, so far such oligomers have not been isolated or identified. The hIAPP monomers also bind to lipid bilayers, but the mode of interaction is not clear. Here, we performed two types of experiments that, to our knowledge, have not been done before. We used x-ray diffraction, in conjunction with circular dichroism measurement, to reveal the location of the peptide bound to a lipid bilayer. We also investigated the effects of hIAPP on giant unilamellar vesicles at various peptide concentrations. We obtained the following qualitative results. Monomeric hIAPP binds within the headgroup region and expands the membrane area of a lipid bilayer. At low concentrations, such binding causes no leakage or damage to the lipid bilayer. At high concentrations, the bound peptides transform to β-aggregates. The aggregates exit the headgroup region and bind to the surface of lipid bilayers. The damage by the surface bound β-aggregates depends on the aggregation size. The initial aggregation extracts lipid molecules, which probably causes ion permeation, but no molecular leakage. However, the initial β-aggregates serve as the seed for larger fibrils, in the manner of the Jarrett-Lansbury seeded-polymerization model, that eventually disintegrate lipid bilayers by electrostatic and hydrophobic interactions.     

Direct imaging of salt effects on lipid bilayer ordering at sub-molecular resolution

The interactions of salts with lipid bilayers are known to alter the properties of membranes and therefore influence their structure and dynamics. Sodium and calcium cations penetrate deeply into the headgroup region and bind to the lipids, whereas potassium ions only loosely associate with lipid molecules and mostly remain outside of the headgroup region. We investigated a dipalmitoylphosphatidylcholine (DPPC) bilayer in the gel phase in the presence of all three cations with a concentration of Ca²⁺ ions an order of magnitude smaller than the Na⁺ and K⁺ ions. Our findings indicate that the area per unit cell does not significantly change in these three salt solutions. However the lipid molecules do re-order non-isotropically under the influence of the three different cations. We attribute this reordering to a change in the highly directional intermolecular interactions caused by a variation in the dipole-dipole bonding arising from a tilt of the headgroup out of the membrane plane. Measurements in different NaCl concentrations also show a non-isotropic re-ordering of the lipid molecules.     

Molecular simulation study of structural and dynamic properties of mixed DPPC/DPPE bilayers

Molecular dynamics simulations have been used to study structural and dynamic properties of fully hydrated mixed 1,2-dipalmitoyl-sn-glycero-3-phosphocholine (DPPC) and 1,2-dipalmitoyl-sn-glycero-3-phosphoethanolamine (DPPE) bilayers at 0, 25, 50, 75, and 100 mol % DPPE. Simulations were performed for 50 ns at 350 K and 1 bar for the liquid-crystalline state of the mixtures. Results show that the average area per headgroup reduces from 0.65 +/- 0.01 nm(2) in pure DPPC to 0.52 +/- 0.01 nm(2) in pure DPPE systems. The lipid tails become more ordered with increasing DPPE concentration, resulting in a slight increase in membrane thickness (3.43 +/- 0.01 nm in pure DPPC to 4.00 +/- 0.01 nm in pure DPPE). The calculated area per headgroup and order parameter for pure DPPE deviates significantly from available experimental measurements, suggesting that the force field employed requires further refinement. In-depth analysis of the hydrogen-bond distribution in DPPE molecules shows that the amine groups strongly interact with the phosphate and carbonyl groups through inter/intramolecular hydrogen bonds. This yields a bilayer structure with DPPE headgroups preferentially located near the lipid phosphate and ester oxygens. It is observed that increasing DPPE concentrations causes competitive hydrogen bonding between the amine groups (hydrogen-donor) and the phosphate/carbonyl groups or water (hydrogen-acceptor). Due to the increasing number of hydrogen-donors from DPPE molecules with increasing concentration, DPPE becomes more hydrated. Trajectory analysis shows that DPPE molecules in the lipid mixtures move laterally and randomly around the membrane surface and the movement becomes more localized with increasing DPPE concentrations. For the conditions and simulation time considered, no aggregation or phase separation was observed between DPPC and DPPE.     

Interactions of liquid crystal-forming molecules with phospholipid bilayers studied by molecular dynamics simulations

Recent experiments have shown that liquid crystals can be used to image mammalian cell membranes and to amplify structural reorganization in phospholipid-laden liquid crystal-aqueous interfaces. In this work, molecular dynamics simulations were employed to explore the interactions between commonly used liquid crystal-forming molecules and phospholipid bilayers. In particular, umbrella sampling was used to obtain the potential of mean force of 4-cyano-4'-pentylbiphenyl (5CB) and 4'-(3,4-difluor-phenyl)-4-pentyl-bicylohexyl (5CF) molecules partitioning into a dipalmitoylphosphatidylcholine bilayer. In addition, results of simulations are presented for systems consisting of a fully hydrated bilayer with 5CB or 5CF molecules at the lowest (4.5 mol %) and highest (20 mol %) concentrations used in recent laboratory experiments. It is found that mesogens preferentially partition from the aqueous phase into the membrane; the potential of mean force exhibits highly favorable free energy differences for partitioning (-18 k(B)T for 5CB and -26 k(B)T for 5CF). The location and orientation of mesogens associated with the most stable free energies in umbrella sampling simulations of dilute systems were found to be consistent with those observed in liquid-crystal-rich bilayers. It is found that the presence of mesogens in the bilayer enhances the order of lipid acyl tails, and changes the spatial and orientational arrangement of lipid headgroup atoms. These effects are more pronounced at higher liquid-crystal concentrations. In comparing the behavior of 5CB and 5CF, a stronger spatial correlation (i.e., possibly leading to aggregation) is observed between 5CB molecules within a bilayer than between 5CF molecules. Also, the range of molecular orientations and positions along the bilayer normal is larger for 5CB molecules. At the same time, 5CF molecules were found to bind more strongly to lipid headgroups, thereby slowing the lateral motion of lipid molecules.     

Atomic view of calcium-induced clustering of phosphatidylserine in mixed lipid bilayers

Membranes play key regulatory roles in biological processes, with bilayer composition exerting marked effects on binding affinities and catalytic activities of a number of membrane-associated proteins. In particular, proteins involved in diverse processes such as vesicle fusion, intracellular signaling cascades, and blood coagulation interact specifically with anionic lipids such as phosphatidylserine (PS) in the presence of Ca(2+) ions. While Ca(2+) is suspected to induce PS clustering in mixed phospholipid bilayers, the detailed structural effects of this ion on anionic lipids are not established. In this study, combining magic angle spinning (MAS) solid-state NMR (SSNMR) measurements of isotopically labeled serine headgroups in mixed lipid bilayers with molecular dynamics (MD) simulations of PS lipid bilayers in the presence of different counterions, we provide site-resolved insights into the effects of Ca(2+) on the structure and dynamics of lipid bilayers. Ca(2+)-induced conformational changes of PS in mixed bilayers are observed in both liposomes and Nanodiscs, a nanoscale membrane mimetic of bilayer patches. Site-resolved multidimensional correlation SSNMR spectra of bilayers containing (13)C,(15)N-labeled PS demonstrate that Ca(2+) ions promote two major PS headgroup conformations, which are well resolved in two-dimensional (13)C-(13)C, (15)N-(13)C, and (31)P-(13)C spectra. The results of MD simulations performed on PS lipid bilayers in the presence or absence of Ca(2+) provide an atomic view of the conformational effects underlying the observed spectra.     

Solute effects on the colloidal and phase behavior of lipid bilayer membranes: ethanol-dipalmitoylphosphatidylcholine mixtures.

By means of the scanning differential calorimetry, x-ray diffractometry, and the dynamic light scattering, we have systematically studied the phase and packing properties of dipalmitoylphosphatidylcholine vesicles or multibilayers in the presence of ethanol. We have also determined the partial ternary phase diagram of such dipalmitoylphosphatidylcholine/water/ethanol mixtures. The directly measured variability of the structural bilayer parameters implies that ethanol binding to the phospholipid bilayers increases the lateral as well as the transverse repulsion between the lipid molecules. This enlarges the hydrocarbon tilt (by up to 23 degrees) and molecular area (by < or = 40%). Ethanol-phospholid association also broadens the interface and, thus, promotes lipid headgroup solvation. This results in excessive swelling (by 130%) of the phosphatidylcholine bilayers in aqueous ethanol solutions. Lateral bilayer expansion, moreover, provokes a successive interdigitation of the hydrocarbon chains in the systems with bulk ethanol concentrations of 0.4-1.2 M. The hydrocarbon packing density as well as the propensity for the formation of lamellar gel phases simultaneously increase. The pretransition temperature of phosphatidylcholine bilayers is more sensitive to the addition of alcohol (initial shift: delta Tp = 22 degrees C/mol) than the subtransition temperature (delta Ts reversible 5 degrees C/mol), whereas the chain-melting phase transition temperature is even less affected (delta Tm = 1.8 degrees C/mol). After an initial decrease of 3 degrees for the bulk ethanol concentrations below 1.2 M, the Tm value increases by 2.5 degrees above this limiting concentration. The gel-phase phosphatidylcholine membranes below Tm are fully interdigitated above this limiting concentration. The chain tilt on the fringe of full chain interdigitation is zero and increases with higher ethanol concentrations. Above Tm, some of the lipid molecules are solubilized by the bound ethanol molecules. More highly concentrated ethanol solutions (> 7 M) solubilize the phosphatidylcholine bilayers with fluid chains fully and result in the formation of mixed lipid-alcohol micelles.     

Reparameterization of all-atom dipalmitoylphosphatidylcholine lipid parameters enables simulation of fluid bilayers at zero tension

Molecular dynamics simulations of dipalmitoylphosphatidylcholine (DPPC) lipid bilayers using the CHARMM27 force field in the tensionless isothermal-isobaric (NPT) ensemble give highly ordered, gel-like bilayers with an area per lipid of approximately 48 A(2). To obtain fluid (L(alpha)) phase properties of DPPC bilayers represented by the CHARMM energy function in this ensemble, we reparameterized the atomic partial charges in the lipid headgroup and upper parts of the acyl chains. The new charges were determined from the electron structure using both the Mulliken method and the restricted electrostatic potential fitting method. We tested the derived charges in molecular dynamics simulations of a fully hydrated DPPC bilayer. Only the simulation with the new restricted electrostatic potential charges shows significant improvements compared with simulations using the original CHARMM27 force field resulting in an area per lipid of 60.4 +/- 0.1 A(2). Compared to the 48 A(2), the new value of 60.4 A(2) is in fair agreement with the experimental value of 64 A(2). In addition, the simulated order parameter profile and electron density profile are in satisfactory agreement with experimental data. Thus, the biologically more interesting fluid phase of DPPC bilayers can now be simulated in all-atom simulations in the NPT ensemble by employing our modified CHARMM27 force field.     

Phase behavior and nanoscale structure of phospholipid membranes incorporated with acylated C14-peptides

The thermotropic phase behavior and lateral structure of dipalmitoylphosphatidylcholine (DPPC) lipid bilayers containing an acylated peptide has been characterized by differential scanning calorimetry (DSC) on vesicles and atomic force microscopy (AFM) on mica-supported bilayers. The acylated peptide, which is a synthetic decapeptide N-terminally linked to a C14 acyl chain (C14-peptide), is incorporated into DPPC bilayers in amounts ranging from 0-20 mol %. The calorimetric scans of the two-component system demonstrate a distinct influence of the C14-peptide on the lipid bilayer thermodynamics. This is manifested as a concentration-dependent downshift of both the main phase transition and the pretransition. In addition, the main phase transition peak is significantly broadened, indicating phase coexistence. In the AFM imaging scans we found that the C14-peptide, when added to supported gel phase DPPC bilayers, inserts preferentially into preexisting defect regions and has a noticeable influence on the organization of the surrounding lipids. The presence of the C14-peptide gives rise to a laterally heterogeneous bilayer structure with coexisting lipid domains characterized by a 10 A height difference. The AFM images also show that the appearance of the ripple phase of the DPPC lipid bilayers is unaffected by the C14-peptide. The experimental results are supported by molecular dynamics simulations, which show that the C14-peptide has a disordering effect on the lipid acyl chains and causes a lateral expansion of the lipid bilayer. These effects are most pronounced for gel-like bilayer structures and support the observed downshift in the phase-transition temperature. Moreover, the molecular dynamics data indicate a tendency of a tryptophan residue in the peptide sequence to position itself in the bilayer headgroup region.     

1-Alkanols and membranes: A story of attraction

Although 1-alkanols have long been known to act as penetration enhancers and anesthetics, the mode of operation is not yet understood. In this study, long-time molecular dynamics simulations have been performed to investigate the effect of 1-alkanols of various carbon chain lengths onto the structure and dynamics of dimyristoylphosphatidylcholine bilayers. The simulations were complemented by microcalorimetry, continuous bleaching and film balance experiments. In the simulations, all investigated 1-alkanols assembled inside the lipid bilayer within tens of nanoseconds. Their hydroxyl groups bound preferentially to the lipid carbonyl group and the hydrocarbon chains stretched into the hydrophobic core of the bilayer. Both molecular dynamics simulations and experiments showed that all 1-alkanols drastically affected the bilayer properties. Insertion of long-chain 1-alkanols decreased the area per lipid while increasing the thickness of the bilayer and the order of the lipids. The bilayer elasticity was reduced and the diffusive motion of the lipids within the bilayer plane was suppressed. On the other hand, integration of ethanol into the bilayer enlarged the area per lipid. The bilayer became softer and lipid diffusion was enhanced.