Specific binding of chloride ions to lipid vesicles and implications at molecular scale

Effect of cholesterol on the ion-membrane interaction: Zeta potential and dynamic light scattering study

Cholesterol in a bio-membrane plays a significant role in many cellular event and is known to regulate the functional activity of protein and ion channel. In this study we report a significant effect of cholesterol on the ion-membrane interaction. We prepare large unilamellar vesicles, composed of zwitterionic lipid DOPC and anionic lipid DOPG with different cholesterol concentration. Electrostatics of anionic membranes containing cholesterol in the presence of NaCl has systematically been explored using dynamic light scattering and zeta potential. Negative zeta potential of the membrane decreases its negative value with increasing ion concentration for all cholesterol concentrations. However, zeta potential itself decreases with increasing cholesterol content even in the absence of monovalent ions. Electrostatic behaviour of the membrane is determined from well-known Gouy Chapmann model. Negative surface charge density of the membrane decreases with increasing cholesterol content. Binding constant, estimated from the electrostatic double layer theory, is found to increase significantly in the presence of cholesterol. Comparison of electrostatic parameters of the membrane in the presence and absence of cholesterol suggests that cholesterol significantly alter the electrostatic behaviour of the membrane.

Temperature-Promoted Giant Unilamellar Vesicle (GUV) Aggregation: A Way of Multicellular Formation

The evolution of unicellular to multicellular life is considered to be an important step in the origin of life, and it is crucial to study the influence of environmental factors on this process through cell models in the laboratory. In this paper, we used giant unilamellar vesicles (GUVs) as a cell model to investigate the relationship between environmental temperature changes and the evolution of unicellular to multicellular life. The zeta potential of GUVs and the conformation of the headgroup of phospholipid molecules at different temperatures were examined using phase analysis light scattering (PALS) and attenuated total reflection-Fourier transform infrared spectroscopy (ATR-FTIR), respectively. In addition, the effect of increasing temperature on the aggregation of GUVs was further investigated in ionic solutions, and the possible mechanisms involved were explored. The results showed that increasing temperature reduced the repulsive forces between cells models and promoted their aggregation. This study could effectively contribute to our understanding of the evolution of primitive unicellular to multicellular life.

A practical guide to biologically relevant molecular simulations with charge scaling for electronic polarization

Molecular simulations can elucidate atomistic-level mechanisms of key biological processes, which are often hardly accessible to experiment. However, the results of the simulations can only be as trustworthy as the underlying simulation model. In many of these processes, interactions between charged moieties play a critical role. Current empirical force fields tend to overestimate such interactions, often in a dramatic way, when polyvalent ions are involved. The source of this shortcoming is the missing electronic polarization in these models. Given the importance of such biomolecular systems, there is great interest in fixing this deficiency in a computationally inexpensive way without employing explicitly polarizable force fields. Here, we review the electronic continuum correction approach, which accounts for electronic polarization in a mean-field way, focusing on its charge scaling variant. We show that by pragmatically scaling only the charged molecular groups, we qualitatively improve the charge-charge interactions without extra computational costs and benefit from decades of force field development on biomolecular force fields.

Induced fusion and aggregation of bacterial outer membrane vesicles: Experimental and theoretical analysis

Recombinantly engineered bacterial outer membrane vesicles (OMVs) are promising vaccine delivery vehicles. The diversity of exogenous antigens delivered by OMVs can be enhanced by induced fusion of OMV populations. To date there are no reports of induced fusion of bacterial OMVs. Here we measure the pH and salt-induced aggregation and fusion of OMVs and analyze the processes against the Derjaguin-Landau-Verwey-Overbeek (DLVO) colloidal stability model. Vesicle aggregation and fusion kinetics were investigated for OMVs isolated from native E. coli (Nissle 1917) and lipopolysaccharide (LPS) modified E. coli (ClearColi) strains to evaluate the effect of lipid type on vesicle aggregation and fusion. Electrolytes and low pHs induced OMV aggregation for both native and modified LPS constructs, approaching a calculated fusion efficiency of ~25% (i.e. ~1/4 of collision events lead to fusion). However, high fusion efficiency was achieved for Nissle OMVs solely with decreased pH as opposed to a combination of low pH and increased divalent counterion concentration for ClearColi OMVs. The lipid composition of the OMVs from Nissle negatively impacted fusion in the presence of electrolytes, causing higher deviations from DLVO-predicted critical coagulation concentrations with monovalent counterions. The outcome of the work is a defined set of conditions under which investigators can induce OMVs to fuse and make various combinations of vesicle compositions.

Influence of External NaCl Salt on Membrane Rigidity of Neutral DOPC Vesicles

Sodium chloride (NaCl) is a very common molecule in biotic and abiotic aqueous environments. In both cases, variation of ionic strength is inevitable. In addition to the osmotic variation posed by such perturbations, the question of whether the interactions of monovalent ions Na⁺ and Cl^-, especially with the neutral head groups of phospholipid membranes are impactful enough to change the membrane rigidity, is still not entirely understood. We investigated the dynamics of 1,2-di-(octadecenoyl)-sn-glycero-3-phosphocholine (DOPC) vesicles with zwitterionic neutral head groups in the fluid phase with increasing external salt concentration. At higher salt concentrations, we observe an increase in bending rigidity from neutron spin echo (NSE) spectroscopy and an increase in bilayer thickness from small-angle X-ray scattering (SAXS). We compared different models to distinguish membrane undulations, lipid tail motions, and the translational diffusion of the vesicles. All of the models indicate an increase in bending rigidity by a factor of 1.3-3.6. We demonstrate that even down to t > 10 ns and for Q > 0.07 Å^-1, the observed NSE relaxation spectra are influenced by translational diffusion of the vesicles. For t < 5 ns, the lipid tail motion dominates the intermediate dynamic structure factor. As the salt concentration increases, this contribution diminishes. We introduced a time-dependent analysis for the bending rigidity that highlights only a limited Zilman-Granek time window in which the rigidity is physically meaningful.

Effect of Zwitterionic Phospholipid on the Interaction of Cationic Membranes with Monovalent Sodium Salts

Cationic lipids have attracted much attention because of their potential for biomedical applications, such as gene delivery. The gene transfection efficiency of cationic lipids is greatly influenced by the counterions as well as salt ions. We have systematically investigated the interaction of different monovalent sodium salts with positively charged membrane, composed of 1,2-dioleoyl- sn-glycero-3-phosphocholine (DOPC)/1,2-dioleoyl-3-trimethylammonium-propane (DOTAP) and DOTAP, using dynamic light scattering, zeta potential, isothermal titration calorimetry (ITC), and fluorescence spectroscopy techniques. Our results reveal that the affinity of anions with cationic membranes follows the sequence I^- ≫ Br^- > Cl^- according to descending order of their sizes and is consistent with the Hofmeister series. Interestingly, the electrostatic behavior of the DOTAP membrane in the presence of monovalent anions differs significantly from the DOPC/DOTAP membrane. This difference is due to the strong interplay between phosphocholine and trimethylammonium-propane (TAP) headgroups leading to the reorientation of the TAP group in the membrane. The binding constant of anions, derived from zeta potential and ITC is in agreement with the affinity of anions mentioned above. Among all anions, I^- shows strongest affinity, as evidenced from the rapid increase in hydrodynamic radius which eventually leads to the formation of large aggregates. The fluorescence spectroscopy of a lypophilic probe Nile red in the presence of cationic vesicles containing ions complements the I^- adsorption onto the membrane. Nonlinear Stern-Volmer plot, consisting of accessible and inaccessible Nile red to I^- is consistent with the zeta potential as well as ITC results.

Role of Counterions in Constant-pH Molecular Dynamics Simulations of PAMAM Dendrimers

Electrostatic interactions play a pivotal role in the structure and mechanism of action of most biomolecules. There are several conceptually different methods to deal with electrostatics in molecular dynamics simulations. Ionic strength effects are usually introduced using such methodologies and can have a significant impact on the quality of the final conformation space obtained. We have previously shown that full system neutralization can lead to wrong lipidic phases in the 25% PA/PC bilayer (J. Chem. Theory Comput. 2014,10, 5483-5492). In this work, we investigate how two limit approaches to the ionic strength treatment (implicitly with GRF or using full system neutralization with either GRF or PME) can influence the conformational space of the second-generation PAMAM dendrimer. Constant-pH MD simulations were used to map PAMAM's conformational space at its full pH range (from 2.5 to 12.5). Our simulations clearly captured the coupling between protonation and conformation in PAMAM. Interestingly, the dendrimer conformational distribution was almost independent of the ionic strength treatment methods, which is in contrast to what we have observed in charged lipid bilayers. Overall, our results confirm that both GRF with implicit ionic strength and a fully neutralized system with PME are valid approaches to model charged globular systems, using the GROMOS 54A7 force field.

Salinity induced effects on the growth rates and mycelia composition of basidiomycete and zygomycete fungi

Soil salinization, as the combination of primary and secondary events, can adversely affect organisms inhabiting this compartment. In the present study, the effects of increased salinity were assessed in four species of terrestrial fungi: Lentinus sajor caju, Phanerochaete chrysosporium, Rhizopus oryzae and Trametes versicolor. The mycelial growth and biochemical composition of the four fungi were determined under three exposure scenarios: 1) exposure to serial dilutions of natural seawater (SW), 2) exposure to serial concentrations of NaCl (potential surrogate of SW); and 3) exposure to serial concentrations of NaCl after a period of pre-exposure to low levels of NaCl. The toxicity of NaCl was slightly higher than that of SW, for all fungi species: the conductivities causing 50% of growth inhibition (EC₅₀) were within 14.9 and 22.0 mScm^-1 for NaCl and within 20.2 and 34.1 mScm^-1 for SW. Phanerochaete chrysosporium showed to be the less sensitive species, both for NaCl and SW. Exposure to NaCl caused changes in the biochemical composition of fungi, mainly increasing the production of polysaccharides. When fungi were exposed to SW this pattern of biochemical response was not observed. Fungi pre-exposed to low levels of salinity presented higher EC₅₀ than fungi non-pre-exposed, though 95% confidence limits overlapped, with the exception of P. chrysosporium. Pre-exposure to low levels of NaCl also induced changes in the biochemical composition of the mycelia of L. sajor caju and R. oryzae, relatively to the respective control. These results suggest that some terrestrial fungi may acquire an increased tolerance to NaCl after being pre-exposed to low levels of this salt, thus, suggesting their capacity to persist in environments that will undergo salinization. Furthermore, NaCl could be used as a protective surrogate of SW to derive safe salinity levels for soils, since it induced toxicity similar or higher than that of SW.

The Charge Properties of Phospholipid Nanodiscs

Phospholipids (PLs) are a major, diverse constituent of cell membranes. PL diversity arises from the nature of the fatty acid chains, as well as the headgroup structure. The headgroup charge is thought to contribute to both the strength and specificity of protein-membrane interactions. Because it has been difficult to measure membrane charge, ascertaining the role charge plays in these interactions has been challenging. Presented here are charge measurements on lipid Nanodiscs at 20°C in 100 mM NaCl, 50 mM Tris, at pH 7.4. Values are also reported for measurements made in the presence of Ca(2+) and Mg(2+) as a function of NaCl concentration, pH, and temperature, and in solvents containing other types of cations and anions. Measurements were made for neutral (phosphatidylcholine and phosphatidylethanolamine) and anionic (phosphatidylserine, phosphatidic acid, cardiolipin, and phosphatidylinositol 4,5-bisphosphate (PIP2)) PLs containing palmitoyl-oleoyl and dimyristoyl fatty acid chains. In addition, charge measurements were made on Nanodiscs containing an Escherichia coli lipid extract. The data collected reveal that 1) POPE is anionic and not neutral at pH 7.4; 2) high-anionic-content Nanodiscs exhibit polyelectrolyte behavior; 3) 3 mM Ca(2+) neutralizes a constant fraction of the charge, but not a constant amount of charge, for POPS and POPC Nanodiscs; 4) in contrast to some previous work, POPC only interacts weakly with Ca(2+); 5) divalent cations interact with lipids in a lipid- and ion-specific manner for POPA and PIP2 lipids; and 6) the monovalent anion type has little influence on the lipid charge. These results should help eliminate inconsistencies among data obtained using different techniques, membrane systems, and experimental conditions, and they provide foundational data for developing an accurate view of membranes and membrane-protein interactions.

Effect of Sodium and Chloride Binding on a Lecithin Bilayer. A Molecular Dynamics Study

The effect of ion binding on the structural, mechanical, dynamic and electrostatic properties of a 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) bilayer in a 0.5 M aqueous NaCl solution is investigated using classical atomistic molecular dynamics simulation with different force-field descriptions for ion-ion and ion-lipid interactions. Most importantly, the repulsive Lennard-Jones parameters for the latter were modified, such that approximately similar binding of cations and anions to the lipid membrane is achieved. This was done to qualitatively improve the apparent ion-lipid binding constants obtained from simulations with the original force field (Berger lipids and GROMOS87 ions in combination with the SPC water model) in comparison to experimental data. Furthermore, various parameters characterizing membrane structure, elasticity, order and dynamics are analyzed. It is found that ion binding as observed in simulations involving the modified in comparison to the original force-field description leads to: (i) a smaller salt-induced change in the area per lipid, which is in closer agreement with the experiment; (ii) a decrease in the area compressibility and bilayer thickness to values comparable to a bilayer in pure water; (iii) lipid deuterium order parameters and lipid diffusion coefficients on nanosecond timescales that are very similar to the values for a membrane in pure water. In general, salt effects on the structural properties of a POPC bilayer in an aqueous sodium-chloride solution appear to be reproduced reasonably well by the new force-field description. An analysis of membrane-membrane disjoining pressure suggests that the smaller salt-induced change in area per lipid induced by the new force-field description is not due to the alteration of membrane-associated net charge, but must rather be understood as a consequence of ion-specific effects on the arrangement of lipid molecules.

Ion-Pairing Contribution to the Liposomal Transport of Topotecan as Revealed by Mechanistic Modeling

Actively loaded liposomal formulations of anticancer agents have been widely explored due to their high drug encapsulation efficiencies and prolonged drug retention. Mathematical models to predict and optimize drug loading and release kinetics from these nanoparticle formulations would be useful in their development and may allow researchers to tune release profiles. Such models must account for the driving forces as influenced by the physicochemical properties of the drug and the microenvironment, and the liposomal barrier properties. This study employed mechanistic modeling to describe the active liposomal loading and release kinetics of the anticancer agent topotecan (TPT). The model incorporates ammonia transport resulting in generation of a pH gradient, TPT dimerization, TPT lactone ring-opening and -closing interconversion kinetics, chloride transport, and transport of TPT-chloride ion-pairs to describe the active loading and release kinetics of TPT in the presence of varying chloride concentrations. Model-based predictions of the kinetics of active loading at varying loading concentrations of TPT and release under dynamic dialysis conditions were in reasonable agreement with experiments. These findings identify key attributes to consider in optimizing and predicting loading and release of liposomal TPT that may also be applicable to liposomal formulations of other weakly basic pharmaceuticals.

Estimation of Nanodiamond Surface Charge Density from Zeta Potential and Molecular Dynamics Simulations

Molecular dynamics simulations of nanoparticles (NPs) are increasingly used to study their interactions with various biological macromolecules. Such simulations generally require detailed knowledge of the surface composition of the NP under investigation. Even for some well-characterized nanoparticles, however, this knowledge is not always available. An example is nanodiamond, a nanoscale diamond particle with surface dominated by oxygen-containing functional groups. In this work, we explore using the harmonic restraint method developed by Venable et al., to estimate the surface charge density (σ) of nanodiamonds. Based on the Gouy-Chapman theory, we convert the experimentally determined zeta potential of a nanodiamond to an effective charge density (σ_eff), and then use the latter to estimate σ via molecular dynamics simulations. Through scanning a series of nanodiamond models, we show that the above method provides a straightforward protocol to determine the surface charge density of relatively large (> ∼100 nm) NPs. Overall, our results suggest that despite certain limitation, the above protocol can be readily employed to guide the model construction for MD simulations, which is particularly useful when only limited experimental information on the NP surface composition is available to a modeler.

Diphenylhexatriene membrane probes DPH and TMA-DPH: A comparative molecular dynamics simulation study

Fluorescence spectroscopy and microscopy have been utilized as tools in membrane biophysics for decades now. Because phospholipids are non-fluorescent, the use of extrinsic membrane probes in this context is commonplace. Among the latter, 1,6-diphenylhexatriene (DPH) and its trimethylammonium derivative (TMA-DPH) have been extensively used. It is widely believed that, owing to its additional charged group, TMA-DPH is anchored at the lipid/water interface and reports on a bilayer region that is distinct from that of the hydrophobic DPH. In this study, we employ atomistic MD simulations to characterize the behavior of DPH and TMA-DPH in 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) and POPC/cholesterol (4:1) bilayers. We show that although the dynamics of TMA-DPH in these membranes is noticeably more hindered than that of DPH, the location of the average fluorophore of TMA-DPH is only ~3-4Å more shallow than that of DPH. The hindrance observed in the translational and rotational motions of TMA-DPH compared to DPH is mainly not due to significant differences in depth, but to the favorable electrostatic interactions of the former with electronegative lipid atoms instead. By revealing detailed insights on the behavior of these two probes, our results are useful both in the interpretation of past work and in the planning of future experiments using them as membrane reporters.

Computer modelling studies of the bilayer/water interface

This review summarises high resolution studies on the interface of lamellar lipid bilayers composed of the most typical lipid molecules which constitute the lipid matrix of biomembranes. The presented results were obtained predominantly by computer modelling methods. Whenever possible, the results were compared with experimental results obtained for similar systems. The first and main section of the review is concerned with the bilayer-water interface and is divided into four subsections. The first describes the simplest case, where the interface consists only of lipid head groups and water molecules and focuses on interactions between the lipid heads and water molecules; the second describes the interface containing also mono- and divalent ions and concentrates on lipid-ion interactions; the third describes direct inter-lipid interactions. These three subsections are followed by a discussion on the network of direct and indirect inter-lipid interactions at the bilayer interface. The second section summarises recent computer simulation studies on the interactions of antibacterial membrane active compounds with various models of the bacterial outer membrane. This article is part of a Special Issue entitled: Biosimulations edited by Ilpo Vattulainen and Tomasz Róg.

Effect of counterions on the binding affinity of Na+ ions with phospholipid membranes

We have systematically investigated the effect of counterions on the interaction of the Na⁺ ion with phospholipid membranes using dynamic light scattering, zeta potential, isothermal titration calorimetry and fluorescence spectroscopy techniques.

Molecular electrometer and binding of cations to phospholipid bilayers

The molecular electrometer – reorientation of lipid head due to bound charge – allows direct quantitative vetting of simulations against noninvasive NMR experiments; most simulation models overestimated lipid–cation affinities.

Interaction of Spin-Labeled Lipid Membranes with Transition Metal Ions

The large values of spin relaxation enhancement (RE) for PC spin-labels in the phospholipid membrane induced by paramagnetic metal salts dissolved in the aqueous phase can be explained by Heisenberg spin exchange due to conformational fluctuations of the nitroxide group as a result of membrane fluidity, flexibility of lipid chains, and, possibly, amphiphilic nature of the nitroxide label. Whether the magnetic interaction occurs predominantly via Heisenberg spin exchange (Ni) or by the dipole–dipole (Gd) mechanism, it is essential for the paramagnetic ion to get into close proximity to the nitroxide moiety for efficient RE. For different salts of Ni the RE in phosphatidylcholine membranes follows the anionic Hofmeister series and reflects anion adsorption followed by anion-driven attraction of paramagnetic cations on the choline groups. This adsorption is higher for chaotropic ions, e.g., perchlorate. (A chaotropic agent is a molecule in water solution that can disrupt the hydrogen bonding network between water molecules.) However, there is no anionic dependence of RE for model membranes made from negatively charged lipids devoid of choline groups. We used Ni-induced RE to study the thermodynamics and electrostatics of ion/membrane interactions. We also studied the effect of membrane composition and the phase state on the RE values. In membranes with cholesterol a significant difference is observed between PC labels with nitroxide tethers long enough vs not long enough to reach deep into the membrane hydrophobic core behind the area of fused cholesterol rings. This study indicates one must be cautious in interpreting data obtained by PC labels in fluid membranes in terms of probing membrane properties at different immersion depths when it can be affected by paramagnetic species at the membrane surface.

Membrane permeation of a peptide: it is better to be positive

A joint experimental and computational study investigates the translocation of a tryptophan molecule through a phospholipid membrane. Time dependent spectroscopy of the tryptophan side chain determines the rate of permeation into 150 nm phospholipid vesicles. Atomically detailed simulations are conducted to calculate the free energy profiles and the permeation coefficient. Different charging conditions of the peptide (positive, negative, or zwitterion) are considered. Both experiment and simulation reproduce the qualitative trend and suggest that the fastest permeation is when the tryptophan is positively charged. The permeation mechanism, which is revealed by molecular dynamics simulations, is of a translocation assisted by a local defect. The influence of long-range electrostatic interactions, such as the membrane dipole potential on the permeation process, is not significant.

Structural properties of archaeal lipid bilayers: small-angle X-ray scattering and molecular dynamics simulation study

Aeropyrum pernix is an aerobic hyperthermophilic archaeon that grows in harsh environmental conditions and as such possesses unique structural and metabolic features. Its membrane interfaces with the extreme environment and is the first line of defense from external factors. Therefore, lipids composing this membrane have special moieties that increase its stability. The membrane of A. pernix is composed predominantly of two polar lipids 2,3-di-O-sesterterpanyl-sn-glicerol-1-phospho-1'(2'-O-α-D-glucosyl)-myo-inositol (AGI) and 2,3-di-O-sesterterpanyl-sn-glicerol-1-phospho-myo-inositol (AI). Both have methyl branches in their lipid tails and ether linkages and carbohydrates in their headgroup. These moieties significantly affect the structure and dynamics of the bilayer. To provide a molecular level insight into these characteristics, we used here Molecular Dynamics (MD) simulations of lipid bilayers of composition similar to those of the archaeal membranes. First, we show that the electron density profiles along the normal to the bilayers derived from the simulations are in good agreement with the profiles obtained by the small-angle X-ray scattering (SAXS) technique, which provides confidence in the force fields used. Analyses of the simulation data show that the archaeal lipid bilayers are less hydrated than conventional phosphatidylcholine (PC) lipids and that their structure is not affected by the salt present in the surrounding solution. Furthermore, the lateral pressure in their hydrophobic core, due to the presence of the branched tails, is much higher than that at PC-based lipid bilayers. Both the methyl branched tails and the special headgroup moieties contribute to slow drastically the lateral diffusion of the lipids. Furthermore, we found that the lipid head groups associate via hydrogen bonding, which affects their reorientational dynamics. All together, our data provide links between the microscopic properties of these membranes and their overall stability in harsh environments.

Continuum approaches to understanding ion and peptide interactions with the membrane

Experimental and computational studies have shown that cellular membranes deform to stabilize the inclusion of transmembrane (TM) proteins harboring charge. Recent analysis suggests that membrane bending helps to expose charged and polar residues to the aqueous environment and polar head groups. We previously used elasticity theory to identify membrane distortions that minimize the insertion of charged TM peptides into the membrane. Here, we extend our work by showing that it also provides a novel, computationally efficient method for exploring the energetics of ion and small peptide penetration into membranes. First, we show that the continuum method accurately reproduces energy profiles and membrane shapes generated from molecular simulations of bare ion permeation at a fraction of the computational cost. Next, we demonstrate that the dependence of the ion insertion energy on the membrane thickness arises primarily from the elastic properties of the membrane. Moreover, the continuum model readily provides a free energy decomposition into components not easily determined from molecular dynamics. Finally, we show that the energetics of membrane deformation strongly depend on membrane patch size both for ions and peptides. This dependence is particularly strong for peptides based on simulations of a known amphipathic, membrane binding peptide from the human pathogen Toxoplasma gondii. In total, we address shortcomings and advantages that arise from using a variety of computational methods in distinct biological contexts.

Structural properties of POPC monolayers under lateral compression: computer simulations analysis

1-Palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC), a lipid comprising a saturated and an unsaturated acyl chain, belongs to the class of glycerophosphatidylcholines, major lipids in eukaryotic cell membranes. To get insight into the structural properties of this lipid within monolayers as membrane models, we performed molecular dynamics (MD) simulations of POPC monolayers under compression at the air/water interface. MD simulations were carried out at 300 K and at different surface pressures using the all-atom general Amber force field (GAFF). A good agreement was found between the simulated data and experimental isotherms. At surface pressures greater than 15 mN/m, two orientations of the head groups clearly appear: one nearly parallel to the monolayer interface and another one pointing toward the water. On the basis of the analysis of headgroup dihedral angles, we propose that the conformational variations around the bonds connecting the phosphorus atom to the adjacent oxygens are involved in these two orientations of the headgroup. The glycerol group orientation is characterized by a large distribution centered around 50° with respect to the monolayer normal. The acyl chains are predominantly in trans configuration from 7.5 to 43 mN/m surface pressures. Moreover, the calculated order parameter profiles of both chains suggest an independent behavior of the saturated and unsaturated chains that could be correlated with the formation of chain-type clusters observed along the simulated trajectories.

Macro- versus microscopic view on the electrokinetics of a water-membrane interface

Electrophoresis is an experimental method widely used to study electrostatic properties of interfaces. Here, we question the validity of the macroscopic theory for the planar geometry by Helmholtz and Smoluchowski by considering a POPC bilayer in an aqueous solution with 500 mM NaCl, using molecular dynamics simulations. We find that POPC shows positive electrophoretic mobility due to adsorption of sodium ions at the lipid headgroups. The theory assumes that the region in which the water density undergoes a transition from the bulk value to zero (interfacial width) is small compared to the Debye screening length. This separation of length scale is not fullfilled in the present case. Hence, contrasting the theory, we observe that the surface is not sharply defined, continuum hydrodynamics is not applicable, the effective viscosity in the double layer is increased compared to the bulk, and the zeta potential is dominated by the dipole potential. Our results might have widespread implications for interpretation of electrokinetic studies in general.