Influence of pyrene-labeling on fluid lipid membranes

We elucidate the influence of pyrene-labeled phospholipids on the structural properties of a fluid dipalmitoylphosphatidylcholine lipid membrane. To this end, we employ extensive atomic-scale molecular dynamics simulations with varying concentrations of pyrene-linked lipids. We find pyrene labeling to perturb the membrane structure significantly in the vicinity of the probe, the correlation length in the bilayer plane being about 1.0-1.5 nm. The local perturbations lead to enhanced ordering and packing of lipid acyl chains located in the vicinity of the probe. Surprisingly, this holds true not only for lipids that reside in the same leaflet as the pyrene-labeled probe but also for lipids in the opposite monolayer. The latter is due to substantial interdigitation of the pyrene moiety into the opposite leaflet, suggesting that occasional excimer formation may take place for probes in different leaflets. As a related issue, we also discuss the location and conformational orientation of the pyrene moieties. In particular, the orientational distribution of pyrene turns out to be more broad and diverse than the distribution of the corresponding acyl tails of nonlabeled lipids.     

Effect of double bond position on lipid bilayer properties: insight through atomistic simulations

Phosphatidylcholines (PCs) are among the most common phospholipids in plasma membranes. Their structural and dynamic properties are known to be strongly affected by unsaturation of lipid hydrocarbon chains, yet the role of the exact positions of the double bonds is poorly understood. In this work, we shed light on this matter through atomic-scale molecular dynamics simulations of eight different one-component lipid bilayers comprised of PCs with 18 carbons in their acyl chains. By introducing a single double bond in each acyl chain and varying its position in a systematic manner, we elucidate the effects of a double bond on various membrane properties. Studies in the fluid phase show that a number of membrane properties depend on the double bond position. In particular, when the double bond in an acyl chain is located close to the membrane-water interface, the area per lipid is considerably larger than that found for a saturated lipid. Further, when the double bond is shifted from the interfacial region toward membrane center, the area per lipid is observed to increase and have a maximum when the double bond is in the middle of the chain. Beyond this point, the surface area decreases systematically as the double bond approaches membrane center. These changes in area per lipid are accompanied by corresponding changes in membrane thickness and ordering of the chains. Further changes are observed in the tilt angles of the chains, membrane hydration together with changes in the number of gauche conformations, and direct head group interactions. All of these effects can be associated with changes in acyl chain conformations and local effects of the double bond on the packing of the surrounding atoms.     

Determination of the line tension of giant vesicles from pore-closing dynamics

Giant vesicles generated from synthetic and natural lipids such as phosphatidylcholines are useful models for understanding mechanical properties of cell membranes. Line tension is the one-dimensional force enabling the closing of transient pores on cell membranes. Transient pores were repeatedly and reproducibly formed on the membrane edge of giant vesicles generated from synthetic and natural phosphatidylcholines employing a nitrogen-pumped coumarin dye laser (440 nm). Line tension was determined at room temperature from closing of these pores that occurred over several seconds when the radius of the vesicle could be considered to be constant. The value of line tension depends on the nature of the lipid for single lipid systems, which, at room temperature, yielded a vesicle bilayer region in the gel, fluid, or mixed gel and fluid phases. The line tension for vesicles generated from phosphatidylcholines with saturated acyl chains of lengths of 12-18 carbon atoms ranges from 1 to 12 pN, exhibiting an increase with chain length. Vesicles generated from the natural Egg-PC, which is a mixture of lipids, are devoid of phase transition and exhibited the largest value of line tension (32 pN). This value is much larger than that estimated from the line tensions of vesicles obtained from lipids with homologous acyl chains. This study, to our knowledge, is the first to employ laser ablation to generate transient pores and determine line tension from the rate of pore closure and demonstrate a relationship between line tension and acyl chain length.     

Coarse-grain model for lipid bilayer self-assembly and dynamics: Multiparticle collision description of the solvent

A mesoscopic coarse-grain model for computationally efficient simulations of biomembranes is presented. It combines molecular dynamics simulations for the lipids, modeled as elastic chains of beads, with multiparticle collision dynamics for the solvent. Self-assembly of a membrane from a uniform mixture of lipids is observed. Simulations at different temperatures demonstrate that it reproduces the gel and liquid phases of lipid bilayers. Investigations of lipid diffusion in different phases reveals a crossover from subdiffusion to normal diffusion at long times. Macroscopic membrane properties, such as stretching and bending elastic moduli, are determined directly from the mesoscopic simulations. Velocity correlation functions for membrane flows are determined and analyzed.     

Structural properties of docosahexaenoyl phospholipid bilayers investigated by solid-state 2H NMR spectroscopy.

Polyunsaturated lipids in cellular membranes are known to play key roles in such diverse biological processes as vision, neuronal signaling, and apoptosis. One hypothesis is that polyunsaturated lipids are involved in second messenger functions in biological signaling. Another current hypothesis affirms that the functional role of polyunsaturated lipids relies on their ability to modulate physical properties of the lipid bilayer. The present research has employed solid-state 2H NMR spectroscopy to acquire knowledge of the molecular organization and material properties of polyunsaturated lipid bilayers. We report measurements for a homologous series of mixed-chain phosphatidylcholines containing a perdeuterated, saturated acyl chain (n:0) at the sn-1 position, adjacent to docosahexaenoic acid (DHA, 22:6omega3) at the sn-2 position. Measurements have been performed on fluid (L(alpha))-state multilamellar dispersions as a function of temperature for saturated acyl chain lengths of n = 12, 14, 16, and 18 carbons. The saturated sn-1 chains are therefore used as an intrinsic probe with site-specific resolution of the polyunsaturated bilayer structure. The 2H NMR order parameters as a function of acyl position (order profiles) have been analyzed using a mean-torque potential model for the chain segments, and the results are discussed in comparison with the homologous series of disaturated lipid bilayers. At a given absolute temperature, as the sn-1 acyl length adjacent to the sn-2 DHA chain is greater, the order of the initial chain segments increases, whereas that of the end segments decreases, in marked contrast with the corresponding disaturated series. For the latter, the order of the end segments is practically constant with acyl length, thus revealing a universal chain packing profile. We find that the DHA-containing series, while more complex, is still characterized by a universal chain packing profile, which is shifted relative to the homologous saturated series. Moreover, we show how introduction of DHA chains translates the order profile along the saturated chains, making more disordered states accessible within the bilayer central region. As a result, the area per lipid headgroup is increased as compared to disaturated bilayers. The systematic analysis of the 2H NMR data provides a basis for studies of lipid interactions with integral membrane proteins, for instance in relation to characteristic biological functions of highly unsaturated lipid membranes.     

Structure and thermotropic phase behavior of fluorinated phospholipid bilayers: a combined attenuated total reflection FTIR spectroscopy and imaging ellipsometry study

Lipid bilayers consisting of lipids with terminally perfluoroalkylated chains have remarkable properties. They exhibit increased stability and phase-separated nanoscale patterns in mixtures with nonfluorinated lipids. In order to understand the bilayer properties that are responsible for this behavior, we have analyzed the structure of solid-supported bilayers composed of 1,2-dipalmitoyl- sn-glycero-3-phosphocholine (DPPC) and of a DPPC analogue with 6 terminal perfluorinated methylene units (F6-DPPC). Polarized attenuated total reflection Fourier-transform infrared spectroscopy indicates that for F6-DPPC, the tilt of the lipid acyl chains to the bilayer normal is increased to 39 degrees as compared to 21 degrees for native DPPC, for both lipids in the gel phase. This substantial increase of the tilt angle is responsible for a decrease of the bilayer thickness from 5.4 nm for DPPC to 4.5 nm for F6-DPPC, as revealed by temperature-controlled imaging ellipsometry on microstructured lipid bilayers and solution atomic force microscopy. During the main phase transition from the gel to the fluid phase, both the relative bilayer thickness change and the relative area change are substantially smaller for F6-DPPC than for DPPC. In light of these structural and thermotropic data, we propose a model in which the higher acyl-chain tilt angle in F6-DPPC is the result of a conformational rearrangement to minimize unfavorable fluorocarbon-hydrocarbon interactions in the center of the bilayer due to chain staggering.     

Interleaflet interaction and asymmetry in phase separated lipid bilayers: molecular dynamics simulations

In order to investigate experimentally inaccessible, molecular-level detail regarding interleaflet interaction in membranes, we have run an extensive series of coarse-grained molecular dynamics simulations of phase separated lipid bilayers. The simulations are motivated by differences in lipid and cholesterol composition in the inner and outer leaflets of biological membranes. Over the past several years, this phenomenon has inspired a series of experiments in model membrane systems which have explored the effects of lipid compositional asymmetry in the two leaflets. The simulations are directed at understanding one potential consequence of compositional asymmetry, that being regions of bilayers where liquid-ordered (L(o)) domains in one leaflet are opposite liquid-disordered (L(d)) domains in the other leaflet (phase asymmetry). The simulated bilayers are of two sorts: 1) Compositionally symmetric leaflets where each of the two leaflets contains an identical, phase separated (L(o)/L(d)) mixture of cholesterol, saturated and unsaturated phospholipid; and 2) Compositionally asymmetric leaflets, where one leaflet contains a phase separated (L(o)/L(d)) mixture while the other contains only unsaturated lipid, which on its own would be in the L(d) phase. In addition, we have run simulations where the lengths of the saturated lipid chains as well as the mole ratios of the three lipid components are varied. Collectively, we report on three types of interleaflet coupling within a bilayer. First, we show the effects of compositional asymmetry on acyl chain tilt and order, lipid rotational dynamics, and lateral diffusion in regions of leaflets that are opposite L(o) domains. Second, we show substantial effects of compositional asymmetry on local bilayer curvature, with the conclusion that phase separated leaflets resist curvature, while inducing large degrees of curvature in an opposing L(d) leaflet. Finally, in compositionally symmetric, phase separated bilayers, we find phase asymmetry (domain antiregistration) between the two leaflets occurs as a consequence of mismatched acyl chain-lengths in the saturated and unsaturated lipids.     

Determination of lipid phase transition temperatures in hybrid bilayer membranes

The main gel-to-liquid-crystal (LC) phase transition temperature, T(m), of the lipid monolayer in hybrid bilayer membranes (HBMs) was investigated using vibrational sum frequency spectroscopy (VSFS). In the gel phase, the acyl chains of the lipid molecules assume an ordered, all-trans configuration, whereas in the LC phase, the acyl chains exhibit a significant number of disordered gauche conformers. VSFS has unique sensitivity to the order/disorder transitions in the acyl chains and was used to determine T(m) for a series of saturated phosphatidylcholine lipids on octadecanethiolate self-assembled monolayers (SAMs). The values obtained for T(m) for all lipids studied are significantly higher than for the corresponding lipids in vesicles in solution. Additionally, the transition widths are broader for the lipids in HBMs. The underlying SAM clearly influences the phase behavior of the overlying lipid monolayer.     

CHARMM‐GUI Membrane Builder toward realistic biological membrane simulations

CHARMM‐GUI Membrane Builder, http://www.charmm‐gui.org/input/membrane , is a web‐based user interface designed to interactively build all‐atom protein/membrane or membrane‐only systems for molecular dynamics simulations through an automated optimized process. In this work, we describe the new features and major improvements in Membrane Builder that allow users to robustly build realistic biological membrane systems, including (1) addition of new lipid types, such as phosphoinositides, cardiolipin (CL), sphingolipids, bacterial lipids, and ergosterol, yielding more than 180 lipid types, (2) enhanced building procedure for lipid packing around protein, (3) reliable algorithm to detect lipid tail penetration to ring structures and protein surface, (4) distance‐based algorithm for faster initial ion displacement, (5) CHARMM inputs for P₂₁ image transformation, and (6) NAMD equilibration and production inputs. The robustness of these new features is illustrated by building and simulating a membrane model of the polar and septal regions of E. coli membrane, which contains five lipid types: CL lipids with two types of acyl chains and phosphatidylethanolamine lipids with three types of acyl chains. It is our hope that CHARMM‐GUI Membrane Builder becomes a useful tool for simulation studies to better understand the structure and dynamics of proteins and lipids in realistic biological membrane environments. © 2014 Wiley Periodicals, Inc.     

Atomistic simulation of lipid and DiI dynamics in membrane bilayers under tension

Membrane tension modulates cellular processes by initiating changes in the dynamics of its molecular constituents. To quantify the precise relationship between tension, structural properties of the membrane, and the dynamics of lipids and a lipophilic reporter dye, we performed atomistic molecular dynamics (MD) simulations of DiI-labeled dipalmitoylphosphatidylcholine (DPPC) lipid bilayers under physiological lateral tensions ranging from -2.6 mN m(-1) to 15.9 mN m(-1). Simulations showed that the bilayer thickness decreased linearly with tension consistent with volume-incompressibility, and this thinning was facilitated by a significant increase in acyl chain interdigitation at the bilayer midplane and spreading of the acyl chains. Tension caused a significant drop in the bilayer's peak electrostatic potential, which correlated with the strong reordering of water and lipid dipoles. For the low tension regime, the DPPC lateral diffusion coefficient increased with increasing tension in accordance with free-area theory. For larger tensions, free area theory broke down due to tension-induced changes in molecular shape and friction. Simulated DiI rotational and lateral diffusion coefficients were lower than those of DPPC but increased with tension in a manner similar to DPPC. Direct correlation of membrane order and viscosity near the DiI chromophore, which was just under the DPPC headgroup, indicated that measured DiI fluorescence lifetime, which is reported to decrease with decreasing lipid order, is likely to be a good reporter of tension-induced decreases in lipid headgroup viscosity. Together, these results offer new molecular-level insights into membrane tension-related mechanotransduction and into the utility of DiI in characterizing tension-induced changes in lipid packing.     

Perturbation of phospholipid bilayer properties by ethanol at a high concentration

Atomistic molecular dynamics (MD) simulations have been carried out at 30 degrees C on a fully hydrated liquid crystalline lamellar phase of dimyrystoylphosphatidylcholine (DMPC) lipid bilayer with embedded ethanol molecules at 1:1 composition, as well as on the pure bilayer phase. The ethanol molecules are found to exhibit a preference to occupy regions near the upper part of the lipid acyl chains and the phosphocholine headgroups. The calculations revealed that the phosphocholine headgroup dipoles (P- --> N+) of the lipids prefer to orient more toward the aqueous layer in the presence of ethanol. It is noticed that the ethanol molecules modify the dynamic properties of both lipids as well as the water molecules in the hydration layer of the lipid headgroups. Both the in-plane "rattling" and out-of-plane "protrusion" motions of the lipids have been found to increase in the presence of ethanol. Most importantly, it is observed that the water molecules within the hydration layer of the lipid headgroups exhibit faster translational and rotational motions in the presence of ethanol. This arises due to faster dynamics of hydrogen bonds between lipid headgroups and water in the presence of ethanol.     

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.     

Reconstitution of rhodopsin into polymerizable planar supported lipid bilayers: influence of dienoyl monomer structure on photoactivation

G-protein-coupled receptors (GPCRs) play key roles in cellular signal transduction and many are pharmacologically important targets for drug discovery. GPCRs can be reconstituted in planar supported lipid bilayers (PSLBs) with retention of activity, which has led to development of GPCR-based biosensors and biochips. However, PSLBs composed of natural lipids lack the high stability desired for many technological applications. One strategy is to use synthetic lipid monomers that can be polymerized to form robust bilayers. A key question is how lipid polymerization affects GPCR structure and activity. Here we have investigated the photochemical activity of bovine rhodopsin (Rho), a model GPCR, reconstituted into PSLBs composed of lipids having one or two polymerizable dienoyl moieties located in different regions of the acyl chains. Plasmon waveguide resonance spectroscopy was used to compare the degree of Rho photoactivation in fluid and poly(lipid) PSLBs. The position of the dienoyl moiety was found to have a significant effect: polymerization near the glycerol backbone significantly attenuates Rho activity whereas polymerization near the acyl chain termini does not. Differences in cross-link density near the acyl chain termini also do not affect Rho activity. In unpolymerized PSLBs, an equimolar mixture of phosphatidylethanolamine and phosphatidylcholine (PC) lipids enhances activity relative to pure PC; however after polymerization, the enhancement is eliminated which is attributed to stabilization of the membrane lamellar phase. These results should provide guidance for the design of robust lipid bilayers functionalized with transmembrane proteins for use in membrane-based biochips and biosensors.     

Modulating membrane properties: the effect of trehalose and cholesterol on a phospholipid bilayer

The protective properties of trehalose on cholesterol-containing lipid dipalmitoylphosphatidylcholine (DPPC) bilayers are studied through molecular simulations. The ability of the disaccharide to interact with the phospholipid headgroups and stabilize the membrane persists even at high cholesterol concentrations and restricts some of the changes to the structure that would otherwise be imposed by cholesterol molecules. Predictions of bilayer properties such as area per lipid, tail ordering, and chain conformation support the notion that the disaccharide decreases the main melting transition in these multicomponent model membranes, which correspond more closely to common biological systems than pure bilayers. Molecular simulations indicate that the membrane dynamics are slowed considerably by the presence of trehalose, indicating that high sugar concentrations would serve to avert possible phase separations that could arise in mixed phospholipid systems. Various time correlation functions suggest that the character of the modifications in lipid dynamics induced by trehalose and cholesterol is different in the hydrophilic and hydrophobic regions of the membrane.     

Molecular dynamics simulations of cardiolipin bilayers

Cardiolipin is a key lipid component in the inner mitochondrial membrane, where the lipid is involved in energy production, cristae structure, and mechanisms in the apoptotic pathway. In this article we used molecular dynamics computer simulations to investigate cardiolipin and its effect on the structure of lipid bilayers. Three cardiolipin/POPC bilayers with different lipid compositions were simulated: 100, 9.2, and 0% cardiolipin. We found strong association of sodium counterions to the carbonyl groups of both lipid types, leaving in the case of 9.2% cardiolipin virtually no ions in the aqueous compartment. Although binding occurred primarily at the carbonyl position, there was a preference to bind to the carbonyl groups of cardiolipin. Ion binding and the small headgroup of cardiolipin gave a strong ordering of the hydrocarbon chains. We found significant effects in the water dipole orientation and water dipole potential which can compensate for the electrostatic repulsion that otherwise should force charged lipids apart. Several parameters relevant for the molecular structure of cardiolipin were calculated and compared with results from analyses of coarse-grained simulations and available X-ray structural data.     

Implicit-solvent mesoscale model based on soft-core potentials for self-assembled lipid membranes

An efficient implicit-solvent model for self-assembled lipid bilayers is presented and analyzed using Langevin molecular dynamics simulations. The model is based on soft interactions between particles and short-range attractive interaction between lipid tails, leading for the self-assembly of a lipid bilayer without an explicit solvent. This allows for efficient simulations of large membranes over long times. The model exhibits a fluid phase at high temperatures and a gel phase at low temperatures, identified with the Lbeta-phase. The melting transition is investigated via analysis of the diffusivity of the lipid molecules, the chain-orientational order parameter, the sixfold bond-orientational order parameter, and the positional and bond-orientational correlation functions. The analysis suggests the existence of a hexatic phase over a narrow range of temperatures around the melting transition. The elastic properties of the membrane in the fluid phase are also investigated.     

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.     

Lipid lateral mobility and membrane phase structure modulation by protein binding

Using a combination of fluorescence correlation and infrared absorption spectroscopies, we characterize lipid lateral diffusion and membrane phase structure as a function of protein binding to the membrane surface. In a supported membrane configuration, cholera toxin binding to the pentasaccharaide headgroup of membrane-incorporated GM1 lipid alters the long-range lateral diffusion of fluorescently labeled probe lipids, which are not involved in the binding interaction. This effect is prominently amplified near the gel-fluid transition temperature, Tm, of the majority lipid component. At temperatures near Tm, large changes in probe lipid diffusion are measured at average protein coverage densities as low as 0.02 area fraction. Spectral shifts of the methylene symmetric and asymmetric stretching modes in the lipid acyl chain confirm that protein binding alters the fraction of lipid in the gel phase.     

Influence of polymyxins on the structural dynamics of Escherichia coli lipid membranes

Polymyxins are a family of nonribosomic cationic peptide antibiotics highly effective against Gram-negative bacteria. Two members of this family, Polymyxins B and E (PxB, PxE), form molecular vesicle-vesicle contacts and promote a selective exchange of phospholipids at very low concentrations in the membrane, a biophysical phenomenon that can be the basis of their antibiotic mode of action. To get more insight into the interaction of these antibiotics with the lipid membrane, their effect on the structural dynamics of bilayers prepared with lipids extracted from the membrane of Escherichia coli was determined using fluorescently labeled phopholipids. Steady-state anisotropy measurements with probes that localize at different positions in the membrane give information on the effects of polymyxins on the mobility of the phospholipids. Results with PxB, PxE, colymycin M and polymyxin B nonapeptide (PxB-NP), a deacylated derivative with no antibiotic properties, are compared. At low peptide concentrations (<2 mol%) PxB and PxE bind to the membranes superficially, affecting very slightly the ordering of the lipids at the outermost part of the bilayer. Above this concentration, PxB and PxE insert more deeply in the bilayer, increasing lipid order both in the gel and liquid-crystal states and modifying phase transitions. Fluorescence experiments with pyrene labeled phospholipids indicate that the increase in lipid packing is accompanied by an enrichment of phospholipids in the bilayers. In contrast, colymycin M and PxB-NP did not modify lipid packing or phase transition, nor did they induce microdomain formation. The possible significance of these results in the antibiotic mode of action of PxB and PxE is discussed. The combination of spectroscopic techniques described here can be useful as part of a general method of screening for new antibiotics that act on the membrane by the same mechanism as polymyxins.