Intramembrane electrostatic interactions destabilize lipid vesicles

Natural Size Variation Amongst Protocells Leads to Survival and Growth Under Hypoosmotic Conditions

Membrane growth is vital to the evolution of cellular life. For model protocells, this is typically achieved through competition between different protocell populations or by adding extra amphiphiles. This work demonstrates an alternative mechanism for protocell membrane growth: hypoosmotic shocks, which could have occurred naturally in the protocell environment, leading to the redistribution of lipids within a single population of vesicles. Here we report that, even without an additional lipid supply, nanoscale and giant fatty acid vesicles can withstand substantial osmotic pressures through membrane growth, whilst also retaining a significant portion of their contents. This stands in contrast to phospholipid systems, which burst and release their contents under the same conditions. Notably, the fatty acid giant vesicles retained contents following hypoosmotic shocks ten times higher than the predicted tolerance levels. The observed robustness is likely enabled by the membrane's incorporation of additional amphiphiles from smaller and less tense vesicles within the same population. The dynamic nature of these fatty acid-based model protocells not only supports membrane growth but also enables protocell survival in hypoosmotic environments.

Investigation into temporal changes in the human bloodstain lipidome

Bloodstains are crucial pieces of physical evidences found at violent crime scenes, providing valuable information for reconstructing forensic cases. However, there is limited data on how bloodstain lipidomes change over time after deposition. Hence, we deployed a high-throughput high-performance liquid chromatography-mass spectrometry (HPLC-MS) approach to construct lipidomic atlases of bloodstains, whole blood, plasma, and blood cells from 15 healthy adults. A time-course analysis was also performed on bloodstains deposited for up to 6 months at room temperature (~ 25°C). The molecular levels of 60 out of 400 detected lipid species differed dramatically between bloodstain and whole blood samples, with major disturbances observed in membrane glycerophospholipids. More than half of these lipids were prevalent in the cellular and plasmic fractions; approximately 27% and 10% of the identified lipids were uniquely derived from blood cells and plasma, respectively. Furthermore, a subset of 65 temporally dynamic lipid species arose across the 6-month room-temperature deposition period, with decreased triacylglycerols (TAGs) and increased lysophosphatidylcholines (LPCs) as representatives, accounting for approximately 8% of the total investigated lipids. The instability of lipids increased linearly with time, with the most variability observed in the first 10 days. This study sheds light on the impact of air-drying bloodstains on blood components at room temperature and provides a list of potential bloodstain lipid markers for determining the age of bloodstains.

Revealing local molecular distribution, orientation, phase separation, and formation of domains in artificial lipid layers: Towards comprehensive characterization of biological membranes

Lipids, together with molecules such as DNA and proteins, are one of the most relevant systems responsible for the existence of life. Selected lipids are able to assembly into various organized structures, such as lipid membranes. The unique properties of lipid membranes determine their complex functions, not only to separate biological environments, but also to participate in regulatory functions, absorption of nutrients, cell-cell communication, endocytosis, cell signaling, and many others. Despite numerous scientific efforts, still little is known about the reason underlying the variability within lipid membranes, and its biochemical significance. In this review, we discuss the structural complexity of lipid membranes, as well as the importance to simplify studied systems in order to understand phenomena occurring in natural, complex membranes. Such systems require a model interface to be analyzed. Therefore, here we focused on analytical studies of artificial systems at various interfaces. The molecular structure of lipid membranes, specifically the nanometric thickens of molecular bilayer, limits in a major extent the choice of highly sensitive methods suitable to study such structures. Therefore, we focused on methods that combine high sensitivity, and/or chemical selectivity, and/or nanometric spatial resolution, such as atomic force microscopy, nanospectroscopy (tip-enhanced Raman spectroscopy, infrared nanospectroscopy), phase modulation infrared reflection-absorption spectroscopy, sum-frequency generation spectroscopy. We summarized experimental and theoretical approaches providing information about molecular structure and composition, lipid spatial distribution (phase separation), organization (domain shape, molecular orientation) of lipid membranes, and real-time visualization of the influence of various molecules (proteins, drugs) on their integrity. An integral part of this review discusses the latest achievements in the field of lipid layer-based biosensors.

Quantification of pulsed electric field for the rupture of giant vesicles with various surface charges, cholesterols and osmotic pressures

Electropermeabilization is a promising phenomenon that occurs when pulsed electric field with high frequency is applied to cells/vesicles. We quantify the required values of pulsed electric fields for the rupture of cell-sized giant unilamellar vesicles (GUVs) which are prepared under various surface charges, cholesterol contents and osmotic pressures. The probability of rupture and the average time of rupture are evaluated under these conditions. The electric field changes from 500 to 410 Vcm^-1 by varying the anionic lipid mole fraction from 0 to 0.60 for getting the maximum probability of rupture (i.e., 1.0). In contrast, the same probability of rupture is obtained for changing the electric field from 410 to 630 Vcm^-1 by varying the cholesterol mole fraction in the membranes from 0 to 0.40. These results suggest that the required electric field for the rupture decreases with the increase of surface charge density but increases with the increase of cholesterol. We also quantify the electric field for the rupture of GUVs containing anionic mole fraction of 0.40 under various osmotic pressures. In the absence of osmotic pressure, the electric field for the rupture is obtained 430 Vcm^-1, whereas the field is 300 Vcm^-1 in the presence of 17 mOsmL^-1, indicating the instability of GUVs at higher osmotic pressures. These investigations open an avenue of possibilities for finding the electric field dependent rupture of cell-like vesicles along with the insight of biophysical and biochemical processes.

Recent developments in the kinetics of ruptures of giant vesicles under constant tension

External tension in membranes plays a vital role in numerous physiological and physicochemical phenomena. In this review, recent developments in the constant electric- and mechanical-tension-induced rupture of giant unilamellar vesicles (GUVs) are considered. We summarize the results relating to the kinetics of GUV rupture as a function of membrane surface charge, ions in the bathing solution, lipid composition, cholesterol content in the membrane, and osmotic pressure. The mechanical stability and line tension of the membrane under these conditions are discussed. The membrane tension due to osmotic pressure and the critical tension of rupture for various membrane compositions are also discussed. The results and their analysis provide a biophysical description of the kinetics of rupture, along with insight into biological processes. Future directions and possible developments in this research area are included.

To Close or to Collapse: The Role of Charges on Membrane Stability upon Pore Formation

Resealing of membrane pores is crucial for cell survival. Membrane surface charge and medium composition are studied as defining regulators of membrane stability. Pores are generated by electric field or detergents. Giant vesicles composed of zwitterionic and negatively charged lipids mixed at varying ratios are subjected to a strong electric pulse. Interestingly, charged vesicles appear prone to catastrophic collapse transforming them into tubular structures. The spectrum of destabilization responses includes the generation of long-living submicroscopic pores and partial vesicle bursting. The origin of these phenomena is related to the membrane edge tension, which governs pore closure. This edge tension significantly decreases as a function of the fraction of charged lipids. Destabilization of charged vesicles upon pore formation is universal-it is also observed with other poration stimuli. Disruption propensity is enhanced for membranes made of lipids with higher degree of unsaturation. It can be reversed by screening membrane charge in the presence of calcium ions. The observed findings in light of theories of stability and curvature generation are interpreted and mechanisms acting in cells to prevent total membrane collapse upon poration are discussed. Enhanced membrane stability is crucial for the success of electroporation-based technologies for cancer treatment and gene transfer.

Analysis of purification of charged giant vesicles in a buffer using their size distribution

We have analyzed the purification of charged giant unilamellar vesicles (GUVs) prepared in a buffer containing various concentrations of salt using their size distribution. The membranes of GUVs were synthesized by a mixture of dioleoylphosphocholine (DOPC) and dioleoylphosphatidylglycerol (DOPG) lipids. The DOPG mole fractions (X) in the membranes of GUVs were 0.10, 0.25, 0.40, 0.55, 0.70, 0.90 in a physiological buffer containing 162 mM salt. In addition, for a fixed value of X the concentrations of salt (C) in the buffer were 12, 62, 112, 162, 212, 312, 362 mM. The size distribution histograms of experimentally investigated unpurified and purified GUVs were fitted with the lognormal distribution and obtained the multiplication factor [Formula: see text] for mean ([Formula: see text]) and [Formula: see text] for standard deviation ([Formula: see text]) of the lognormal distribution. The key parameters [Formula: see text] and [Formula: see text] were responsible for changing the average size and size distribution of unpurified GUVs to purified ones. The theoretically fitting equation of experimentally obtained X- and C-dependent values of [Formula: see text] and [Formula: see text] provided the calibration equation for estimating the average size of purified GUVs theoretically for any values of X and C. The estimated size of purified GUVs increased with the increase in electrostatic effect (i.e., increase in vesicle surface charge density or decrease in salt concentration in buffer). The estimated size of purified GUVs varied with X and C, which supported the previous report qualitatively. These investigations might be helpful in the field of cell/chemical biology for understanding the process of purification of vesicles/cells investigated by any other techniques.

An investigation into the critical tension of electroporation in anionic lipid vesicles

Irreversible electroporation (IRE) is a technique for the disruption of localized cells or vesicles by a series of short and high-frequency electric pulses which has been used for tissue ablation and treatment in certain diseases. It is well reported that IRE induces lateral tension in the membranes of giant unilamellar vesicles (GUVs). The GUVs are prepared by a mixture of anionic lipid dioleoylphosphatidylglycerol (DOPG) and neutral lipid dioleoylphosphatidylcholine (DOPC) using the natural swelling method. Here the influence of DOPG mole fraction, XDOPG, on the critical tension of electroporation in GUVs has been investigated in sodium chloride-containing PIPES buffer. The critical tension decreases from 9.0 ± 0.3 to 6.0 ± 0.2 mN/m with the increase of XDOPG from 0.0 to 0.60 in the membranes of GUVs. Hence an increase in XDOPG greatly decreases the mechanical stability of membranes. We develop a theoretical equation that fits the XDOPG dependent normalized critical tension, and obtain a binding constant for the lipid-ion interaction of 0.75 M-1. The decrease in the energy barrier for formation of the nano-size nascent or prepore state, due to the increase in XDOPG, is the main factor explaining the decrease in critical tension of electroporation in vesicles.

Electrostatic effects on the electrical tension-induced irreversible pore formation in giant unilamellar vesicles

Irreversible electroporation (IRE) is a new technique in which a series of short pulses with high frequency electrical energy is applied on the targeted regions of cells or vesicles for their destruction or rupture formation. IRE induces lateral tension in the membranes of vesicles. We have investigated the electrostatic interaction effects on the constant electrical tension-induced rate constant of irreversible pore formation in the membranes of giant unilamellar vesicles (GUVs). The electrostatic interaction has been varied by changing the salt concentration in buffer and the surface charge density of membranes. The membranes of GUVs are synthesized by a mixture of negatively charged lipid dioleoylphosphatidylglycerol (DOPG) and neutral lipid dioleoylphosphatidylcholine (DOPC) using the natural swelling method. The rate constant of pore formation increases with the decrease of salt concentration in buffer along with the increase of surface charge density of membranes. The tension dependent probability of pore formation and the rate constant of pore formation are fitted to the theoretical equation, and obtained the line tension of membranes. The decrease in energy barrier of a prepore due to electrostatic interaction is the key factor causing an increase of rate constant of pore formation.

Kinetics of irreversible pore formation under constant electrical tension in giant unilamellar vesicles

Stretching in the plasma membranes of cells and lipid membranes of vesicles plays important roles in various physiological and physicochemical phenomena. Irreversible electroporation (IRE) is a minimally invasive non-thermal tumor ablation technique where a series of short electrical energy pulses with high frequency is applied to destabilize the cell membranes. IRE also induces lateral tension due to stretching in the membranes of giant unilamellar vesicles (GUVs). Here, the kinetics of irreversible pore formation under constant electrical tension in GUVs has been investigated. The GUVs are prepared by a mixture of dioleoylphosphatidylglycerol and dioleoylphosphatidylcholine using the natural swelling method. An IRE signal of frequency 1.1 kHz is applied to the GUVs through a gold-coated electrode system. Stochastic pore formation is observed for several 'single GUVs' at a particular constant tension. The time course of the fraction of intact GUVs among all the examined GUVs is fitted with a single-exponential decay function from which the rate constant of pore formation in the vesicle, k_p, is calculated. The value of k_p increases with an increase of membrane tension. An increase in the proportion of negatively charged lipids in a membrane gives a higher k_p. Theoretical equations are fitted to the tension-dependent k_p and to the probability of pore formation, which allows us to obtain the line tension of the membranes. The decrease in the energy barrier for formation of the nano-size nascent or prepore state, due to the increase in electrical tension, is the main factor explaining the increase of k_p.

Electrostatic interaction effects on the size distribution of self-assembled giant unilamellar vesicles

The influence of electrostatic conditions (salt concentration of the solution and vesicle surface charge density) on the size distribution of self-assembled giant unilamellar vesicles (GUVs) is considered. The membranes of GUVs are synthesized by a mixture of dioleoylphosphatidylglycerol and dioleoylphosphatidylcholine in a physiological buffer using the natural swelling method. The experimental results are presented in the form of a set of histograms. The log-normal distribution is used for statistical treatment of results. It is obtained that the decrease of salt concentration and the increase of vesicle surface charge density of the membranes increase the average size of the GUV population. To explain the experimental results, a theory using the Helmholtz free energy of the system describing the GUV vesiculation is developed. The size distribution histograms and average size of GUVs under various conditions are fitted with the proposed theory. It is shown that the variation of the bending modulus due to changing of electrostatic parameters of the system is the main factor causing a change in the average size of GUVs.

Membrane interaction of off-pathway prion oligomers and lipid-induced on-pathway intermediates during prion conversion: A clue for neurotoxicity

Soluble oligomers of prion proteins (PrP), produced during amyloid aggregation, have emerged as the primary neurotoxic species, instead of the fibrillar end-products, in transmissible spongiform encephalopathies. However, whether the membrane is among their direct targets, that mediate the downstream adverse effects, remains a question of debate. Recently, questions arise from the formation of membrane-active oligomeric species generated during the β-aggregation pathway, either in solution, or in lipid environment. In the present study, we characterized membrane interaction of off-pathway oligomers from recombinant prion protein generated along the amyloid aggregation and compared to lipid-induced intermediates produced during lipid-accelerated fibrillation. Using calcein-leakage assay, we show that the soluble prion oligomers are the most potent in producing leakage with negatively charged vesicles. Binding affinities, conformational states, mode of action of the different PrP assemblies were determined by thioflavin T binding-static light scattering experiments on DOPC/DOPS vesicles, as well as by FTIR-ATR spectroscopy and specular neutron reflectivity onto the corresponding supported lipid bilayers. Our results indicate that the off-pathway PrP oligomers interact with lipid membrane via a distinct mechanism, compared to the inserted lipid-induced intermediates. Thus, separate neurotoxic mechanisms could exist following the puzzling intermediates generated in the different cell compartments. These results not only reveal an important regulation of lipid membrane on PrP behavior but may also provide clues for designing stage-specific and prion-targeted therapy.

Accelerated molecular dynamics simulation analysis of MSI-594 in a lipid bilayer

Multidrug resistance against the existing antibiotics is one of the most challenging threats across the globe. Antimicrobial peptides (AMPs), in this regard, are considered to be one of the effective alternatives that can overcome bacterial resistance. MSI-594, a 24-residue linear alpha-helical cationic AMP, has been shown to function via the carpet mechanism to disrupt bacterial membrane systems. To better understand the role of lipid composition in the function of MSI-594, in the present study, eight different model membrane systems have been studied using accelerated molecular dynamics (aMD) simulations. The simulated results are helpful in discriminating the particular effects of cationic MSI-594 against zwitterionic POPC, anionic POPG and POPS, and neutral POPE lipid moieties. Additionally, the effects of various heterogeneous POPC/POPG (7 : 3), POPC/POPS (7 : 3), and POPG/POPE (1 : 3 and 3 : 1) bilayer systems on the dynamic interaction of MSI-594 have also been investigated. The effect on the lipid bilayer due to the interaction with the peptide is characterized by lipid acyl-chain order, membrane thickness, and acyl-chain dynamics. Our simulation results show that the lipid composition affects the membrane interaction of MSI-594, suggesting that membrane selectivity is crucial to its mechanism of action. The results reported in this study are helpful to obtain accurate atomistic-level information governing MSI-594 and its membrane disruptive antimicrobial mechanism of action, and to design next generation potent antimicrobial peptides.

Flow Cytometric Analysis To Evaluate Morphological Changes in Giant Liposomes As Observed in Electrofusion Experiments

Liposome fusion is a way of supplying additional components for in-liposome biochemical reactions. Electrofusion is a method that does not require the addition of fusogens, which often alter the liposome dispersion, and is therefore useful for repetitive liposome fusion. However, the details of electrofusion have not been elucidated because of the limitations surrounding observing liposomes using a microscope. Therefore, we introduced fluorescent markers and high-throughput flow cytometry to analyze the morphological changes that occur in liposome electrofusion. (i) The content mixing was evaluated by a calcein-Co²⁺-EDTA system, in which green fluorescence from dequenched free calcein is detected when the quenched calcein-Co²⁺ complex and EDTA are mixed together. (ii) Liposome destruction was evaluated from the decrease in the total membrane volume of giant liposomes. (iii) Liposome fission was evaluated from the increase in the number of giant liposomes. By applying the flow cytometric analysis, we investigated the effect of three parameters (DC pulse, AC field, and lipid composition) on liposome electrofusion. The larger numbers or higher voltages of DC pulses induced liposome fusion and destruction with higher probability. The longer application time of the AC field induced liposome fusion, fission, and destruction with higher probability. Higher content of negatively charged POPG (≥19%) strongly inhibited liposome electrofusion.

Permeability modes in fluctuating lipid membranes with DNA-translocating pores

Membrane pores can significantly alter not only the permeation dynamics of biological membranes but also their elasticity. Large membrane pores able to transport macromolecular contents represent an interesting model to test theoretical predictions that assign active-like (non-equilibrium) behavior to the permeability contributions to the enhanced membrane fluctuations existing in permeable membranes [Maneville et al. Phys. Rev. Lett. 82, 4356 (1999)]. Such high-amplitude active contributions arise from the forced transport of solvent and solutes through the open pores, which becomes even dominant at large permeability. In this paper, we present a detailed experimental analysis of the active shape fluctuations that appear in highly permeable lipid vesicles with large macromolecular pores inserted in the lipid membrane, which are a consequence of transport permeability events occurred in an osmotic gradient. The experimental results are found in quantitative agreement with theory, showing a remarkable dependence with the density of membrane pores and giving account of mechanical compliances and permeability rates that are compatible with the large size of the membrane pore considered. The presence of individual permeation events has been detected in the fluctuation time-series, from which a stochastic distribution of the permeation events compatible with a shot-noise has been deduced. The non-equilibrium character of the membrane fluctuations in a permeation field, even if the membrane pores are mere passive transporters, is clearly demonstrated. Finally, a bio-nano-technology outlook of the proposed synthetic concept is given on the context of prospective uses as active membrane DNA-pores exploitable in gen-delivery applications based on lipid vesicles.

Fluorescence monitoring of the effect of oxidized lipids on the process of protein fibrillization

The kinetics of lysozyme and insulin amyloid formation in the presence of the oxidized phospholipids (oxPLs) was investigated using Thioflavin T fluorescence assay. The kinetic parameters of fibrillization process (lag time and apparent rate constant) have been determined upon varying the following experimental parameters: the type of lipid assemblies (premicellar aggregates and lipid bilayer vesicles), pH, temperature and lipid-to-protein molar ratio. It was found that oxPLs premicellar aggregates induced the more pronounced increase of the maximum Thioflavin T fluorescence, which is proportional to the extent of fibril formation, compared to the vesicles composed of the oxidized and unoxidized lipids. In contrast, the oxPLs, used as dispersions or included into vesicles, inhibited fibril nucleation and elongation under near-physiological conditions in vitro compared to liposomes containing unoxidized lipids. The results obtained provide deeper insight into the molecular mechanisms of the oxidative stress-modulated conformational diseases, and could be employed for the anti-amyloid drug development.

Osmotically-induced tension and the binding of N-BAR protein to lipid vesicles

The binding affinity of a curvature-sensing protein domain (N-BAR) is measured as a function of applied osmotic stress while the membrane curvature is nearly constant. Varying the osmotic stress allows us to control membrane tension, which provides a probe of the mechanism of binding. We study the N-BAR domain of the Drosophila amphiphysin and monitor its binding on 50 nm-radius vesicles composed of 90 mol% DOPC and 10 mol% PIP. We find that the bound fraction of N-BAR is enhanced by a factor of approximately 6.5 when the tension increases from zero to 2.6 mN m(-1). This tension-induced response can be explained by the hydrophobic insertion mechanism. From the data we extract a hydrophobic domain area that is consistent with known structure. These results indicate that membrane stress and strain could play a major role in the previously reported curvature-affinity of N-BAR.

A New View of Ras Isoforms in Cancers

Does small GTPase K-Ras4A have a single state or two states, one resembling K-Ras4B and the other N-Ras? A recent study of K-Ras4A made the remarkable observation that even in the absence of the palmitoyl, K-Ras4A can be active at the plasma membrane. Importantly, this suggests that K-Ras4A may exist in two distinct signaling states. In state 1, K-Ras4A is only farnesylated, like K-Ras4B; in state 2, farnesylated and palmitoylated, like N-Ras. The K-Ras4A hypervariable region sequence is positively charged, in between K-Ras4B and N-Ras. Taken together, this raises the possibility that the farnesylated but nonpalmitoylated state 1, like K-Ras4B, binds calmodulin and is associated with colorectal and other adenocarcinomas like lung cancer and pancreatic ductal adenocarcinoma. On the other hand, state 2 may be associated with melanoma and other cancers where N-Ras is a major contributor, such as acute myeloid leukemia. Importantly, H-Ras has two, singly and doubly, palmitoylated states that may also serve distinct functional roles. The multiple signaling states of palmitoylated Ras isoforms question the completeness of small GTPase Ras isoform statistics in different cancer types and call for reevaluation of concepts and protocols. They may also call for reconsideration of oncogenic Ras therapeutics.

Structural Significance of Lipid Diversity as Studied by Small Angle Neutron and X-ray Scattering

We review recent developments in the rapidly growing field of membrane biophysics, with a focus on the structural properties of single lipid bilayers determined by different scattering techniques, namely neutron and X-ray scattering. The need for accurate lipid structural properties is emphasized by the sometimes conflicting results found in the literature, even in the case of the most studied lipid bilayers. Increasingly, accurate and detailed structural models require more experimental data, such as those from contrast varied neutron scattering and X-ray scattering experiments that are jointly refined with molecular dynamics simulations. This experimental and computational approach produces robust bilayer structural parameters that enable insights, for example, into the interplay between collective membrane properties and its components (e.g., hydrocarbon chain length and unsaturation, and lipid headgroup composition). From model studies such as these, one is better able to appreciate how a real biological membrane can be tuned by balancing the contributions from the lipid's different moieties (e.g., acyl chains, headgroups, backbones, etc.).

Electrostatic interaction effects on tension-induced pore formation in lipid membranes

We investigated the effects of electrostatic interactions on the rate constant (k(p)) for tension-induced pore formation in lipid membranes of giant unilamellar vesicles under constant applied tension. A decrease in salt concentration in solution as well as an increase in surface charge density of the membranes increased k(p). These data indicate that k(p) increases as the extent of electrostatic interaction increases. We developed a theory on the effect of the electrostatic interactions on the free energy profile of the membrane containing a prepore and also on the values of k(p); this theory explains the experimental results and fits the experimental data reasonably well in the presence of weak electrostatic interactions. Based on these results, we conclude that a decrease in the free energy barrier of the prepore state due to electrostatic interactions is the main factor causing an increase in k(p).

pH-Driven Ordering Transitions in Liquid Crystal Induced by Conformational Changes of Cardiolipin

We report an investigation of interfacial phenomena occurring at aqueous-liquid crystal (LC) interfaces that triggers an orientational ordering transition of the LC in the presence of cardiolipin (CL) by varying pH, salt concentration and valence. In particular, the effects of three different conformational isomeric forms of the CL are observed to cause the response of the LC ordering to vary significantly from one to another at those interfaces. An ordering transition of the LC was observed when the CL is mostly in undissociated (at pH 2) and/or in bicyclic (at pH 4) conformation in which LC shows changes in the optical appearance from bright to dark. By contrast, no change in the optical appearance of the LC was observed when the pH of the system increases to 8 or higher in which the CL mostly exists in the open conformation. Fluorescence microscopy measurements further suggest that pH-dependent conformational forms of the CL have different ability to self-assemble (thus different packing efficiency) at aqueous-LC interfaces leading to dissimilar orientational behavior of the LC. Specifically, we found that change in headgroup-headgroup repulsion of the central phosphatidyl groups of the CL plays a key role in tuning the lipid packing efficiency and thus responses to interfacial phenomena. Orientational ordering transition of the LC was also observed as a function of increasing the ionic strength (buffer capacity) and strongly influenced in the presence of mono and divalent cations. Langmuir-Blodgett (LB) and polarization modulation infrared reflection absorption spectroscopy (PM-IRRAS) measurements provide further insight in modulation of the lipid packing efficiency and alkyl chain conformation of the CL at different pH and ionic conditions. Overall, the results presented in this paper establish that LCs offer a promising approach to differentiate different conformations (label free detection) of the CL through ordering transition of the LC at aqueous-LC interfaces.

The free energy of nanopores in tense membranes

Membrane nanopores are central players for a range of important cellular membrane remodeling processes as well as membrane rupture. Understanding pore formation in tense membranes requires comprehension of the molecular mechanism of pore formation and the associated free energy change as a function of the membrane tension. Here we propose a scheme to calculate the free energy change associated with the formation of a nanometer sized pore in molecular dynamics simulations as a function of membrane tension, which requires the calculation of only one computationally expensive potential of mean force. We show that membrane elastic theory can be used to estimate the pore formation free energy at different tension values from the free energy change in a relaxed membrane and the area expansion curves of the membranes. We have computed the pore formation free energy for a dipalmitoyl-phosphatidylcholine (DPPC) membrane at two different lateral pressure values, 1 bar and -40 bar, by calculating the potential of mean force acting on the head group of a single lipid molecule. Unrestrained simulations of the closing process confirm that the intermediate states along this reaction coordinate are reasonable and show that hydrophilic indentations spanning half the bilayer connected by a hydrophobic pore segment represent the corresponding high energy transition state. A comparison of the stability of simulated membranes to experiment at high loading rates show that, contrary to expectation, pores form too easily in small simulated membrane patches. This discrepancy originates from a combination of the absence of ions in the simulations and the small membrane size.

Lysylated phospholipids stabilize models of bacterial lipid bilayers and protect against antimicrobial peptides

Aminoacylated phosphatidylglycerols are common lipids in bacterial cytoplasmic membranes. Their presence in Staphylococcus aureus has been linked to increased resistance to a number of antibacterial agents, including antimicrobial peptides. Most commonly, the phosphatidylglycerol headgroup is esterified to lysine, which converts anionic phosphatidylglycerol into a cationic lipid with a considerably increased headgroup size. In the present work, we investigated the interactions of two well-studied antimicrobial peptides, cecropin A and mastoparan X, with lipid vesicles composed of 1-palmitoyl-2-oleoyl-phosphatidylcholine (POPC) and 1-palmitoyl-2-oleoyl-phosphatidylglycerol (POPG), containing varying fractions of an aminoacylated phosphatidylethanolamine, a stable analog of the corresponding phosphatidylglycerol-derivative. To differentiate between the effects of headgroup size and charge on peptide-lipid interactions, we synthesized two different derivatives. In one, the headgroup was modified by the addition of lysine, and in the other, by glutamine. The modification by glutamine results in a phospholipid with a headgroup size comparable to that of the lysylated version. However, whereas lysyl-phosphatidylethanolamine (Lys-PE) is cationic, glutaminyl-phosphatidylethanolamine (Gln-PE) is zwitterionic. We found that binding of mastoparan X and cecropin A was not significantly altered if the content of aminoacylated phosphatidylethanolamines did not exceed 20mol.%, which is the concentration found in bacterial membranes. However, a lysyl-phosphatidylethanolamine content of 20mol% significantly inhibits dye release from lipid vesicles, to a degree that depends on the peptide. In the case of mastoparan X, dye release is essentially abolished at 20mol.% lysyl-phosphatidylethanolamine, whereas cecropin A is less sensitive to the presence of lysyl-phosphatidylethanolamine. These observations are understood through the complex interplay between peptide binding and membrane stabilization as a function of the aminoacylated lipid content. This article is part of a Special Issue entitled: Interfacially Active Peptides and Proteins. Guest Editors: William C. Wimley and Kalina Hristova.

Cardiolipin effects on membrane structure and dynamics

Cardiolipin (CL) is a lipid with unique properties solely found in membranes generating electrochemical potential. It contains four acyl chains and tends to form nonlamellar structures, which are believed to play a key role in membrane structure and function. Indeed, CL alterations have been linked to disorders such as Barth syndrome and Parkinson's disease. However, the molecular effects of CL on membrane organization remain poorly understood. Here, we investigated the structure and physical properties of CL-containing membranes using confocal microscopy, fluorescence correlation spectroscopy, and atomic force microscopy. We found that the fluidity of the lipid bilayer increased and its mechanical stability decreased with CL concentration, indicating that CL decreases the packing of the membrane. Although the presence of up to 20% CL gave rise to flat, stable bilayers, the inclusion of 5% CL promoted the formation of flowerlike domains that grew with time. Surprisingly, we often observed two membrane-piercing events in atomic force spectroscopy experiments with CL-containing membranes. Similar behavior was observed with a lipid mixture mimicking the mitochondrial outer membrane composition. This suggests that CL promotes the formation of membrane areas with apposed double bilayers or nonlamellar structures, similar to those proposed for mitochondrial contact sites. All together, we show that CL induces membrane alterations that support the role of CL in facilitating bilayer structure remodeling, deformation, and permeabilization.

Direct measurements of effect of counterion concentration on mechanical properties of cationic vesicles

Theoretical analyses of charged membranes in aqueous solutions have long predicted that the electric double layer surrounding them contributes significantly to their mechanical properties. Here we report the first, direct experimental measurements of the effect of counterion concentration on the bending and area expansion modulus of cationic surfactant vesicles. Using the classical technique of micropipet aspiration coupled with a modified experimental protocol that is better suited for cationic vesicles, we successfully measure the mechanical properties of a double-tailed cationic surfactant, diethylesterdimethyl ammonium chloride (diC18:1 DEEDMAC) in CaCl2 solutions. It is observed that the area expansion modulus of the charged membrane exhibits no measurable dependence on the counterion concentration, in accordance with existing models of bilayer elasticity. The measured bending modulus, however, is found to vary nonmonotonically and exhibits a minimum in its variation with counterion concentration. The experimental results are interpreted based on theoretical calculations of charged and bare membrane mechanics. It is determined that the initial decrease in bending modulus with increasing counterion concentration may be attributed to a decreasing double layer thickness, while the subsequent increase is likely due to an increasing membrane thickness. These mechanical moduli measurements qualitatively confirm, for the first time, theoretical predictions of a nonmonotonic behavior and the opposing effects of ionic strength on the bending rigidity of charged bilayers.

Novicidin's membrane permeabilizing activity is driven by membrane partitioning but not by helicity: a biophysical study of the impact of lipid charge and cholesterol

We have investigated the interactions between the antimicrobial peptide Novicidin (Nc) and vesicles containing the phospholipid DOPC, with various amounts of DOPG and cholesterol using circular dichroism spectroscopy, calcein release, equilibrium dialysis and isothermal titration calorimetry. Nc adopts a random coil structure in the absence of lipids and in the presence of vesicles containing 100% DOPC. Lipids with 25-40% DOPG induce the highest level of helicity in Nc; higher DOPG levels lead to lower helicity levels and an altered tertiary arrangement of the peptide. However, the ability of Nc to permeabilize vesicles correlates not with helicity but rather with its overall membrane affinity, which is enthalpically favorable but opposed by entropy. Permeabilization declines with increasing mole percentage PG. Changes in helicity correlate with changes in enthalpy, reflecting the enthalpy of helix formation, but not with affinity. There is also a large favorable enthalpic interaction between Nc and lipids in the absence of negative charge and structural changes. Cholesterol slightly reduces membrane permeabilization but has little effect on Nc affinity and secondary structure, and probably protects the membrane by inducing the liquid ordered state. We conclude that helicity is not a prerequisite for activity, and charge-charge interactions are not the only major driving force for AMP interactions with membranes. Our data are compatible with a model in which a superficial binding mode with a large membrane surface binding area per peptide is more efficient than a more intimate embedding within the membrane environment.

Model-based approaches for the determination of lipid bilayer structure from small-angle neutron and X-ray scattering data

Some of our recent work has resulted in the detailed structures of fully hydrated, fluid phase phosphatidylcholine (PC) and phosphatidylglycerol (PG) bilayers. These structures were obtained from the joint refinement of small-angle neutron and X-ray data using the scattering density profile (SDP) models developed by Kučerka et al. (Biophys J 95:2356-2367, 2008; J Phys Chem B 116:232-239, 2012). In this review, we first discuss models for the standalone analysis of neutron or X-ray scattering data from bilayers, and assess the strengths and weaknesses inherent to these models. In particular, it is recognized that standalone data do not contain enough information to fully resolve the structure of naturally disordered fluid bilayers, and therefore may not provide a robust determination of bilayer structure parameters, including the much-sought-after area per lipid. We then discuss the development of matter density-based models (including the SDP model) that allow for the joint refinement of different contrast neutron and X-ray data, as well as the implementation of local volume conservation within the unit cell (i.e., ideal packing). Such models provide natural definitions of bilayer thicknesses (most importantly the hydrophobic and Luzzati thicknesses) in terms of Gibbs dividing surfaces, and thus allow for the robust determination of lipid areas through equivalent slab relationships between bilayer thickness and lipid volume. In the final section of this review, we discuss some of the significant findings/features pertaining to structures of PC and PG bilayers as determined from SDP model analyses.

Molecular structures of fluid phase phosphatidylglycerol bilayers as determined by small angle neutron and X-ray scattering

We have determined the molecular structures of commonly used phosphatidylglycerols (PGs) in the commonly accepted biologically relevant fluid phase. This was done by simultaneously analyzing small angle neutron and X-ray scattering data, with the constraint of measured lipid volumes. We report the temperature dependence of bilayer parameters obtained using the one-dimensional scattering density profile model - which was derived from molecular dynamics simulations - including the area per lipid, the overall bilayer thickness, as well as other intrabilayer parameters (e.g., hydrocarbon thickness). Lipid areas are found to be larger than their phosphatidylcholine (PC) counterparts, a result likely due to repulsive electrostatic interactions taking place between the charged PG headgroups even in the presence of sodium counterions. In general, PG and PC bilayers show a similar response to changes in temperature and chain length, but differ in their response to chain unsaturation. For example, compared to PC bilayers, the inclusion of a first double bond in PG lipids results in a smaller incremental change to the area per lipid and bilayer thickness. However, the extrapolated lipid area of saturated PG lipids to infinite chain length is found to be similar to that of PCs, an indication of the glycerol-carbonyl backbone's pivotal role in influencing the lipid-water interface.

Phosphatidylserine is a critical modulator for Akt activation

Akt activation relies on the binding of Akt to phosphatidylinositol-3,4,5-trisphosphate (PIP(3)) in the membrane. Here, we demonstrate that Akt activation requires not only PIP(3) but also membrane phosphatidylserine (PS). The extent of insulin-like growth factor-induced Akt activation and downstream signaling as well as cell survival under serum starvation conditions positively correlates with plasma membrane PS levels in living cells. PS promotes Akt-PIP(3) binding, participates in PIP(3)-induced Akt interdomain conformational changes for T308 phosphorylation, and causes an open conformation that allows for S473 phosphorylation by mTORC2. PS interacts with specific residues in the pleckstrin homology (PH) and regulatory (RD) domains of Akt. Disruption of PS-Akt interaction by mutation impairs Akt signaling and increases susceptibility to cell death. These data identify a critical function of PS for Akt activation and cell survival, particularly in conditions with limited PIP(3) availability. The novel molecular interaction mechanism for Akt activation suggests potential new targets for controlling Akt-dependent cell survival and proliferation.

Membrane Permeabilization by Oligomeric α-Synuclein: In Search of the Mechanism

Background

The question of how the aggregation of the neuronal protein α-synuclein contributes to neuronal toxicity in Parkinson's disease has been the subject of intensive research over the past decade. Recently, attention has shifted from the amyloid fibrils to soluble oligomeric intermediates in the α-synuclein aggregation process. These oligomers are hypothesized to be cytotoxic and to permeabilize cellular membranes, possibly by forming pore-like complexes in the bilayer. Although the subject of α-synuclein oligomer-membrane interactions has attracted much attention, there is only limited evidence that supports the pore formation by α-synuclein oligomers. In addition the existing data are contradictory.

Methodology/Principal Findings

Here we have studied the mechanism of lipid bilayer disruption by a well-characterized α-synuclein oligomer species in detail using a number of in vitro bilayer systems and assays. Dye efflux from vesicles induced by oligomeric α-synuclein was found to be a fast all-or-none process. Individual vesicles swiftly lose their contents but overall vesicle morphology remains unaltered. A newly developed assay based on a dextran-coupled dye showed that non-equilibrium processes dominate the disruption of the vesicles. The membrane is highly permeable to solute influx directly after oligomer addition, after which membrane integrity is partly restored. The permeabilization of the membrane is possibly related to the intrinsic instability of the bilayer. Vesicles composed of negatively charged lipids, which are generally used for measuring α-synuclein-lipid interactions, were unstable to protein adsorption in general.

Conclusions/Significance

The dye efflux from negatively charged vesicles upon addition of α-synuclein has been hypothesized to occur through the formation of oligomeric membrane pores. However, our results show that the dye efflux characteristics are consistent with bilayer defects caused by membrane instability. These data shed new insights into potential mechanisms of toxicity of oligomeric α-synuclein species.

Dilution technique to determine the hydrodynamic volume fraction of a vesicle suspension

A simple dilution method to determine the hydrodynamic volume fraction of vesicle suspensions is presented. The vesicle suspension is diluted with a solution containing a tracer Y, which is similar to a component X already present in the suspending fluid and which does not bind to or permeate through the vesicles. The concentrations of X and Y in the suspending fluid measured after dilution are used to determine the volume fraction. Using this technique, the volume fractions of vesicle suspensions comprising cationic vesicles prepared in solutions of CaCl(2) (X) were measured by dilution with MgCl(2) (Y) solutions. Various experimental parameters such as the concentration of the MgCl(2) diluents and the dilution volume ratio were studied and their effects optimized to arrive at a robust recipe for measuring the volume fraction. It is demonstrated that the technique can be applied to concentrated suspensions containing multilamellar and polydisperse vesicles.

The SCHOOL of nature: II. Protein order, disorder and oligomericity in transmembrane signaling

Recent reports have revealed that many proteins that do not adopt globular structures under native conditions, thus termed intrinsically disordered proteins (IDPs), are involved in cell signaling. Intriguingly, physiologically relevant oligomerization of IDPs has been recently observed and shown to exhibit unique biophysical characteristics, including the lack of significant changes in chemical shift and peak intensity upon binding. In this work, I summarize several distinct features of protein disorder that are especially important as related to receptor-mediated transmembrane signal transduction. I also hypothesize that interactions of IDPs with their protein or lipid partners represent a general biphasic process with the "no disorder-to-order" fast interaction which, depending on the interacting partner, may or may not be accompanied by the slow formation of a secondary structure. Further, I suggest signaling-related functional connections between protein order, disorder, and oligomericity and hypothesize that receptor oligomerization induced or tuned upon ligand binding outside the cell is translated across the membrane into protein oligomerization inside the cell, thus providing a general platform, the Signaling Chain HOmoOLigomerization (SCHOOL) platform, for receptor-mediated signaling. This structures our current multidisciplinary knowledge and views of the mechanisms governing the coupling of recognition to signal transduction and cell response. Importantly, this approach not only reveals previously unrecognized striking similarities in the basic mechanistic principles of function of numerous functionally diverse and unrelated surface membrane receptors, but also suggests the similarity between therapeutic targets, thus opening new horizons for both fundamental and clinically relevant studies.

Effect of hydrogen bonding on the rotational and translational dynamics of a headgroup-bound chromophore in bilayer lipid membranes

We have studied the interactions of the chromophore 1,2-dioleoyl-sn-glycero-3-phosphoethanolamine-N-7-nitro-2-1,3-benzoxadiazol-4-yl (18:1 NBD-PE) imbedded in the headgroup region of bilayer lipid membranes consisting of 1,2-dioleoyl-sn-glycero-3-phosphocholine (DOPC) and 1,2-dioleoyl-sn-glycero-3-[phospho-rac-(1-glycerol)] (DOPG). We have examined the molecular and mesoscale dynamics of the chromophore using time-correlated single photon counting (TCSPC) to measure rotational diffusion dynamics in lipid vesicles and fluorescence recovery after pattern photobleaching (FRAPP) to determine translational diffusion coefficients and mobile fractions in supported lipid bilayers. TCSPC data reveal that chromophore rotational diffusion rates in DOPG vesicles are statistically the same as in DOPC and mixed DOPC/DOPG vesicles, suggesting that the NBD-PE chromophore does not interact strongly with the headgroup region of these bilayers; however, FRAPP experiments show that lateral diffusion is statistically lower in mixed DOPC/DOPG-supported bilayers than in DOPC-supported bilayers. These results suggest that bilayers containing DOPG likely undergo interlipid headgroup hydrogen bonding interactions that suppress translational diffusion.

Ion dynamics in cationic lipid bilayer systems in saline solutions

Positively charged lipid bilayer systems are a promising class of nonviral vectors for safe and efficient gene and drug delivery. Detailed understanding of these systems is therefore not only of fundamental but also of practical biomedical interest. Here, we study bilayers comprising a binary mixture of cationic dimyristoyltrimethylammoniumpropane (DMTAP) and zwitterionic (neutral) dimyristoylphosphatidylcholine (DMPC) lipids. Using atomistic molecular dynamics simulations, we address the effects of bilayer composition (cationic to zwitterionic lipid fraction) and of NaCl electrolyte concentration on the dynamical properties of these cationic lipid bilayer systems. We find that, despite the fact that DMPCs form complexes via Na(+) ions that bind to the lipid carbonyl oxygens, NaCl concentration has a rather minute effect on lipid diffusion. We also find the dynamics of Cl(-) and Na(+) ions at the water-membrane interface to differ qualitatively. Cl(-) ions have well-defined characteristic residence times of nanosecond scale. In contrast, the binding of Na(+) ions to the carbonyl region appears to lack a characteristic time scale, as the residence time distributions displayed power-law features. As to lateral dynamics, the diffusion of Na(+) ions within the water-membrane interface consists of two qualitatively different modes of motion: very slow diffusion when ions are bound to DMPC, punctuated by fast rapid jumps when detached from the lipids. Overall, the prolonged dynamics of the Na(+) ions are concluded to be interesting for the physics of the whole membrane, especially considering its interaction dynamics with charged macromolecular surfaces.

Simulations of lipid transfer between a supported lipid bilayer and adsorbing vesicles

Recent experiments demonstrate transfer of lipid molecules between a charged, supported lipid membrane (SLB) and vesicles of opposite charge when the latter adsorb on the SLB. A simple phenomenological bead model has been developed to simulate this process. Beads were defined to be of three types, 'n', 'p', and '0', representing POPS (negatively charged), POEPC (positively charged), and POPC (neutral but zwitterionic) lipids, respectively. Phenomenological bead-bead interaction potentials and lipid transfer rate constants were used to account for the overall interaction and transfer kinetics. Using different bead mixtures in both the adsorbing vesicle and in the SLB (representing differently composed/charged vesicles and SLBs as in the reported experiments), we clarify under which circumstances a vesicle adsorbs to the SLB, and whether it, after lipid transfer and changed composition of the SLB and vesicle, desorbs back to the bulk again or not. With this model we can reproduce and provide a conceptual picture for the experimental findings.

Membrane binding mode of intrinsically disordered cytoplasmic domains of T cell receptor signaling subunits depends on lipid composition

Intrinsically disordered cytoplasmic domains of T cell receptor (TCR) signaling subunits including zeta(cyt) and CD3epsilon(cyt) all contain one or more copies of an immunoreceptor tyrosine-based activation motif (ITAM), tyrosine residues of which are phosphorylated upon receptor triggering. Membrane binding-induced helical folding of zeta(cyt) and CD3epsilon(cyt) ITAMs is thought to control TCR activation. However, the question whether or not lipid binding of zeta(cyt) and CD3epsilon(cyt) is necessarily accompanied by a folding transition of ITAMs remains open. In this study, we investigate whether the membrane binding mechanisms of zeta(cyt) and CD3epsilon(cyt) depend on the membrane model used. Circular dichroic and fluorescence data indicate that binding of zeta(cyt) and CD3epsilon(cyt) to detergent micelles and unstable vesicles is accompanied by a disorder-to-order transition, whereas upon binding to stable vesicles these proteins remain unfolded. Using electron microscopy and dynamic light scattering, we show that upon protein binding, unstable vesicles fuse and rupture. In contrast, stable vesicles remain intact under these conditions. This suggests different membrane binding modes for zeta(cyt) and CD3epsilon(cyt) depending on the bilayer stability: (1) coupled binding and folding, and (2) binding without folding. These findings explain the long-standing puzzle in the literature and highlight the importance of the choice of an appropriate membrane model for protein-lipid interactions studies.

Extensive bilayer perforation coupled with the phase transition region of an anionic phospholipid

At low ionic strength dimyristoylphosphatidylglycerol (DMPG) exhibits a broad phase transition region characterized by several superimposed calorimetric peaks. Peculiar properties, such as sample transparency, are observed only in the transition region. In this work we use differential scanning calorimetry (DSC), turbidity, and optical microscopy to study the narrowing of the transition region with the increase of ionic strength (0-500 mM NaCl). Upon addition of salt, the temperature extension of the transition region is reduced, and the number of calorimetric peaks decreases until a single cooperative event at T(m) = 23 degrees C is observed in the presence of 500 mM NaCl. The transition region is always coupled with a decrease in turbidity, but a transparent region is detected within the melting process only in the presence of up to 20 mM NaCl. The vanishing of the transparent region is associated with one of the calorimetric peaks. Optical microscopy of giant vesicles shows that bilayers first rupture when the transition region is reached and subsequently lose optical contrast. Fluorescence microscopy reveals a blurry and undefined image in the transparent region, suggesting a different lipid self-assembly. Overall sample turbidity can be directly related to the bilayer optical contrast. Our observations are discussed in terms of the bilayer being perforated along the transition region. In the narrower temperature interval of the transparent region, dependent on the ionic strength, the perforation is extensive and the bilayer completely loses the optical contrast.

Lipid transfer between charged supported lipid bilayers and oppositely charged vesicles

The bidirectional transfer of phospholipids between a charged, supported lipid bilayer (SLB) on SiO(2) and oppositely charged, unilamellar vesicles was studied by means of quartz crystal microbalance with dissipation (QCM-D) and optical reflectometry techniques. SLBs and vesicles were prepared from binary mixtures of 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) mixed with different fractions of either 1-palmitoyl-2-oleoyl-sn-glycero-3-[phospho-l-serine] (POPS) (negatively charged) or 1-palmitoyl-2-oleoyl-sn-glycero-3-ethylphosphocholine (POEPC) (positively charged). The interaction process consists of an attachment-transfer-detachment (ATD) sequence, where added vesicles first attach to and interact with the SLB, after which they detach, leaving behind a compositionally modified SLB and ditto vesicles. When the process is complete, there is no net addition or reduction of total lipid mass in the SLB, but lipid exchange has occurred. The time scale of the process varies from a few to many tens of minutes depending on the type of charged lipid molecule and the relative concentration of charged lipids in the two membranes. Electrostatically symmetric cases, where only the charge sign (but not the fraction of charged lipid) was reversed between the SLB and the vesicles, produce qualitatively similar but quantitatively different kinetics. The time scale of the interaction varies significantly between the two cases, which is attributed to a combination of the differences in the molecular structure of the lipid headgroup for the positively and the negatively charged lipids used, and to nonsymmetric distribution of charged lipids in the lipid membranes. The maximum amounts of attached vesicles during the ATD process were estimated to be 25-40% of a full monolayer of vesicles, with the precise amount depending on the actual charge fractions in the vesicles and the SLB. Interrupted vesicle exposure experiments, and experiments where the bulk concentration of vesicles was varied, show that vesicles in some cases may be trapped irreversibly on the SLB, when only partial transfer of lipid molecules has occurred. Additional supply of vesicles and further transfer induces detachment, when a sufficient amount of oppositely charged lipids has been transferred to the SLB, so that the latter becomes repulsive to the attached vesicles. Possible mechanistic scenarios, including monomer insertion and hemifusion models, are discussed. The observed phenomena and the actual SLB preparation process form a platform both for studies of various intermembrane molecular transfer processes and for modifying the composition of SLBs in a controlled way, for example, for biosensor and cell culture applications.

Vesicle diffusion close to a membrane: intermembrane interactions measured with fluorescence correlation spectroscopy

The protein machinery controlling membrane fusion (or fission) has been well studied; however, the role of vesicle diffusion near membranes in these critical processes remains unclear. We experimentally and theoretically investigated the dynamics of small vesicles (approximately 50 nm in diameter) that are diffusing near supported planar bilayers acting as "target" membranes. Using total internal reflection-fluorescence correlation spectroscopy, we examined the validity of theoretical analyses of vesicle-membrane interactions. Vesicles were hindered by hydrodynamic drag as a function of their proximity to the planar bilayer. The population distributions and diffusion kinetics of the vesicles were further affected by changing the ionic strength and pH of the buffer, as well as the lipid composition of the planar membrane. Effective surface charges on neutral bilayers were also analyzed by comparing experimental and theoretical data, and we show the possibility that vesicle dynamics can be modified by surface charge redistribution of the planar bilayer. Based on these results, we hypothesize that the dynamics of small vesicles, diffusing close to biomembranes, may be spatially restricted by altering local physiological conditions (e.g., salt concentration, lipid composition, and pH), which may represent an additional mechanism for controlling fusion (or fission) dynamics.

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.

Examining the contributions of lipid shape and headgroup charge on bilayer behavior

To better understand bilayer property dependency on lipid electrostatics and headgroup size, we use atomistic molecular dynamics simulations to study negatively charged and neutral lipid membranes. We compare the negatively charged phosphatidic acid (PA), which at physiological pH and salt concentration has a negative spontaneous curvature, with the negatively charged phosphatidylglycerol (PG) and neutrally charged phosphatidylcholine (PC), both of which have zero spontaneous curvature. The PA lipids are simulated using two different sets of partial charges for the headgroup and the varied charge distribution between the two PA systems results in significantly different locations for the Na(+) ions relative to the water/membrane interface. For one PA system, the Na(+) ions are localized around the phosphate group. In the second PA system, the Na(+) ions are located near the ester carbonyl atoms, which coincides with the preferred location site for the PG Na(+) ions. We find that the Na(+) ion location has a larger effect on bilayer fluidity properties than lipid headgroup size, where the A(lipid) and acyl chain order parameter values are more similar between the PA and PG bilayers that have Na(+) ions located near the ester groups than between the two PA bilayers.