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.     

Study of the interaction of myoglobin with lipid bilayer membranes by potentiodynamic method

For modeling the interaction of myoglobin with mitochondrial membranes, the adsorption of different ligand forms, the physiologically active reduced MbO₂ and inactive oxidized met-Mb, on one of the surfaces of artificial bilayer lipid membrane (BLM) was studied using a potentiodynamic technique known as the “capacity minimization” method. As mitochondrial membranes are negatively charged, BLM of the negatively charged palmitoyl-2-oleyl-phosphatidyl glycerol (POPG) and neutral soybean phosphatidylcholine (lecithin) were used. It is shown that both myoglobins strongly interact with BLM in the pH range 6–8. The dependence of the potential difference between cis-and trans-surfaces of the lipid membrane (ΔE, mV) on the protein concentration is characteristic of the Langmuir adsorption isotherm, and the saturation level (ΔE _max, mV) corresponds to monolayer of myoglobin. The protein adsorption is essentially electrostatic in nature, as adsorption activity increases sharply in the case of the membrane from POPG: ∼15-fold in the case of MbO₂ and ∼2.5 times for met-Mb. The parameters of the MbO₂ and met-Mb adsorption on BLM of lecithin and POPG do not change in the pH 6–8 range. It can be assumed that the anionic groups of phospholipids associate with the cationic groups of the protein, the charge state of those does not change in the pH 6–8 range. The most likely candidates for interaction with phospholipids of BLM are invariant lysines and arginines in the environment of the myoglobin heme cavity.     

Effect of membrane tension on the electric field and dipole potential of lipid bilayer membrane

The dipole potential of lipid bilayer membrane controls the difference in permeability of the membrane to oppositely charged ions. We have combined molecular dynamics (MD) simulations and experimental studies to determine changes in electric field and electrostatic potential of 1,2-dioleoyl-sn-glycero-3-phosphocholine (DOPC) lipid bilayer in response to applied membrane tension. MD simulations based on CHARMM36 force field showed that electrostatic potential of DOPC bilayer decreases by ~45mV in the physiologically relevant range of membrane tension values (0 to 15dyn/cm). The electrostatic field exhibits a peak (~0.8×10(9)V/m) near the water/lipid interface which shifts by 0.9? towards the bilayer center at 15dyn/cm. Maximum membrane tension of 15dyn/cm caused 6.4% increase in area per lipid, 4.7% decrease in bilayer thickness and 1.4% increase in the volume of the bilayer. Dipole-potential sensitive fluorescent probes were used to detect membrane tension induced changes in DOPC vesicles exposed to osmotic stress. Experiments confirmed that dipole potential of DOPC bilayer decreases at higher membrane tensions. These results are suggestive of a potentially new mechanosensing mechanism by which mechanically induced structural changes in the lipid bilayer membrane could modulate the function of membrane proteins by altering electrostatic interactions and energetics of protein conformational states.     

Modelling of proteins in membranes

This review describes some recent theories and simulations of mesoscopic and microscopic models of lipid membranes with embedded or attached proteins. We summarize results supporting our understanding of phenomena for which the activities of proteins in membranes are expected to be significantly affected by the lipid environment. Theoretical predictions are pointed out, and compared to experimental findings, if available. Among others, the following phenomena are discussed: interactions of interfacially adsorbed peptides, pore-forming amphipathic peptides, adsorption of charged proteins onto oppositely charged lipid membranes, lipid-induced tilting of proteins embedded in lipid bilayers, protein-induced bilayer deformations, protein insertion and assembly, and lipid-controlled functioning of membrane proteins.     

Lipid bilayers cushioned with polyelectrolyte-based films on doped silicon surfaces

Silicon semiconductors with a thin surface layer of silica were first modified with polyelectrolytes (polyethyleneimine, polystyrene sulfonate and poly(allylamine)) via a facile layer-by-layer deposition approach. Subsequently, lipid vesicles were added to the preformed polymeric cushion, resulting in the adsorption of intact vesicles or fusion and lipid bilayer formation. To study involved interactions we employed optical reflectometry, electrochemical impedance spectroscopy and fluorescent recovery after photobleaching. Three phospholipids with different charge of polar head groups, i.e. 1,2-dioleoyl-sn-glycero-3-phosphocholine (DOPC), 1,2-dioleoyl-sn-glycero-3-phospho-l-serine (DOPS) and 1,2-dioleoyl-3-trimethylammonium-propane (DOTAP) were used to prepare vesicles with varying surface charge. We observed that only lipid vesicles composed from 1:1 (mole:mole) mixture of DOPC/DOPS have the ability to fuse onto an oppositely charged terminal layer of polyelectrolyte giving a lipid bilayer with a resistance of >100 kΩ. With optical reflectometry we found that the vesicle surface charge is directly related to the amount of mass adsorbed onto the surface. An interesting observation was that zwitterionic polar head groups of DOPC allow the adsorption on both positively and negatively charged surfaces. As found with fluorescent recovery after photobleaching, positively charged surface governed by the presence of poly(allylamine) as the terminal layer resulted in intact DOPC lipid vesicles adsorption whereas in the case of a negatively charged silica surface formation of lipid bilayers was observed, as expected from literature.     

Lipid demixing and protein-protein interactions in the adsorption of charged proteins on mixed membranes

The adsorption free energy of charged proteins on mixed membranes, containing varying amounts of (oppositely) charged lipids, is calculated based on a mean-field free energy expression that accounts explicitly for the ability of the lipids to demix locally, and for lateral interactions between the adsorbed proteins. Minimization of this free energy functional yields the familiar nonlinear Poisson-Boltzmann equation and the boundary condition at the membrane surface that allows for lipid charge rearrangement. These two self-consistent equations are solved simultaneously. The proteins are modeled as uniformly charged spheres and the (bare) membrane as an ideal two-dimensional binary mixture of charged and neutral lipids. Substantial variations in the lipid charge density profiles are found when highly charged proteins adsorb on weakly charged membranes; the lipids, at a certain demixing entropy penalty, adjust their concentration in the vicinity of the adsorbed protein to achieve optimal charge matching. Lateral repulsive interactions between the adsorbed proteins affect the lipid modulation profile and, at high densities, result in substantial lowering of the binding energy. Adsorption isotherms demonstrating the importance of lipid mobility and protein-protein interactions are calculated using an adsorption equation with a coverage-dependent binding constant. Typically, at bulk-surface equilibrium (i.e., when the membrane surface is "saturated" by adsorbed proteins), the membrane charges are "overcompensated" by the protein charges, because only about half of the protein charges (those on the hemispheres facing the membrane) are involved in charge neutralization. Finally, it is argued that the formation of lipid-protein domains may be enhanced by electrostatic adsorption of proteins, but its origin (e.g., elastic deformations associated with lipid demixing) is not purely electrostatic.     

Balance of Electrostatic and Hydrophobic Interactions in the Lysis of Model Membranes by E. coliα-Haemolysin

The relative weight of electrostatic interactions and hydrophobic forces in the process of membrane disruption caused by E. coliα-haemolysin (HlyA) has been studied with a purified protein preparation and a model system consisting of large unilamellar vesicles loaded with water-soluble fluorescent probes. Vesicles were prepared in buffers of different ionic strengths, or pHs, and the net surface charge of the bilayers was also modified by addition of negatively (e.g., phosphatidylinositol) or positively (e.g., stearylamine) charged lipids. The results can be interpreted in terms of a multiple equilibrium in which α-haemolysin may exist: aggregated HlyA ⇄ monomeric HlyA ⇄ membrane-bound HlyA. In these equilibria both electrostatic and hydrophobic forces are significant. Electrostatic forces become substantial under certain circumstances, e.g., membrane binding when bilayer and protein have opposite electric charges. Protein adsorption to the bilayer is more sensitive to electrostatic forces than membrane disruption itself. In the latter case, the irreversible nature of protein insertion may overcome electrostatic repulsions. Also of interest is the complex effect of pH on the degree of aggregation of an amphipathic toxin like α-haemolysin, since pH changes are not only influencing the net protein charge but may also be inducing protein conformational transitions shown by changes in the protein intrinsic fluorescence and in its susceptibility to protease digestion, that appear to regulate the presence of hydrophobic patches at the surface of the molecule, thus modifying the ability of the toxin to either aggregate or become inserted in membranes. Received: 29 October 1996/Revised: 4 February 1997     

On the adsorption of ions at charged sites in an electric field: II

A study is made of the adsorption of one kind of monovalent positive ion at a long chain of alternating monovalent negative fixed charged (“lattice”) and uncharged (“interstitial”) sites both of one type in an electric field. Considering only nearest neighbor interactions an expression is obtained for the grand partition function. The fractions of sites of both types which are occupied and unoccupied are determined. It is shown that an equilibrium constant can be defined for the adsorption of ions at oppositely charged sites.     

Interaction of low-molecular-weight chitosan with mimic membrane studied by electrochemical methods and surface plasmon resonance

Chitosan has shown its potential as a non-viral gene carrier and an adsorption enhancer for subsequent drug delivery to cells. These results showed that chitosan acted as a membrane perturbant. However, there is currently a lack of direct experimental evidence of this membrane perturbance effect, especially for chitosans with low molecular weight (LMW). In this report, the interaction between a lipid (didodecyl dimethylammonium bromide; DDAB) bilayer and chitosan with molecular weight (MW) of 4200 Da was studied with cyclic voltammetry (CV), electrochemical impedance spectroscopy and surface plasmon resonance (SPR). A lipid bilayer was formed by fusion of oppositely charged lipid vesicles on a mercaptopropionic acid (MPA)-modified gold surface to mimic a cell membrane. The results showed that the LMW chitosan could disrupt the lipid bilayer, and the effect seemed to be in a concentration-dependent manner.     

The Role of Extramembranous Cytoplasmic Termini in Assembly and Stability of the Tetrameric K+-Channel KcsA

Membrane-active alcohol 2,2,2-trifluoroethanol has been proven to be an attractive tool in the investigation of the intrinsic stability of integral membrane protein complexes by taking K⁺-channel KcsA as a suitable and representative ion channel. In the present study, the roles of both cytoplasmic N and C termini in channel assembly and stability of KcsA were determined. The N terminus (1–18 residues) slightly increased tetramer stability via electrostatic interactions in the presence of 30 mol.% acidic phosphatidylglycerol (PG) in phosphatidylcholine lipid bilayer. Furthermore, the N terminus was found to be potentially required for efficient channel (re)assembly. In contrast, truncation of the C terminus (125–160 residues) greatly facilitated channel reversibility from either a partially or a completely unfolded state, and this domain was substantially involved in stabilizing the tetramer in either the presence or absence of PG in lipid bilayer. These studies provide new insights into how extramembranous parts play their crucial roles in the assembly and stability of integral membrane protein complexes.     

Intermolecular interactions of influenza M1 proteins on the model lipid membrane surface: A study using the inner field compensation method

M1 protein binding to the lipid bilayer membrane (BLM) was recorded by the inner field compensation technique as a change of the boundary potential. After the protein was added to the bulk solution, the M1 adsorption produced a slow increase in boundary potential to a stationary value that was reached within the time period dependent on the quantity of the added protein. The stationary value of the potential grew with the decrease of pH or KCl concentration in the medium and was higher in the presence of negatively charged lipids in the BLM. It was shown that the potential growth with the decrease of pH is due to an increase of M1 molecule charge and not due to the increase of the M1 surface concentration or to the change of lipid charge. As the potential did not change after the removal of the protein from the bulk solution, we consider the protein adsorption on the BLM irreversible. The obtained results suggest that the protein adsorption is influenced by both electrostatic and hydrophobic interactions of M1 molecules with each other and with lipid membrane. We offer a mechanism of dissociation of the viral shell formed by M1 matrix protein. The protein shell is destabilized due to electrostatic repulsion of protein molecules caused by the increase of their positive charge.     

Molecular dynamics simulations of local anesthetic articaine in a lipid bilayer

In order to investigate structural and dynamical properties of local anesthetic articaine in a model lipid bilayer, a series of molecular dynamics simulations have been performed. Simulations were carried out for neutral and charged (protonated) forms of articaine inserted in fully hydrated dimyristoylphosphatidylcholine (DMPC) lipid bilayer. For comparison purpose, a fully hydrated DMPC bilayer without articaine was also simulated. The length of each simulation was 200 ns. Various properties of the lipid bilayer systems in the presence of both charged and uncharged forms of articaine taken at two different concentrations have been examined: membrane area per lipid, mass density distributions, order parameters, radial distribution functions, head group tilt, diffusion coefficients, electrostatic potential, etc, and compared with results of previous simulations of DMPC bilayer in the presence of lidocaine. It was shown that addition of both charged and neutral forms of articaine causes increase of the dipole electrostatic potential in the membrane interior.     

Interaction between K+ channel gate modifier hanatoxin and lipid bilayer membranes analyzed by molecular dynamics simulation

Hanatoxin (HaTx) is an ellipsoidal-shaped peptide that binds to the voltage sensor of voltage-dependent channels. Of physicochemical interest, HaTx has a “ring” of charged residues around its periphery and a hydrophobic protrusion. It has previously been postulated that HaTx binds to and functions on the surface of membranes, but a recent fluorescent-quenching study has implied a fairly deep positioning of HaTx in the lipid bilayer membrane. We carried out numerous molecular dynamic simulations of HaTx1, a well-studied variant of HaTx, in fully hydrated phospholipid bilayers. The system reproduced the surface-binding mode of HaTx1, in which HaTx1 resided in the extracellular side (outer) of the water/membrane interface with the hydrophobic patch of HaTx1 facing the membrane interior. On the other hand, analyses with various parameter settings suggested that the surface-binding mode was unstable because of the substantial attractive electrostatic force between HaTx1 and the lipid head groups of the inner (opposite) leaflet. Compared with this electrostatic force, the energetic cost for membrane deformation involving meniscus formation appeared to be small. In an attempt to interpret the quenching data, we consider the possibility of dimpling (meniscus formation) that brings HaTx1 inward (only ~0.7–0.8 nm above the bilayer center), while accounting for the flexibility of both leaflets of the membrane and the long-range interaction between positively charged residues of the membrane-bound peptide and the polar head groups of the opposite leaflet of the membrane. It is suggested that molecular dynamics simulations taking into account the flexibility of the membrane surface is potentially useful in interpreting the fluorescence-quenching data.Electronic Supplementary Material Supplementary material is available for this article at and is accessible for authorized users.     

DNA-lipid complexes: stability of honeycomb-like and spaghetti-like structures.

A molecular level theory is presented for the thermodynamic stability of two (similar) types of structural complexes formed by (either single strand or supercoiled) DNA and cationic liposomes, both involving a monolayer-coated DNA as the central structural unit. In the "spaghetti" complex the central unit is surrounded by another, oppositely curved, monolayer, thus forming a bilayer mantle. The "honeycomb" complex is a bundle of hexagonally packed DNA-monolayer units. The formation free energy of these complexes, starting from a planar cationic/neutral lipid bilayer and bare DNA, is expressed as a sum of electrostatic, bending, mixing, and (for the honeycomb) chain frustration contributions. The electrostatic free energy is calculated using the Poisson-Boltzmann equation. The bending energy of the mixed lipid layers is treated in the quadratic curvature approximation with composition-dependent bending rigidity and spontaneous curvature. Ideal lipid mixing is assumed within each lipid monolayer. We found that the most stable monolayer-coated DNA units are formed when the charged/neutral lipid composition corresponds (nearly) to charge neutralization; the optimal monolayer radius corresponds to close DNA-monolayer contact. These conclusions are also valid for the honeycomb complex, as the chain frustration energy is found to be negligible. Typically, the stabilization energies for these structures are on the order of 1 k(B)T/A of DNA length, reflecting mainly the balance between the electrostatic and bending energies. The spaghetti complexes are less stable due to the additional bending energy of the external monolayer. A thermodynamic analysis is presented for calculating the equilibrium lipid compositions when the complexes coexist with excess bilayer.     

Protein adsorption onto auto-assembled polyelectrolyte films

Surface modification by deposition of ordered protein systems constitutes one of the major objectives of bio-related chemistry and biotechnology. In this respect a concept has recently been reported aimed at fabricating multilayers by the consecutive adsorption of positively and negatively charged polyelectrolytes. We investigate the adsorption processes between polyelectrolyte multilayers and a series of positively and negatively charged proteins. The film buildup and adsorption experiments were followed by Scanning Angle Reflectometry (SAR). We find that proteins strongly interact with the polyelectrolyte film whatever the sign of the charge of both the multilayer and the protein. When charges of the multilayer and the protein are similar, one usually observes the formation of protein monolayers, which can become dense. We also show that when the protein and the multilayer become oppositely charged, the adsorbed amounts are usually larger and the formation of thick protein layers extending up to several times the largest dimension of the protein can be observed. Our results confirm that electrostatic interactions dominate protein/polyelectrolyte multilayer interactions.     

Mathematical model of the interaction of cations and anions in a channel

A mathematical model of the ionic channel permeable both to anions and cations is considered. The model takes into account the electrostatic interaction between oppositely charged ions and does not suppose single-file movement. An equation for zero-current potential is derived, which leads to the Goldman equation in the limit of low ion concentrations. The model is used to describe concentration relationships of zero-current potentials on a lipid bilayer with amphotericin B channels which cannot be described on the basis of the independence principle.     

Evidence that Electrostatic Interactions between Vesicle-associated Membrane Protein 2 and Acidic Phospholipids May Modulate the Fusion of Transport Vesicles with the Plasma Membrane

The juxtamembrane domain of vesicle-associated membrane protein (VAMP) 2 (also known as synaptobrevin2) contains a conserved cluster of basic/hydrophobic residues that may play an important role in membrane fusion. Our measurements on peptides corresponding to this domain determine the electrostatic and hydrophobic energies by which this domain of VAMP2 could bind to the adjacent lipid bilayer in an insulin granule or other transport vesicle. Mutation of residues within the juxtamembrane domain that reduce the VAMP2 net positive charge, and thus its interaction with membranes, inhibits secretion of insulin granules in β cells. Increasing salt concentration in permeabilized cells, which reduces electrostatic interactions, also results in an inhibition of insulin secretion. Similarly, amphipathic weak bases (e.g., sphingosine) that reverse the negative electrostatic surface potential of a bilayer reverse membrane binding of the positively charged juxtamembrane domain of a reconstituted VAMP2 protein and inhibit membrane fusion. We propose a model in which the positively charged VAMP and syntaxin juxtamembrane regions facilitate fusion by bridging the negatively charged vesicle and plasma membrane leaflets.     

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     

Cholesterol promotes the interaction of Alzheimer β-amyloid monomer with lipid bilayer

Cholesterol in the plasma membrane plays an important role in the pathogenesis of Alzheimer's disease, but the exact function of cholesterol in the regulation of amyloid-β (Aβ) generation, aggregation, and toxicity remains elusive. To gain insight into the bioactivity of cholesterol, we investigate the effect of cholesterol levels on the interaction of Aβ(1-42) monomer with the zwitterionic 1-palmitoyl-2-oleoylphosphatidylcholine (POPC) bilayer containing different mole fractions of cholesterol from χ=0, 0.2, to 0.4 using all-atom molecular dynamics simulations. Simulation results show that an increased cholesterol level alters the structure, dynamics, and surface chemistry of the lipid bilayer, leading to increased bilayer thickness, hydrophobic chain order, surface hydrophobicity, and decreased lipid mobility. All these effects significantly promote the binding of Aβ to the lipid bilayer. Mechanistically, the adsorption of Aβ on the bilayer is a cooperative process. First, charged residues act as anchors to establish the initial binding of Aβ to phosphate headgroups of the bilayer driven by electrostatic interactions, which further facilitates hydrophobic residues to reside on the bilayer. Once hydrophobic residues especially from C-terminus are locked on the bilayer, the interactions among charged residues, lipid bilayer, and calcium ions are optimized to provide additional attractive forces to stabilize Aβ adsorbed on or inserted into the lipid bilayer. Inclusion of cholesterol makes this binding process more energetically favorable. Upon adsorption on the bilayer, Aβ appears to preferentially adopt α-helical or unstructured conformation. This work supports that cholesterol acts as a promoter for Aβ--membrane interactions, which would facilitate Aβ aggregation and membrane insertion.     

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.