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.     

Phosphatidylethanolamine-induced cholesterol domains chemically identified with mass spectrometric imaging

Coexisting liquid phases of model membrane systems are chemically identified using imaging time-of-flight secondary ion mass spectrometry (TOF-SIMS). The systems studied were Langmuir-Blodgett (LB) model membranes of cholesterol (CH) with two different phospholipids, one a major component in the outer plasma membrane bilayer leaflet (dipalmitoylphosphatidylcholine (PC)) and the other a major component in the inner leaflet (dipalmitoylphosphatidylethanolamine (PE)). Binary mixtures of CH with each of the phospholipids were investigated, as well as a ternary system. A single homogeneous phase is evident for PC/CH, whereas both systems containing PE show lateral heterogeneity with phospholipid-rich and CH-rich regions. The interaction between CH and the two phospholipids differs due to the disparity between the phospholipid headgroups. Imaging TOF-SIMS offers a novel opportunity to chemically identify and differentiate the specific membrane locations of CH and phospholipid in membrane regions without the use of fluorescent dyes. This unique imaging method has been used to demonstrate the formation of micrometer-size CH domains in phosphatidylethanolamine-rich systems and is further evidence suggesting that CH may facilitate transport and signaling across the two leaflets of the plasma membrane.     

Asymmetry of lipid bilayers induced by monovalent salt: atomistic molecular-dynamics study

Interactions between salt ions and lipid components of biological membranes are essential for the structure, stability, and functions of the membranes. The specific ionic composition of aqueous buffers inside and outside of the cell is known to differ considerably. To model such a situation we perform atomistic molecular-dynamics (MD) simulations of a single-component phosphatidylcholine lipid bilayer which separates two aqueous reservoirs with and without NaCl salt. To implement the difference in electrolyte composition near two membrane sides, a double bilayer setup (i.e., two bilayers in a simulation box) is employed. It turns out that monovalent salt, being in contact with one leaflet only, induces a pronounced asymmetry in the structural, electrostatic, and dynamical properties of bilayer leaflets after 50 ns of MD simulations. Binding of sodium ions to the carbonyl region of the leaflet which is in contact with salt results in the formation of "Na-lipids" complexes and, correspondingly, reduces mobility of lipids of this leaflet. In turn, attractive interactions of chloride ions (mainly located in the aqueous phase close to the water-lipid interface) with choline lipid groups lead to a substantial (more vertical) reorientation of postphatidylcholine headgroups of the leaflet adjoined to salt. The difference in headgroup orientation on two sides of a bilayer, being coupled with salt-induced reorientation of water dipoles, leads to a notable asymmetry in the charge-density profiles and electrostatic potentials of bilayer constitutes of the two leaflets. Although the overall charge density of the bilayer is found to be almost insensitive to the presence of salt, a slight asymmetry in the charge distribution between the two bilayer leaflets results in a nonzero potential difference of about 85 mV between the two water phases. Thus, a transmembrane potential of the order of the membrane potential in a cell can arise without ionic charge imbalance between two aqueous compartments.     

Electrostatic Localization of RNA to Protocell Membranes by Cationic Hydrophobic Peptides

Cooperative interactions between RNA and vesicle membranes on the prebiotic earth may have led to the emergence of primitive cells. The membrane surface offers a potential platform for the catalysis of reactions involving RNA, but this scenario relies upon the existence of a simple mechanism by which RNA could become associated with protocell membranes. Here, we show that electrostatic interactions provided by short, basic, amphipathic peptides can be harnessed to drive RNA binding to both zwitterionic phospholipid and anionic fatty acid membranes. We show that the association of cationic molecules with phospholipid vesicles can enhance the local positive charge on a membrane and attract RNA polynucleotides. This phenomenon can be reproduced with amphipathic peptides as short as three amino acids. Finally, we show that peptides can cross bilayer membranes to localize encapsulated RNA. This mechanism of polynucleotide confinement could have been important for primitive cellular evolution.     

Structural and dynamic characterization of the interaction of the putative fusion peptide of the S2 SARS-CoV virus protein with lipid membranes

The SARS coronavirus (SARS-CoV) envelope spike (S) glycoprotein, a Class I viral fusion protein, is responsible for the fusion between the membranes of the virus and the target cell. In the present work, we report a study of the binding and interaction with model membranes of a peptide pertaining to the putative fusion domain of SARS-CoV, SARS FP, as well as the structural changes that take place in both the phospholipid and the peptide molecules upon this interaction. From fluorescence and infrared spectroscopies, the peptide ability to induce membrane leakage, aggregation and fusion, as well as its affinity toward specific phospholipids, was assessed. We demonstrate that SARS FP strongly partitions into phospholipid membranes, more specifically with those containing negatively charged phospholipids, increasing the water penetration depth and displaying membrane-activity modulated by the lipid composition of the membrane. Interestingly, peptide organization is different depending if SARS FP is in water or bound to the membrane. These data suggest that SARS FP could be involved in the merging of the viral and target cell membranes by perturbing the membrane outer leaflet phospholipids and specifically interacting with negatively charged phospholipids located in the inner leaflet.     

Spontaneous formation of asymmetric lipid bilayers by adsorption of vesicles

We have investigated the adsorption of phospholipid mixtures using neutron reflection. Small sonicated unilamellar vesicles (SUV) composed of DOPC and d(62)-DPPC were incubated at 50 degrees C in contact with a silica surface using a method commonly employed to form supported model membranes. The composition of the mixed supported bilayer was found to be substantially different from that of the bulk vesicles in a direction indicating a higher affinity of DPPC for the silica surface. Formation of an asymmetric bilayer arrangement was also discovered in all the cases studied. DPPC tended to dominate the composition of the leaflet next to silica, while the outer leaflet was generally closer to the bulk composition. The supported bilayers also exhibited increasing interfacial roughness in the outer membrane leaflet in the region of the DOPC-DPPC gel-liquid immiscibility region. To our knowledge, this is the first time that both the structure and the absolute composition of a mixed-lipid supported bilayer have been resolved, and the results raise a number of questions regarding the adsorption of vesicles and the properties of supported bilayers, which are discussed in terms of the bulk phase diagram of DOPC and DPPC.     

制菌霉素与固体支撑磷脂膜的相互作用

为了更好地了解制菌霉素的作用机理, 本文利用表面等离子体共振(SPR)和交流阻抗两种技术, 考察了制菌霉素与不含固醇的固体支撑纯磷脂膜的相互作用, 结果发现, 制菌霉素可与纯磷脂膜相互作用, 并可能在膜上形成微孔.     

Asymmetric droplet interface bilayers

In cell membranes, the lipid compositions of the inner and outer leaflets differ. Therefore, a robust model system that enables single-channel electrical recording with asymmetric bilayers would be very useful. We and others recently developed the droplet interface bilayer (DIB), which is formed by connecting lipid monolayer-encased aqueous droplets submerged in an oil-lipid mixture. Here, we incorporate lipid vesicles of different compositions into aqueous droplets and immerse them in an oil bath to form asymmetric DIBs (a-DIBs). Both alpha-helical and beta-barrel membrane proteins insert readily into a-DIBs, and their activity can be measured by single-channel electrical recording. We show that the gating behavior of outer membrane protein G (OmpG) from Escherichia coli differs depending on the side of insertion in an asymmetric DIB with a positively charged leaflet opposing a negatively charged leaflet. The a-DIB system provides a general platform for studying the effects of bilayer leaflet composition on the behavior of ion channels and pores.     

A class of supported membranes: formation of fluid phospholipid bilayers on photonic band gap colloidal crystals

We report the formation of a new class of supported membranes consisting of a fluid phospholipid bilayer coupled directly to a broadly tunable colloidal crystal with a well-defined photonic band gap. For nanoscale colloidal crystals exhibiting a band gap at the optical frequencies, substrate-induced vesicle fusion gives rise to a surface bilayer riding onto the crystal surface. The bilayer is two-dimensionally continuous, spanning multiple beads with lateral mobilities which reflect the coupling between the bilayer topography and the curvature of the supporting colloidal surface. In contrast, the spreading of vesicles on micrometer scale colloidal crystals results in the formation of bilayers wrapping individual colloidal beads. We show that simple UV photolithography of colloidal crystals produces binary patterns of crystal wettabilities, photonic stopbands, and corresponding patterns of lipid mono- and bilayer morphologies. We envisage that these approaches will be exploitable for the development of optical transduction assays and microarrays for many membrane-mediated processes, including transport and receptor-ligand interactions.     

Molecular level investigation of organization in ternary lipid bilayer: a computational approach

The differential organization of lipid components in a multicomponent membrane leads to formation of domains having diverse composition and size. Cholesterol and glycosphingolipids are known to be important components of such lateral assembly. We report here the ordering of cholesterol around ganglioside GM1 and the nature of the cluster from an all-atom simulation of a ternary lipid system. The results are compared with a binary bilayer and a pure phospholipid bilayer. The difference in molecular rearrangements in ternary and binary lipid mixture shows the role of GM1 in the rearrangement of cholesterol. Calculation of the radial distribution function, rotational reorientation, and residence time analysis of cholesterol shows that cholesterol is preferentially accumulating near gangliosides, while the lateral translational motion, rotational diffusion, and order parameter of phospholipids characterize the amount of rigidity imparted on the phospholipid bilayer.     

Transbilayer complementarity of phospholipids. A look beyond the fluid mosaic model

Lipid-lipid interactions across a phospholipid bilayer were probed by measuring the nearest-neighbor preferences of exchangeable phospholipid monomers derived from 1,2-dimyristoyl-sn-glycero-3-phosphoethanolamine (DMPE) and 1,2-distearoyl-sn-glycero-3-phosphoethanolamine (DSPE) in the presence of nonexchangeable DMPE- or DSPE-based dimers. Each of these permanent dimers promoted homophospholipid association to the same extent, whereas the corresponding nonexchangeable monomers were without effect. These results support a model in which the longer phospholipids in one monolayer leaflet preferentially associate with shorter ones in the adjoining monolayer. Such transbilayer complementarity is likely to play an important role in stabilizing biological membranes and also in promoting a compositional interdependence of their two lipid leaflets.     

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.     

Influence of organotin compounds on phosphatidylserine membranes

Organotin compounds are widely distributed toxicants. They are membrane‐active molecules with broad biological toxicity. We have studied the interaction of tributyltin and triphenyltin with phosphatidylserine model membranes using differential scanning calorimetry, infrared spectroscopy and X‐ray diffraction techniques. Organotin compounds produced a broadening of the gel to the liquid‐crystalline phase transition of the phospholipid and a shifting of the phase transition temperature to lower values. Infrared spectroscopy experiments showed that tributyltin exerted a fluidizing effect on the apolar part of the bilayer, and that both tributyl‐ and triphenyltin interact with the interfacial region of the bilayer, making the carbonyl groups less accessible to water. As seen by X‐ray diffraction experiments, organotin compounds were unable to change the bilayer macroscopic organization of the phospholipid, but they were able to reduce the long‐range order of the multibilayer system and to disorder the packing of the phospholipid molecules. The observed interaction between organotin compounds and phosphatidylserine membranes promotes physical perturbations that could affect membrane function and may mediate some of their toxic effects. Copyright © 2004 John Wiley & Sons, Ltd.     

Bias-dependent admittance in hybrid bilayer membranes

Artificial bilayer membranes provide a platform for bioelectronic devices based on their structural, sensing, and transport functions. In this letter, we report on the impedance response of an engineered membrane with a lower leaflet of octadecanethiol on gold and an outer leaflet of dioleoylphosphatidylcholine with the monomeric channel protein gramicidin. This hybrid bilayer exhibits an electrical response analogous to a solid-state diode: the admittance is very low (<10(-)(7) Omega(-)(1) cm(-)(2)) over a wide potential range but increases exponentially at negative potentials.     

Interrogating interfacial organization in planar bilayer structures

The field of biomimetic planar lipid membranes is finding increased importance as the need to devise sensing systems for biologically important species increases. We approach this area with an eye toward understanding how to interrogate local organization in these complex media. Our primary tools for this purpose are spectroscopy and electrochemistry, where imbedded reporter molecules serve as the information transducers. We use Langmuir-Blodgett (LB) and Langmuir-Schaefer (LS) methods to construct planar lipid membranes on hydrophilic solid substrates (Au for electrochemistry, SiO x for spectroscopy). Pyrene tethered to the substrate acts as our probe and AC voltammetry was used to evaluate its redox kinetics, showing slow, distance independent electron transfer between the pyrene moieties and the electrode for both monolayer and bilayer systems. Time-resolved fluorescence data indicate that tethered pyrene resides in a highly rigid environment and that the addition of the top lipid leaflet improves the organization of the bottom lipid leaflet. These data point to the cooperative effect of the bilayer leaflets in creating a system that is comparatively rigid on short length scales, and capable of mediating motion and accessibility of imbedded species.     

Exploiting Benzophenone Photoreactivity To Probe the Phospholipid Environment and Insertion Depth of the Cell‐Penetrating Peptide Penetratin in Model Membranes

Penetratin (RQIKIWFQNRRMKWKK) enters cells by different mechanisms, including membrane translocation, thus implying that the peptide interacts with the lipid bilayer. Penetratin also crosses the membrane of artificial vesicles, depending on their phospholipid content. To evaluate the phospholipid preference of penetratin, as the first step of translocation, we exploited the benzophenone triplet kinetics of hydrogen abstraction, which is slower for secondary than for allylic hydrogen atoms. By using multilamellar vesicles of varying phospholipid content, we identified and characterized the cross‐linked products by MALDI‐TOF mass spectrometry. Penetratin showed a preference for negatively charged (vs. zwitterionic) polar heads, and for unsaturated (vs. saturated) and short (vs. long) saturated phospholipids. Our study highlights the potential of using benzophenone to probe the environment and insertion depth of membranotropic peptides in membranes.     

Layer-by-layer PMIRRAS characterization of DMPC bilayers deposited on a Au111 electrode surface

A combination of Langmuir-Blodgett and Langmuir-Schaefer techniques was employed to deposit 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC) bilayers at a gold electrode surface. One leaflet consisted of hydrogen-substituted acyl chains, and the second leaflet was composed of molecules with deuterium-substituted acyl chains. This architecture allowed for layer-by-layer analysis of the structure of the bilayer. Photon polarization modulation infrared reflection absorption spectroscopy (PM-IRRAS) was used to determine the conformation and orientation of the acyl chains of DMPC molecules in the individual leaflets as a function of the potential applied to the gold electrode. The bilayer is adsorbed onto the metal surface when the field applied to the membrane does not exceed approximately 108 V/m. When adsorbed, the bottom leaflet is in contact with a hydrophobic metal surface, and the top leaflet is interacting with the aqueous solution. The asymmetry of the environment has an effect on the orientation of the DMPC molecules in each leaflet. The tilt angle of the acyl chains of the DMPC molecules in the bottom leaflet that is in contact with the gold is approximately 10 degrees smaller than that observed for the top leaflet that is exposed to the solution. These studies provide direct evidence that the structure of a phospholipid bilayer deposited at an electrode surface is affected by interaction with the metal.     

Organization and dynamics of NBD-labeled lipids in membranes analyzed by fluorescence recovery after photobleaching

Lateral diffusion of membrane constituents plays an important role in membrane organization and represents a central theme in current models describing the structure and function of biological membranes. Fluorescence recovery after photobleaching (FRAP) is a widely used approach that provides information regarding dynamic properties and spatial distribution of membrane constituents. On the basis of the unique concentration-dependent fluorescence emission properties of a fluorescently labeled cholesterol analogue modified at the tail region, 25-[N-[(7-nitrobenz-2-oxa-1,3-diazol-4-yl)-methyl]amino]-27-norcholesterol (25-NBD-cholesterol), we have previously shown that it exhibits local organization even at very low concentrations in membranes. In this paper, we address aspects regarding the molecular size and dynamics of such an organized assembly of 25-NBD-cholesterol by monitoring its lateral diffusion characteristics using FRAP. To obtain a comprehensive understanding of the organization and dynamics of 25-NBD-cholesterol in the membrane, we compare its diffusion properties to that of a fluorescent phospholipid analogue 1,2-dipalmitoyl-sn-glycero-3-phosphoethanolamine-N-(7-nitro-2-(1,3-benzoxadiazol-4-yl)) (NBD-PE). Our results indicate significant differences in the membrane dynamics of these NBD-labeled lipids. Importantly, on the basis of a novel wavelength-selective FRAP approach, our results show that the organization of 25-NBD-cholesterol is heterogeneous, with the presence of fast- and slow-diffusing species which could correspond to predominant populations of monomers and dimers of 25-NBD-cholesterol. The potential application of the wavelength-selective FRAP approach to monitor the organization and dynamics of molecules in membranes therefore represents an exciting possibility.     

Cation and anion transport through hydrophilic pores in lipid bilayers

To understand the origin of transmembrane potentials, formation of transient pores, and the movement of anions and cations across lipid membranes, we have performed systematic atomistic molecular dynamics simulations of palmitoyl-oleoyl-phosphatidylcholine (POPC) lipids. A double bilayer setup was employed and different transmembrane potentials were generated by varying the anion (Cl-) and cation (Na+) concentrations in the two water compartments. A transmembrane potential of approximately 350 mV was thereby generated per bilayer for a unit charge imbalance. For transmembrane potential differences of up to approximately 1.4 V, the bilayers were stable, over the time scale of the simulations (10-50 ns). At larger imposed potential differences, one of the two bilayers breaks down through formation of a water pore, leading to both anion and cation translocations through the pore. The anions typically have a short residence time inside the pore, while the cations show a wider range of residence times depending on whether they bind to a lipid molecule or not. Over the time scale of the simulations, we do not observe the discharge of the entire potential difference, nor do we observe pore closing, although we observe that the size of the pore decreases as more ions translocate. We also observed a rare lipid flip-flop, in which a lipid molecule translocated from one bilayer leaflet to the opposite leaflet, assisted by the water pore.     

Shape fluctuations of large unilamellar lipid vesicles observed by laser light scattering: influence of the small-scale structure

In the present paper, we apply the dynamic laser light scattering technique to investigate the dependence of the characteristic times of thermally induced shape fluctuation of large unilamellar vesicles (LUVs) on bilayer composition. After addressing single-component LUVs made of two common phospholipids, dimyristoylphosphatidylcholine (DMPC) and dipalmitoylphosphatidylcholine (DPPC), we investigate the changes in vesicle shape fluctuation times due to the presence of cholesterol and gangliosides (GM1), added in small amounts. The experimental results show that the addition of a second component, even in small amount, to DMPC vesicles induces a change in membrane fluctuation times. Moreover, in the case of ganglioside addition, also the disposition of GM1 within the bilayer is of importance. Quite unexpectedly, the symmetric or asymmetric disposition of GM1 has opposite effects on bilayer dynamics, the first resulting in a "hardening" and the second in a "softening" of the membrane. Those results support that the small-scale structure of the bilayer is important in determining the overall dynamics of the vesicle. They also suggest that the physiological disposition of GM1 in the outer leaflet of real cells has a significative result in mechanical terms, positively affecting the dynamics of the membrane.