Chain ordering in liquid crystals: II. Structure of bilayer membranes

The ordering of the hydrocarbon chain interior of bilayer membranes has been calculated using the molecular field approximation developed in previous work on liquid crystals. Different statistical averages are evaluated by exact summation over all conformations of a single chain in the field due to neighboring molecules. The internal energy of each conformation, as well as contributions arising from interaction with the molecular field and from a lateral pressure on the chain have been included.The results describe properties of both lipid monolayers and bilayers. For monolayers, the calculated pressure-area relationships are in good agreement with experimental observations. The order parameter for hydrocarbon chains in bilayers (or monolayers) as a function of temperature, lateral pressure and position along the chain, is shown and compared with the available NMR data. Combining the results of calculation and NMR measurements we obtain the value for intrinsic lateral pressure within bilayer membranes, in excellent agreement with direct measurements on surface monolayers.The calculation also gives average length of hydrocarbon chains, thermal expansion coefficient and fraction of bonds in gauche conformations. The effect of cholesterol and proteins within the bilayer is qualitatively described, and the contribution of the bilayer interior to membrane elasticity is determined.     

The water permeability of the monoolein/triolein bilayer membrane

The water permeability of the lipid bilayer can be used as a probe of membrane structure. A simple model of the bilayer, the liquid hydrocarbon model, views the membrane as a thin slice of bulk hydrocarbon liquid. A previous study (Petersen, D. (1980) Biochim. Biophys. Acta 600, 666–677) showed that this model does not accurately predict the water permeability of the monoolein/n-hexadecane bilayer: the measured activation energy for water permeation is 50% above the predicted value. From this it was inferred that the hydrocarbon chains in the lipid bilayer are more ordered than in the bulk hydrocarbon liquid. The present study tests the liquid hydrocarbon model for the monoolein/triolein bilayer, which has been shown to contain very little triolein in the plane of the membrane (Waldbillig, R.C. and Szabo, G. (1979) Biochim. Biophys. Acta 557, 295–305). Measurements of the water permeability coefficient of the bilayer are compared with predictions of the liquid hydrocarbon model based on measurements of the water permeability coefficient of bulk 8-heptadecene. The predicted and measured values agree quite closely over the temperature range studied (15–35°C): the predicted activation energy is 11.1±0.2 kcal/mol, whereas the measured activation energy for the bilayer is 9.8±0.7 kcal/mol. This close agreement is in contrast with the monoolein/n-hexadecane results and suggests that, insofar as water permeation is concerned, the liquid hydrocarbon model quite closely represents the monoolein/triolein bilayer.     

The molecular mechanism of action of the proton ionophore FCCP (carbonylcyanide p-trifluoromethoxyphenylhydrazone).

We propose a simple model that accounts for the ability of the weak acid FCCP (Carbonylcyanide-p-trifluoromethoxyphenylhydrazone) to both transport protons across phospholipid bilayer membranes and uncouple oxidation from phosphorylation in mitochondria. Four parameters are required to characterize this model: the rate constant for the movement of A- across the membrane, kA, the rate constant for the movement of HA across the membrane, kHA, the adsorption coefficient of A- onto the membrane-solution interface, beta A, and the surface pK. These four parameters were determined from kinetic measurements on planar bilayer membranes using the charge-pulse and voltage-clamp techniques. We confirmed the adequacy of the model by determining each of these parameters independently, utilizing equilibrium dialysis, zeta potential, membrane potential, spectrophotometric, and conductance measurements. For a phosphatidylethanolamine bilayer the values of the parameters are kHA = 10(4)S-1, beta A = 3 10(-3) cm, and 6.0 less than pK less than 6.4. As predicted theoretically, the value of KA depends on both the applied voltage, V, and dielectric constant of the membrane, epsilon r; when V approaches zero and the membrane contains chlorodecane (epsilon r congruent to 2.7) kA = 700 s-1. If oxidation is coupled to phosphorylation by means of a delta microH+, and V er congruent to 2.7 for the inner membrane of the mitochondrion, the model predicts that FCCP should exert maximal uncoupling activity at a pH congruent to pK. This prediction agrees with the published experimental results.     

Compression of lipid membranes as observed at varying membrane positions.

We have measured the microscopic isothermal compressibility of dioleoyl- and dimyristyl-phosphatidylcholine multilayers and bilayers as a function of membrane depth by the pressure dependence of the polarization of a series of anthroyloxy fatty acids. In both systems, within experimental error, the compressibility did not change with membrane depth. The magnitudes of the compressibilities matched those of organic solids and those reported for dipalmitoylphosphatidylcholine multilayers from neutron diffraction measurements (Braganza, L. F., and D. L. Worcester. 1986. Biochemistry, 25:7484-7488). The bilayer compressibility decreased with temperature and this decrease was similar with membrane depth consistent with the isotropic thermal expansion of membranes previously observed (Scarlata, S. 1989. Biophys. J. 55:1215-1223). The vertical compressibility in the z direction is much lower than the horizontal (xy planes) for probes that lie parallel to the hydrocarbon chains which is consistent with an increase in bilayer thickness. The compressibility for probes that lie perpendicular to the hydrocarbon chains is more isotropic due to their limited spatial access to the z plane.     

Evaluation of the thickness of lipid bilayers of biological membranes by measurement of specific capacity

The estimate of the thickness of the hydrophobic part of biological membranes by means of its formula for plane capacitor leads to a necessity of very high values of protein dielectric permeability. The current estimate can be performed if the membrane possesses a wavy cross-section. Due to the existence of integral proteins which occupy about 1/3 of the bilayer surface and have dielectric permeability epsilon = 6 divided by 8 and also due to some peculiarities of intermembrane electric field the thickness of hydrophobic bilayer region can be estimated from membrane specific capacity as 3 divided by 3.5 nm, this value coincides with X-ray data.     

Evidence of partial rotational order in gel phase DPPC.

It is shown how the dichroic ratio of the symmetric methylene stretching modes depends upon both the rotational order of hydrocarbon chains about their long axis and the tilting of the long chains with respect to the bilayer normal. Use of a recent determination of the tilt angle from x-ray measurements together with recent dichroic infrared data yields a rotational order parameter g = -0.30 compared to g = 0 for complete disorder and g = +/- 1 for complete order. The negative value of g corresponds to a preference for the plane defined by the chain carbons to be more perpendicular than parallel to the plane defined by the tilt direction and the bilayer normal.     

The noneffect of a large linear hydrocarbon, squalene, on the phosphatidylcholine packing structure.

The interaction of squalene with liposomes and monolayers of dipalmitoyl phosphatidylcholine (DPL) has been studied by differential scanning calorimetry, Raman spectroscopy, and surface potential measurements. Mole ratios of squalene to DPL up to 9 to 1 were studied. In contrast to small, nonpolar molecules, which profoundly influence the structure of lipid bilayers as detected by changes in both their thermodynamic phase transition parameters and membrane fluidity, this large, nonpolar, linear hydrocarbon is devoid of such influences. It is clear from our data that a large nonpolar molecule such as squalene, having no polar group that might anchor it to the aqueous interface, cannot intercalate between the acyl chains either below or above the phase transition of DPL. This behavior is not compatible with models that treat the bilayer interior as a bulk hydrocarbon, and suggests that great caution should be exercised in extrapolating partition coefficients based on bulk hydrocarbon measurements to lipid bilayers.     

Cholesterol dynamics in membranes.

Time-resolved fluorescence anisotropy of the sterol analogue, cholestatrienol, and 13C nuclear magnetic resonance (NMR) spin lattice relaxation time (T1c) measurements of [13C4] labeled cholesterol were exploited to determine the correlation times characterizing the major modes of motion of cholesterol in unsonicated phospholipid multilamellar liposomes. Two modes of motion were found to be important: (a) rotational diffusion and (b) time dependence of the orientation of the director for axial diffusion, or "wobble." From the time-resolved fluorescence anisotropy decays of cholestatrienol in egg phosphatidylcholine (PC) bilayers, a value for tau perpendicular, the correlation time for wobble, of 0.9 x 10(-9) s and a value for S perpendicular, the order parameter characterizing the same motion, of 0.45 s were calculated. Both tau perpendicular and S perpendicular were relatively insensitive to temperature and cholesterol content of the membranes. The T1c measurements of [13C4] labeled cholesterol did not provide a quantitative determination of tau parallel, the correlation time for axial diffusion. T1c from the lipid hydrocarbon chains suggested a value for tau perpendicular similar to that for cholesterol. Steady-state anisotropy measurements and time-resolved anisotropy measurements of cholestatrienol were used to probe sterol behavior in a variety of pure and mixed lipid multilamellar liposomes. Both the lipid headgroups and the lipid hydrocarbons chains contributed to the determination of the sterol environment in the membrane, as revealed by these fluorescence measurements. In particular, effects of the phosphatidylethanolamine (PE) headgroup and of multiple unsaturation in the lipid hydrocarbon chains were observed. However, while the steady-state anisotropy was sensitive to these factors, the time-resolved fluorescence analysis indicated that tau perpendicular was not strongly affected by the lipid composition of the membrane. S perpendicular may be increased by the presence of PE. Both steady-state anisotropy measurements and time-resolved anisotropy measurements of cholestatrienol were used to probe sterol behavior in three biological membranes: bovine rod outer segment (ROS) disk membranes, human erythrocyte plasma membranes, and light rabbit muscle sarcoplasmic reticulum membranes. In the ROS disk membranes the value for S perpendicular was marginally higher than in the PC membranes, perhaps reflecting the influence of PE. The dramatic difference noted was in the value for tau perpendicular. In both the ROS disk membranes and the erythrocyte membranes, tau perpendicular was one-third to one-fifth of tau perpendicular in the phospholipid bilayers. This result may reveal an influence of membrane proteins on sterol behavior.     

Dehydration of model membranes induced by lectins from Ricinus communis and Viscum album.

The effects of ribosome-inactivating proteins (RIPs) from Ricinus communis and from Viscum album on the water permeability, Pf, and the surface dielectric constant, epsilon, of model membranes were studied. Pf was calculated from microelectrode measurements of the ion concentration distribution in the immediate vicinity of a planar membrane, and epsilon was obtained from the fluorescence of dansyl phosphatidylethanolamine incorporated into unilamellar vesicles. Pf and epsilon of fully saturated phosphatidylcholine membranes were affected only in the presence of a lectin receptor (monosialoganglioside, GM1) in the bilayer. It is suggested that the membrane area occupied by clustered lectin-receptor complexes is markedly less permeable to water. Protein binding to the receptor was not a prelude for hydrophobic lipid-protein interactions when the membranes were formed from a mixture of natural phospholipids with a high content of unsaturated fatty acids. These membranes, characterized by a high initial water permeability, were found to interact with the RIPs unspecifically. From a decrease of both Pf and epsilon it was concluded that not only water partitioning but also protein adsorption correlates with looser packing of polyunsaturated lipids at the lipid-water interface.     

Interaction among molecules in mixtures of ceramide/stearic acid,ceramide/cholesterol and ceramide/stearic acid/cholesterol

To elucidate the interaction among the molecules which constitute intercellular lipids of stratum corneum, the phase diagrams in the binary mixtures of N-octadecanoyl-phytosphingosine (CER)/stearic acid (SA) and CER/cholesterol (CHOL) were studied by differential scanning calorimetry and by small- and wide-angle X-ray diffraction. These phase diagrams are mostly expressed by a eutectic type one. However, from their detailed analyses, it was revealed that in the phase diagram of CER/SA a new solid structure is formed just above the eutectic temperature. The lamellar spacing of the new structure is nearly equal to the length given by the sum of the two molecules of CER and/or SA, that is, in the lipid bilayer the hydrocarbon chains of CER and SA lie almost perpendicular to the lipid bilayer surface and the two kinds of molecules distribute homogeneously. On the other hand, in the binary mixture of CER/CHOL, CHOL molecules are apt to be isolated from the mixture. In a ternary mixture composed of equimolar lipids of CER, CHOL and SA, it was found that a pseudo-hexagonal structure takes place even in the solid state. This fact indicates that the three components are miscible and the hydrocarbon chains lie perpendicular to the lipid bilayer surface. We can draw the conclusion that the multi-component mixtures containing ceramide are apt to form the lamellar structure where even in the solid state the hydrocarbon chains lie perpendicular to the lipid bilayer surface and the components with hydrocarbon chains distribute homogeneously.     

Electrostatic potential and Born energy of charged molecules interacting with phospholipid membranes: calculation via 3-D numerical solution of the full Poisson equation.

Understanding the physicochemical basis of the interaction of molecules with lipid bilayers is fundamental to membrane biology. In this study, a new, three-dimensional numerical solution of the full Poisson equation including local dielectric variation is developed using finite difference techniques in order to model electrostatic interactions of charged molecules with a non-uniform dielectric. This solution is used to describe the electric field and electrostatic potential profile of a charged molecule interacting with a phospholipid bilayer in a manner consistent with the known composition and structure of the membrane. Furthermore, the Born interaction energy is then calculated by appropriate integration of the electric field over whole space. Numerical computations indicate that the electrostatic potential profile surrounding a charge molecule and its resultant Born interaction energy are a function of molecular position within the membrane and change most significantly within the polar region of the bilayer. The maximum interaction energy is observed when the charge is placed at the center of the hydrophobic core of the membrane and is strongly dependent on the size of the charge and on the thickness of the hydrocarbon core of the bilayer. The numerical results of this continuum model are compared with various analytical approximations for the Born energy including models established for discontinuous slab dielectrics. The calculated energies agree with the well-known Born analytical expression only when the charge is located near the center of a hydrocarbon core of greater than 60 A in thickness. The Born-image model shows excellent agreement with the numerical results only when modified to include an appropriate effective thickness of the low dielectric region. In addition, a newly derived approximation which considers the local mean dielectric provides a simple and continuous solution that also agrees well with the numerical results.     

The effect of the lipid bilayer state on fluorescence intensity of fluorescein-PE in a saturated lipid bilayer

Fluorescein-PE is a fluorescence probe that is used as a membrane label or a sensor of surface associated processes. Fluorescein-PE fluorescence intensity depends not only on bulk pH, but also on the local electrostatic potential, which affects the local membrane interface proton concentration. The pH sensitivity and hydrophilic character of the fluorescein moiety was used to detect conformational changes at the lipid bilayer surface. When located in the dipalmitoylphosphatidylcholine (DPPC) bilayer, probe fluorescence depends on conformational changes that occur during phase transitions. Relative fluorescence intensity changes more at pretransition than at the main phase transition temperature, indicating that interface conformation affects the condition in the vicinity of the membrane. Local electrostatic potential depends on surface charge density, the local dielectric constant, salt concentration and water organisation. Initial increase in fluorescence intensity at temperatures preceding that of pretransition can be explained by the decreased value of the dielectric constant in the lipid polar headgroups region related in turn to decreased water organisation within the membrane interface. The abrupt decrease in fluorescence intensity at temperatures between 25 degrees C and 35 degrees C (DPPC pretransition) is likely to be caused by an increased value of the electrostatic potential, induced by an elevated value of the dielectric constant within the phosphate group region. Further increase in the fluorescence intensity at temperatures above that of the gel-liquid phase transition correlates with the calculated decreased surface electrostatic potential. Above the main phase transition temperature, fluorescence intensity increase at a salt concentration of 140 mM is larger than with 14 mM. This results from a sharp decline of the electrostatic potential induced by the phosphocholine dipole as a function of distance from the membrane surface.     

Synthesis of 6-thio pseudo glycolipids and their orientation on a gold slide studied by IRRAS

We have synthesized four 6-thio pseudo glycolipid analogues and assessed how two of them self-assembled on a gold surface. These structures were designed as candidate tethers molecules to anchor bilayer lipid membranes on gold. 6-Deoxy-6-thiogalactose was chosen to anchor the macromolecule to the gold and define an aqueous zone at the gold surface. A long alkane chain (C-12 or C-18) linked to the anomeric position of the sugar residue was chosen to anchor a bilayer lipid membrane. The linkage between the carbohydrate and the hydrophobic chains is either a glycosidic bond or a 1,4-disubstituted triazole formed by copper(I)-catalysed alkyne-azide cycloaddition (CuAAC) of the propargyl glycoside with azido-dodecane and azido-octadecane. We are expecting that the hydrocarbon chains will orient themselves perpendicular to the gold surface and be incorporated into the first leaflet of the bilayer membrane. We have studied self assembled monolayers of the C-12 aglycone analogues on gold using infrared reflection absorption spectroscopy (IRRAS). We compared the results given by the IRRAS experiments to the IR spectra recorded by attenuated total reflection (ATR) spectroscopy on films of the randomly oriented analogues. Our results demonstrate that the C-12 analogues did bind to gold and did orient themselves perpendicular to the gold slide.     

Interaction of apoprotein from porcine high-density lipoprotein with dimyristoyl lecithin. 1. The structure of the complexes.

The morphology and structural organisation of the complexes formed from the apoprotein of porcine high-density lipoprotein and dimyristoyl phosphatidylcholine (lecithin) have been studied using the technique of small-angle X-ray scattering. Scattering measurements made in solvents of varying electron density were interpreted in terms of a scattering-equivalent model for the structure of the complex. This model is described by an oblate ellipsoidal morphology with dimensions at 20 degrees C: major axis 11.0 nm, minor axis 5.5 nm. Within this overall shape the lipid hydrocarbon chains are organised in an apolar core whilst the lipid polar head groups and protein are located in a outer shell 0.85 nm in thickness. The oblate morphology demonstrates that the structure of the complex is directed by the fundamental bilayer organisation of the lecithin. The dimension of the minor axis (5.5 nm) indicates that phospholipid hydrocarbon chains are orientated perpendicular to the interface.     

The tilting of the hydrocarbon chains in a single bilayer of phospholipid.

The tilt angle of the hydrocarbon chains to the planes of a dipalmitoryl lecithin single bilayer and multilayers were estimated by the asymmetry of the electron diffraction patterns of respective hydrated specimens. The chains in a single bilayer were found to be perpendicular to the bilayer plane, whereas the chains in the multilayers were found to be tilted with respect to the normal of the plane. Thermal analysis data also supported this conclusion.     

Water permeation through the lipid bilayer membrane. Test of the liquid hydrocarbon model

According to the liquid hydrocarbon model, the lipid bilayer is viewed simply as a thin slice of bulk hydrocarbon liquid. This allows the water permeability of the bilayer to be calculated from bulk properties. In this paper the prediction of the liquid hydrocarbon model is compared with the known water permeability coefficient of the glycerol monoolein/n-hexadecane bilayer (Fettiplace, R. (1978) Biochim. Biophys. Acta 513, 1–10). As the alkyl chain of glycerol monoolein is equivalent to 8-heptadecene, the water permeability coefficient of 8-heptadecene/n-hexadecane mixtures was measured for temperatures between 20 and 35°C. The mole fraction of n-hexadecane in the bulk liquid was chosen at each temperature to match the known mole fraction of n-hexadecane in the bilayer (White, S. (1976) Nature 262, 421–422). The predicted water permeability coefficient agrees with the measured value at 32°C but is 40% above the measured value at 20°C. The apparent activation energy predicted by the liquid hydrocarbon model is 9.0 ± 0.3 kcal/mol, while the measured value is 14.2 ± 1.0 kcal/mol. The failure of the liquid hydrocarbon model probably results from a different molecular organization of the hydrocarbon chains in the bilayer and in the bulk liquid.     

Dielectric control of counterion-induced single-chain folding transition of DNA

In the presence of condensing agents, single chains of giant double-stranded DNA undergo a first-order phase transition between an elongated coil state and a folded compact state. To connect this like-charged attraction phenomenon to counterion condensation, we performed a series of single-chain experiments on aqueous solutions of DNA, where we varied the extent of counterion condensation by varying the relative dielectric constant epsilon(r) from 80 to 170. Single-chain observations of changes in the conformation of giant DNA were performed by transmission electron microscopy and fluorescence microscopy, with tetravalent spermine (SPM(4+)) as a condensing agent. At a fixed dielectric constant, single DNA chains fold into a compact state upon the addition of spermine, whereas at a constant spermine concentration single DNA chains unfold with an increase in epsilon(r). In both cases, the transition is largely discrete at the level of single chains. We found that the critical concentration of spermine necessary to induce the single-chain folding transition increases exponentially as the dielectric constant increases, corresponding to 87-88% of the DNA charge neutralized at the onset of the transition. We also observed that the toroidal morphology of compact DNA partially unfolds when epsilon(r) is increased.     

Correlation of electron transfer rate in photosynthetic reaction centers with intraprotein dielectric properties

A number of the electrogenic reactions in photosystem I, photosystem II, and bacterial reaction centers (RC) were comparatively analyzed, and the variation of the dielectric permittivity (epsilon) in the vicinity of electron carriers along the membrane normal was calculated. The value of epsilon was minimal at the core of the complexes and gradually increased towards the periphery. We found that the rate of electron transfer (ET) correlated with the value of the dielectric permittivity: the fastest primary ET reactions occur in the low-polarity core of the complexes within the picosecond time range, whereas slower secondary reactions take place at the high-polarity periphery of the complexes within micro- to millisecond time range. The observed correlation was quantitatively interpreted in the framework of the Marcus theory. We calculated the reorganization energy of ET carriers using their van der Waals volumes and experimentally determined epsilon values. The electronic coupling was calculated by the empirical Moser-Dutton rule for the distance-dependent electron tunneling rate in nonadiabatic ET reactions. We concluded that the local dielectric permittivity inferred from the electrometric measurements could be quantitatively used to estimate the rate constant of ET reactions in membrane proteins with resolved atomic structure with the accuracy of less than one order of magnitude.     

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.