Curvature elasticity of multilamellar lipid bilayers close to the chain-melting transition

We directly measured curvature elasticity of dipalmitoylphosphatidylcholine multilamellar bilayers close to the chain-melting transition using the method of electric-field-induced bending deformation of the cylindrical tubes. The result shows that the bending modulus, kappa(c), decreases remarkably at temperatures close to the melting transition temperature. This reflects a softening of the bilayer resulted from the area fluctuations as predicted theoretically. However, the decrease of kappa(c) near the transition is far smaller than that predicted. This is due to the experimental method and the narrow transition width of the multilamellar bilayers. Nevertheless, the result obtained gives direct evidence of the kappa(c) reduction predicted for multilamellar membranes in the transition regime. Below about 41 degrees C, almost of all cylindrical tubes cannot response to the electric field, indicating a very large bending rigidity.     

Effects of the anesthetic dibucaine on the kinetics of the gel-liquid crystalline transition of dipalmitoylphosphatidylcholine multilamellar vesicles.

The effects of the anesthetic dibucaine on the relaxation kinetics of the gel-liquid crystalline transition of dipalmitoylphosphatidylcholine (DC16PC) multilamellar vesicles have been investigated using volume-perturbation calorimetry. The temperature and pressure responses to a periodic volume perturbation were measured in real time. Data collected in the time domain were subsequently converted into and analyzed in the frequency domain using Fourier series representations of the perturbation and response functions. The Laplace transform of the classical Kolmogorov-Avrami kinetic relation was employed to describe the relaxation dynamics in the frequency domain. The relaxation time of anesthetic-lipid mixtures, as a function of the fractional degree of melting, appears to be qualitatively similar to that of pure lipid systems, with a pronounced maximum, tau max, observed at a temperature corresponding to greater than 75% melting. The tau max decreases by a factor of approximately 2 as the nominal anesthetic/lipid mole ratio increases from 0 to 0.013 and exhibits no further change as the nominal anesthetic/lipid mole ratio is increased. However, the fractional dimensionality of the relaxation process decreases monotonically from slightly less than two to approximately one as the anesthetic/lipid mole ratio increases from 0 to 0.027. At higher ratios, the dimensionality appears to be less than one. These results are interpreted in terms of the classical kinetic theory and related to those obtained from Monte Carlo simulations. Specifically, low concentrations of dibucaine appear to reduce the average cluster size and cause the fluctuating lipid clusters to become more ramified. At the highest concentration of dibucaine, where n < 1, the system must be kinetically heterogeneous.     

The kinetics of interactions of bilirubin with lipid bilayers and with serum albumin.

Rate constants for the hydration of bilirubin bound to unilamellar bilayers of dioleoylphosphatidylcholine and albumin were measured by stopped-flow methods. Rate constants for association of bilirubin with these vesicles and albumin were calculated from measured rate constants for dissociation and the equilibrium binding constants of bilirubin and lipids or albumin. Rate constants for hydration (dissociation) for bilirubin bound to dioleoylphosphatidylcholine and albumin were 71 s-1 and 1.8 s-1 respectively. Rate constants for association were 4.0 10(7) s-1 and 1.1 10(9) M-1 s-1, respectively. Both rates for interactions of bilirubin with bilayers were essentially independent of temperature in the range 0-40 degrees C, indicating that barriers to entry and exit of bilirubin from bilayers were entropic. Rates of transbilayer movement of bilirubin in dioleoylphosphatidylcholine were too fast to resolve by measuring rates of hydration of bilirubin. Rate constants for hydration of bilirubin bound to bilayers with less avidity for bilirubin as compared with dioleoylphosphatidylcholine also were too fast to measure with stopped-flow methods. In addition to providing details of the energetic basis for interactions between bilirubin and membranes, the data allow for calculating the maximal rates at which bilirubin could transfer spontaneously from sites on albumin in blood to the interior of cells. The data show, in this regard, that this rate is 10-50 fold faster than measured rates of uptake of bilirubin by intact liver.     

Interaction of small peptides with lipid bilayers.

Molecular dynamics simulations of the tripeptide Ala-Phe-Ala-O-tert-butyl interacting with dimyristoylphosphatidylcholine lipid bilayers have been carried out. The lipid and aqueous environments of the peptide, the alkyl chain order, and the lipid and peptide dynamics have been investigated with use of density profiles, radial distribution functions, alkyl chain order parameter profiles, and time correlation functions. It appears that the alkyl chain region accommodates the peptides in the bilayer with minimal perturbation to this region. The peptide dynamics in the bilayer bound form has been compared with that of the free peptide in water. The peptide structure does not vary on the simulation time scale (of the order of hundreds of picoseconds) compared with the solution structure in which a random structure is observed.     

Segregation of chlorophyll a incorporated into lipid bilayers.

Absorption and fluorescence spectra are reported for chlorophyll a incorporated into a number of aqueous phospholipid dispersions. Absorption spectra show that in dipalmitoylphosphatidylcholine bilayers, monomeric and oligomeric forms of chlorophyll a are present in both the gel and liquid crystalline phases. The formation of aggregates of chlorophyll a is reflected in the fluorescence spectra by a marked concentration quenching. In bilayers conatining small proportions of chlorophyll a, a marked increase in aggregation occurs at the transition temperatures that can be detected calorimetrically. At higher concentrations (greater than 1 chlorophyll:100 lipid), the "pretransition" is abolished in the phosphatidylcholines, and the main transition is broadened, consistent with an orientation for the chlorophyll a with the chlorine ring in the head group region and the phytol chain in the fatty acid chain region of the bilayer. In mixtures of saturated and unsaturated lipids, there is no preferential segregation of the chlorophyll a into the unsaturated lipid.     

Transport kinetics of uncoupling proteins. Analysis of UCP1 reconstituted in planar lipid bilayers

According to alternative hypotheses, mitochondrial uncoupling protein 1 (UCP1) is either a proton channel ("buffering model") or a fatty acid anion carrier ("fatty acid cycling"). Transport across the proton channel along a chain of hydrogen bonds (Grotthus mechanism) may include fatty acid carboxyl groups or occur in the absence of fatty acids. In this work, we demonstrate that planar bilayers reconstituted with UCP1 exhibit an increase in membrane conductivity exclusively in the presence of fatty acids. Hence, we can exclude the hypothesis considering a preexisting H+ channel in UCP1, which does not require fatty acid for function. The augmented conductivity is nearly completely blocked by ATP. Direct application of transmembrane voltage and precise current measurements allowed determination of ATP-sensitive conductances at 0 and 150 mV as 11.5 and 54.3 pS, respectively, by reconstituting nearly 3 x 10(5) copies of UCP1. The proton conductivity measurements carried out in presence of a pH gradient (0.4 units) allowed estimation of proton turnover numbers per UCP1 molecule. The observed transport rate of 14 s-1 is compatible both with carrier and channel nature of UCP1.     

Iodide penetration into lipid bilayers as a probe of membrane lipid organization.

The quenching efficiency of iodide as a penetrating fluorescence quencher for a membrane-associated fluorophore was utilized to measure the molecular packing of lipid bilayers. The KI quenching efficiency of tryptophan-fluorescence from melittin incorporated in DMPC bilayer vesicles peaks at the phase transition temperature (24 degrees C) of DMPC, whereas acrylamide quenching efficiency does not depend on temperature. The ability of iodide to penetrate the hydrocarbon region of the bilayer was examined by measuring the fluorescence quenching of the pyrene-phosphatidylcholine incorporated into DMPC vesicles (pyrene was attached to the 10th carbon of the sn-2 chain). The quenching efficiency of pyrene by iodide again shows a maximum at the lipid phase transition. We conclude that iodide penetrates the membrane hydrocarbon region at phase transition through an increased number of bilayer defects. The magnitude of change in quenching efficiency of iodide during lipid phase transition provides a sensitive technique to probe the lipid organization in membranes.     

Kinetics of carrier-mediated ion transport in two new types of solvent-free lipid bilayers.

In contrast with the usual glyceryl-monooleate/decane (GMO-D) bilayer lipid membranes, new membranes, formed from a mixture of GMO in squalene (GMO-S) or from a mixture of GMO in triolein (GMO-T), seem to be almost solvent free. Our results from voltage-jump relaxation studies, using these "solvent-free" membranes with the homologue carriers, nonactin, monactin, dinactin, trinactin, and tetranactin, are compared with the corresponding ones for GMO-D membranes. With all homologues, solvent-free membranes show an increase of the free carrier translocation rate, ks, by a factor of 2.5, a decrease in the dissociation rate constant of the complex, kDi, by a factor of 1.5 and no significant change in its formation rate constant, kRi. However, the principal effect of the absence of solvent in these membranes is an increase by a factor of approximately 10 of the translocation rate constant for moving the complex across the membrane, kis. This increase varies regularly from a factor of 7-15 with decreasing carrier size, and is always larger for GMO-T than for GMO-S membranes. These solvent-free effects are interpreted in terms of modifications of electrostatic and hydrophobic energy profiles in the membrane.     

Concentration of oxygen in lipid bilayers using a spin-label method.

The concentration of oxygen in the hydrocarbon region of lipid bilayer has been determined using a novel electron spin resonance (ESR) nitroxide-radical spin-probe method. For dimyristoylphosphatidylcholine (DMPC), the partition coefficient above the main transition temperature is approximately 3. Rapid decrease to 0.2 occurs below the pretransition temperature indicating exclusion of oxygen in the crystalline phase. The differences of molar free energy, enthalpy, and entropy of mixing between water and lipid have been determined for each phase.     

Sodium channels in planar lipid bilayers. Channel gating kinetics of purified sodium channels modified by batrachotoxin

Single channel currents of sodium channels purified from rat brain and reconstituted into planar lipid bilayers were recorded. The kinetics of channel gating were investigated in the presence of batrachotoxin to eliminate inactivation and an analysis was conducted on membranes with a single active channel at any given time. Channel opening is favored by depolarization and is strongly voltage dependent. Probability density analysis of dwell times in the closed and open states of the channel indicates the occurrence of one open state and several distinct closed states in the voltage (V) range-120 mV less than or equal to V less than or equal to +120 mV. For V less than or equal to 0, the transition rates between stages are exponentially dependent on the applied voltage, as described in mouse neuroblastoma cells (Huang, L. M., N. Moran, and G. Ehrenstein. 1984. Biophysical Journal. 45:313-322). In contrast, for V greater than or equal to 0, the transition rates are virtually voltage independent. Autocorrelation analysis (Labarca, P., J. Rice, D. Fredkin, and M. Montal. 1985. Biophysical Journal. 47:469-478) shows that there is no correlation in the durations of successive open or closing events. Several kinetic schemes that are consistent with the experimental data are considered. This approach may provide information about the mechanism underlying the voltage dependence of channel activation.     

Interaction of fluorescently labeled pardaxin and its analogues with lipid bilayers.

Fluorescence measurements were used to monitor the interaction of the neurotoxin pardaxin and its analogues with membranes. Eight peptides were selectively labeled with the fluorophore 7-nitrobenz-2-oxa-1,3-diazole-4-yl, either at their N-terminal or at their C-terminal. No detectable changes in membrane permeability or hemolytic activity were observed upon modification. Upon the titration of solutions containing the different peptides with small unilamellar vesicles, the fluorescent emission spectra of 7-nitrobenz-2-oxa-1,3-diazole-4-yl-labeled pardaxin and its analogues, but not those of control peptides, displayed blue shifts in addition to enhanced intensities upon relocation of the probe to a more apolar environment. The results revealed that the N terminus of pardaxin is buried within the lipid bilayer while the C terminus is located at the bilayer's surface. Binding isotherms were obtained from the observed increases in the fluorescence emission yields, from which surface partition constants, in the range of 10(4) M-1, were in turn derived. The existence of an aggregation process was suggested by the shape of the binding isotherms. Furthermore, the results show good correlation between the incidence of aggregation and the ability of the different analogues to induce the release of relatively large molecules from vesicles. As such, our results suggest that the mechanism of pore formation employed by pardaxin and its analogues could be described by the "barrel stave" model.     

Structure stability of lytic peptides during their interactions with lipid bilayers.

In this work, molecular dynamics simulations were used to examine the consequences of a variety of analogs of cecropin A on lipid bilayers. Analog sequences were constructed by replacing either the N- or C-terminal helix with the other helix in native or reverse sequence order, by making palindromic peptides based on both the N- and C-terminal helices, and by deleting the hinge region. The structure of the peptides was monitored throughout the simulation. The hinge region appeared not to assist in maintaining helical structure but help in motion flexibility. In general, the N-terminal helix of peptides was less stable than the C-terminal one during the interaction with anionic lipid bilayers. Sequences with hydrophobic helices tended to regain helical structure after an initial loss while sequences with amphipathic helices were less able to do this. The results suggests that hydrophobic design peptides have a high structural stability in an anionic membrane and are the candidates for experimental investigation.     

Interactions of alpha-helices with lipid bilayers: a review of simulation studies

Membrane proteins, of which the majority seem to contain one or more alpha-helix, constitute approx. 30% of most genomes. A complete understanding of the nature of helix/bilayer interactions is necessary for an understanding of the structural principles underlying membrane proteins. This review describes computer simulation studies of helix/bilayer interactions. Key experimental studies of the interactions of alpha-helices and lipid bilayers are briefly reviewed. Surface associated helices are found in some membrane-bound enzymes (e.g. prostaglandin synthase), and as stages in the mechanisms of antimicrobial peptides and of pore-forming bacterial toxins. Transmembrane alpha-helices are found in most integral membrane proteins, and also in channels formed by amphipathic peptides or by bacterial toxins. Mean field simulations, in which the lipid bilayer is approximated as a hydrophobic continuum, have been used in studies of membrane-active peptides (e.g. alamethicin, melittin, magainin and dermaseptin) and of simple membrane proteins (e.g. phage Pf1 coat protein). All atom molecular dynamics simulations of fully solvated bilayers with transmembrane helices have been applied to: the constituent helices of bacteriorhodopsin; peptide-16 (a simple model TM helix); and a number of pore-lining helices from ion channels. Surface associated helices (e.g. melittin and dermaseptin) have been simulated, as have alpha-helical bundles such as bacteriorhodopsin and alamethicin. From comparison of the results from the two classes of simulation, it emerges that a major theoretical challenge is to exploit the results of all atom simulations in order to improve the mean field approach.     

Interaction of sporidesmin,a mycotoxin fromPithomyces chartarum,with lipid bilayers

Sporidesmin, a mycotoxin fromPithomyces chartarum is a hydrophobic molecule. It can therefore be easily incorporated in the cell membrane, where it is likely to cause changes in the bilayer organization and the properties of membrane proteins. In order to understand the redox behaviour of sporidesmin in a hydrophobic environment, we have investigated the effects of oxidized and reduced sporidesmin on the phase transition properties of bilayers and on the susceptibility of bilayers to pancreatic phospholipase A₂ (PLA₂). The changes induced by sporidesmin in the thermotropic phase transition profiles of dimyristoyl-sn-3-phosphatidyl choline (DMPC) bilayers were similar to those caused by solutes known to localize in the glycerol-backbone region of the lipid bilayer, suggesting a similar localization for oxidized and reduced sporidesmin. Neither form of toxin disrupt the bilayer or membrane organization even at relatively high mole fractions. At concentrations <10 mole% both forms partitioned equally well in the gel and liquid-crystalline phases, whereas at higher concentrations (30 mole%) reduced sporidesmin is preferentially localized in the liquid-crystalline phase. These effects of sporidesmin on the phase properties of DMPC vesicles were also reported by the fluorescence behavior of 10-pyrenedecanoic acid (PDA). The effects of oxidized and reduced sporidesmins on PLA₂ kinetics are consistent with their ability to perturb bilayer organisation.     

Interaction of a two-transmembrane-helix peptide with lipid bilayers and dodecyl sulfate micelles

To probe structural changes that occur when a membrane protein is transferred from lipid bilayers to SDS micelles, a fragment of bacteriorhodopsin containing transmembrane helical segments A and B was studied by fluorescence spectroscopy, molecular dynamics (MD) simulation, and stopped flow kinetics. In lipid bilayers, F?rster resonance energy transfer (FRET) was observed between tyrosine 57 on helix B and tryptophans 10 and 12 on helix A. FRET efficiency decreased substantially when the peptide was transferred to SDS. MD simulation showed no evidence for significant disruption of helix-helix interactions in SDS micelles. However, a cluster of water molecules was observed to form a hydrogen-bonded network with the phenolic hydroxyl group of tyrosine 57, which probably causes the disappearance of tyrosine-to-tryptophan FRET in SDS. The tryptophan quantum yield decreased in SDS, and the change occurred at nearly the same rate as membrane solubilization. The results provide a clear example of the importance of corroborating distance changes inferred from FRET by using complementary methods.     

Detection of transient capacitance increase associated with channel formation in lipid bilayers.

The insertion of alpha- and beta-latrotoxins and sea anemone (Radianthus macrodactilus) toxin into bilayer lipid membranes (BLMs) was investigated using the method of simultaneous conductance/capacitance measurement. All the toxins investigated induced capacitance changes which preceded toxin-induced conductance increases. The processes that may underlie the observed effect are discussed.     

Interactions of cholesterol with lipid bilayers: the preferred configuration and fluctuations.

The free energy difference associated with the transfer of a single cholesterol molecule from the aqueous phase into a lipid bilayer depends on its final location, namely on its insertion depth and orientation within the bilayer. We calculated desolvation and lipid bilayer perturbation contributions to the water-to-membrane transfer free energy, thus allowing us to determine the most favorable location of cholesterol in the membrane and the extent of fluctuations around it. The electrostatic and nonpolar contributions to the solvation free energy were calculated using continuum solvent models. Lipid layer perturbations, resulting from both conformational restrictions of the lipid chains in the vicinity of the (rigid) cholesterol backbone and from cholesterol-induced elastic deformations, were calculated using a simple director model and elasticity theory, respectively. As expected from the amphipathic nature of cholesterol and in agreement with the available experimental data, our results show that at the energetically favorable state, cholesterol's hydrophobic core is buried within the hydrocarbon region of the bilayer. At this state, cholesterol spans approximately one leaflet of the membrane, with its OH group protruding into the polar (headgroup) region of the bilayer, thus avoiding an electrostatic desolvation penalty. We found that the transfer of cholesterol into a membrane is mainly driven by the favorable nonpolar contributions to the solvation free energy, whereas only a small opposing contribution is caused by conformational restrictions of the lipid chains. Our calculations also predict a strong tendency of the lipid layer to elastically respond to (thermally excited) vertical fluctuations of cholesterol so as to fully match the hydrophobic height of the solute. However, orientational fluctuations of cholesterol were found to be accompanied by both an elastic adjustment of the surrounding lipids and by a partial exposure of the hydrophobic cholesterol backbone to the polar (headgroup) environment. Our calculations of the molecular order parameter, which reflects the extent of orientational fluctuations of cholesterol in the membrane, are in good agreement with available experimental data.     

Effect of the lipid phase transition on the kinetics of H+/OH− diffusion across phosphatidic acid bilayers

The kinetics of H⁺/OH^? diffusion across dimyristoyl phosphatidic acid bilayer membranes was measured by following the absorbance of the pH-sensitive indicator Cresol red ($\text{o-}\text{cresolsulfonphthalein}$) entrapped in single lamellar vesicles after rapidly changing the external pH in a stopped-flow apparatus. The H⁺/OH^?-permeability coefficient was found to be in the 10^?5 to 10^?3 cm·s^?1 range. The lipid phase transition has a strong influence on the permeation kinetics as the permeability coefficients in the liquid-crystalline phase are drastically higher. The permeability shows no maximum at the phase transition temperature as is the case for other ions, but displays a similar temperature dependence as water permeation. This is also reflected in the high activation energy of approx. 20 kcal/mol and supports the hypothesis (Nichols, J.W. and Deamer, D.W. (1980) Proc. Natl. Acad. Sci. U.S.A. 77, 2038–2042) of H⁺/OH^? permeation via hydrogen bonded water molecules. A second slower kinetic phase is also observed, where the permeation is obviously controlled by counterion diffusion. The temperature dependence of this slow process displays the for ion diffusion characteristic maximum in the permeability at the phase-transition temperature.