Electrostatic interactions at charged lipid membranes. Measurement of surface pH with fluorescent lipoid pH indicators.

The 5-dimethylaminonapthalene-1-sulfonyl (dansyl) chromophore attached to the polar head groups of lipids has been used as a fluorescent lipoid pH indicator to evaluate the interfacial pH in lipid-water lamellar systems prepared from negatively charged lipids. The pH in the vicinity of the charged lipid bilayers is different from the pH of the bulk aqueous phase and the difference is a function of the electrolyte concentration in the aqueous phase and of the lipid packing in the bilayer. At a fixed electrolyte concentration in the aqueous phase, the observed interfacial pH is 0.6 to 0.7 pH units lower above the thermal phase transition of the lipid than it is below this temperature. A quantitative interpretation of the results is given on the basis of the Gouy-Chapman theory. The results indicate that the dansyl chromophore is located in front of the charged surface and its distance from this surface increases with a decrease in lipid packing.     

Binding of small basic peptides to membranes containing acidic lipids: theoretical models and experimental results.

We measured directly the binding of Lys3, Lys5, and Lys7 to vesicles containing acidic phospholipids. When the vesicles contain 33% acidic lipids and the aqueous solution contains 100 mM monovalent salt, the standard Gibbs free energy for the binding of these peptides is 3, 5, and 7 kcal/mol, respectively. The binding energies decrease as the mol% of acidic lipids in the membrane decreases and/or as the salt concentration increases. Several lines of evidence suggest that these hydrophilic peptides do not penetrate the polar headgroup region of the membrane and that the binding is mainly due to electrostatic interactions. To calculate the binding energies from classical electrostatics, we applied the nonlinear Poisson-Boltzmann equation to atomic models of the phospholipid bilayers and the basic peptides in aqueous solution. The electrostatic free energy of interaction, which arises from both a long-range coulombic attraction between the positively charged peptide and the negatively charged lipid bilayer, and a short-range Born or image charge repulsion, is a minimum when approximately 2.5 A (i.e., one layer of water) exists between the van der Waals surfaces of the peptide and the lipid bilayer. The calculated molar association constants, K, agree well with the measured values: K is typically about 10-fold smaller than the experimental value (i.e., a difference of about 1.5 kcal/mol in the free energy of binding). The predicted dependence of K (or the binding free energies) on the ionic strength of the solution, the mol% of acidic lipids in the membrane, and the number of basic residues in the peptide agree very well with the experimental measurements. These calculations are relevant to the membrane binding of a number of important proteins that contain clusters of basic residues.     

Calorimetric and spectroscopic studies of the polymorphic phase behavior of a homologous series of n-saturated 1,2-diacyl phosphatidylethanolamines.

The polymorphic phase behavior of a homologous series of n-saturated 1,2-diacyl phosphatidylethanolamines was investigated by differential scanning calorimetry, 31P-nuclear magnetic resonance, and Fourier transform infrared spectroscopy. Upon heating, aqueous dispersions of dried samples of the short- and medium-chain homologues (n < or = 17) exhibit single, highly energetic transitions from a dry, crystalline form to the fully hydrated, liquid-crystalline bilayer at temperatures higher than the lamellar gel-liquid-crystalline phase transition exhibited by fully hydrated samples. In contrast, the longer chain homologues (n > or = 18) first exhibit a transition from a dehydrated solid form to the hydrated L beta gel phase followed by the gel-liquid-crystalline phase transition normally observed with fully hydrated samples. The fully hydrated, aqueous dispersions of these lipids all exhibit reversible, fairly energetic gel-liquid-crystalline transitions at temperatures that are significantly higher than those of the corresponding phosphatidylcholines. In addition, at still higher temperatures, the longer chain members of this series (n > or = 16) exhibit weakly energetic transitions from the lamellar phase to an inverted nonlamellar phase. Upon appropriate incubation at low temperatures, aqueous dispersions of the shorter chain members of this homologous series (n < or = 16) form a highly ordered crystal-like phase that, upon heating, converts directly to the liquid-crystalline phase at the same temperature as do the aqueous dispersions of the dried lipid. The spectroscopic data indicate that unlike the n-saturated diacyl phosphatidylcholines, the stable crystal-like phases of this series of phosphatidylethanolamines describe an isostructural series in which the hydrocarbon chains are packed in an orthorhombic subcell and the headgroup and polar/apolar interfacial regions of the bilayer are effectively immobilized and substantially dehydrated. Our results suggest that many of the differences between the properties of these phosphatidylethanolamine bilayers and their phosphatidylcholine counterparts can be rationalized on the basis of stronger intermolecular interactions in the headgroup and interfacial regions of the phosphatidylethanolamine bilayers. These are probably the result of differences in the hydration and hydrogen bonding interactions involving the phosphorylethanolamine headgroup and moieties in the polar/apolar interfacial regions of phosphatidylethanolamine bilayers.     

Phosphorylation reverses the membrane association of peptides that correspond to the basic domains of MARCKS and neuromodulin.

Several groups have observed that phosphorylation causes the MARCKS (Myristoylated Alanine-Rich C Kinase Substrate) protein to move off cell membranes and phospholipid vesicles. Our working hypothesis is that significant membrane binding of MARCKS requires both hydrophobic insertion of the N-terminal myristate into the bilayer and electrostatic association of the single cluster of basic residues in the protein with acidic lipids and that phosphorylation reverses this electrostatic association. Membrane binding measurements with myristoylated peptides and phospholipid vesicles show this hydrophobic moiety could, at best, barely attach proteins to plasma membranes. We report here membrane binding measurements with basic peptides that correspond to the phosphorylation domains of MARCKS and neuromodulin. Binding of these peptides increases sigmoidally with the percent acidic lipid in the phospholipid vesicle and can be described by a Gouy-Chapman/mass action theory that explains how electrostatics and reduction of dimensionality produce apparent cooperativity. The electrostatic affinity of the MARCKS peptide for membranes containing 10% acidic phospholipids (10(4) M-1 = chi/[P], where chi is the mole ratio of peptide bound to the outer monolayer of the vesicles and [P] is the concentration of peptide in the aqueous phase) is the same as the hydrophobic affinity of the myristate moiety for bilayer membranes. Phosphorylation decreases the affinity of the MARCKS peptide for membranes containing 15% acidic lipid about 1000-fold and produces a rapid (t1/2 < 30 s) dissociation of the peptide from phospholipid vesicles.     

Large divalent cations and electrostatic potentials adjacent to membranes. A theoretical calculation

We have extended the Gouy-Chapman theory of the electrostatic diffuse double layer by considering the finite size of divalent cations in the aqueous phase adjacent to a charged surface. The divalent cations are modeled as either two point charges connected by an infinitely thin, rigid "rod" or two noninteracting point charges connected by an infinitely thin, flexible "string." We use the extended theory to predict the effects of a cation of length 10 A (1 nm) on the zeta and surface potentials of phospholipid bilayer membranes. The predictions of the rod and string models are similar to one another but differ markedly from the predictions of the Gouy-Chapman theory. Specifically, the extended model predicts that a large divalent cation will have a smaller effect on the potential adjacent to a negatively charged bilayer membrane than a point divalent cation, that the magnitude of this discrepancy will decrease as the Debye length increases, and that a large divalent cation will produce a negative zeta potential on a membrane formed from zwitterionic lipids. These predictions agree qualitatively with the experimental results obtained with the large divalent cation hexamethonium. We discuss the biological relevance of our calculations in the context of the interaction of cationic drugs with receptor sites on cell membranes.     

The Mechanism of Collapse of Heterogeneous Lipid Monolayers

Collapse of homogeneous lipid monolayers is known to proceed via wrinkling/buckling, followed by folding into bilayers in water. For heterogeneous monolayers with phase coexistence, the mechanism of collapse remains unclear. Here, we investigated collapse of lipid monolayers with coexisting liquid-liquid and liquid-solid domains using molecular dynamics simulations. The MARTINI coarse-grained model was employed to simulate monolayers of ∼80 nm in lateral dimension for 10–25 μs. The monolayer minimum surface tension decreased in the presence of solid domains, especially if they percolated. Liquid-ordered domains facilitated monolayer collapse due to the spontaneous curvature induced at a high cholesterol concentration. Upon collapse, bilayer folds formed in the liquid (disordered) phase; curved domains shifted the nucleation sites toward the phase boundary. The liquid (disordered) phase was preferentially transferred into bilayers, in agreement with the squeeze-out hypothesis. As a result, the composition and phase distribution were altered in the monolayer in equilibrium with bilayers compared to a flat monolayer at the same surface tension. The composition and phase behavior of the bilayers depended on the degree of monolayer compression. The monolayer-bilayer connection region was enriched in unsaturated lipids. Percolation of solid domains slowed down monolayer collapse by several orders of magnitude. These results are important for understanding the mechanism of two-to-three-dimensional transformations in heterogeneous thin films and the role of lateral organization in biological membranes. The study is directly relevant for the function of lung surfactant, and can explain the role of nanodomains in its surface activity and inhibition by an increased cholesterol concentration.     

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.     

Role of cholesterol in the formation and nature of lipid rafts in planar and spherical model membranes

Sterols play a crucial regulatory and structural role in the lateral organization of eukaryotic cell membranes. Cholesterol has been connected to the possible formation of ordered lipid domains (rafts) in mammalian cell membranes. Lipid rafts are composed of lipids in the liquid-ordered (l(o)) phase and are surrounded with lipids in the liquid-disordered (l(d)) phase. Cholesterol and sphingomyelin are thought to be the principal components of lipid rafts in cell and model membranes. We have used fluorescence microscopy and fluorescence recovery after photobleaching in planar supported lipid bilayers composed of porcine brain phosphatidylcholine (bPC), porcine brain sphingomyelin (bSM), and cholesterol to map the composition-dependence of l(d)/l(o) phase coexistence. Cholesterol decreases the fluidity of bPC bilayers, but disrupts the highly ordered gel phase of bSM, leading to a more fluid membrane. When mixed with bPC/bSM (1:1) or bPC/bSM (2:1), cholesterol induces the formation of l(o) phase domains. The fraction of the membrane in the l(o) phase was found to be directly proportional to the cholesterol concentration in both phospholipid mixtures, which implies that a significant fraction of bPC cosegregates into l(o) phase domains. Images reveal a percolation threshold, i.e., the point where rafts become connected and fluid domains disconnected, when 45-50% of the total membrane is converted to the l(o) phase. This happens between 20 and 25 mol % cholesterol in 1:1 bPC/bSM bilayers and between 25 and 30 mol % cholesterol in 2:1 bPC/bSM bilayers at room temperature, and at approximately 35 mol % cholesterol in 1:1 bPC/bSM bilayers at 37 degrees C. Area fractions of l(o) phase lipids obtained in multilamellar liposomes by a fluorescence resonance energy transfer method confirm and support the results obtained in planar lipid bilayers.     

Effects of membrane surface charge and calcium on the gating of rat brain sodium channels in planar bilayers [published erratum appears in J Gen Physiol 1989 Apr;93(4):following 760]

The voltage-dependent gating of single, batrachotoxin-activated Na channels from rat brain was studied in planar lipid bilayers composed of negatively charged or neutral phospholipids. The relationship between the probability of finding the Na channel in the open state and the membrane potential (Po vs. Vm) was determined in symmetrical NaCl, both in the absence of free Ca2+ and after the addition of Ca2+ to the extracellular side of the channel, the intracellular side, or both. In the absence of Ca2+, neither the midpoint (V0.5) of the Po vs. Vm relation, nor the steepness of the gating curve, was affected by the charge on the bilayer lipid. The addition of 7.5 mM Ca2+ to the external side caused a depolarizing shift in V0.5. This depolarizing shift was approximately 17 mV in neutral bilayers and approximately 25 mV in negatively charged bilayers. The addition of the same concentration of Ca2+ to only the intracellular side caused hyperpolarizing shifts in V0.5 of approximately 7 mV (neutral bilayers) and approximately 14 mV (negatively charged bilayers). The symmetrical addition of Ca2+ caused a small depolarizing shift in Po vs. Vm. We conclude that: (a) the Na channel protein possesses negatively charged groups on both its inner and outer surfaces. Charges on both surfaces affect channel gating but those on the outer surface exert a stronger influence. (b) Negative surface charges on the membrane phospholipid are close enough to the channel's gating machinery to substantially affect its operation. Charges on the inner and outer surfaces of the membrane lipid affect gating symmetrically. (c) Effects on steady-state Na channel activation are consistent with a simple superposition of contributions to the local electrostatic potential from charges on the channel protein and the membrane lipid.     

Structure and dynamics of the phosphatidylcholine and the phosphatidylethanolamine head group in L-M fibroblasts as studied by deuterium nuclear magnetic resonance.

Mouse fibroblast L-M cells were grown in tissue culture medium containing selectively deuterated choline or ethanolamine. Both compounds were incorporated into the corresponding phospholipids at levels greater than 50% thus leading to a selective deuteration of these phospholipid head groups. Choline and ethanolamine were labeled at either the alpha- or the beta-carbon atom and well-resolved deuterium and phosphorus n.m.r. spectra were obtained from intact cells, crude plasma membranes and lipid extracts, leading to the following conclusions. (i) A large fraction, if not all, of the phospholipids in the intact L-M cell membranes were organized in a liquid crystalline bilayer. (ii) The phosphoethanolamine and the phosphocholine head group conformation were found to be remarkably similar in pure lipid bilayers and in intact L-M cell membranes with the head group dipoles being oriented parallel to the membrane surface. (iii) The deuterium T1 spin lattice relaxation times fell in the range of 7-25 ms and were similar in intact L-M cells and in pure lipid model membranes, suggesting that the two head groups are not involved in strong interactions with membrane proteins. The rotational diffusion rate of the two head groups was reduced by at least a factor of 10 compared to molecules of the same size in aqueous solution. (iv) The phosphocholine head group was sensitive to the size and sign of membrane surface charges as verified in mixing experiments with charged lipids. In L-M cell membranes the phosphocholine appeared to sense an electrically neutral environment in spite of the fact that L-M cell membranes contain 10-20% negatively charged lipids.     

Molecular Dynamics Simulations of PIP2 and PIP3 in Lipid Bilayers: Determination of Ring Orientation, and the Effects of Surface Roughness on a Poisson-Boltzmann Description

Molecular dynamics (MD) simulations of phosphatidylinositol (4,5)-bisphosphate (PIP₂) and phosphatidylinositol (3,4,5)-trisphosphate (PIP₃) in 1-palmitoyl 2-oleoyl phosphatidylcholine (POPC) bilayers indicate that the inositol rings are tilted ∼40° with respect to the bilayer surface, as compared with 17° for the P-N vector of POPC. Multiple minima were obtained for the ring twist (analogous to roll for an airplane). The phosphates at position 1 of PIP₂ and PIP₃ are within an Ångström of the plane formed by the phosphates of POPC; lipids in the surrounding shell are depressed by 0.5-0.8 Å, but otherwise the phosphoinositides do not substantially perturb the bilayer. Finite size artifacts for ion distributions are apparent for systems of ∼26 waters/lipid, but, based on simulations with a fourfold increase of the aqueous phase, the phosphoinositide positions and orientations do not show significant size effects. Electrostatic potentials evaluated from Poisson-Boltzmann (PB) calculations show a strong dependence of potential height and ring orientation, with the maxima on the −25 mV surfaces (17.1 ± 0.1 Å for PIP₂ and 19.4 ± 0.3 Å for PIP₃) occurring near the most populated orientations from MD. These surfaces are well above the background height of 10 Å estimated for negatively charged cell membranes, as would be expected for lipids involved in cellular signaling. PB calculations on microscopically flat bilayers yield similar maxima as the MD-based (microscopically rough) systems, but show less fine structure and do not clearly indicate the most probable regions. Electrostatic free energies of interaction with pentalysine are also similar for the rough and flat systems. These results support the utility of a rigid/flat bilayer model for PB-based studies of PIP₂ and PIP₃ as long as the orientations are judiciously chosen.     

Bilayer interactions of pHLIP, a peptide that can deliver drugs and target tumors

The pH-dependent insertion of pHLIP across membranes is proving to be a useful property for targeting acidic tissues or tumors and delivering drugs attached to its C-terminus. It also serves as a model peptide for studies of protein insertion into membranes, so further elucidation of the insertion mechanism of pHLIP and its features is desirable. We examine how the peptide perturbs a model phosphatidylcholine membrane and how it associates with the lipid bilayer using an array of fluorescence techniques, including fluorescence anisotropy measurements of TMA-DPH anchored in bilayers, quenching of pHLIP fluorescence by brominated lipids and acrylamide, and measurements of energy transfer between aromatic residues of pHLIP and TMA-DPH. When pHLIP is bound to the surface of bilayers near neutral pH, the membrane integrity is preserved whereas the elastic properties of bilayers are changed as reported by an increase of membrane viscosity. When it is inserted, there is little perturbation of the lipids. The results also suggest that pHLIP can bind to the membrane surface in a shallow or a deep mode depending on the phase state of the lipids. Using parallax analysis, the change of the penetration depth of pHLIP was estimated to be 0.4 Å from the bilayer center and 2.8 Å from the membrane surface after the liquid-to-gel phase transition.     

The Basic Protein of CNS Myelin: Its Structure and Ligand Binding

Consideration of the evidence presented in this review leads to the following conclusions: (a) Isolated MBP in aqueous solution has little ordered secondary or tertiary structure. (b) In this state, the protein can associate with a wide range of hydrophobic and amphiphilic compounds, these interactions involving limited sections of the protein. (c) The strength of binding to bilayers and the accompanying conformational changes in the protein are greatest for systems containing acidic lipids, presumably because of the involvement of ionic interactions. (d) When bound to bilayers of acidic lipids, MBP will have substantially more ordered secondary structure than it manifests in aqueous solution, and it is likely to be oligomeric (possibly hexameric). (e) MBP does affect the organization of lipid aggregates. It influences strongly the separation of bilayers in multilayers of purified lipids, and at present this must be viewed as its prime role within myelin. The greatest impediment to our understanding of MBP is the lack of an assayable biological activity. In contrast to the situation with enzymes, for example, we have no functional test for changes in protein structure or changes accompanying interactions with other molecules. Current evidence suggests that the protein has a structural role within myelin and that its own three-dimensional structure is strongly dependent on the molecules with which it is associated. If this picture is correct, studies of the isolated protein or of the protein in reconstituted lipid systems may yield, at best, a rough guide to the structure within its biological environment. Further clarification of the structure and function of MBP may have to await development of more powerful techniques for studying proteins bound to large molecular aggregates, such as lipid bilayers. The paucity of generally applicable methods is reflected in the fact that even low resolution structures are known for only a handful of intrinsic membrane proteins, and even more limited information exists for proteins associated with membrane surfaces. However, the increasing use of a combination of electron microscopy and diffraction on two-dimensional arrays of proteins formed on lipid bilayers (Henderson et al., 1990) offers the hope that it may not be too long before it will be possible to study at moderate resolution the three-dimensional structure of MBP bound to a lipid membrane.     

Structure, stability, and thermodynamics of lamellar DNA-lipid complexes.

We develop a statistical thermodynamic model for the phase evolution of DNA-cationic lipid complexes in aqueous solution, as a function of the ratios of charged to neutral lipid and charged lipid to DNA. The complexes consist of parallel strands of DNA intercalated in the water layers of lamellar stacks of mixed lipid bilayers, as determined by recent synchrotron x-ray measurements. Elastic deformations of the DNA and the lipid bilayers are neglected, but DNA-induced spatial inhomogeneities in the bilayer charge densities are included. The relevant nonlinear Poisson-Boltzmann equation is solved numerically, including self-consistent treatment of the boundary conditions at the polarized membrane surfaces. For a wide range of lipid compositions, the phase evolution is characterized by three regions of lipid to DNA charge ratio, rho: 1) for low rho, the complexes coexist with excess DNA, and the DNA-DNA spacing in the complex, d, is constant; 2) for intermediate rho, including the isoelectric point rho = 1, all of the lipid and DNA in solution is incorporated into the complex, whose inter-DNA distance d increases linearly with rho; and 3) for high rho, the complexes coexist with excess liposomes (whose lipid composition is different from that in the complex), and their spacing d is nearly, but not completely, independent of rho. These results can be understood in terms of a simple charging model that reflects the competition between counterion entropy and inter-DNA (rho < 1) and interbilayer (rho > 1) repulsions. Finally, our approach and conclusions are compared with theoretical work by others, and with relevant experiments.     

Solvent induced conformational changes in renin inhibitor polypeptide

Conformation of the renin inhibitor peptide, Pro-His-Pro-Phe-His-Phe-Phe-Val-Tyr-Lys (RIP) has been studied in aqueous solution and in lipid bilayers using 500 MHz 1H NMR spectroscopy. Analysis of the NMR parameters indicates that in aqueous solution, RIP exists as a random coil. On incorporation into lipid bilayers, the peptide adopts a rigid and well defined conformation. The N-terminal end is stabilized by the hydrophobic environment of the lipid bilayer. The C-terminal end is located near the lipid-water interface and attains rigidity due to interaction with the phosphate groups of lipids. The observations emphasize the role of environment in stabilizing significantly different conformations of RIP in three different media--D2O, DMSO and lipid bilayers.     

Asymmetric charge distributions in planar bilayer systems.

Using the simple argument based on irreversible thermodynamics and the Gouy-Chapman theory of the double layer, we show that the equilibrium distribution of charged lipid molecules between the two surfaces of a bilayer is asymmetric if the two solutions bathing the surfaces have the same ionic strength but contain ions of different valencies. For example, if one bathing solution contains 0.10 M NaCl and the other contains 0.70 M NaCl and 0.10 M CaCl2, the ratio of charged lipid molecules of the two surfaces in a membrane that contains 50% total negative lipids is 1.46, leading to a transbilayer potential of 18 mV. A complete set of such numerical results is presented in four figures.     

Ether- versus Ester-Linked Phospholipid Bilayers Containing either Linear or Branched Apolar Chains

We studied the properties of bilayers formed by ether-and ester-containing phospholipids, whose hydrocarbon chains can be either linear or branched, using sn-1,2 dipalmitoyl, dihexadecyl, diphytanoyl, and diphytanyl phosphatidylcholines (DPPC, DHPC, DPhoPC, and DPhPC, respectively) either pure or in binary mixtures. Differential scanning calorimetry and confocal fluorescence microscopy of giant unilamellar vesicles concurred in showing that equimolar mixtures of linear and branched lipids gave rise to gel/fluid phase coexistence at room temperature. Mixtures containing DHPC evolved in time (0.5 h) from initial reticulated domains to extended solid ones when an equilibrium was achieved. The nanomechanical properties of supported planar bilayers formed by each of the four lipids studied by atomic force microscopy revealed average breakdown forces F_b decreasing in the order DHPC ≥ DPPC > DPhoPC >> DPhPC. Moreover, except for DPPC, two different F_b values were found for each lipid. Atomic force microscopy imaging of DHPC was peculiar in showing two coexisting phases of different heights, probably corresponding to an interdigitated gel phase that gradually transformed, over a period of 0.5 h, into a regular tilted gel phase. Permeability to nonelectrolytes showed that linear-chain phospholipids allowed a higher rate of solute + water diffusion than branched-chain phospholipids, yet the former supported a smaller extent of swelling of the corresponding vesicles. Ether or ester bonds appeared to have only a minor effect on permeability.     

The structure and thermotropic phase behaviour of dipalmitoylphosphatidylcholine codispersed with a branched-chain phosphatidylcholine

The structure and thermotropic phase behaviour of a fully hydrated binary mixture of dipalmitoylphosphatidylcholine and a branched-chain phosphatidylcholine, 1, 2-di(4-dodecyl-palmitoyl)-sn-glycero-3-phosphocholine, were examined using differential scanning calorimetry, synchrotron X-ray diffraction and freeze-fracture electron microscopy. The branched-chain lipid forms a nonlamellar phase when dispersed alone in aqueous medium. Mixed aqueous dispersions of the two phospholipids containing less than 33 mol% of the branched-chain lipid form lamellar phases over the whole temperature range were studied (4 degrees C to 60 degrees C). When present in proportions greater than 33 mol% it induces a hexagonal phase in mixed aqueous dispersions with dipalmitoylphosphatidylcholine at temperatures above the fluid phase transition. At temperatures below 35 degrees C a hexagonal phase coexists with a gel bilayer phase. The lamellar<-->nonlamellar transition can be explained satisfactorily on the basis of the shape of the molecule expressed in terms of headgroup and chain cross-sectional areas. At temperatures below 35 degrees C macroscopic phase separation of two gel phases takes place. Freeze-fracture electron microscopy revealed that one gel phase consists of bilayers with a highly regular, periodic superstructure (macro-ripples) whereas the other phase forms flat, planar bilayers. The macro-ripple phase appears to represent a relaxation structure required to adapt to the packing constraints imposed by the incorporation of the branched-chain lipid into the dipalmitoylphosphatidylcholine host bilayer. The data suggest that structural changes that take place on cooling the mixed dispersion below the lamellar<-->nonlamellar phase transition temperature cannot be adequately described using the molecular form concept. Instead it is necessary to take into account the detailed molecular form of the guest lipid as well as its physical properties.     

Membrane surface-charge titration probed by gramicidin A channel conductance.

We manipulate lipid bilayer surface charge and gauge its influence on gramicidin A channel conductance by two strategies: titration of the lipid charge through bulk solution pH and dilution of a charged lipid by neutral. Using diphytanoyl phosphatidylserine (PS) bilayers with CsCl aqueous solutions, we show that the effects of lipid charge titration on channel conductance are masked 1) by conductance saturation with Cs+ ions in the neutral pH range and 2) by increased proton concentration when the bathing solution pH is less than 3. A smeared charge model permits us to separate different contributions to the channel conductance and to introduce a new method for "bilayer pKa" determination. We use the Gouy-Chapman expression for the charged surface potential to obtain equilibria of protons and cations with lipid charges. To calculate cation concentration at the channel mouth, we compare different models for the ion distribution, exact and linearized forms of the planar Poisson-Boltzmann equation, as well as the construction of a "Gibbs dividing surface" between salt bath and charged membrane. All approximations yield the intrinsic pKain of PS lipid in 0.1 M CsCl to be in the range 2.5-3.0. By diluting PS surface charge at a fixed pH with admixed neutral diphytanoyl phosphatidylcholine (PC), we obtain a conductance decrease in magnitude greater than expected from the electrostatic model. This observation is in accord with the different conductance saturation values for PS and PC lipids reported earlier (, Biochim. Biophys. Acta. 552:369-378) and verified in the present work for solvent-free membranes. In addition to electrostatic effects of surface charge, gramicidin A channel conductance is also influenced by lipid-dependent structural factors.     

The ELBA force field for coarse-grain modeling of lipid membranes

A new coarse-grain model for molecular dynamics simulation of lipid membranes is presented. Following a simple and conventional approach, lipid molecules are modeled by spherical sites, each representing a group of several atoms. In contrast to common coarse-grain methods, two original (interdependent) features are here adopted. First, the main electrostatics are modeled explicitly by charges and dipoles, which interact realistically through a relative dielectric constant of unity (ε(r) = 1). Second, water molecules are represented individually through a new parametrization of the simple Stockmayer potential for polar fluids; each water molecule is therefore described by a single spherical site embedded with a point dipole. The force field is shown to accurately reproduce the main physical properties of single-species phospholipid bilayers comprising dioleoylphosphatidylcholine (DOPC) and dioleoylphosphatidylethanolamine (DOPE) in the liquid crystal phase, as well as distearoylphosphatidylcholine (DSPC) in the liquid crystal and gel phases. Insights are presented into fundamental properties and phenomena that can be difficult or impossible to study with alternative computational or experimental methods. For example, we investigate the internal pressure distribution, dipole potential, lipid diffusion, and spontaneous self-assembly. Simulations lasting up to 1.5 microseconds were conducted for systems of different sizes (128, 512 and 1058 lipids); this also allowed us to identify size-dependent artifacts that are expected to affect membrane simulations in general. Future extensions and applications are discussed, particularly in relation to the methodology's inherent multiscale capabilities.