The effect of calcium on the properties of charged phospholipid bilayers

We have performed molecular dynamics simulations to investigate the structure and dynamics of charged bilayers as well as the distribution of counterions at the bilayer interface. For this, we have considered the negatively charged di-myristoyl-phosphatidyl-glycerol (DMPG) and di-myristoyl-phosphatidyl-serine (DMPS) bilayers as well as a protonated di-myristoyl-phosphatidyl-serine (DMPSH) bilayer. We were particularly interested in calcium ions due to their important role in biological systems. Simulations performed in the presence of calcium ions (DMPG, DMPS) or sodium ions (DMPS) were run for 45-60 ns. Simulation results for DMPG are compared with fluorescence measurements. The average areas per molecule were 47.4 ± 0.5 Å² (DMPG with calcium), 47.3 ± 0.5 Å² (DMPS with calcium), 51.3 ± 1.0 Å² (DMPS with sodium) and 45.3 ± 0.5 Å² (DMPSH). The structure of the negatively charged lipids is significantly affected by the counterions, where calcium ions have a more pronounced effect than sodium ions. Calcium ions were found to be tightly bound to the anionic groups of the lipid molecules and as such appear to constitute an integral part of the membrane interface on nanoseconds time scales. In contrast to sodium ions, calcium ions are localised in a narrow (∼ 10 Å) band around the phosphate group. The interaction of calcium with the lipid molecules enhances the molecular packing of the PG and PS lipids. This observation is in good agreement with emission spectra of the membrane partitioning probe Laurdan in DMPG multilamellar vesicles that indicate an increase in the ordering of the DMPG bilayer due to the presence of calcium. Our results indicate that calcium ions, which often function as a second messengers in living cells have a pronounced effect on membrane structures, which may have implications during signal transduction events.     

The effect of calcium on the properties of charged phospholipid bilayers

We have performed molecular dynamics simulations to investigate the structure and dynamics of charged bilayers as well as the distribution of counterions at the bilayer interface. For this, we have considered the negatively charged di-myristoyl-phosphatidyl-glycerol (DMPG) and di-myristoyl-phosphatidyl-serine (DMPS) bilayers as well as a protonated di-myristoyl-phosphatidyl-serine (DMPSH) bilayer. We were particularly interested in calcium ions due to their important role in biological systems. Simulations performed in the presence of calcium ions (DMPG, DMPS) or sodium ions (DMPS) were run for 45-60 ns. Simulation results for DMPG are compared with fluorescence measurements. The average areas per molecule were 47.4+/-0.5 A2 (DMPG with calcium), 47.3+/-0.5 A2 (DMPS with calcium), 51.3+/-1.0 A2 (DMPS with sodium) and 45.3+/-0.5 A2 (DMPSH). The structure of the negatively charged lipids is significantly affected by the counterions, where calcium ions have a more pronounced effect than sodium ions. Calcium ions were found to be tightly bound to the anionic groups of the lipid molecules and as such appear to constitute an integral part of the membrane interface on nanoseconds time scales. In contrast to sodium ions, calcium ions are localised in a narrow (approximately 10 A) band around the phosphate group. The interaction of calcium with the lipid molecules enhances the molecular packing of the PG and PS lipids. This observation is in good agreement with emission spectra of the membrane partitioning probe Laurdan in DMPG multilamellar vesicles that indicate an increase in the ordering of the DMPG bilayer due to the presence of calcium. Our results indicate that calcium ions, which often function as a second messengers in living cells have a pronounced effect on membrane structures, which may have implications during signal transduction events.     

Variation in hydration forces between neutral phospholipid bilayers: evidence for hydration attraction

It is now generally recognized that hydration forces dominate close interactions of lipid hydrophilic surfaces. The commonality of their characteristics has been reasonably established. However, differences in measured net repulsion, particularly evident when phosphatidylethanolamine (PE) and phosphatidylcholine (PC) bilayers are compared, suggest there exists a variety of behavior wider than expected from earlier models of hydration and fluctuation repulsion balanced by van der Waals attraction. To find a basis for this diverse behavior, we have looked more closely at measured structural parameters, degrees of hydration, and interbilayer repulsive forces for the lamellar phases of the following lipids: 1-palmitoyl-2-oleoyl-PE (POPE), egg PE, transphosphatidylated egg PE (egg PE-T), mono- and dimethylated egg PE-T (MMPE and DMPE), 1-stearoyl-2-oleoyl-PC (SOPC), and mixtures of POPE and SOPC. POPE and SOPC bilayers differ not only in their maximum degrees of hydration but also in the empirical hydration force coefficients and decay lengths that characterize their interaction. When mixed with POPE, SOPC effects sudden and disproportionate increases in hydration. POPE, egg PE, and egg PE-T differ in their degree of hydration, molecular area, and hydration repulsion. A single methylation of egg PE-T almost completely converts its hydration and bilayer repulsive properties to those of egg PC; little progression of hydration is seen with successive methylations. In order to reconcile these observations with the conventional scheme of balancing interbilayer hydration and fluctuation-enhanced repulsion with van der Waals attraction, it is necessary to relinquish the fundamental idea that the decay of hydration forces is a constant determined by the properties of the aqueous medium. Alternatively, one can retain that fundamental idea if one recognizes the possibility that polar group hydration has an attractive component to it. In the latter view, that attractive component originates from interbilayer hydrogen-bonded water bridges between apposing bilayer surfaces, arising from correlation of zwitterionic or other complementary polar groups or from factors that affect polar group solubility. The same Marcelja and Radic formalism that accounts so well for the repulsive component also leads to an estimate of the attractive one. We suggest that the full range of degrees of hydration and of interbilayer spacings observed for different neutral bilayers results in part from variable contributions of the attractive and repulsive hydration components.(ABSTRACT TRUNCATED AT 400 WORDS)     

Effects of protein perturbants on phospholipid bilayers

Series of alcohols, amides, ureas, and sulfoxides with increasingly longer hydrocarbon chains have been shown to lower progressively the thermal denaturation temperature of proteins. This effect is presumably due to a hydrophobic interaction between the solute and nonpolar domains of the protein. Theoretically, these interactions should occur between the solute and any macromolecular structure having a nonpolar region to which the solute has access. A recent review by Arakawa et al. has summarized evidence for such an interaction between organic solutes and proteins and suggested that these interactions are favored at higher temperatures. The present study investigates the effects of several classes of compounds on the stability of phospholipid vesicles. The results show that many compounds that are known to perturb protein function also destabilize phospholipid bilayers as reflected by solute-induced loss of vesicle contents.     

Calculations of the electrostatic potential adjacent to model phospholipid bilayers.

We used the nonlinear Poisson-Boltzmann equation to calculate electrostatic potentials in the aqueous phase adjacent to model phospholipid bilayers containing mixtures of zwitterionic lipids (phosphatidylcholine) and acidic lipids (phosphatidylserine or phosphatidylglycerol). The aqueous phase (relative permittivity, epsilon r = 80) contains 0.1 M monovalent salt. When the bilayers contain < 11% acidic lipid, the -25 mV equipotential surfaces are discrete domes centered over the negatively charged lipids and are approximately twice the value calculated using Debye-Hückel theory. When the bilayers contain > 25% acidic lipid, the -25 mV equipotential profiles are essentially flat and agree well with the values calculated using Gouy-Chapman theory. When the bilayers contain 100% acidic lipid, all of the equipotential surfaces are flat and agree with Gouy-Chapman predictions (including the -100 mV surface, which is located only 1 A from the outermost atoms). Even our model bilayers are not simple systems: the charge on each lipid is distributed over several atoms, these partial charges are non-coplanar, there is a 2 A ion-exclusion region (epsilon r = 80) adjacent to the polar headgroups, and the molecular surface is rough. We investigated the effect of these four factors using smooth (or bumpy) epsilon r = 2 slabs with embedded point charges: these factors had only minor effects on the potential in the aqueous phase.     

The binding of organic ions to phospholipid bilayers

The binding of organic anions and cations, mainly tetraphenylboride and tetraphenylarsonium, to phospholipid membranes has been studied using an NMR method. Binding is appreciable and is affected by cholesterol in the membrane and counterions in solution. The passage of the organic anions through the membrane has also been followed. These measurements indicate that it is naive to use organic anions to measure membrane potentials in a simple manner.     

Molecular dynamics simulations of charged and neutral lipid bilayers: treatment of electrostatic interactions

Molecular dynamics (MD) simulations complement experimental methods in studies of the structure and dynamics of lipid bilayers. The choice of algorithms employed in this computational method represents a trade-off between the accuracy and real calculation time. The largest portion of the simulation time is devoted to calculation of long-range electrostatic interactions. To speed-up evaluation of these interactions, various approximations have been used. The most common ones are the truncation of long-range interactions with the use of cut-offs, and the particle-mesh Ewald (PME) method. In this study, several multi-nanosecond cut-off and PME simulations were performed to establish the influence of the simulation protocol on the bilayer properties. Two bilayers were used. One consisted of neutral phosphatidylcholine molecules. The other was a mixed lipid bilayer consisting of neutral phosphatidylethanolamine and negatively charged phosphatidylglycerol molecules. The study shows that the cut-off simulation of a bilayer containing charge molecules generates artefacts; in particular the mobility and order of the charged molecules are vastly different from those determined experimentally. In the PME simulation, the bilayer properties are in general agreement with experimental data. The cut-off simulation of bilayers containing only uncharged molecules does not generate artefacts, nevertheless, the PME simulation gives generally better agreement with experimental data.     

Binding of imipramine to phospholipid bilayers using radioligand binding assay

Binding of the tricyclic antidepressant imipramine (IMI) to neutral and negatively charged lipid membranes was investigated using a radioligand binding assay combined with centrifugation or filtration. Lipid bilayers were composed of brain phosphatidylcholine (PC) and phosphatidylserine (PS). IMI binding isotherms were measured up to IMI concentration of 0.5 mmol/l. Due to electrostatic attraction, binding between the positively charged IMI and the negatively charged surfaces of PS membranes was augmented compared to binding to neutral PC membranes. After correction for electrostatic effects by means of the Gouy-Chapman theory, the binding isotherms were described both by surface partition coefficients and by binding parameters (association constants and binding capacities). It was confirmed that binding of IMI to model membranes is strongly affected by negatively charged phospholipids and that the binding is heterogeneous; in fact, weak surface adsorption and incorporation of the drug into the hydrophobic core of lipid bilayer can be seen and characterized. These results support the hypothesis suggesting that the lipid part of biological membranes plays a role in the mechanism of antidepressant action.     

Spontaneous interbilayer transfer of hexosylceramides between phospholipid bilayers

The kinetics of spontaneous transfer of various glucosyl- and galactosylceramides between 1-palmitoyl-2-oleoylphosphatidylcholine vesicles have been examined at 45 degrees C. Bovine brain galactosylceramides, kerasin and phrenosin, were found to transfer with biexponential kinetics. The kerasin fast pool is approximately 17% with a half-time of 29 h and the slow pool approximately 83% with a half-time of 2700 h. In contrast, semisynthetic N-palmitoylgalactosylceramide at the same temperature transfers with single-exponential kinetics with a half-time of 32 h. The half-time for N-lignoceroylgalactosylceramide under the same conditions proved to be greater than 3500 h. No concentration dependence for these half-times was found in the concentration range studied (0-10 mol%). Similar results were obtained for semisynthetic glucosylceramides. The biexponential kinetics observed for the two bovine brain ceramides, both of which are mixtures of short and long acyl chain molecules, are most probably a reflection of the strong dependence of transfer rate on acyl chain length. The very slow transfer rates of the long acyl chain hexosylceramides ensure that these molecules are lost very slowly, if at all, by spontaneous transfer from the external surface of plasma membranes; a result that is consistent with the putative biological role of glycosphingolipids.     

Measurement and modification of forces between lecithin bilayers.

We probe in two different ways the competing attractive and repulsive forces that create lamellar arrays of the phospholipid lecithin when in equilibrium with pure water. The first probe involves the addition of low molecular weight solutes, glucose and sucrose, to a system where the phospholipid is immersed in a large excess of water. Small solutes can enter the aqueous region between bilayers. Their effect is first to increase and then to decrease the separation between bilayers as sugar concentration increases. We interpret this waxing and waning of the lattice spacing in terms of the successive weakening and strengthening of the attractive van der Waals forces originally responsible for creation of a stable lattice. The second probe is an "osmotic stress method," in which very high molecular weight neutral polymer is added to the pure water phase but is unable to enter the multilayers. The polymer competes for water with the lamellar lattice, and thereby compresses it. From the resulting spacing (determined by X-ray diffraction) and the directly measured osmotic pressure, we find a force vs. distance curve for compressing the lattice (or, equivalently, the free energy of transfer to bulk water of water between bilayers. This method reveals a very strong, exponentially varying "hydration force" with a decay distance of about 2 A.     

Formation of supported phospholipid bilayers on molecular surfaces: role of surface charge density and electrostatic interaction

Electrostatic interaction is known to play important roles in the adsorption of charged lipids on oppositely charged surfaces. Here we show that, even for charge neutral (zwitterionic) lipids, electrostatic interaction is critical in controlling the adsorption and fusion of lipid vesicles to form supported phospholipid bilayers (SPBs) on surfaces. We use terminally functionalized alkanethiol self-assembled monolayers (SAMs) to systematically control the surface charge density. Charge neutral egg phophatidylcholine (eggPC) vesicles readily fuse into SPBs on either a positively charged 11-aminino-1-undecanethiol SAM or a negatively charged 10-carboxy-1-decanethiol SAM when the density of surface charge groups is > or = 80%. These processes depend critically on the buffer environment: fusion of adsorbed vesicles to form SPBs on each charged molecular surface does not occur when the molecular ion of the buffer used is of the opposite charge type. We attribute this to the high entropic repulsion (electric double layer repulsion) due to the large size of molecular counterions. On the other hand, such a critical dependence on buffer type is not observed when charged lipids are used. This study suggests the general importance of controlling electrostatic interaction in the formation of stable SPBs.     

Role of electrostatics in the binding of charged metallophthalocyanines to neutral and charged phospholipid membranes

Binding of the cationic tetra(tributylammoniomethyl)-substituted hydroxoaluminum phthalocyanine (AlPcN(4)) to bilayer lipid membranes was studied by fluorescence correlation spectroscopy (FCS) and intramembrane field compensation (IFC) methods. With neutral phosphatidylcholine membranes, AlPcN(4) appeared to bind more effectively than the negatively charged tetrasulfonated aluminum phthalocyanine (AlPcS(4)), which was attributed to the enhancement of the coordination interaction of aluminum with the phosphate moiety of phosphatidylcholine by the electric field created by positively charged groups of AlPcN(4). The inhibitory effect of fluoride ions on the membrane binding of both AlPcN(4) and AlPcS(4) supported the essential role of aluminum-phosphate coordination in the interaction of these phthalocyanines with phospholipids. The presence of negative or positive charges on the surface of lipid membranes modulated the binding of AlPcN(4) and AlPcS(4) in accord with the character (attraction or repulsion) of the electrostatic interaction, thus showing the significant contribution of the latter to the phthalocyanine adsorption on lipid bilayers. The data on the photodynamic activity of AlPcN(4) and AlPcS(4) as measured by sensitized photoinactivation of gramicidin channels in bilayer lipid membranes correlated well with the binding data obtained by FCS and IFC techniques. The reduced photodynamic activity of AlPcN(4) with neutral membranes violating this correlation was attributed to the concentration quenching of singlet excited states as proved by the data on the AlPcN(4) fluorescence quenching.     

Role of electrostatics in the binding of charged metallophthalocyanines to neutral and charged phospholipid membranes

Binding of the cationic tetra(tributylammoniomethyl)-substituted hydroxoaluminum phthalocyanine (AlPcN₄) to bilayer lipid membranes was studied by fluorescence correlation spectroscopy (FCS) and intramembrane field compensation (IFC) methods. With neutral phosphatidylcholine membranes, AlPcN₄ appeared to bind more effectively than the negatively charged tetrasulfonated aluminum phthalocyanine (AlPcS₄), which was attributed to the enhancement of the coordination interaction of aluminum with the phosphate moiety of phosphatidylcholine by the electric field created by positively charged groups of AlPcN₄. The inhibitory effect of fluoride ions on the membrane binding of both AlPcN₄ and AlPcS₄ supported the essential role of aluminum-phosphate coordination in the interaction of these phthalocyanines with phospholipids. The presence of negative or positive charges on the surface of lipid membranes modulated the binding of AlPcN₄ and AlPcS₄ in accord with the character (attraction or repulsion) of the electrostatic interaction, thus showing the significant contribution of the latter to the phthalocyanine adsorption on lipid bilayers. The data on the photodynamic activity of AlPcN₄ and AlPcS₄ as measured by sensitized photoinactivation of gramicidin channels in bilayer lipid membranes correlated well with the binding data obtained by FCS and IFC techniques. The reduced photodynamic activity of AlPcN₄ with neutral membranes violating this correlation was attributed to the concentration quenching of singlet excited states as proved by the data on the AlPcN₄ fluorescence quenching.