Direct measurement of the interaction between phosphatidylglycerol bilayers in aqueous electrolyte solutions

A closer look at calcium-induced interactions between phosphatidylserine-(PS) doped liposomes and the structural effects caused by inclusion of gangliosides or polyethylene glycol- (PEG) modified lipids

The effects of polyethylene glycol- (PEG) modified lipids and gangliosides on the Ca2+ induced interaction between liposomes composed of palmitoyl-oleoyl phosphatidylethanolamine (POPE) and palmitoyl-oleoyl phosphatidylserine (POPS) was investigated at physiological ionic strength. Förster resonance energy transfer (FRET) studies complemented with dynamic light scattering (DLS) and cryo-transmission electron microscopy (Cryo-EM) show that naked liposomes tend to adhere, rupture, and collapse on each other's surfaces upon addition of Ca2+, eventually resulting in the formation of large multilamellar aggregates and bilayer sheets. Noteworthy, the presence of gangliosides or PEGylated lipids does not prevent the adhesion-rupture process, but leads to the formation of small, long-lived bilayer fragments/disks. PEGylated lipids seem to be more effective than gangliosides at stabilizing these structures. Attractive interactions arising from ion correlation are proposed to be a driving force for the liposome-liposome adhesion and rupture processes. The results suggest that, in contrast with the conclusions drawn from previous solely FRET-based studies, direct liposome-liposome fusion is not the dominating process triggered by Ca2+ in the systems studied.

Small-Angle Neutron Scattering for Studying Lipid Bilayer Membranes

Small-angle neutron scattering (SANS) is a powerful tool for studying biological membranes and model lipid bilayer membranes. The length scales probed by SANS, being from 1 nm to over 100 nm, are well-matched to the relevant length scales of the bilayer, particularly when it is in the form of a vesicle. However, it is the ability of SANS to differentiate between isotopes of hydrogen as well as the availability of deuterium labeled lipids that truly enable SANS to reveal details of membranes that are not accessible with the use of other techniques, such as small-angle X-ray scattering. In this work, an overview of the use of SANS for studying unilamellar lipid bilayer vesicles is presented. The technique is briefly presented, and the power of selective deuteration and contrast variation methods is discussed. Approaches to modeling SANS data from unilamellar lipid bilayer vesicles are presented. Finally, recent examples are discussed. While the emphasis is on studies of unilamellar vesicles, examples of the use of SANS to study intact cells are also presented.

How Well Can You Tailor the Charge of Lipid Vesicles?

Knowledge and control of surface charge or potential is important for tailoring colloidal interactions. In this work, we compare widely used zeta potential (ζ) measurements of charged lipid vesicle surface potential to direct measurements using the surface force apparatus (SFA). Our measurements show good agreement between the two techniques. On varying the fraction of anionic lipids dimyristoylphosphatidylserine (DMPS) or dimyristoylphosphatidylglycerol (DMPG) mixed with zwitterionic dimyristoylphosphatidylcholine (DMPC) from 0 to 100 mol % we observed a near-linear increase in membrane surface charge or potential up to 20-30 mol % charged lipids beyond which charge saturation occurred in physiological (high) salt conditions. Similarly, in low salt concentrations, a linear increase in charge/potential was found but only up to ∼5-10 mol % charged lipids beyond which the surface charge or potential leveled off. While a lower degree of ionization is expected due to the lower dielectric constant (ε ∼ 4) of the lipid acyl chain environment, increasing intramembrane electrostatic repulsion between neighboring charged lipid head groups at higher charge loading contributes to charge suppression. Measured potentials in physiological salt solutions were consistent with predictions using the Gouy-Chapman-Stern-Grahame (GCSG) model of the electrical double layer with Langmuir binding of counterions, but in low salt conditions, the model significantly overestimated the surface charge/potential. The much lower ionization in low salt (maximum ∼1-2% of total lipids ionized) instead was consistent with counterion condensation at the bilayer surface which limited the charge that could be obtained. The strong interplay between membrane composition, lipid headgroup ionization, electrolyte concentration, and solution pH complicates exact prediction and tuning of membrane surface charge for applications. However, the theoretical frameworks used here can provide guidelines to understand this interplay and establish a range of achievable potentials for a system and predict the response to triggers like pH and salt concentration changes.

Attraction between Like-Charged Macroions Mediated by Specific Counterion Configurations

Attraction between like-charged macroions is fundamental to many processes in biology, chemistry, and physics. It also plays an important role in industrial applications such as ion-extraction processes or catalysis. In this work, we report a novel mechanism by which attraction can be realized between spherical macroions at high ionic strength. It consists of specific configurations of two, three, and more counterions that appear between macroions with high statistical probability. The attraction is manifested in a minimum in the potential of mean force between the macroions at short distances. Its depth increases with increasing charge of the macroion, demonstrating that the attraction is electrostatic in nature. It is shown that the implicit solvent model with a distance-dependent dielectric constant can capture both the geometry and thermodynamics of charge-stabilized macroion dimers on the qualitative level. The results obtained for a model colloid with a smooth surface are extrapolated to more realistic systems. Evidence is found that the reported mechanism can be observed in small chemical compounds with encapsulated ions such as fullerenes.

Sculpting and fusing biomimetic vesicle networks using optical tweezers

Constructing higher-order vesicle assemblies has discipline-spanning potential from responsive soft-matter materials to artificial cell networks in synthetic biology. This potential is ultimately derived from the ability to compartmentalise and order chemical species in space. To unlock such applications, spatial organisation of vesicles in relation to one another must be controlled, and techniques to deliver cargo to compartments developed. Herein, we use optical tweezers to assemble, reconfigure and dismantle networks of cell-sized vesicles that, in different experimental scenarios, we engineer to exhibit several interesting properties. Vesicles are connected through double-bilayer junctions formed via electrostatically controlled adhesion. Chemically distinct vesicles are linked across length scales, from several nanometres to hundreds of micrometres, by axon-like tethers. In the former regime, patterning membranes with proteins and nanoparticles facilitates material exchange between compartments and enables laser-triggered vesicle merging. This allows us to mix and dilute content, and to initiate protein expression by delivering biomolecular reaction components.

Assembly of higher-order artificial vesicles can unlock new applications. Here, the authors use optical tweezers to construct user-defined 2D and 3D architectures of chemically distinct vesicles and demonstrate inter-vesicle communication and light-enabled compartment merging.

Cospreading of Anionic Phospholipids with Peptides of the Structure (KX)4K at the Air-Water Interface: Influence of Lipid Headgroup Structure and Hydrophobicity of the Peptide on Monolayer Behavior

Mixtures of anionic phospholipids (PG, PA, PS, and CL) with cationic peptides were cospread from a common organic solvent at the air-water interface. The compression of the mixed film was combined with epifluorescence microscopy or infrared reflection adsorption spectroscopy (IRRAS) to gain information on the interactions of the peptide with the different lipids. To evaluate the influence of the amino acid X of peptides with the sequence (KX)₄K on the binding, 1,2-dipalmitoyl-sn-glycero-3-phosphoglycerol (DPPG) was mixed with different peptides with increasing hydrophobicity of the uncharged amino acid X. The monolayer isotherms of DPPG/(KX)₄K mixtures show an increased area for the lift-off due to incorporation of the peptide into the liquid-expanded (LE) state of the lipid. The surface pressure for the transition from LE to the liquid-condensed (LC) state is slightly increased for peptides with amino acids X with moderate hydrophobicity. For the most hydrophobic peptide (KL)₄K two plateaus are seen at a charge ratio PG to K of 5:1, and a strongly increased transition pressure is observed for a charge ratio of 1:1. Epifluorescence microscopy images and infrared spectroscopy show that the lower plateau corresponds to the LE-LC phase transition of the lipid. The upper plateau is connected with a squeeze-out of the peptide into the subphase. To test the influence of the lipid headgroup structure on peptide binding (KL)₄K was cospread with different anionic phospholipids. The shift of the isotherm to larger areas for lift-off and to higher surface pressure for the LE-LC phase transition was observed for all tested anionic lipids. Epifluorescence microscopy reveals the formation of LC domains with extended filaments indicating a decrease in line tension due to accumulation of the peptides at the LC-domain boundaries. This effect depends on the size of the headgroup of the anionic phospholipid.

Interaction of the Antimicrobial Peptide Aurein 1.2 and Charged Lipid Bilayer

Aurein 1.2 is a potent antimicrobial peptide secreted by frog Litoria aurea. As a short membrane-active peptide with only 13 amino acids in sequence, it has been found to be residing on the surface of lipid bilayer and permeabilizing bacterial membranes at high concentration. However, the detail at the molecular level is largely unknown. In this study, we investigated the action of Aurein 1.2 in charged lipid bilayers composed of DMPC/DMPG. Oriented Circular Dichroism results showed that the peptide was on the surface of lipid bilayer regardless of the charged lipid ratio. Only at a very high peptide-to-lipid ratio (~1/10), the peptide became perpendicular to the bilayer, however no pore was detected by neutron in-plane scattering. To further understand how it interacted with charged lipid bilayers, we employed Small Angle Neutron Scattering to probe lipid distribution across bilayer leaflets in lipid vesicles. The results showed that Aurein 1.2 interacted strongly with negatively charged DMPG, causing strong asymmetry in lipid bilayer. At high concentration, while the vesicles were intact, we found additional structure feature on the bilayer. Our study provides a glimpse into how Aurein 1.2 disturbs anionic lipid-containing membranes without pore formation.

Important roles for membrane lipids in haloarchaeal bioenergetics

Recent advances in lipidomic analysis in combination with various physiological experiments set the stage for deciphering the structure-function of haloarchaeal membrane lipids. Here we focused primarily on changes in lipid composition of Haloferax volcanii, but also performed a comparative analysis with four other haloarchaeal species (Halobacterium salinarum, Halorubrum lacusprofundi, Halorubrum sodomense and Haloplanus natans) all representing distinctive cell morphologies and behaviors (i.e., rod shape vs. pleomorphic behavior). Common to all five haloarchaea, our data reveal an extraordinary high level of menaquinone, reaching up to 72% of the total lipids. This ubiquity suggests that menaquinones may function beyond their ordinary role as electron and proton transporter, acting simultaneously as ion permeability barriers and as powerful shield against oxidative stress. In addition, we aimed at understanding the role of cations interacting with the characteristic negatively charged surface of haloarchaeal membranes. We propose for instance that by bridging the negative charges of adjacent anionic phospholipids, Mg²⁺ acts as surrogate for cardiolipin, a molecule that is known to control curvature stress of membranes. This study further provides a bioenergetic perspective as to how haloarchaea evolved following oxygenation of Earth's atmosphere. The success of the aerobic lifestyle of haloarchaea includes multiple membrane-based strategies that successfully balance the need for a robust bilayer structure with the need for high rates of electron transport - collectively representing the molecular basis to inhabit hypersaline water bodies around the planet.

Aggregation of Calcium Silicate Hydrate Nanoplatelets

We study the aggregation of calcium silicate hydrate nanoplatelets on a surface by means of Monte Carlo and molecular dynamics simulations at thermodynamic equilibrium. Calcium silicate hydrate (C-S-H) is the main component formed in cement and is responsible for the strength of the material. The hydrate is formed in early cement paste and grows to form platelets on the nanoscale, which aggregate either on dissolving cement particles or on auxiliary particles. The general result is that the experimentally observed variations in these dynamic processes generically called growth can be rationalized from interaction free energies, that is, from pure thermodynamic arguments. We further show that the surface charge density of the particles determines the aggregate structures formed by C-S-H and thus their growth modes.

Coarse-Grained Model of SNARE-Mediated Docking

Synaptic transmission requires that vesicles filled with neurotransmitter molecules be docked to the plasma membrane by the SNARE protein complex. The SNARE complex applies attractive forces to overcome the long-range repulsion between the vesicle and membrane. To understand how the balance between the attractive and repulsive forces defines the equilibrium docked state we have developed a model that combines the mechanics of vesicle/membrane deformation with an apparently new coarse-grained model of the SNARE complex. The coarse-grained model of the SNARE complex is calibrated by comparison with all-atom molecular dynamics simulations as well as by force measurements in laser tweezer experiments. The model for vesicle/membrane interactions includes the forces produced by membrane deformation and hydration or electrostatic repulsion. Combining these two parts, the coarse-grained model of the SNARE complex with membrane mechanics, we study how the equilibrium docked state varies with the number of SNARE complexes. We find that a single SNARE complex is able to bring a typical synaptic vesicle to within a distance of ∼ 3 nm from the membrane. Further addition of SNARE complexes shortens this distance, but an overdocked state of >4-6 SNAREs actually increases the equilibrium distance.

Ionizable Nitroxides for Studying Local Electrostatic Properties of Lipid Bilayers and Protein Systems by EPR

Electrostatic interactions are known to play a major role in the myriad of biochemical and biophysical processes. Here, we describe biophysical methods to probe local electrostatic potentials of proteins and lipid bilayer systems that are based on an observation of reversible protonation of nitroxides by electron paramagnetic resonance (EPR). Two types of probes are described: (1) methanethiosulfonate derivatives of protonatable nitroxides for highly specific covalent modification of the cysteine's sulfhydryl groups and (2) spin-labeled phospholipids with a protonatable nitroxide tethered to the polar head group. The probes of both types report on their ionization state through changes in magnetic parameters and degree of rotational averaging, thus, allowing the electrostatic contribution to the interfacial pKa of the nitroxide, and, therefore, the local electrostatic potential to be determined. Due to their small molecular volume, these probes cause a minimal perturbation to the protein or lipid system. Covalent attachment secures the position of the reporter nitroxides. Experimental procedures to characterize and calibrate these probes by EPR, and also the methods to analyze the EPR spectra by simulations are outlined. The ionizable nitroxide labels and the nitroxide-labeled phospholipids described so far cover an exceptionally wide range of ca. 2.5-7.0 pH units, making them suitable to study a broad range of biophysical phenomena, especially at the negatively charged lipid bilayer surfaces. The rationale for selecting proper electrostatically neutral interface for probe calibration, and examples of lipid bilayer surface potential studies, are also described.

Study of curcumin behavior in two different lipid bilayer models of liposomal curcumin using molecular dynamics simulation

Liposomal formulation of curcumin is an important therapeutic agent for the treatment of various cancers. Despite extensive studies on the biological effects of this formulation in cancer treatment, much remains unknown about curcumin-liposome interactions. Understanding how different lipid bilayers respond to curcumin molecule may help us to design more effective liposomal curcumin. Here, we used molecular dynamics simulation method to investigate the behavior of curcumin in two lipid bilayers commonly used in preparation of liposomal curcumin, namely dipalmitoylphosphatidylcholine (DPPC) and dimyristoylphosphatidylglycerol (DMPG). First, the free energy barriers for translocation of one curcumin molecule from water to the lipid bilayer were determined by using the potential of mean force (PMF). The computed free energy profile exhibits a global minimum at the solvent-headgroup interface (LH region) for both lipid membranes. We also evaluated the free energy difference between the equilibrium position of curcumin in the lipid bilayer and bulk water as the excess chemical potential. Our results show that curcumin has the higher affinity in DMPG compared to DPPC lipid bilayer (-8.39 vs. -1.69 kBT) and this is related to more hydrogen bond possibility for curcumin in DMPG lipid membrane. Next, using an unconstrained molecular dynamic simulation with curcumin initially positioned at the center of lipid bilayer, we studied various properties of each lipid bilayer system in the presence of curcumin molecule that was in full agreement with PMF and experimental data. The results of these simulation studies suggest that membrane composition could have a large effect on interaction of curcumin-lipid bilayer.

Surface charge effects on optical trapping of nanometer-sized lipid vesicles

Optical trapping of nanometer-sized lipid vesicles has been challenging due to the low refractive index contrast of the thin lipid bilayer to the aqueous medium. Using an "optical bottle", a recently developed technique to measure interactions of nanoparticles trapped by an infrared laser, we report, for the first time, quantitative measurements of the trapping energy of charged lipid vesicles. We found that the trapping energy increases with the relative amount of anionic lipids (DOPG) to neutral lipids (DOPC) in vesicles. Moreover, as monovalent salt is added into the exterior solution of vesicles, the trapping energy rapidly approaches zero, and this decrease in trapping energy strongly depends on the amount of anionic lipids in vesicles. A simple model with our experimental observations explains that the trapping energy of charged lipid vesicles is highly correlated with the surface charge density and electric double layer. In addition, we demonstrated selective trapping of a binary mixture of vesicles in different mole fractions of charged lipids, a strategy that has potential implications on charge selective vesicle sorting for engineering applications.

Experimental and theoretical study of the silica particle interactions in the presence of multivalent rod-like ions

The silica particle interactions in the presence of spermidine were systematically investigated both from experimental and theoretical points of view. The hydrodynamic radii and the corresponding polydispersity indices of the colloidal silica particles were determined by dynamic light scattering (DLS) as a function of spermidine concentration. Whereas the effective size of the silica particles increases with increasing spermidine concentration (pointing to the particle aggregation), the polydispersity index first increases reaches a maximum and then further decreases with the increasing spermidine concentration. From the mobility measurements it was concluded that the increase in spermidine concentration causes less negative values of zeta potential, meaning that the adsorption of spermidine leads to the less negative silica surface. Moreover, Monte Carlo (MC) simulations also confirmed that the addition of spermidine reduces the repulsion between silica particles. The MC concentration profiles of spermidine close to the charged silica particle are in a very good agreement with the results obtained by theory. An important motivation for our study is the effectiveness of multivalent ions to coagulate colloidal suspensions; e.g., the multivalent ions are exploited in the water purification process.

Simulations of anionic lipid membranes: development of interaction-specific ion parameters and validation using NMR data

Overbinding of ions to lipid head groups is a potentially serious artifact in simulations of charged lipid bilayers. In this study, the Lennard-Jones radii in the CHARMM force field for interactions of Na(+) and lipid oxygen atoms of carboxyl, phosphate, and ester groups were revised to match osmotic pressure data on sodium acetate and electrophoresis data on palmitoyloleoyl phosphatidylcholine (POPC) vesicles. The new parameters were then validated by successfully reproducing previously published experimental NMR deuterium order parameters for dimyristoyl phosphatidylglycerol (DMPG) and newly obtained values for palmitoyloleoyl phosphatidylserine (POPS). Although the increases in Lennard-Jones diameters are only 0.02-0.12 Å, they are sufficient to reduce Na+ binding, and thereby increase surface areas per lipid by 5-10% compared with the unmodified parameters.

Surface electrostatics of lipid bilayers by EPR of a pH-sensitive spin-labeled lipid

Many biophysical processes such as insertion of proteins into membranes and membrane fusion are governed by bilayer electrostatic potential. At the time of this writing, the arsenal of biophysical methods for such measurements is limited to a few techniques. Here we describe a, to our knowledge, new spin-probe electron paramagnetic resonance (EPR) approach for assessing the electrostatic surface potential of lipid bilayers that is based on a recently synthesized EPR probe (IMTSL-PTE) containing a reversibly ionizable nitroxide tag attached to the lipids' polar headgroup. EPR spectra of the probe directly report on its ionization state and, therefore, on electrostatic potential through changes in nitroxide magnetic parameters and the degree of rotational averaging. Further, the lipid nature of the probe provides its full integration into lipid bilayers. Tethering the nitroxide moiety directly to the lipid polar headgroup defines the location of the measured potential with respect to the lipid bilayer interface. Electrostatic surface potentials measured by EPR of IMTSL-PTE show a remarkable (within ±2%) agreement with the Gouy-Chapman theory for anionic DMPG bilayers in fluid (48°C) phase at low electrolyte concentration (50 mM) and in gel (17°C) phase at 150-mM electrolyte concentration. This agreement begins to diminish for DMPG vesicles in gel phase (17°C) upon varying electrolyte concentration and fluid phase bilayers formed from DMPG/DMPC and POPG/POPC mixtures. Possible reasons for such deviations, as well as the proper choice of an electrostatically neutral reference interface, have been discussed. Described EPR method is expected to be fully applicable to more-complex models of cellular membranes.

Interactions between charged surfaces mediated by molecules with spatially distributed charges

A short review of recent theoretical advances in studies of the interaction between highly charged systems is presented. Such a system could not be described by the mean field theory. More advanced methods have to be used in order to introduce the correlations between highly charged particles. In this work I focus on the system of highly charged surfaces, separated by a solution of molecules with spatially distributed charge. Two different representations of the molecular shape will be considered: rod-like and spherical. The system will be theoretically described by the density functional theory. For sufficiently long molecules and large surface charge densities, an attractive force between like-charged surfaces arises due to the spatially distributed charges within the molecules. The added salt has influence on the condition for the attractive force between like-charged surfaces. The theoretical results will be compared with Monte Carlo (MC) simulations. Recent measurements with multivalent rigid rod-like particles will be discussed.

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.

Phase transitions and spatially ordered counterion association in ionic-lipid membranes: theory versus experiment

Aqueous dispersions of phosphatidylglycerol (PG) lipids may present an anomalous chain-melting transition at low ionic strengths, as seen by different experimental techniques such as calorimetry or light scattering. The anomaly disappears at high ionic strengths or for longer acyl-chain lengths. In this article, we use a statistical model for the bilayer that distinguishes both lipid chain and headgroup states in order to compare model and experimental thermotropic and electrical properties. The effective van der Waals interactions among hydrophobic chains compete with the electrostatic repulsions between polar headgroups, which may be ionized (counterion dissociated) or electrically neutral (associated with counterions). Electric degrees of freedom introduce new thermotropic charge-ordered phases in which headgroup charges may be spatially ordered, depending on the electrolyte ionic strength, introducing a new rationale for experimental data on PGs. The thermal phases presented by the model for different chain lengths, at fixed ionic strength, compare well with an experimental phase diagram constructed on the basis of differential scanning calorimetry profiles. In the case of dispersions of DMPG (dimyristoyl phosphatidylglycerol) with added monovalent salt, the model properties reproduce the main features displayed by data from differential scanning calorimetry as well as the characteristic profile for the degree of ionization of the bilayer surface across the anomalous transition region, obtained from the theoretical interpretation of electrokinetic (conductivity and electrophoretic mobility) measurements.

Phase transitions and spatially ordered counterion association in ionic-lipid membranes: a statistical model

We propose a statistical model to account for the gel-fluid anomalous phase transitions in charged bilayer- or lamellae-forming ionic lipids. The model Hamiltonian comprises effective attractive interactions to describe neutral-lipid membranes as well as the effect of electrostatic repulsions of the discrete ionic charges on the lipid headgroups. The latter can be counterion dissociated (charged) or counterion associated (neutral), while the lipid acyl chains may be in gel (low-temperature or high-lateral-pressure) or fluid (high-temperature or low-lateral-pressure) states. The system is modeled as a lattice gas with two distinct particle types--each one associated, respectively, with the polar-headgroup and the acyl-chain states--which can be mapped onto an Ashkin-Teller model with the inclusion of cubic terms. The model displays a rich thermodynamic behavior in terms of the chemical potential of counterions (related to added salt concentration) and lateral pressure. In particular, we show the existence of semidissociated thermodynamic phases related to the onset of charge order in the system. This type of order stems from spatially ordered counterion association to the lipid headgroups, in which charged and neutral lipids alternate in a checkerboard-like order. Within the mean-field approximation, we predict that the acyl-chain order-disorder transition is discontinuous, with the first-order line ending at a critical point, as in the neutral case. Moreover, the charge order gives rise to continuous transitions, with the associated second-order lines joining the aforementioned first-order line at critical end points. We explore the thermodynamic behavior of some physical quantities, like the specific heat at constant lateral pressure and the degree of ionization, associated with the fraction of charged lipid headgroups.

Peptide-induced asymmetric distribution of charged lipids in a vesicle bilayer revealed by small-angle neutron scattering

Cellular membranes are complex mixtures of lipids, proteins, and other small molecules that provide functional, dynamic barriers between the cell and its environment, as well as between environments within the cell. The lipid composition of the membrane is highly specific and controlled in terms of both content and lipid localization. The membrane structure results from the complex interplay between the wide varieties of molecules present. Here, small-angle neutron scattering and selective deuterium labeling were used to probe the impact of the membrane-active peptides melittin and alamethicin on the structure of lipid bilayers composed of a mixture of the lipids dimyristoyl phosphatidylglycerol (DMPG) and chain-perdeuterated dimyristoyl phosphatidylcholine (DMPC). We found that both peptides enriched the outer leaflet of the bilayer with the negatively charged DMPG, creating an asymmetric distribution of lipids. The level of enrichment is peptide concentration-dependent and is stronger for melittin than it is for alamethicin. The enrichment between the inner and outer bilayer leaflets occurs at very low peptide concentrations and increases with peptide concentration, including when the peptide adopts a membrane-spanning, pore-forming state. The results suggest that these membrane-active peptides may have a secondary stressful effect on target cells at low concentrations that results from a disruption of the lipid distribution between the inner and outer leaflets of the bilayer that is independent of the formation of transmembrane pores.

Molecular dynamics simulations of mixed acidic/zwitterionic phospholipid bilayers

Anionic lipids are key components in the cell membranes. Many cell-regulatory and signaling mechanisms depend upon a complicated interplay between them and membrane-bound proteins. Phospholipid bilayers are commonly used as model systems in experimental or theoretical studies to gain insight into the structure and dynamics of biological membranes. We report here 200-ns-long MD simulations of pure (DMPC and DMPG) and mixed equimolar (DMPC/DMPG, DMPC/DMPS, and DMPC/DMPA) bilayers that each contain 256 lipids. The intra- and intermolecular interaction patterns in pure and mixed bilayers are analyzed and compared. The effect of monovalent ions (Na+) on the formation of salt-bridges is investigated. In particular, the number of Na(+)-mediated clusters in the presence of DMPS is higher than with DMPG and DMPA. We observe a preferential clustering of DMPS (and to some extent DMPA) lipids together rather than with DMPC molecules, which can explain the phase separation observed experimentally for DMPC/DMPS and DMPC/DMPA bilayers.

Forces between surfaces in liquids

Recent developments in the direct measurements of forces between surfaces in liquids at the ångstrom resolution level are reviewed. The results reveal a rich variety of interactions and interaction potentials that depend on the nature of the surfaces and intervening liquids. These results also shed new insights into liquid structure adjacent to surfaces and the interactions occurrig in complex systems, with implications in many different areas of chemical physics, biology, and technology. The origin of some important fundamental interactions, such as repulsive "hydration" forces and attractive "hydrophobic" forces, are still not understood and offer a challenge for experimental and theoretical work in this area.

Electrostatic interaction in atomic force microscopy

In atomic force microscopy, the stylus experiences an electrostatic force when imaging in aqueous medium above a charged surface. This force has been calculated numerically with continuum theory for a silicon nitrite or silicon oxide stylus. For comparison, the Van der Waals force was also calculated. In contrast to the Van der Waals attraction, the electrostatic force is repulsive. At a distance of 0.5 nm the electrostatic force is typically 10(-12)-10(-10)N and thus comparable in strength to the Van der Waals force. The electrostatic force increases with increasing surface charge density and decreases roughly exponentially with distance. It can be reduced by imaging in high salt concentrations. Below surface potentials of approximately 50 mV, a simple analytical approximation of the electrostatic force is described.

Measuring electrostatic, van der Waals, and hydration forces in electrolyte solutions with an atomic force microscope

In atomic force microscopy, the tip experiences electrostatic, van der Waals, and hydration forces when imaging in electrolyte solution above a charged surface. To study the electrostatic interaction force vs distance, curves were recorded at different salt concentrations and pH values. This was done with tips bearing surface charges of different sign and magnitude (silicon nitride, Al(2)O(3), glass, and diamond) on negatively charged surfaces (mica and glass). In addition to the van der Waals attraction, neutral and negatively charged tips experienced a repulsive force. This repulsive force depended on the salt concentration. It decayed exponentially with distance having a decay length similar to the Debye length. Typical forces were about 0.1 nN strong. With positively charged tips, purely attractive forces were observed. Comparing these results with calculations showed the electrostatic origin of this force.In the presence of high concentrations (> 3 M) of divalent cations, where the electrostatic force can be completely ignored, another repulsive force was observed with silicon nitride tips on mica. This force decayed roughly exponentially with a decay length of 3 nm and was approximately 0.07-nN strong. This repulsion is attributed to the hydration force.

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.

Models for phosphatidylglycerol lipids put to a structural test

Three atomistic empirical models for phosphatidylglycerol (PG) lipids are tested against structural data in the crystal and liquid crystal states. Simulations of the anhydrous crystal of dimyristoyl-phosphatidylglycerol (DMPG) show that only the CHARMM force field describes the conformation and interactions of PG head groups accurately. The other two models do not reproduce the native network of hydrogen bonds, suggesting the presence of biases in their conformational and nonbonded interaction properties. The CHARMM model is further validated in the biologically relevant liquid crystal phase by comparing experimental small-angle X-ray scattering spectra from DMPG unilamellar vesicles with data calculated from fluid bilayer simulations. The good agreement found in this model-free comparison implies that liquid crystal PG bilayers as described by CHARMM exhibit realistic bilayer thickness and lateral packing. Last, this model is used to simulate a fluid bilayer of palmitoyl-oleoyl-phosphatidylglycerol (POPG). The resulting view of the POPG bilayer structure is at variance with that proposed previously based on simulations, in particular, with respect to lateral packing of head groups and the role of counterions.

Role of nanomechanical properties in the tribological performance of phospholipid biomimetic surfaces

The role of phospholipid bilayers in controlling and reducing frictional forces between biological surfaces is investigated by three complementary experiments: friction forces are measured using a homemade tribometer, mechanical resistance to indentation is measured by AFM, and lipid bilayer degradation is controlled in situ during friction testing using fluorescence microscopy. DPPC lipid bilayers in the solid phase generate friction coefficients as low as 0.002 (comparable to that found for cartilage) that are stable through time. DOPC bilayers formed by the vesicle fusion method or the adsorption of mixed micelles generate higher friction coefficients. These coefficients increased through time, during which the bilayers degraded. The friction coefficient is correlated with the force needed to penetrate the bilayer with the AFM tip. With only one bilayer in the contact region, the friction increased to a similar value of about 0.08 for the DPPC and DOPC. Our study therefore shows that good mechanical stability of the bilayers is essential and suggests that the low friction coefficient is ensured by the hydration layers between adjacent lipid bilayers.

Thermal phase behavior of DMPG: the exclusion of continuous network and dense aggregates

1,2-Dimyristoyl-sn-glycero-3-phospho-rac-glycerol has been suggested to form at intermediate temperatures and at high concentrations in low-salt solutions as a continuous sponge phase (Heimburg, T.; Biltonen, R. L. Biochemistry 1994, 33, 9477-9488). In the present study, the changes in signals seen for a range of fluorescent probes during phase transformations of this phospholipid indicate continuous melting and a change in lipid packing, in accordance with previous reports. However, in accordance with Lamy-Freund and Riske (Lamy-Freund, M. T.; Riske, K. A. Chem. Phys. Lipids 2003, 122, 19-32), no enhancement of lipid mixing within the putative sponge phase region was seen, suggesting a lack of a connected lipid surface. Accordingly, a typical sponge phase cannot account for the properties of the intermediate phase. The low scattering intensities of the latter have also been taken as evidence for disaggregation. While dynamic light scattering and data for membranes containing poly(ethylene glycol)-ylated lipids could lend credence to disaggregation, the most likely explanation for the scattering data would appear to be a shape transition without significant changes in neither vesicle aggregation nor bilayer connectivity. An abrupt change in light scattering and signals from some of the fluorescent probes used reveals a new transition at Tt approximately 43 degrees C, with the formation of a more ordered interface.

Properties of mixed lipid monolayers assembled on hydrophobic surfaces through vesicle adsorption

Supported lipid films are becoming increasingly important tools for the study of membrane protein function because of the availability of high-sensitivity surface analytical and patterning techniques. In this study, we have characterized the physical chemical properties of lipid films assembled on hydrophobic surfaces through the spontaneous adsorption of large unilamellar lipid vesicles composed of dioleoylphosphatidylglycerol (DOPG) and dioleoylphosphatidylcholine (DOPC). The density of the lipid films was measured with surface plasmon resonance spectroscopy as the lipid composition of the vesicles and ionic concentration were varied. As expected, monolayer films were formed, but the density of the monolayers was found to be weakly dependent on the lipid composition of the vesicles and strongly dependent on the ionic concentration of the solution in contact with the monolayer. Atomic force microscopy (AFM) images of the lipid films indicate that they are composed of a homogeneous monolayer. Surface force measurements were used to determine the surface charge and DOPG density of the monolayers. The DOPG content of the films was found to be weakly dependent on the DOPG composition of the vesicles and strongly dependent on the salt concentration of the environment. A model has been developed to describe the behavior of the lipid composition of the films in terms of the hydrophobic, electrostatic, and steric forces acting on the lipid monolayer on the hydrophobic surface.

1,2-Dimyristoyl-sn-glycero-3-phosphoglycerol (DMPG) monolayers: influence of temperature, pH, ionic strength and binding of alkaline earth cations

Ion binding and lipid ionization of the acidic phospholipid 1,2-dimyristoyl-sn-glycero-3-phosphoglycerol (DMPG) in monolayers was studied by measuring the lateral pressure Pi as a function of the molecular area A at the air/water interface at different temperatures. The pH of the subphase (pH 2 and 7) and the ionic strength (NaCl) was varied. In addition, different divalent cations (1mM MgCl2, CaCl2 and SrCl2, pH 7) were added. DMPG is partly protonated on pure water at pH 7. An increase in the NaCl concentration in the subphase leads to film expansion. This effect is caused by an ionization of the headgroup of DMPG, i.e. a shift of the apparent pK. More condensed films are obtained on pure water at pH 2, due to the reduction of electrostatic repulsion by headgroup protonation and the possibility for the formation of a hydrogen bonding network. The divalent cations Mg2+, Ca2+ and Sr2+ interact differently with a DMPG monolayer in pure water at pH 7. In the presence of 1mM CaCl2 a condensation of the DMPG film is induced, whereas an expansion of the monolayer is observed in the presence of Mg2+ and Sr2+. Two counteracting effects are operative: (a) ionization of the headgroup due to electrostatic screening leads to film expansion and (b) binding of the divalent cations to the lipid headgroups leads to condensation. The latter effect is more pronounced in the case of Ca2+, whereas the binding of Mg2+ and Sr2+ to DMPG is weaker. Site-specific cation binding has to be assumed in addition to electrostatic effects.

Molecular dynamics simulations and 2H NMR study of the GalCer/DPPG lipid bilayer

Molecular dynamics simulations were performed on a two-component lipid bilayer system in the liquid crystalline phase at constant pressure and constant temperature. The lipid bilayers were composed of a mixture of neutral galactosylceramide (GalCer) and charged dipalmitoylphosphatidylglycerol (DPPG) lipid molecules. Two lipid bilayer systems were prepared with GalCer:DPPG ratio 9:1 (10%-DPPG system) and 3:1 (25%-DPPG system). The 10%-DPPG system represents a collapsed state lipid bilayer, with a narrow water space between the bilayers, and the 25%-DPPG system represents an expanded state with a fluid space of approximately 10 nm. The number of lipid molecules used in each simulation was 1024, and the length of the production run simulation was 10 ns. The simulations were validated by comparing the results with experimental data for several important aspects of the bilayer structure and dynamics. Deuterium order parameters obtained from (2)H NMR experiments for DPPG chains are in a very good agreement with those obtained from molecular dynamics simulations. The surface area per GalCer lipid molecule was estimated to be 0.608 +/- 0.011 nm(2). From the simulated electron density profiles, the bilayer thickness defined as the distance between the phosphorus peaks across the bilayer was calculated to be 4.21 nm. Both simulation systems revealed a tendency for cooperative bilayer undulations, as expected in the liquid crystalline phase. The interaction of water with the GalCer and DPPG oxygen atoms results in a strong water ordering in a spherical hydration shell and the formation of hydrogen bonds (H-bonds). Each GalCer lipid molecule makes 8.6 +/- 0.1 H-bonds with the surrounding water, whereas each DPPG lipid molecule makes 8.3 +/- 0.1 H-bonds. The number of water molecules per GalCer or DPPG in the hydration shell was estimated to be 10-11 from an analysis of the radial distribution functions. The formation of the intermolecular hydrogen bonds was observed between hydroxyl groups from the opposing GalCer sugar headgroups, giving an energy of adhesion in the range between -1.0 and -3.4 erg/cm(2). We suggest that this value is the contribution of the hydrogen-bond component to the net adhesion energy between GalCer bilayers in the liquid crystalline phase.

Electrostatic interactions between model mitochondrial membranes

Lipids are very diverse in both their respective structures and functions; and cells exquisitely control membrane composition. One intriguing issue is the specific role of lipids in modulating the physical properties of membranes. Cardiolipin (CL) is a unique four-tailed, doubly negatively charged lipid found predominately within the inner mitochondrial membrane, and is thought to be influential in determining the inner mitochondrial membrane potential and permeability. To determine the role of cardiolipin in modulating the charge properties of membranes, this study investigated the electrostatic interactions between mixed cardiolipin and phosphatidylcholine bilayers as a function of cardiolipin concentration. For physiologically relevant concentrations of cardiolipin, the surface charge density of the membrane was found to increase linearly with increasing concentration of cardiolipin. However, only a fraction of the cardiolipin molecules predicted to carry a charge from pK-values were ionized. Clearly environmental factors, beyond that of pH, play a role in determining the charge of bilayers containing cardiolipin.

Thermodynamics of cationic lipid-DNA complex formation as studied by isothermal titration calorimetry

The detailed analysis of the cationic lipid-DNA complex formation by means of isothermal titration calorimetry is presented. Most experiments were done using 1,2-dioleyl-sn-glycero-3-ethylphosphocholine (EDOPC), but basic titrations were also done using DOTAP, DOTAP:DOPC, and DOTAP:DOPE mixtures. Complex formation was endothermic with less than 1 kcal absorbed per mole of lipid or DNA charge. This enthalpy change was attributed to DNA-DNA mutual repulsion within the lamellar complex. The exception was DOTAP:DOPE-containing lipoplex for which the enthalpy of formation was exothermic, presumably because of DOPE amine group protonation. Experimental conditions, namely, direction and titration increment as well as concentration of titrant, which dictate the structure of resulting lipoplex (whether lamellar complex or DNA-coated vesicle), were found to affect the apparent thermodynamics of complex formation. The structure, in turn, influences the biological properties of the lipoplex. If the titration of lipid into DNA was carried out in large increments, the DeltaH was larger than when the injection increments were smaller, a finding that is consistent with increased vesicle disruption under large increments and which is expected theoretically. Cationic lipid-DNA binding was weak in high ionic strength solutions, however, the effective binding constant is within micromolar range because of macromolecular nature of the interaction.

Lateral diffusion coefficients in membranes measured by resonance energy transfer and a new algorithm for diffusion in two dimensions

We describe measurements of lateral diffusion in membranes using resonance energy transfer. The donor was a rhenium (Re) metal-ligand complex lipid, which displays a donor decay time near 3 micros. The long donor lifetime resulted in an ability to measure lateral diffusion coefficient below 10(-8) cm(2)/s. The donor decay data were analyzed using a new numerical algorithm for calculation of resonance energy transfer for donors and acceptors randomly distributed in two dimensions. An analytical solution to the diffusion equation in two dimensions is not known, so the equation was solved by the relaxation method in Laplace space. This algorithm allows the donor decay in the absence of energy transfer to be multiexponential. The simulations show that mutual lateral diffusion coefficients of the donor and acceptor on the order of 10(-8) cm(2)/s are readily recovered from the frequency-domain data with donor decay times on the microsecond timescale. Importantly, the lateral diffusion coefficients and acceptor concentrations can be recovered independently despite correlation between these parameters. This algorithm was tested and verified using the donor decays of a long lifetime rhenium lipid donor and a Texas red-lipid acceptor. Lateral diffusion coefficients ranged from 4.4 x 10(-9) cm(2)/s in 1,2-dimyristoyl-sn-glycero-3-[phospho-rac-(1-glycerol)] (DMPG) at 10 degrees C to 1.7 x 10(-7) cm(2)/s in 1,2-dioleoyl-sn-glycero-3-phosphocholine (DOPC) at 35 degrees C. These results demonstrated the possibility of direct measurements of lateral diffusion coefficients using microsecond decay time luminophores.