Modification of the erythrocyte membrane dielectric constant by alcohols

Summary Aliphatic alcohols are found to stimulate the transmembrane fluxes of a hydrophobic cation (tetraphenylarsonium, TPA) and anion (AN-12) 5–20 times in red blood cells. The results are analyzed using the Born-Parsegian equation (Parsegian, A., 1969,Nature (London) 221:844–846), together with the Clausius-Mossotti equation to calculate membrane dielectric energy barriers. Using established literature values of membrane thickness, native membrane dielectric constant, TPA ionic radius, and alcohol properties (partition coefficient, molar volume, dielectric constant), the TPA permeability data is predicted remarkably well by theory. If the radius of AN-12 is taken as 1.9 Å, its permeability in the presence of butanol is also described by our analysis. Further, the theory quantitatively accounts for the data of Gutknecht and Tosteson (Gutknecht, J., Tosteson, D.C., 1970,J. Gen. Physiol. 55:359–374) covering alcohol-induced conductivity changes of 3 orders of magnitude in artificial bilayers. Other explanations including perturbations of membrane fluidity, surface charge, membrane thickness, and dipole potential are discussed. However, the large magnitude of the stimulation, the more pronounced effect on smaller ions, and the acceleration of both anions and cations suggest membrane dielectric constant change as the primary basis of alcohol effects.     

Oestradiol changes the dielectric structure of bilayer membranes

The addition of the hormone Oestradiol to Phosphatidylcholine-Cholesterol membrane changes the frequency dependence of the membrane impedance. It increases severalfold the electrical admittance of the polar regions and consequently provides a conducting shunt from the hydrocarbon region to the aqueous phase.     

Molecular dynamics simulation of melittin in a dimyristoylphosphatidylcholine bilayer membrane.

Molecular dynamics trajectories of melittin in an explicit dimyristoyl phosphatidylcholine (DMPC) bilayer are generated to study the details of lipid-protein interactions at the microscopic level. Melittin, a small amphipathic peptide found in bee venom, is known to have a pronounced effect on the lysis of membranes. The peptide is initially set parallel to the membrane-solution interfacial region in an alpha-helical conformation with unprotonated N-terminus. Solid-state nuclear magnetic resonance (NMR) and polarized attenuated total internal reflectance Fourier transform infrared (PATIR-FTIR) properties of melittin are calculated from the trajectory to characterize the orientation of the peptide relative to the bilayer. The residue Lys7 located in the hydrophobic moiety of the helix and residues Lys23, Arg24, Gln25, and Gln26 at the C-terminus hydrophilic form hydrogen bonds with water molecules and with the ester carbonyl groups of the lipids, suggesting their important contribution to the stability of the helix in the bilayer. Lipid acyl chains are closely packed around melittin, contributing to the stable association with the membrane. Calculated density profiles and order parameters of the lipid acyl chains averaged over the molecular dynamics trajectory indicate that melittin has effects on both layers of the membrane. The presence of melittin in the upper layer causes a local thinning of the bilayer that favors the penetration of water through the lower layer. The energetic factors involved in the association of melittin at the membrane surface are characterized using an implicit mean-field model in which the membrane and the surrounding solvent are represented as structureless continuum dielectric material. The results obtained by solving the Poisson-Bolztmann equation numerically are in qualitative agreement with the detailed dynamics. The influence of the protonation state of the N-terminus of melittin is examined. After 600 ps, the N-terminus of melittin is protonated and the trajectory is continued for 400 ps, which leads to an important penetration of water molecules into the bilayer. These observations provide insights into how melittin interacts with membranes and the mechanism by which it enhances their lysis.     

Effect of different drugs on a cytochrome c--phospholipid bilayer membrane.

Liposomes were prepared from dipalmitoyl phosphatidylcholine and dicetylphosphate and their interaction with the extrinsic membrane protein cytochrome c examined in terms of changes in 22Na permeability, electrophoretic mobility, protein binding, and motion of an incorporated spin label. The amount of cytochrome c bound displays no significant temperature dependence over the temperature range studied (from 30 to 55 degrees C) whereas cytochrome c causes an increase in 22Na efflux only above the phospholipid phase transition temperature. Interaction of the protein with the lipid vesicles causes no significant disturbance in the bilayer interior as monitored by the motion of the incorporated spin probe. The drugs 2,4-dinitrophenol and ethacrynic acid, both of which increase the magnitude of the vesicle negative charge, enhance both cytochrome c binding and its effect on 22Na permeability. In contrast, the local anesthetic dibucaine, which induces a positive surface charge on these liposomes, reduces both protein binding and the protein-induced increase in 22Na efflux. Finally, the chemicals butylated hydroxytoluene, 2-tert-butylphenol, and tert-butylbenzene, all of which cause early 'melting' of the phospholipid fatty acyl chains, block the capacity of cytochrome c to enhance 22Na permeability while having no effect on its binding to the vesicles.     

Phospholipid exchange between bilayer membrane vesicles.

The turbidity of lipid vesicles, freshly prepared by sonicating purified dimyristoyllecithin (DML) in dilute KCl solutions, was measured as a function of time at various temperatures. A sharp maximum in the rate of increase of turbidity is found just above the crystal:liquid-crystal phase transition temperature (Tm). The initial rate of turbidity increase is first order with respect to DML concentration. Electron and light microscopy reveal large vesicles which are not present before incubation or after incubation at temperatures far from the Tm. When temperature, rather than time, is the independent variable, a sharp drop in turbidity is seen at the Tm. The magnitude of this drop and the temperature at which it occurs were used to measure the rate of lipid transfer between vesicles composed of different lipids. A mixture of DML vesicles and dipalmitoyllecithin (DPL) vesicles exhibits sharp drops in turbidity at 24 and 41 degrees, the corresponding Tm's. With time, the magnitude of the transition at 24 degrees decreases while that which was originally at 41 degrees moves to lower temperatures and increases in magnitude. At equilibrium there is a single transition at 32.5 degrees characteristic of vesicles composed of equimolar DPL and DML. The rate at which equilibrium is approached increases at around 24 degrees and again around 41 degrees. These observations indicate that vesicles are in equilibrium with monomolecular lipid, the concentration of the latter being higher the shorter the lipid acyl group or the smaller the vesicle. DML molecules are therefore lost from small vesicles to large vesicles (DML system) or lost from DML vesicles to DML-DPL vesicles (mixed system). When DML vesicles containing a few percent brain gangliosides were studied, different behavior was observed; the initial rate of increase of turbidity becomes second order in lipid concentration, and the rate constant increases with increasing concentrations of KCl. The kinetic order, coupled with the fact that electrolyte reduces intervesicle electrostatic repulsion, argues that in this situation the mechanism of vesicle growth requires vesicle collision.     

Molecular dynamics simulations of a fluid bilayer of dipalmitoylphosphatidylcholine at full hydration, constant pressure, and constant temperature.

Molecular dynamics simulations of 500 ps were performed on a system consisting of a bilayer of 64 molecules of the lipid dipalmitoylphosphatidylcholine and 23 water molecules per lipid at an isotropic pressure of 1 atm and 50 degrees C. Special attention was devoted to reproduce the correct density of the lipid, because this quantity is known experimentally with a precision better than 1%. For this purpose, the Lennard-Jones parameters of the hydrocarbon chains were adjusted by simulating a system consisting of 128 pentadecane molecules and varying the Lennard-Jones parameters until the experimental density and heat of vaporization were obtained. With these parameters the lipid density resulted in perfect agreement with the experimental density. The orientational order parameter of the hydrocarbon chains agreed perfectly well with the experimental values, which, because of its correlation with the area per lipid, makes it possible to give a proper estimate of the area per lipid of 0.61 +/- 0.01 nm2.     

Thickness dependence of monoglyceride bilayer membrane conductance.

We have studied the conductance properties of unmodified monoglyceride membranes as a function of monoglyceride chain length. As membrane thickness decreases from 31 to 20 nm, the steepness of the current-voltage (I-V) curve increases from 80 mV per e-fold current increase to 52 mV per e-fold current increase. The zero-voltage conductance increases more than 1,000-fold and the apparent activation energy of conductance decreases from 18.4 to 14.2 kcal/mol. We have analyzed our results using both the Nernst-Planck equation and absolute rate theory. Both approaches are consistent with our results and give consistent values for the parameters describing the I-V curves. We conclude that both the surface ion concentration and the distance from the surface of the membrane at which the energy of an ion rises appreciably above its value in solution (position of the barrier) are invariant with thickness.     

Phosphatidylethanolamine-phosphatidylglycerol bilayer as a model of the inner bacterial membrane

Phosphatidylethanolamine (PE) and phosphatidylglycerol (PG) are the main lipid components of the inner bacterial membrane. A computer model for such a membrane was built of palmitoyloleoyl PE (POPE) and palmitoyloleoyl PG (POPG) in the proportion 3:1, and sodium ions (Na+) to neutralize the net negative charge on each POPG (POPE-POPG bilayer). The bilayer was simulated for 25 ns. A final 10-ns trajectory fragment was used for analyses. In the bilayer interfacial region, POPEs and POPGs interact readily with one another via intermolecular hydrogen (H) bonds and water bridges. POPE is the main H-bond donor in either PEPE or PEPG H-bonds; PGPG H-bonds are rarely formed. Almost all POPEs are H-bonded and/or water bridged to either POPE or POPG but PE-PG links are favored. In effect, the atom packing in the near-the-interface regions of the bilayer core is tight. Na+ does not bind readily to lipids, and interlipid links via Na+ are not numerous. Although POPG and POPE comprise one bilayer, their bilayer properties differ. The average surface area per POPG is larger and the average vertical location of the POPG phosphate group is lower than those of POPE. Also, the alkyl chains of POPG are more ordered and less densely packed than the POPE chains. The main conclusion of this study is that in the PE-PG bilayer PE interacts more strongly with PG than with PE. This is a likely molecular-level event behind a regulating mechanism developed by the bacteria to control its membrane permeability and stability consisting in changes of the relative PG/PE concentration in the membrane.     

Electrostatics of phosphoinositide bilayer membranes. Theoretical and experimental results.

We made fluorescence, electron paramagnetic resonance (EPR), electrophoretic mobility, and ionizing electrode measurements to study the effect of the monovalent lipid phosphatidylinositol (PI) and the trivalent lipid phosphatidylinositol 4,5-bisphosphate (PIP2) on the electrostatic potential adjacent to bilayer membranes. When the membranes were formed from mixtures of PI and the zwitterionic lipid phosphatidylcholine (PC), the Gouy-Chapman-Stern (GCS) theory described adequately the dependence of potential on distance (0, 1, 2 nm) from the membrane, mole % negative lipid, and [KCI]. Furthermore, all EPR and fluorescence probes reported identical surface potentials with a PC/PI membrane. With PC/PIP2 membranes, however, the anionic (coion) probes reported less negative potentials than the cationic (counterion) probes; the deviations from the GCS theory were greater for the coions than the counterions. Discreteness-of-charge theories based on the Poisson-Boltzmann equation incorrectly predict that deviations from the GCS theory should be greater for counterions than for coions. We discuss a consistent statistical mechanical theory that takes into account three effects ignored in the GCS theory: the finite size of the ions in the double layer, the electrical interaction between pairs of ions (correlation effects), and the mobile discrete nature of the surface charges. This theory correctly predicts that deviations from GCS theory should be negligible for monovalent lipids, significant for trivalent lipids, and greater for coions than for counterions.     

Weak acid transport across bilayer lipid membrane in the presence of buffers. Theoretical and experimental pH profiles in the unstirred layers.

This paper presents a simple model to describe experimental data on weak acid transport across planar bilayer lipid membrane separating two buffered solutions. The model takes into account multiple proton-transfer reactions occurring in the unstirred layers (ULs) adjacent to the membrane. Differential equations of the model are shown to be reduced to a set of nonlinear algebraic equations. Since the latter equations depend monotonically on unknown variables, they can be easily solved numerically, using bisection method. For the particular system studied experimentally (with acetate as the weak acid and TRIS+MES as the buffer mixture) pH profiles in the ULs are calculated from the model. These results are compared with experimental data obtained using pH microelectrode. The agreement between theoretical and experimental pH profiles is found to be satisfactory. The most pronounced deviations are observed at the UL/bulk solution boundary. To obtain a better correlation between the theoretical and experimental results, two other, less idealized models are considered. They take into account, respectively, (a) the electric field arising in the ULs from ion diffusion and (b) finiteness of the rates of proton-transfer reactions. However, both acetate membrane fluxes and pH profiles in the ULs computed from these models are found to be close to those of the simple model. One can thus conclude that the difference between experimental and theoretical pH profiles is due to the inconsistency of the generally accepted model of the "unstirred layer", assuming the existence of a strict boundary between the regions of "pure diffusion" and "ideal stirring".     

Amphiphile-induced tubular budding of the bilayer membrane

Amphiphile-induced tubular budding of the erythrocyte membrane was studied using transmission electron microscopy. No chiral patterns of the intramembraneous particles were found, either on the cylindrical buds, or on the tubular nanoexovesicles. In agreement with these observations, the tubular budding may be explained by in-plane ordering of anisotropic membrane inclusions in the buds where the difference between the principal membrane curvatures is very large. In contrast to previously reported theories, no direct external mechanical force is needed to explain tubular budding of the bilayer membrane.     

Alignment of cholesterol in the membrane bilayer

Freeze-fracture replicas of filipin-treated samples of guinea pig colon mucosa reveal areas in the membrane of the goblet cell granules labeled by filipin-cholesterol complexes (FCC) intermingled with regions patterned by lines. The FCC and lines are arranged in an approximately rhombic pattern. Other membranes of the same cell or of other cells display either FCC only, aligned and occasionally ordered in rhombs, lines only, with a similar pattern, or randomly distributed FCC. Optical diffraction was used to analyze and compare replicas of membranes with ordered FCC and lines, as well as randomly distributed FCC. The results demonstrate that all these structures are reciprocally related through a common distribution pattern in the membrane. This observation supports the assumption that cholesterol has a preferential ordered distribution within the membrane bilayer.     

On stability of circular hole in membrane bilayer.

It was observed recently that nonionic surfactant octaethyleneglycol dodecylether (Cl2E8) decreases threshold for irreversible electroporation in membrane bilayers. In accordance, it is shown theoretically in this work that anisotropic C12E8 membrane inclusions may stabilize circular hole in a flat membrane segment.     

A theoretical study of the dielectric constant of protein

The dielectric properties of a protein molecule were investigated by calculating a 'local dielectric constant' with the aid of normal mode analysis. This local dielectric constant was calculated from the electronic polarization of atoms and the orientational polarization of local dipoles. The former was obtained from atomic polarizations of the whole atoms in a protein molecule. The latter was determined from the fluctuation-dissipation theorem. The degree of dipole fluctuation was calculated from the positional fluctuation of each atom obtained by the normal mode analysis. Assuming a minimum volume for a continuum model, the resulting local dielectric constants ranged from 1 to 20 inside the protein.     

Experimental evaluation of the effective dielectric constant of proteins

Chemical modifications that alter the net charge of residues in reduction-oxidation proteins influence the redox potential of the protein by changing the electrostatic potential at the redox center. If the locations of the modified charges are known, the shift in redox potential may be used to determine the effective dielectric constant for the interactions between the redox center and modified residues. From the shift in redox potential upon charge neutralization of specific lysines in the hemoprotein cytochrome c, an effective dielectric constant of approximately 50 is calculated for the electrostatic interaction between the modified lysines and heme iron in the native protein.     

Fusion of phospholipid vesicles with planar phospholipid bilayer membranes. II. Incorporation of a vesicular membrane marker into the planar membrane

Fusion of multilamellar phospholipid vesicles with planar phospholipid bilayer membranes was monitored by the rate of appearance in the planar membrane of an intrinsic membrane protein present in the vesicle membranes. An essential requirement for fusion is an osmotic gradient across the planar membrane, with the cis side (the side containing the vesicles) hyperosmotic to the opposite (trans) side; for substantial fusion rates, divalent cation must also be present on the cis side. Thus, the low fusion rates obtained with 100 mM excess glucose in the cis compartment are enhanced orders of magnitude by the addition of 5-10 mM CaCl2 to the cis compartment. Conversely, the rapid fusion rates induced by 40 mM CaCl2 in the cis compartment are completely suppressed when the osmotic gradient (created by the 40 mM CaCl2) is abolished by addition of an equivalent amount of either CaCl2, NaCl, urea, or glucose to the trans compartment. We propose that fusion occurs by the osmotic swelling of vesicles in contact with the planar membrane, with subsequent rupture of the vesicular and planar membranes in the region of contact. Divalent cations catalyze this process by increasing the frequency and duration of vesicle-planar membrane contact. We argue that essentially this same osmotic mechanism drives biological fusion processes, such as exocytosis. Our fusion procedure provides a general method for incorporating and reconstituting transport proteins into planar phospholipid bilayer membranes.