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.     

Determining the Local Shear Viscosity of a Lipid Bilayer System by Reverse Non‐Equilibrium Molecular Dynamics Simulations

The parallel shear viscosity of a dipalmitoylphosphatidylcholine (DPPC) bilayer system is studied by reverse non‐equilibrium molecular dynamics simulations (RNEMD) with two different united‐atom force fields. The results are related to diffusion coefficients and structural distributions obtained by equilibrium molecular simulations. We investigate technical issues of the algorithm in the bilayer setup, namely, the dependence of the velocity profiles on the imposed flux and the influence of the thermostat on the calculated shear viscosity. We introduce the concept of local shear viscosity and investigate its dependence on the slip velocity of the monolayers and the particle density at the headgroup–water interface and the tail–tail interface. With this we demonstrate that the lipid bilayer is more viscous than the surrounding water phase, and that slip takes place near the headgroup region and in the centre of the bilayer where the alkyl tails meet. We also quantify the apparent increase in viscosity of the water molecules entangled at the water–headgroup interface.     

Effects of albumin and erythrocyte membranes on spread monolayers of lung surfactant lipids

Dipalmitoyl phosphatidylcholine (DPPC), one of the main constituents of lung surfactant is mainly responsible for reduction of surface tension to near 0 mN/m during expiration, resisting alveolar collapse. Other unsaturated phospholipids like palmitoyloleoyl phosphatidylglycerol (PG), palmitoyloleoyl phosphatidylcholine (POPC) and neutral lipids help in adsorption of lung surfactant to the air-aqueous interface. Lung surfactant lipids may interact with plasma proteins and hematological agents flooding the alveoli in diseased states. In this study, we evaluated the effects of albumin and erythrocyte membranes on spread films of DPPC alone and mixtures of DPPC with each of PG, POPC, palmitoyloleoyl phosphatidylethanolamine (PE), cholesterol (CHOL) and palmitic acid (PA) in 9:1 molar ratios. Surface tension-area isotherms were recorded using a Langmuir-Blodgett (LB) trough at 37 degrees C with 0.9% saline as the sub-phase. In the presence of erythrocyte membranes, DPPC and DPPC+PA monolayers reached minimum surface tensions of 7.3+/-0.9 and 9.6+/-1.4 mN/m, respectively. Other lipid combinations reached significantly higher minimum surface tensions >18 mN/m in presence of membranes (Newman Keul's test, p<0.05). The relative susceptibility to membrane inhibition was [(DPPC+PG, 7:3)=(DPPC+PG, 9:1)=(DPPC+POPC)=(DPPC+PE)=(DPPC+CHOL)]>[(DPPC+PA)=(DPPC)]. The differential response was more pronounced in case of albumin with DPPC and DPPC+PA monolayers reaching minimum surface tensions less than 2.4 mN/m in presence of albumin, whereas DPPC+PG and DPPC+POPC reached minimum surface tensions of around 20 mN/m in presence of albumin. Descending order of susceptibility of the spread monolayers of lipid mixtures to albumin destabilization was as follows: [(DPPC+PG, 7:3)=(DPPC+PG, 9:1)=(DPPC+POPC)]>[(DPPC+PE)=(DPPC+CHOL)]>[(DPPC+PA)=(DPPC)] The increase in minimum surface tension in presence of albumin and erythrocyte membranes was accompanied by sudden increases in compressibility at surface tensions of 15-30 mN/m. This suggests a monolayer destabilization and could be indicative of phase transitions in the mixed lipid films due to the presence of the hydrophobic constituents of erythrocyte membranes.     

Effects of lipid chain unsaturation and headgroup type on molecular interactions between paclitaxel and phospholipid within model biomembrane

Molecular interactions between paclitaxel, an anticancer drug, and phospholipids of various chain unsaturations and headgroup types were investigated in the present study by Langmuir film balance and differential scanning calorimetry. Both the lipid monolayer at the air-water interface and the lipid bilayer vesicles (liposomes) were employed as model cell membranes. It was found that, regardless of the difference in molecular structure of the lipid chains and headgroup, the drug can form nonideal, miscible systems with the lipids at the air-water interface over a wide range of paclitaxel mole fractions. The interaction between paclitaxel and phospholipid within the monolayer was dependent on the molecular area of the lipids at the interface and can be explained by intermolecular forces or geometric accommodation. Paclitaxel is more likely to form thermodynamically stable systems with 1,2-dipalmitoyl-sn-glycerol-3-phosphocholine (DPPC) and 1,2-dielaidoyl-sn-glycero-3-phosphocholine (DEPC) than with 1,2-dipalmitoyl-sn-glycero-3-phosphoethanolamine (DPPE) and 1,2-distearoyl-sn-glycero-3-phosphocholine (DSPC). Investigation of the drug penetration into the lipid monolayer showed that DPPC and DEPC have higher incorporation abilities for the drug than DPPE and DSPC. A similar trend was also evidenced by DSC investigation with liposomes. While little change of DSC profiles was observed for the DPPE/paclitaxel and DSPC/paclitaxel liposomes, paclitaxel caused noticeable changes in the thermographs of DPPC and DEPC liposomes. Paclitaxel was found to cause broadening of the main phase transition without significant change in the peak melting temperature of the DPPC bilayers, which demonstrates that paclitaxel was localized in the outer hydrophobic cooperative zone of the bilayer, i.e., in the region of the C1-C8 carbon atoms of the acyl chain or binding at the polar headgroup site of the lipids. However, it may penetrate into the deeper hydrophobic zone of the DEPC bilayers. These findings provide useful information for liposomal formulation of anticancer drugs as well as for understanding drug-cell membrane interactions.     

New protocols for preparing dipalmitoylphosphatidylcholine dispersions and controlling surface tension and competitive adsorption with albumin at the air/aqueous interface

The adsorption behavior of dipalmitoylphosphatidylcholine (DPPC), which is the major component of lung surfactant, at the air/aqueous interface and the competitive adsorption with bovine serum albumin (BSA) were studied with tensiometry, infrared reflection absorption spectroscopy (IRRAS), and ellipsometry. Dynamic surface tensions lower than 1 mN/m were observed for DPPC dispersions, with mostly vesicles, prepared with new protocols, involving extensive sonication above 50 °C. The lipid adsorbs faster and more extensively for DPPC dispersions with vesicles than with liposomes. For DPPC dispersions by a certain preparation procedure at T > T_c, when lipid particles were observed on the surface, dynamic surface tensions as low as 1 mN/m were measured. Moreover, IRRAS intensities and ellipsometric δΔ values were found to be much higher than the values for other DPPC dispersions or spread DPPC monolayers, suggesting that a larger amount of liposomes or vesicles adsorb on the surface. For DPPC/BSA mixtures, the tension behavior is controlled primarily by BSA, which prevents the formation of a dense DPPC monolayer. When BSA is injected into the subphase with a spread DPPC monolayer or into a DPPC dispersion with preadsorbed layers, little or no BSA adsorbs and the DPPC layer remains on the surface. When a DPPC monolayer is spread on a BSA solution at 0.1 wt% at 25 °C, then DPPC lipid can displace the adsorbed BSA molecules. The lack of BSA adsorption, and the expulsion of BSA by DPPC monolayer is probably due to the strong hydrophilicity of the lipid polar headgroup. When a DPPC dispersion is introduced with Trurnit's method or when dispersion drops are sprayed onto the surface of a DPPC/BSA mixture, the surface tension becomes lower and is controlled by DPPC, which can prevent the adsorption of BSA. The results may be important in understanding inhibition of lung surfactants by serum proteins and in designing efficient protocols of surfactant preparation and administration.     

Performance of the general amber force field in modeling aqueous POPC membrane bilayers

The aim of this work was to answer the question of whether the general amber force field (GAFF) is good enough to simulate fully hydrated POPC membrane bilayers. The test system contained 128 POPC and 2985 TIP3P water molecules. The equilibration was carried out in a nonarbitrary manner to reach the stable liquid-crystalline phase. The simulations were performed by using particle mesh Ewald electrostatics implemented in molecular dynamics packages Amber8 (NPT ensembles) and NAMD2 (NPgammaT ensembles). The computational results were assessed against the following experimental membrane properties: (i) area per lipid, (ii) area compressibility modulus, (iii) order parameter, (iv) gauche conformations per acyl chain, (v) lateral diffusion coefficients, (vi) electron density profile, and (vii) bound water at the lipid/water interface. The analyses revealed that the tested force field combination approximates the experimental values at an unexpectedly good level when the NPgammaT ensemble is applied with a surface tension of 60 mN m(-1) per bilayer. It is concluded that the GAFF/TIP3P combination can be utilized for aqueous membrane bilayer simulations, as it provides acceptable accuracy for biomolecular modeling.     

Equilibrium distributions of dipalmitoyl phosphatidylcholine and dilauroyl phosphatidylcholine in a mixed lipid bilayer: atomistic semigrand canonical ensemble simulations

Conventional molecular dynamics (MD) simulations are seriously limited by the slow rate of diffusive mixing in their ability to predict lateral distributions of different lipid types within mixed-lipid bilayers using atomistic models. A method to overcome this limitation, using configuration-bias Monte Carlo (MC) "mutation" moves to transform lipids from one type to another in dynamic equilibrium, is demonstrated in binary fluid-phase mixtures of lipids whose tails differ in length by four carbons. The hybrid MC-MD method operates within a semigrand canonical ensemble, so that an equilibrium composition of the mixture is determined by a constant difference in chemical potential (Delta(mu)) chosen for the components. Within several nanoseconds, bilayer structures initiated as pure dipalmitoyl phosphatidylcholine (DPPC) or pure dilauroyl phosphatidylcholine (DLPC) converge to a common composition and structure in independent simulations conducted at the same Delta(mu). Trends in bilayer thickness, area per lipid, density distributions across the bilayer, and order parameters have been investigated at three mixture compositions and compared with results from the pure bilayers at 323 K. The mixtures showed a moderate increase in DPPC acyl tail sites crossing the bilayer midplane relative to pure DPPC. Correlations between lateral positions of the two lipid types within or across the bilayer were found to be weak or absent. While the lateral distribution is consistent with nearly ideal mixing, the dependence of composition on Delta(mu) indicates a positive excess free energy of mixing.     

Molecular dynamics simulations of phospholipid bilayers: Influence of artificial periodicity, system size, and simulation time

This article investigates the convergence of structural and dynamical properties with system size and with time in molecular dynamics simulations of solvated phospholipid bilayers performed at constant volume under periodic boundary conditions using lattice-sum electrostatics. The electron density profile across the bilayer, the carbon-deuterium order parameters, and the surface tension are shown to be converged for a bilayer containing 36 lipids per leaflet and simulated over a period of 3-4 ns. Reasonable estimates for these properties can already be obtained from a system containing 16 lipids per leaflet. The convergence limit of 36 lipids per leaflet and the investigation of the correlation between lipid headgroup dipoles suggest a correlation length of about 3-5 nm in the lateral directions for a hydrated DPPC bilayer in the liquid-crystalline phase. Although these (relatively small) system sizes and (relatively short) time scales appear sufficient to obtain converged collective structural properties at constant volume, two restrictions should be kept in mind: (i) the relaxation times associated with the motion of individual lipids may be much longer and (ii) simulated properties converge significantly faster under constant volume conditions as compared to constant pressure conditions. Therefore, an accurate assessment of the dynamical properties of the system or of the relaxation of the bilayer under constant pressure conditions may require longer simulation time scales.     

Surface tension effect on transmembrane channel stability in a model membrane

The effect of surface tension on the lipid bilayer membrane is a question that has drawn considerable research effort. This interest has been driven both by the desire to determine the surface tension effects on the lipid bilayer and from the suggestion that adding finite surface tension to a small membrane system may provide more realistic lipid properties in molecular dynamics simulations. Here, the effect of surface tension on a palmitololelylphosphatidylcholine (POPC) bilayer membrane containing a four-helix transmembrane alamethicin peptide bundle is investigated. Simulations of 10 ns were undertaken for two different ensembles, NPT and NP(z)gammaT with a surface tension, gamma, of 20 mN m(-1) per interface, which is near the pore-forming region. The significance of differences between the tension-free and surface tension simulations was determined using nonparametric statistical analysis on replicate simulations with different initial conditions. The results suggest that, when the membrane is under surface tension, the peptide helical structure is perturbed from that in the tension-free state but that the bundle conformation is more stable than that in the tension-free state, with hydrogen bonding playing an important stabilizing role. Surface tension counteracts the influence of the transmembrane helix bundle on nearby lipid order, making the lipid order more uniform throughout the membrane in the tension state. Conversely, the lipid mobility was less uniform in the tension state, with lipids far from the bundle being significantly more mobile than those near the bundle. One general implication of the results is that surface tension can affect the membrane nonuniformly, in that the properties of lipids near the peptide are different from those further away.     

Differential effects of cholesterol, ergosterol and lanosterol on a dipalmitoyl phosphatidylcholine membrane: a molecular dynamics simulation study

Lipid raft/domain formation may arise as a result of the effects of specific sterols on the physical properties of membranes. Here, using molecular dynamics simulation, we examine the effects of three closely-related sterols, ergosterol, cholesterol, and lanosterol, at a biologically relevant concentration (40 mol %) on the structural properties of a model dipalmitoyl phosphatidylcholine (DPPC) membrane at 309 and 323 K. All three sterols are found to order the DPPC acyl tails and condense the membrane relative to the DPPC liquid-phase membrane, but each one does this to a significantly different degree. The smooth alpha-face of ergosterol, together with the presence of tail unsaturation in this sterol, leads to closer interaction of ergosterol with the lipids and closer packing of the lipids with each other, so ergosterol has a higher condensing effect on the membrane, as reflected by the area per lipid. Moreover, ergosterol induces a higher proportion of trans lipid conformers, a thicker membrane, and higher lipid order parameters and is aligned more closely with the membrane normal. Ergosterol also positions itself closer to the bilayer/water interface. In contrast, the rough alpha-face of lanosterol leads to a less close interaction of the steroid ring system with the phospholipid acyl chains, and so lanosterol orders, straightens, and packs the lipid acyl chains less well and is less closely aligned with the membrane normal. Furthermore, lanosterol lies closer to the relatively disordered membrane center than do the other sterols. The behavior of cholesterol in all the above respects is intermediate between that of lanosterol and ergosterol. The findings here may explain why ergosterol is the most efficient of the three sterols at promoting the liquid-ordered phase and lipid domain formation and may also furnish part of the explanation as to why cholesterol is evolutionarily preferred over lanosterol in higher-vertebrate plasma membranes.     

Effect of D-α-tocopheryl polyethylene glycol 1000 succinate (TPGS) on surfactant monolayers

In the present study, the effects of an amphiphilic polymer, d-α-tocopheryl polyethylene glycol 1000 succinate (TPGS) on model surfactant monolayers dipalmitoylphosphatidylcholine (DPPC), a binary mixture of DPPC with palmitoyloleoyl phosphatidylglycerol (DPPC-POPG) 9:1 (w/w) and binary mixture of DPPC and oleic acid (DPPC-OA) were evaluated. The ability of TPGS to act as an antioxidant adjuvant for pulmonary surfactants was also evaluated. Compression isotherms of surfactant monolayers at 37 °C in a Langmuir-Blodgett trough showed that DPPC and DPPC:TPGS mixed monolayers (1:0.25-1:1, w/w) exhibited low minimum surface tensions (MST) of 1-2 mN/m. Similarly [DPPC:POPG (9:1, w/w)]:TPGS mixed films of 1:0.25-1:1 weight ratios reached 1-2 mN/m MST. DPPC:POPG:TPGS liposomes adsorbed to surface tensions of 29-31 mN/m within 1s. While monolayers of DPPC:OA (1:1, w/w) reached high MST of ～11 mN/m, DPPC:OA:TPGS (1:1:0.25, w/w) film reached near zero MST suggesting that low concentrations of TPGS reverses the effect of OA on DPPC monolayer. Capillary surfactometer studies showed DPPC:TPGS and [DPPC:POPG (9:1, w/w)]:TPGS liposomes maintained 84-95% airway patency. Fluorescence spectroscopy of Laurdan loaded DPPC:TPGS and DPPC:POPG:TPGS liposomes revealed no segregation of lipid domains in the lipid bilayer. Addition of TPGS to soybean liposome significantly reduced thiobarbituric acid reactive substance (TBARS) by 29-39% confirming its antioxidant nature. The results suggest a potential use of TPGS as an adjuvant to improve the surfactant activity as well as act as an antioxidant by scavenging free radicals.     

Phase transition of a DPPC bilayer induced by an external surface pressure: from bilayer to monolayer behavior. a molecular dynamics simulation study

Understanding the lipid phase transition of lipid bilayers is of great interest from biophysical, physicochemical, and technological points of view. With the aim of elucidating the structural changes that take place in a DPPC phospholipid bilayer induced by an external isotropic surface pressure, five computer simulations were carried out in a range from 0.1 to 40 mN/m. Molecular dynamics simulations provided insight into the structural changes that took place in the lipid structure. It was seen that low pressures ranging from 0.1 to 1 mN/m had hardly any effect on the structure, electrical properties, or hydration of the lipid bilayer. However, for pressures above 40 mN/m, there was a sharp change in the lipid-lipid interactions, hydrocarbon lipid fluidity, and electrostatic potential, corresponding to the mesomorphic transition from a liquid crystalline state (L(alpha)) to its gel state (P'(beta)). The head lipid orientation remained almost unaltered, parallel to the lipid layer, as the surface pressure was increased, although a noticeable change in its angular distribution function was evident with the phase transition.     

Structural studies of the monolayers and bilayers formed by a novel cholesterol-phospholipid chimera

Langmuir isotherm, neutron reflectivity, and small angle neutron scattering studies have been conducted to characterize the monolayers and vesicular bilayers formed by a novel chimeric phospholipid, ChemPPC, that incorporates a cholesteryl moeity and a C-16 aliphatic chain, each covalently linked via a glycerol backbone to phosphatidylcholine. The structures of the ChemPPC monolayers and bilayers are compared against those formed from pure dipalmitoylphoshatidylcholine (DPPC) and those formed from a 60:40 mol % mixture of DPPC and cholesterol. In accord with previous findings showing that very similar macroscopic properties were exhibited by ChemPPC and 60:40 mol % DPPC/cholesterol vesicles, it is found here that the chimeric lipid and lipid/sterol mixture have very similar monolayer structures (each having a monolayer thickness of ～26 ?), and they also form vesicles with similar lamellar structure, each having a bilayer thickness of ～50 ? and exhibiting a repeat spacing of ～65 ?. The interfacial area of ChemPPC, however, is around 10 ?(2) greater than that of the combined DPPC/cholesterol unit in the mixed lipid monolayer (viz., 57 ± 1 vs 46 ± 1 ?(2), at 35 mN·m(-1)), and this difference in area is attributed to the succinyl linkage which joins the ChemPPC steroid and glyceryl moieties. The larger area of the ChemPPC is reflected in a slightly thicker monolayer solvent distribution width (9.5 vs 9 ? for the DPPC/cholesterol system) and by a marginal increase in the level of lipid headgroup hydration (16 vs 13 H(2)O per lipid, at 35 mN·m(-1)).     

Effect of acetone accumulation on structure and dynamics of lipid membranes studied by molecular dynamics simulations

The modulation of the properties and function of cell membranes by small volatile substances is important for many biomedical applications. Despite available experimental results, molecular mechanisms of action of inhalants and organic solvents, such as acetone, on lipid membranes remain not well understood. To gain a better understanding of how acetone interacts with membranes, we have performed a series of molecular dynamics (MD) simulations of a POPC bilayer in aqueous solution in the presence of acetone, whose concentration was varied from 2.8 to 11.2 mol%. The MD simulations of passive distribution of acetone between a bulk water phase and a lipid bilayer show that acetone favors partitioning into the water-free region of the bilayer, located near the carbonyl groups of the phospholipids and at the beginning of the hydrocarbon core of the lipid membrane. Using MD umbrella sampling, we found that the permeability barrier of ∼0.5 kcal/mol exists for acetone partitioning into the membrane. In addition, a Gibbs free energy profile of the acetone penetration across a bilayer demonstrates a favorable potential energy well of −3.6 kcal/mol, located at 15–16 Å from the bilayer center. The analysis of the structural and dynamics properties of the model membrane revealed that the POPC bilayer can tolerate the presence of acetone in the concentration range of 2.8–5.6 mol%. The accumulation of the higher acetone concentration of 11.2 mol% results, however, in drastic disordering of phospholipid packing and the increase in the membrane fluidity. The acetone molecules push the lipid heads apart and, hence, act as spacers in the headgroup region. This effect leads to the increase in the average headgroup area per molecule. In addition, the acyl tail region of the membrane also becomes less dense. We suggest, therefore, that the molecular mechanism of acetone action on the phospholipid bilayer has many common features with the effects of short chain alcohols, DMSO, and chloroform.     

Influence of ethanol on lipid membranes: from lateral pressure profiles to dynamics and partitioning

We have combined experiments with atomic-scale molecular dynamics simulations to consider the influence of ethanol on a variety of lipid membrane properties. We first employed isothermal titration calorimetry together with the solvent-null method to study the partitioning of ethanol molecules into saturated and unsaturated membrane systems. The results show that ethanol partitioning is considerably more favorable in unsaturated bilayers, which are characterized by their more disordered nature compared to their saturated counterparts. Simulation studies at varying ethanol concentrations propose that the partitioning of ethanol depends on its concentration, implying that the partitioning is a nonideal process. To gain further insight into the permeation of alcohols and their influence on lipid dynamics, we also employed molecular dynamics simulations to quantify kinetic events associated with the permeation of alcohols across a membrane, and to characterize the rotational and lateral diffusion of lipids and alcohols in these systems. The simulation results are in agreement with available experimental data and further show that alcohols have a small but non-vanishing effect on the dynamics of lipids in a membrane. The influence of ethanol on the lateral pressure profile of a lipid bilayer is found to be prominent: ethanol reduces the tension at the membrane-water interface and reduces the peaks in the lateral pressure profile close to the membrane-water interface. The changes in the lateral pressure profile are several hundred atmospheres. This supports the hypothesis that anesthetics may act by changing the lateral pressure profile exerted on proteins embedded in membranes.     

Cholesterol slows down the lateral mobility of an oxidized phospholipid in a supported lipid bilayer

We investigated the mobility and phase-partitioning of the fluorescent oxidized phospholipid analogue 1-palmitoyl-2-glutaroyl-sn-glycero-3-phospho-N-Alexa647-ethanolamine (PGPE-Alexa647) in supported lipid bilayers. Compared to the conventional phospholipid dihexadecanoylphosphoethanolamine (DHPE)-Bodipy we found consistently higher diffusion constants. The effect became dramatic when immobile obstacles were inserted into the bilayer, which essentially blocked the diffusion of DHPE-Bodipy but hardly influenced the movements of PGPE-Alexa647. In a supported lipid bilayer made of 1,2-dioleoyl-sn-glycero-3-phosphocholine (DOPC), the differences in probe mobility leveled off with increasing cholesterol content. Using coarse-grained molecular dynamics simulations, we could ascribe this effect to increased interactions between the oxidized phospholipid and the membrane matrix, concomitant with a translation in the headgroup position of the oxidized phospholipid: at zero cholesterol content, its headgroup is shifted to the outside of the DOPC headgroup region, whereas increasing cholesterol concentrations pulls the headgroup into the bilayer plane.     

Effects of antidepressants on the conformation of phospholipid headgroups studied by solid-state NMR

The effect of tricyclic antidepressants (TCA) on phospholipid bilayer structure and dynamics was studied to provide insight into the mechanism of TCA-induced intracellular accumulation of lipids (known as lipidosis). Specifically we asked if the lipid-TCA interaction was TCA or lipid specific and if such physical interactions could contribute to lipidosis. These interactions were probed in multilamellar vesicles and mechanically oriented bilayers of mixed phosphatidylcholine-phosphatidylglycerol (PC-PG) phospholipids using (31)P and (14)N solid-state NMR techniques. Changes in bilayer architecture in the presence of TCAs were observed to be dependent on the TCA's effective charge and steric constraints. The results further show that desipramine and imipramine evoke distinguishable changes on the membrane surface, particularly on the headgroup order, conformation and dynamics of phospholipids. Desipramine increases the disorder of the choline site at the phosphatidylcholine headgroup while leaving the conformation and dynamics of the phosphate region largely unchanged. Incorporation of imipramine changes both lipid headgroup conformation and dynamics. Our results suggest that a correlation between TCA-induced changes in bilayer architecture and the ability of these compounds to induce lipidosis is, however, not straightforward as imipramine was shown to induce more dramatic changes in bilayer conformation and dynamics than desipramine. The use of (14)N as a probe was instrumental in arriving at the presented conclusions.     

Single molecule probes of membrane structure: orientation of BODIPY probes in DPPC as a function of probe structure

Single molecule fluorescence measurements have recently been used to probe the orientation of fluorescent lipid analogs doped into lipid films at trace levels. Using defocused polarized total internal reflection fluorescence microscopy (PTIRF-M), these studies have shown that fluorophore orientation responds to changes in membrane surface pressure and composition, providing a molecular level marker of membrane structure. Here we extend those studies by characterizing the single molecule orientations of six related BODIPY probes doped into monolayers of DPPC. Langmuir-Blodgett monolayers transferred at various surface pressures are used to compare the response from fluorescent lipid analogs in which the location of the BODIPY probe is varied along the length of the acyl chain. For each BODIPY probe location along the chain, comparisons are made between analogs containing phosphocholine and smaller fatty acid headgroups. Together these studies show a general propensity of the BODIPY analogs to insert into membranes with the BODIPY probe aligned along the acyl chains or looped back to interact with the headgroups. For all BODIPY probes studied, a bimodal orientation distribution is observed which is sensitive to surface pressure, with the population of BODIPY probes aligned along the acyl chains increasing with elevated surface pressure. Trends in the single molecule orientations for the six analogs reveal a configuration where optimal placement of the BODIPY probe within the acyl chain maximizes its sensitivity to the surrounding membrane structure. These results are discussed in terms of balancing the effects of headgroup association with acyl chain length in designing the optimal placement of the BODIPY probe.     

Interface water dynamics and porating electric fields for phospholipid bilayers

Lipid bilayers, normally a barrier to charged species and large molecules, are permeabilized by electric fields, a phenomenon exploited by cell biologists and geneticists for porating and transfecting cells and tissues. Recent molecular simulation studies have advanced our understanding of electroporation, but the relative contributions of atomically local details (interface water and headgroup dipole and counterion configurations) and medium- and long-range electrostatic gradients and changes in membrane structural shifts to the initiating conditions and mechanisms of pore formation remain unclear. Molecular dynamics simulations of electroporation in several lipid systems presented here reveal the effects of lipid hydrocarbon tail length and composition on the magnitude of the field required for poration and on the location of the initial sites of field-driven water intrusion into the bilayer. Minimum porating external fields of 260 mV nm(-1), 280 mV nm(-1), 320 mV nm(-1), and 380 mV nm(-1) were found for 1,2-dilauroyl-sn-glycero-3-phosphatidylcholine (DLPC), 1,2-dipalmitoyl-sn-glycero-3-phosphatidylcholine (DPPC), 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphatidylcholine (POPC), and 1,2-dioleoyl- sn-glycero-3-phosphatidylcholine (DOPC) bilayers, respectively, and correlated most strongly with the bilayer thickness. These phospholipid systems share several common features including a wide, dynamic distribution of the headgroup dipole angle with the bilayer normal ranging from 0 to 155 degrees that is only slightly shifted in a porating electric field, and similar electric field-induced shifts in water dipole orientation, although the mean water dipole moment profile at the aqueous-membrane interface is more sensitive to the electric field for DOPC than for the other phospholipids. The location of pore initiation, at the anode- or cathode-facing leaflet, varies with the composition of the bilayer and correlates with a change in the polarity of the localized electric field at the interface.     

Phase separation of saturated and mono-unsaturated lipids as determined from a microscopic model

A molecular model is proposed of a bilayer consisting of fully saturated dipalmitoylphosphatidylcholine (DPPC) and mono-unsaturated dioleoylphosphatidylcholine (DOPC). The model not only encompasses the constant density within the hydrophobic core of the bilayer, but also the tendency of chain segments to align. It is solved within self-consistent field theory. A model bilayer of DPPC undergoes a main-chain transition to a gel phase, while a bilayer of DOPC does not do so above zero degrees centigrade because of the double bond which disrupts order. We examine structural and thermodynamic properties of these membranes and find our results in reasonable accord with experiment. In particular, order-parameter profiles are in good agreement with NMR experiments. A phase diagram is obtained for mixtures of these lipids in a membrane at zero tension. The system undergoes phase separation below the main-chain transition temperature of the saturated lipid. Extensions to the ternary DPPC, DOPC, and cholesterol system are outlined.