Lipid bilayer disruption by oligomeric α-synuclein depends on bilayer charge and accessibility of the hydrophobic core

Soluble oligomeric aggregates of α-synuclein have been implicated to play a central role in the pathogenesis of Parkinson's disease. Disruption and permeabilization of lipid bilayers by α-synuclein oligomers is postulated as a toxic mechanism, but the molecular details controlling the oligomer–membrane interaction are still unknown. Here we show that membrane disruption strongly depends on the accessibility of the hydrophobic membrane core and that charge interactions play an important but complex role. We systematically studied the influence of the physical membrane properties and solution conditions on lipid bilayer disruption by oligomers using a dye release assay. Varying the lipid headgroup composition revealed that membrane disruption only occurs for negatively charged bilayers. Furthermore, the electrostatic repulsion between the negatively charged α-synuclein and the negative surface charge of the bilayer inhibits vesicle disruption at low ionic strength. The disruption of negatively charged vesicles further depends on lipid packing parameters. Bilayer composition changes that result in an increased lipid headgroup spacing make vesicles more prone to disruption, suggesting that the accessibility of the bilayer hydrocarbon core modulates oligomer–membrane interaction. These data shed important new insights into the driving forces governing the highly debated process of oligomer–membrane interactions.     

Electrostatic interactions across a charged lipid bilayer

We present theoretical work in which the degree of electrostatic coupling across a charged lipid bilayer in aqueous solution is analyzed on the basis of nonlinear Poisson–Boltzmann theory. In particular, we consider the electrostatic interaction of a single, large macroion with the two apposed leaflets of an oppositely charged lipid bilayer where the macroion is allowed to optimize its distance to the membrane. Three regimes are identified: a weak and a high macroion charge regime, separated by a regime of close macroion–membrane contact for intermediate charge densities. The corresponding free energies are used to estimate the degree of electrostatic coupling in a lamellar cationic lipid–DNA complex. That is, we calculate to what extent the one-dimensional DNA arrays in a sandwich-like lipoplex interact across the cationic membranes. We find that, in spite of the low dielectric constant inside a lipid membranes, there can be a significant electrostatic contribution to the experimentally observed cross-bilayer orientational ordering of the DNA arrays. Our approximate analytical model is complemented and supported by numerical calculations of the electrostatic potentials and free energies of the lamellar lipoplex geometry. To this end, we solve the nonlinear Poisson–Boltzmann equation within a unit cell of the lamellar lipoplex using a new lattice Boltzmann method. Dedicated to Prof. K. Arnold on the occasion of his 65th birthday.     

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.     

A method to determine dielectric constants in nonhomogeneous systems: application to biological membranes

Continuum electrostatic models have had quantitative success in describing electrostatic-mediated phenomena on atomistic scales; however, there continues to be significant disagreement about how to assign dielectric constants in mixed, nonhomogeneous systems. We introduce a method for determining a position-dependent dielectric profile from molecular dynamics simulations. In this method, the free energy of introducing a test charge is computed two ways: from a free energy perturbation calculation and from a numerical solution to Poisson's Equation. The dielectric profile of the system is then determined by minimizing the discrepancy between these two calculations simultaneously for multiple positions of the test charge. We apply this method to determine the dielectric profile of a lipid bilayer surrounded by water. We find good agreement with dielectric models for lipid bilayers obtained by other approaches. The free energy of transferring an ion from bulk water to the lipid bilayer computed from the atomistic simulations indicates that large errors are introduced when the bilayer is represented as a single slab of low dielectric embedded in the higher-dielectric solvent. Significant improvement results from introducing an additional layer of intermediate dielectric (∼3) on each side of the low dielectric core extending from ∼12 Å to 18 Å. A small dip in transfer free energy just outside the lipid headgroups indicates the presence of a very high dielectric. These results have implications for the design of implicit membrane models and our understanding of protein-membrane interactions.     

Elastic, electrostatic and electrokinetic forces influencing membrane curvature

Many cellular and intracellular processes critically depend on membrane shape, but the shape generating mechanisms are still to be fully understood. In this study we evaluate how electrostatic/electrokinetic forces contribute to membrane curvature. Membrane bilayer had finite thickness and was either elastically anisotropic or anisotropic overall, but isotropic per sections (heads and tails). The physics of the situation was evaluated using a coupled system of elastic and electrostatic/electrokinetic (Poisson-Nernst-Planck) equations. The fixed charges present only on the upper membrane surface lead to the accumulation of counter-ions and depletion of co-ions that decay spatially very rapidly (Debye length<1nm), as does the potential and electric field. Spatially uneven electric field and the permittivity mismatch also induce charges at the membrane-solution interface, which are not fixed but influence the electrostatics nevertheless. Membrane bends due to - Coulomb force (caused by fixed membrane charges in the electric field) and the dielectric force (due to the non-uniform electric field and the permittivity mismatch between the membrane and the solution). Both act as membrane surface forces, and both depend supra-linearly on the fixed charge density. Regardless of sign of the fixed charges, the membrane bends toward the charged (upper) surface owing to the action of the Coulomb force, but this is opposed by the smaller dielectric force. The spontaneous membrane curvature becomes very pronounced at high fixed charge densities, leading to very small spontaneous radii (<50nm). In conclusion the electrostatic/electrokinetic forces contribute significantly to the membrane curvature.     

Electrostatic Exclusion of Neutral Solutes from Condensed DNA and Other Charged Phases

Motivated by experiments on condensed DNA phases in binary mixtures of water and a low-dielectric solute, we develop a theory for the electrostatic contribution to solute exclusion from a highly charged phase, within the continuum approximation of the medium. Because the electric field is maximum at the surface of each ion, the electrostatic energy is dominated by the Born energy; interactions between charges are of secondary importance. Neglecting interactions and considering only the competition between the Born energy and the free energy of mixing, we predict that low dielectric solutes are excluded from condensed DNA phases in water-cosolvent mixtures. This suggests that the traditional continuum electrostatic approach of modeling binary mixtures with a uniform dielectric constant needs to be modified. The linking of solute exclusion to solute dielectric properties also suggests a mechanism for predicting the electrostatic contribution to preferential hydration of polar and charged surfaces.     

Brownian dynamics simulation of the lateral distribution of charged membrane components

Brownian dynamics simulations were performed to study the contribution of electric interactions between charged membrane components to their lateral distribution in a two-dimensional viscous liquid (bilayer lipid membrane). The electrostatic interaction potential was derived from an analytical solution of the linearized Poisson-Boltzmann equation for point charges in an electrolyte solution — membrane — electrolyte solution system. Equilibrium as well as dynamic quantities were investigated. The lateral organization of membrane particles, modelled by mobile cylinders in a homogeneous membrane separating two electrolyte solutions was described by spatial distribution functions, diffusion coefficients and cluster statistics. Disorder, local order and crystal-like arrangements were observed as a function of the particle charge, the closest possible distances between the charges and the particle density. The simulations revealed that the system is very sensitive to the position of the charges with respect to the electrolyte solution — membrane interface. Electrostatic interactions of charges placed directly on the membrane surface were almost negligible, whereas deeper charges demonstrated pronounced interaction. Biologically relevant parameters corresponded at most to local and transient ordering. It was found that lateral electric forces can give rise to a preferred formation of clusters with an even number of constituents provided that the closest possible charge-charge distances are small. It is concluded that lateral electrostatic interactions can account for local particle aggregations, but their impact on the global arrangement and movement of membrane components is limited. Correspondence to: D. Walther     

Photogating of ionic currents across lipid bilayers. Electrostatics of ions and dipoles inside the membrane.

The conductances of the lipophilic ions tetraphenylboride and tetraphenylphosphonium across a lipid bilayer can be increased or decreased, i.e., gated, by the photoformation of closed-shell metalloporphyrin cations within the bilayer. The gating can be effected by pulsed or continuous light or by chemical oxidants. At high concentrations of lipophilic anions where the dark conductance is saturated due to space charge in the bilayer, the photogated conductance can increase 15-fold. The formation of porphyrin cations allows the conductance to increase to its nonspace charge limited value. Conversely, the decrease of conductance in the light of phosphonium cations diminishes toward zero as the dark conductance becomes space charge limited. We present electrostatic models of the space charge limited conductance that accurately fit the data. One model includes an exponentially varying dielectric constant for the polar regions of the bilayer that allows an analytical solution to the electrostatic problem. The exponential variation of the dielectric constant effectively screens the potential and implies that the inside and outside of real dielectric interfaces can be electrically isolated from one another. The charge density, the distance into the membrane of the ions, about one-quarter of its thickness, and the dielectric constant at that position are determined by these models. These calculations indicate that there is insufficient porphyrin charge density to cancel the boride ion space charge and the following article proposes a novel ion chain mechanism to explain these effects. These models indicate that the positive potential arising from oriented carbonyl ester groups, previously used to explain the 10(3)-fold larger conductance of hydrophobic anions over cations, is smaller than previously estimated. However, the synergistic movement of the positive choline group into the membrane can account for the large positive potential.     

Interaction of charged and uncharged calcium channel antagonists with phospholipid membranes. Binding equilibrium, binding enthalpy, and membrane location.

The membrane location and the binding mechanism of two Ca2+ channel antagonists, amlodipine and nimodipine, in pure lipid membranes were investigated with deuterium and phosphorus-31 nuclear magnetic resonance, with thermodynamic methods such as high-sensitivity titration calorimetry, and by measuring the membrane surface charge via the zeta-potential. The two drugs exhibit quite different physical-chemical properties. The noncharged nimodipine is strongly hydrophobic, and selective deuteration of the lipid membrane reveals a homogeneous distribution of nimodipine across the whole hydrocarbon layer, but no interaction at the lipid headgroup level. The membrane behavior of the amiphiphilic amlodipine (electric charge z = +1) is distinctly more complex. Deuterium magnetic resonance demonstrates that amlodipine adopts a well-defined position in the bilayer membrane. In particular, the charged ethanolamine side group of amlodipine is located near the water-lipid interface, interacting with the dipoles of the headgroup region according to a nonspecific, electrostatic mechanism and inducing a reorientation of the phosphocholine dipoles toward the water phase. At the level of the hydrocarbon segment, the nonpolar ring system of amlodipine interacts specifically with the cis double bond of the membrane lipid, forming a weak association complex. With increasing amlodipine concentration the deuterium signal of the cis double bond gradually loses intensity, a phenomenon previously observed only in related studies on protein-lipid interactions. The binding equilibrium of amlodipine to phosphatidylcholine membranes was studied by measuring the electrophoretic mobility of lipid vesicles and with a centrifugation assay. Hydrophobic interactions of the nonpolar ring systems and electrostatic repulsions at the membrane surface contribute to the binding energy.(ABSTRACT TRUNCATED AT 250 WORDS)     

Equilibrium surface charge distribution in phospholipid vesicles. I. Method of calculation

This paper presents a method of calculation of the surface charge equilibrium distribution between the two surfaces of a spherically closed phospholipid bilayer suspended in aqueous electrolyte solution. The net surface charge is supposed to be provided by the ionized polar groups of the phospholipid molecules. Its equilibrium distribution is found by minimization of the free electrostatic energy. The procedure of minimization utilizes the solution of the Poisson-Boltzmann equation which describes the double electric layers of the membrane and an expression for the membrane potential derived under the assumption of absence of charges in the membrane phase. An analytical solution of the problem in the range of validity of the linearized Poisson-Boltzman equation is obtained. It is shown that in this case an equilibrium transmembrane potential exists, and the surface charge density is greater at the outer surface of the vesicle.     

Anchoring of a monotopic membrane protein: the binding of prostaglandin H2 synthase-1 to the surface of a phospholipid bilayer

Prostaglandin H₂ synthases (PGHS-1 and -2) are monotopic peripheral membrane proteins that catalyse the synthesis of prostaglandins in the arachidonate cascade. Picot et al. (1994) proposed that the enzyme is anchored to one leaflet of the bilayer by a membrane anchoring domain consisting of a right-handed spiral of amphipathic helices (residues 73–116) forming a planar motif. Two different computational approaches are used to examine the association of the PGHS-1 membrane anchoring domain with a membrane via the proposed mechanism. The electrostatic contribution to the free energy of solvation is obtained by solving numerically the finite-difference Poisson equation for the protein attached to a membrane represented as a planar slab of low dielectric. The nonpolar cavity formation and van der Waals contributions to the solvation free energy are assumed to be proportional to the water accessible surface area. Based on the optimum position determined from the continuum solvent model, two atomic models of the PGHS-1 anchoring domain associated with an explicit dimyristoylphosphatidylcholine (DMPC) bilayer differing by the thickness of the membrane bilayer were constructed. A total of 2 ns molecular dynamics simulation were performed to study the details of lipid- protein interactions at the microscopic level. In the simulations the lipid hydrocarbon chains interacting with the anchoring domain assume various shapes, suggesting that the plasticity of the membrane is significant. The hydrophobic residues in the membrane side of the helices interact with the hydrophobic membrane core, while the positively charged residues interact with the lipid polar headgroups to stabilize the anchoring of the membrane domain to the upper half of the bilayer. The phosphate headgroup of one DMPC molecule disposed at the center of the spiral formed by helices A, B, C and D interacts strongly with Arg120, a residue on helix D that has previously been identified as being important in the activity of PGHS-1. In the full enzyme structure, this position corresponds to the entrance of a long hydrophobic channel leading to the cyclooxygenase active site. These observations provide insights into the association of the arachidonic acid substrate to the cyclooxygenase active site of PGHS-1. Received: 20 December 1999 / Revised version: 26 March 2000 / Accepted: 26 March 2000     

Optical and electrical properties of thin monoolein lipid bilayers

Summary Monoolein lipid bilayers were formed using a monolayer transfer technique and from dispersions of monoolein in squalene, triolein, 1-chlorodecane and 1-bromodecane. Measurements of optical reflectance and electrical capacitance were used to determine the thickness and dielectric constant of the bilayers. The thickness of the hydrocarbon region of the five bilayer systems ranged from 2.5 to 3.0 nm. Two of the bilayer systems (made from 1-chlorodecane and 1-bromodecane solvents) had a high dielectric constant (2.8 to 2.9) whereas the other bilayer systems had dielectric constants close to that of pure hydrocarbons (2.2). The charge-pulse technique was used to study the transport kinetics of three lipophilic ions and two ion carrier complexes in the bilayers. For the low dielectric constant bilayers, the transport of the lipophilic ions tetraphenylborate, tetraphenylarsonium and dipicrylamine was governed mainly by the thickness of the hydrocarbon region of the bilayer whereas the transport of the ion-carrier complexes proline valinomycin-K⁺ and valinomycin-Rb⁺ was nearly independent of thickness. This is consistent with previous studies on thicker monoolein bilayers. The transport of lipophilic anions across bilayers with a high dielectric constant was 20 to 50 times greater than expected on the basis of thickness alone. This agrees qualitatively with predictions based on Born charging energy calculations. High dielectric constant bilayers were three times more permeable to the proline valinomycin-K⁺ complex than were low dielectric constant bilayers but were just as permeable as low dielectric constant bilayers to the valinomycin-Rb⁺ complex.     

Identification of the hydrophobic thickness of a membrane protein using fluorescence spectroscopy: studies with the mechanosensitive channel MscL

The hydrophobic thickness of a membrane protein is an important parameter, defining how the protein sits within the hydrocarbon core of the lipid bilayer that surrounds it in a membrane. Here we show that Trp scanning mutagenesis combined with fluorescence spectroscopy can be used to define the hydrophobic thickness of a membrane protein. The mechanosensitive channel of large conductance (MscL) contains two transmembrane alpha-helices, of which the second (TM2) is lipid-exposed. The region of TM2 that spans the hydrocarbon core of the bilayer when MscL is reconstituted into bilayers of dioleoylphosphatidylcholine runs from Leu-69 to Leu-92, giving a hydrophobic thickness of ca. 25 A. The results obtained using Trp scanning mutagenesis were confirmed using Cys residues labeled with the N-methyl-amino-7-nitroben-2-oxa-1,3-diazole [NBD] group; both fluorescence emission maxima and fluorescence lifetimes for the NBD group are sensitive to solvent dielectric constant over the range (2-40) thought to span the lipid headgroup region of a lipid bilayer. Changing phospholipid fatty acyl chain lengths from C14 and C24 results in no significant change for the fluorescence of the interfacial residues, suggesting very efficient hydrophobic matching between the protein and the surrounding lipid bilayer.     

Strong electrolyte continuum theory solution for equilibrium profiles, diffusion limitation, and conductance in charged ion channels.

The solution for the ion flux through a membrane channel that incorporates the electrolyte nature of the aqueous solution is a difficult theoretical problem that, until now, has not been properly formulated. The difficulty arises from the complicated electrostatic problem presented by a high dielectric aqueous channel piercing a low dielectric lipid membrane. The problem is greatly simplified by assuming that the ratio of the dielectric constant of the water to that of the lipid is infinite. It is shown that this is a good approximation for most channels of biological interest. This assumption allows one to derive simple analytical expressions for the Born image potential and the potential from a fixed charge in the channel, and it leads to a differential equation for the potential from the background electrolyte. This leads to a rigorous solution for the ion flux or the equilibrium potential based on a combination of the Nernst-Planck equation and strong electrolyte theory (i.e., Gouy-Chapman or Debye-Huckel). This approach is illustrated by solving the system of equations for the specific case of a large channel containing fixed negative charges. The following characteristics of this channels are discussed: anion and mono- and divalent cation conductance, saturation of current with increasing concentration, current-voltage relationship, influence of location and valence of fixed charge, and interaction between ions. The qualitative behavior of this channel is similar to that of the acetylcholine receptor channel.     

Membrane thickness and acyl chain length

It appears reasonable to expect that the primary result of a change in the length of the acyl chains within a lipid bilayer is a similar change in the bilayer thickness. In the present communication we draw attention to the somewhat more complicated effects which are found experimentally for phosphatidylcholine bilayers as the hydrocarbon chain is varied from twelve to eighteen carbons in length. The major change in dimension which occurs with variation in acyl chain length is the area occupied per molecule rather than the bilayer thickness. The same effect is seen with solute hydrocarbon such as hexane which partition into the membrane and cause only a small variation in membrane thickness but a large increase in the molecular area of the lipid. The origin of this effect arises from the almost isotropic distribution of the additional hydrocarbon to the lipid core of the membrane.     

Membrane bending is critical for the stability of voltage sensor segments in the membrane

The interaction between membrane proteins and the surrounding membrane is becoming increasingly appreciated for its role in regulating protein function, protein localization, and membrane morphology. In particular, recent studies have suggested that membrane deformation is needed to stably accommodate proteins harboring charged amino acids in their transmembrane (TM) region, as it is energetically prohibitive to bury charge in the hydrophobic core of the bilayer. Unfortunately, current computational methods are poorly equipped for describing such deformations, as atomistic simulations are often too short to observe large-scale membrane reorganization and most continuum approaches assume a flat membrane. Previously, we developed a method that overcomes these shortcomings by using elasticity theory to characterize equilibrium membrane distortions in the presence of a TM protein, while using traditional continuum electrostatic and nonpolar energy models to determine the energy of the protein in the membrane. Here, we linked the elastostatics, electrostatics, and nonpolar numeric solvers to permit the calculation of energies for nontrivial membrane deformations. We then coupled this procedure to a robust search algorithm that identifies optimal membrane shapes for a TM protein of arbitrary chemical composition. This advance now permits us to explore a host of biological phenomena that were beyond the scope of our original method. We show that the energy required to embed charged residues in the membrane can be highly nonadditive, and our model provides a simple mechanical explanation for this nonadditivity. Our results also predict that isolated voltage sensor segments do not insert into rigid membranes, but membrane bending dramatically stabilizes these proteins in the bilayer despite their high charge content. Additionally, we use the model to explore hydrophobic mismatch with regard to nonpolar peptides and mechanosensitive channels. Our method is in quantitative agreement with molecular dynamics simulations at a tiny fraction of the computational cost.     

Calculation of the electric potential in the active site cleft due to alpha-helix dipoles

A macroscopic dielectric model has been used to set up the electrostatic equation for the protein-solvent system. A numerical method of solution has been applied, enabling calculation of the electric potential outside a protein due to the charges within the protein. The glycolytic enzyme phosphoglycerate mutase, which is an α/β protein binding negatively charged substrates, has been studied. Modelling the helix dipoles with positive and negative charges shows that the α-helical structure could stabilize negatively charged substrates in the active site cleft of an enzyme with an energy of a few kT.     

Electrostatic free energy and spontaneous curvature of spherical charged layered membrane

A biophysical model for the equilibrium curvature of a composite membrane element is derived taking into account the mechanical bilayer properties and the adjacent charged protein layers. The minimum of the total free energy density with respect to the curvature of such a membrane curved was estimated from the sum of the electrostatic free energy density of the charges of the membrane and the elastic surface energy density due to bending the lipid bilayer membrane. It was shown that the equilibrium curvature, i.e. the spontaneous curvature, of such a charged composite sandwich-like membrane depends inversely on the bending stiffness of the lipid membrane itself and directly on the charge amount inside and outside the membrane to the second power. Furthermore the geometric and electrostatic structure of the protein layers and the physico-chemical environment conditions are involved. Corresponding to the model developed a "standard RBC" membrane element has a negative spontaneous curvature, accounting for a discocyte RBC shape. The shape change from a discocyte to a more stomatocytic shape (increase in the negative spontaneous curvature) after reducing the charges in the glycocalyx is also explained within this model.