Thomas Heimburg

Effect of sodium chloride on a lipid bilayer

Electrostatic interactions govern structural and dynamical properties of membranes and can vary considerably with the composition of the aqueous buffer. We studied the influence of sodium chloride on a pure POPC lipid bilayer by fluorescence correlation spectroscopy experiments and molecular dynamics simulations. Increasing sodium chloride concentration was found to decrease the self-diffusion of POPC lipids within the bilayer. Self-diffusion coefficients calculated from the 100 ns simulations agree with those measured on a millisecond timescale, suggesting that most of the relaxation processes relevant for lipid diffusion are faster than the simulation timescale. As the dominant effect, the molecular dynamics simulations revealed a tight binding of sodium ions to the carbonyl oxygens of on average three lipids leading to larger complexes with reduced mobility. Additionally, the bilayer thickens by approximately 2 A, which increases the order parameter of the fatty acyl chains. Sodium binding alters the electrostatic potential, which is largely compensated by a changed polarization of the aqueous medium and a lipid dipole reorientation.

Mechanical aspects of membrane thermodynamics. Estimation of the mechanical properties of lipid membranes close to the chain melting transition from calorimetry.

Changes in the internal energy of lipids with temperature are related to both lipid volume and area changes. Close to the chain melting transition of lipid bilayers volume and enthalpy fluctuations generally follow proportional functions. This makes it possible to calculate the relationship between membrane excess heat capacity with lipid volume, area compressibility and the membrane bending modulus, if the area fluctuations of the two monolayers are assumed to be mainly decoupled. Thus, compressibility and elasticity display pronounced maxima at the chain melting transition. These maxima can also be related to pronounced minima of the sound velocity in the lipid transition range, which were found in ultrasonic experiments. In the present study heat capacity profiles and volume changes were obtained. The compressibilities and the bending modulus were then deduced from the specific heat. The relevance of these findings for structural transitions and for the curvature dependence of heat capacities is discussed.

A Model for the Lipid Pretransition: Coupling of Ripple Formation with the Chain-Melting Transition

Below the thermotropic chain-melting transition, lipid membrane c(P) traces display a transition of low enthalpy called the lipid pretransition. It is linked to the formation of periodic membrane ripples. In the literature, these two transitions are usually regarded as independent events. Here, we present a model that is based on the assumption that both pretransition and main transition are caused by the same physical effect, namely chain melting. The splitting of the melting process into two peaks is found to be a consequence of the coupling of structural changes and chain-melting events. On the basis of this concept, we performed Monte Carlo simulations using two coupled monolayer lattices. In this calculation, ripples are considered to be one-dimensional defects of fluid lipid molecules. Because lipids change their area by approximately 24% upon melting, line defects are the only ones that are topologically possible in a triangular lattice. The formation of a fluid line defect on one monolayer leads to a local bending of the membrane. Geometric constraints result in the formation of periodic patterns of gel and fluid domains. This model, for the first time, is able to predict heat capacity profiles, which are comparable to the experimental c(P) traces that we obtained using calorimetry. The basic assumptions are in agreement with a large number of experimental observations.

Investigation of secondary and tertiary structural changes of cytochrome c in complexes with anionic lipids using amide hydrogen exchange measurements: an FTIR study

The structure of cytochrome c bound to anionic lipid membranes composed of dimyristoyl, dipalmitoyl, or dioleoyl phosphatidylglycerols, or of bovine heart cardiolipin, has been investigated by Fourier transform infrared spectroscopy. Only small changes in secondary structure, as registered by the amide I band of cytochrome c, were observed upon binding at temperatures below that of denaturation of the protein, and these were not coupled to the thermotropic phase transitions of the lipid. The denaturation temperature of the protein decreased by approximately 25-30 degrees upon binding, in a progression which correlated with that of the lipid phase transition temperatures, being approximately 7 degrees lower for complexes with dioleoyl than with dipalmitoyl phosphatidylglycerol. Large changes in the amide proton exchange characteristics, as monitored by the spectral shifts in the amide I band of the protein in D2O, were observed on binding cytochrome c to the lipid membranes. For the slowly exchanging population, the amide deuteration rates of the free protein were nearly independent of temperature, whereas those of the bound protein increased by up to two orders of magnitude over the temperature range from 10 to 40 degrees C. In addition, the extent of exchange differed between the bound and unbound protein. A structural transition in the bound protein was detected as a discontinuous step in Arrhenius plots of the deuterium exchange rates which occurred at a temperature in the region of 22 to 29 degrees C, depending on the lipid, far below that of denaturation. The temperature of this transition was determined by the physical state of the lipid, being 7 degrees lower for the lipids in the fluid state than for those in the gel state, and, for complexes with dimyristoyl phosphatidylglycerol, occurred at an intermediate temperature, being controlled by the lipid chain-melting transition at 27-28 degrees C. These results provide evidence for a coupling of the tertiary structure of the membrane-bound protein with the physical state of the membrane lipids.

Lipid ion channels

The interpretation of electrical phenomena in biomembranes is usually based on the assumption that the experimentally found discrete ion conduction events are due to a particular class of proteins called ion channels while the lipid membrane is considered being an inert electrical insulator. The particular protein structure is thought to be related to ion specificity, specific recognition of drugs by receptors and to macroscopic phenomena as nerve pulse propagation. However, lipid membranes in their chain melting regime are known to be highly permeable to ions, water and small molecules, and are therefore not always inert. In voltage-clamp experiments one finds quantized conduction events through protein-free membranes in their melting regime similar to or even undistinguishable from those attributed to proteins. This constitutes a conceptual problem for the interpretation of electrophysiological data obtained from biological membrane preparations. Here, we review the experimental evidence for lipid ion channels, their properties and the physical chemistry underlying their creation. We introduce into the thermodynamic theory of membrane fluctuations from which the lipid channels originate. Furthermore, we demonstrate how the appearance of lipid channels can be influenced by the alteration of the thermodynamic variables (e.g., temperature, pressure, tension and chemical potentials) in a coherent description that is free of parameters. This description leads to pores that display dwell times closely coupled to the fluctuation lifetime via the fluctuation-dissipation theorem. Drugs as anesthetics and neurotransmitters are shown to influence the channel likelihood and their lifetimes in a predictable manner. We also discuss the role of proteins in influencing the likelihood of lipid channel formation.

The Temperature Dependence of Lipid Membrane Permeability, its Quantized Nature, and the Influence of Anesthetics

We investigate the permeability of lipid membranes for fluorescence dyes and ions. We find that permeability reaches a maximum close to the chain melting transition of the membranes. Close to transitions, fluctuations in area and compressibility are high, leading to an increased likelihood of spontaneous lipid pore formation. Fluorescence correlation spectroscopy reveals the permeability for rhodamine dyes across 100-nm vesicles. Using fluorescence correlation spectroscopy, we find that the permeability of vesicle membranes for fluorescence dyes is within error proportional to the excess heat capacity. To estimate defect size we measure the conductance of solvent-free planar lipid bilayer. Microscopically, we show that permeation events appear as quantized current events very similar to those reported for channel proteins. Further, we demonstrate that anesthetics lead to a change in membrane permeability that can be predicted from their effect on heat capacity profiles. Depending on temperature, the permeability can be enhanced or reduced. We demonstrate that anesthetics decrease channel conductance and ultimately lead to blocking of the lipid pores in experiments performed at or above the chain melting transition. Our data suggest that the macroscopic increase in permeability close to transitions and microscopic lipid ion channel formation are the same physical process.

Protein surface-distribution and protein-protein interactions in the binding of peripheral proteins to charged lipid membranes

The binding of native cytochrome c to negatively charged lipid dispersions of dioleoyl phosphatidylglycerol has been studied over a wide range of ionic strengths. Not only is the strength of protein binding found to decrease rapidly with increasing ionic strength, but also the binding curves reach an apparent saturation level that decreases rapidly with increasing ionic strength. Analysis of the binding isotherms with a general statistical thermodynamic model that takes into account not only the free energy of the electrostatic double layer, but also the free energy of the surface distribution of the protein, demonstrates that the apparent saturation effects could arise from a competition between the out-of-plane binding reaction and the lateral in-plane interactions between proteins at the surface. It is found that association with nonlocalized sites results in binding isotherms that display the apparent saturation effect to a much more pronounced extent than does the Langmuir adsorption isotherm for binding to localized sites. With the model for nonlocalized sites, the binding isotherms of native cytochrome c can be described adequately by taking into account only the entropy of the surface distribution of the protein, without appreciable enthalpic interactions between the bound proteins. The binding of cytochrome c to dioleoyl phosphatidylglycerol dispersions at a temperature at which the bound protein is denatured on the lipid surface, but is nondenatured when free in solution, has also been studied. The binding curves for the surface-denatured protein differ from those for the native protein in that the apparent saturation at high ionic strength is less pronounced. This indicates the tendency of the denatured protein to aggregate on the lipid surface, and can be described by the binding isotherms for nonlocalized sites only if attractive interactions between the surface-bound proteins are included in addition to the distributional entropic terms. Additionally, it is found that the binding capacity for the native protein is increased at low ionic strength to a value that is greater than that for complete surface coverage, and that corresponds more closely to neutralization of the effective charge (determined from the ionic strength dependence), rather than of the total net charge, on the protein. Electron spin resonance experiments with spin-labeled lipids indicate that this different mode of binding arises from a penetration or disturbance of the bilayer surface by the protein that may alleviate the effects of in-plane interactions under conditions of strong binding.

Binding of Peripheral Proteins to Mixed Lipid Membranes: Effect of Lipid Demixing upon Binding

Binding isotherms have been determined for the association of horse heart cytochrome c with dioleoyl phosphatidylglycerol (DOPG)/dioleoyl phosphatidylcholine (DOPC) bilayer membranes over a range of lipid compositions and ionic strengths. In the absence of protein, the DOPG and DOPC lipids mix nearly ideally. The binding isotherms have been analyzed using double layer theory to account for the electrostatics, either the Van der Waals or scaled particle theory equation of state to describe the protein surface distribution, and a statistical thermodynamic formulation consistent with the mass-action law to describe the lipid distribution. Basic parameters governing the electrostatics and intrinsic binding are established from the binding to membranes composed of anionic lipid (DOPG) alone. Both the Van der Waals and scaled particle equations of state can describe the effects of protein distribution on the DOPG binding isotherms equally well, but with different values of the maximum binding stoichiometry (13 lipids/protein for Van der Waals and 8 lipids/protein for scaled particle theory). With these parameters set, it is then possible to derive the association constant, Kr, of DOPG relative to DOPC for surface association with bound cytochrome c by using the binding isotherms obtained with the mixed lipid membranes. A value of Kr (DOPG:DOPC) = 3.3-4.8, depending on the lipid stoichiometry, is determined that consistently describes the binding at different lipid compositions and different ionic strengths. Using the value of Kr obtained it is possible to derive the average in-plane lipid distribution and the enhancement in protein binding induced by lipid redistribution using the statistical thermodynamic theory.

Relaxation kinetics of lipid membranes and its relation to the heat capacity

We investigated the relaxation behavior of lipid membranes close to the chain-melting transition using pressure jump calorimetry with a temperature accuracy of approximately 10(-3) K. We found relaxation times in the range from seconds up to about a minute, depending on vesicular state. The relaxation times are within error proportional to the heat capacity. We provide a statistical thermodynamics theory that rationalizes the close relation between heat capacity and relaxation times. It is based on our recent finding that enthalpy and volume changes close to the melting transition are proportional functions.

ON THE ACTION POTENTIAL AS A PROPAGATING DENSITY PULSE AND THE ROLE OF ANESTHETICS

The Hodgkin-Huxley model of nerve pulse propagation relies on ion currents through specific resistors called ion channels. We discuss a number of classical thermodynamic findings on nerves that are not contained in this classical theory. Particularly striking is the finding of reversible heat changes, thickness and phase changes of the membrane during the action potential. Data on various nerves rather suggest that a reversible density pulse accompanies the action potential of nerves. Here, we attempted to explain these phenomena by propagating solitons that depend on the presence of cooperative phase transitions in the nerve membrane. These transitions are, however, strongly influenced by the presence of anesthetics. Therefore, the thermodynamic theory of nerve pulses suggests a explanation for the famous Meyer-Overton rule that states that the critical anesthetic dose is linearly related to the solubility of the drug in the membranes.