Lipid enrichment and selectivity of integral membrane proteins in two-component lipid bilayers

A model recently used to study lipid-protein interactions in one-component lipid bilayers (Sperotto and Mouritsen, 1991 a, b) has been extended in order to include two different lipid species characterized by different acyl-chain lengths. The model, which is a statistical mechanical lattice model, assumes that hydrophobic matching between lipid-bilayer hydrophobic thickness and hydrophobic length of the integral protein is an important aspect of the interactions. By means of Monte Carlo simulation techniques, the lateral distribution of the two lipid species near the hydrophobic protein-lipid interface in the fluid phase of the bilayer has been derived. The results indicate that there is a very structured and heterogeneous distribution of the two lipid species near the protein and that the protein-lipid interface is enriched in one of the lipid species. Out of equilibrium, the concentration profiles of the two lipid species away from the protein interface are found to develop a long-range oscillatory behavior. Such dynamic membrane heterogeneity may be of relevance for determining the physical factors involved in lipid specificity of protein function.     

Purification of squid rhodopsin and reassembly into lipid bilayers

Rhodopsin from squid photoreceptor membranes was solubilized in octyl glucoside and purified to a single band on SDS-polyacrylamide gels of $\text{M}{}_{\mathrm{<mtext>r</mtext>}}$ 46 000. Purified rhodopsin was recombined with phospholipids to form vesicles by detergent dialysis. Spectroscopic analysis of the rhodopsin-lipid vesicles showed that the interconversion between acid and basic metarhodopsin had a $\text{p}\text{K}$ of 8. Furthermore, rhodopsin in the vesicles could be photoregenerated from metarhodopsin in solutions of either neutral or alkaline pH. These two spectroscopic properties are comparable to those for rhodopsin in photoreceptor membranes. The results indicate that the native conformation of rhodopsin is preserved during purification and after recombination with phospholipids into vesicles. This preparation is, therefore, an active starting point for functional reconstitution studies.     

Computational investigation of pressure profiles in lipid bilayers with embedded proteins

The distribution of surface tension within a lipid bilayer, also referred to as the lateral pressure profile, has been the subject of theoretical scrutiny recently due to its potential to radically alter the function of biomedically important membrane proteins. Experimental measurements of the pressure profile are still hard to come by, leaving first-principles all-atom calculations of the profile as an important investigative tool. We describe and validate an efficient implementation of pressure profile calculations in the molecular dynamics package NAMD, capable of distinguishing between internal, bonded and nonbonded contributions as well as those of selected atom groups. The new implementation can also be used in conjunction with Ewald summation for long-range electrostatics, improving the accuracy and reproducibility of the calculated profiles. We then describe results of the calculation of a pressure profile for a simple protein–lipid system consisting of melittin embedded in a DMPC bilayer. While the lateral pressure in the protein–lipid system is nearly the same as that of the bilayer alone, partitioning of the lateral pressure by atom type revealed substantial perturbation of the pressure profile and surface tension in an asymmetric manner.     

Outer membrane proteins: comparing X-ray and NMR structures by MD simulations in lipid bilayers

The structures of three bacterial outer membrane proteins (OmpA, OmpX and PagP) have been determined by both X-ray diffraction and NMR. We have used multiple (7 × 15 ns) MD simulations to compare the conformational dynamics resulting from the X-ray versus the NMR structures, each protein being simulated in a lipid (DMPC) bilayer. Conformational drift was assessed via calculation of the root mean square deviation as a function of time. On this basis the ‘quality’ of the starting structure seems mainly to influence the simulation stability of the transmembrane β-barrel domain. Root mean square fluctuations were used to compare simulation mobility as a function of residue number. The resultant residue mobility profiles were qualitatively similar for the corresponding X-ray and NMR structure-based simulations. However, all three proteins were generally more mobile in the NMR-based than in the X-ray simulations. Principal components analysis was used to identify the dominant motions within each simulation. The first two eigenvectors (which account for >50% of the protein motion) reveal that such motions are concentrated in the extracellular loops and, in the case of PagP, in the N-terminal α-helix. Residue profiles of the magnitude of motions corresponding to the first two eigenvectors are similar for the corresponding X-ray and NMR simulations, but the directions of these motions correlate poorly reflecting incomplete sampling on a ∼10 ns timescale.     

Gap junction channels reconstituted in two closely apposed lipid bilayers

Intercellular communication mediated by gap junction channels plays an important role in many cellular processes. In contrast to other channels, gap junction channels span two plasma membranes resulting in an intracellular location for both ends of the junctional pore and the regulatory sites for channel gating. This configuration presents unique challenges for detailed experimental studies of junctional channel physiology and ligand-activation in situ. Availability of an appropriate model system would significantly facilitate future studies of gap junction channel function and structure. Here we show that the double-membrane channel can be reconstituted in pairs of closely apposed lipid bilayers, as experienced in cells. We have trapped the calcium-sensitive dye, arsenazo III (AIII), partially calcium-saturated (AIII-Ca), in one population of connexin32 reconstituted-liposomes, and EGTA in a second one. In such mixtures, the interaction of EGTA with AIII-Ca was measured by a large color shift from blue to red (decreased absorbance at 652 nm). The exchange of these compounds through gap junctions was proportional to these decrements. Results indicate that these connexon-mediated interliposomal channels are functional and are inhibited by the addition of alpha-glycyrrhetinic acid and by flufenamic acid, two gap junction communication inhibitors. Future use of this model system has the potential to improve our understanding of the permeability and modulation of junctional channels in its native intercellular assembly.     

Domain coupling in asymmetric lipid bilayers

Biological membranes are heterogeneous assemblies of lipids, proteins, and cholesterol that are organized as asymmetric bimolecular leaflets of lipids with embedded proteins. Modulated by the concentration of cholesterol lipids and proteins may segregate into two or more liquid phases with different physical properties that can coexist in the same membrane. In this review, we summarize recent advances on how this situation can be recreated in a supported bilayer format and how this system has been used to demonstrate the induction of ordered lipid domains in lipid compositions that are typical for the inner leaflet by lipid compositions that are typical for the outer leaflet of mammalian plasma membranes. Proteins are shown to differentially target such induced inner leaflet domains.     

Softening of lipid bilayers

The softening of wet lipid bilayer membranes during their gel-to-fluid first-order phase transition is studied by computer simulation of a family of two-dimensional microscopic interaction models. The models include a variable number, q, of lipid chain conformational states, where 2q10. Results are presented as functions of q and temperature for a number of bulk properties, such as internal energy, specific heat, and lateral compressibility. A quantitative account is given of the statistics of the lipid clusters which are found to form in the neighborhood of the transition. The occurrence of these clusters is related to the softening and the strong thermal density fluctuations which dominate the specific heat and the lateral compressibility for the high-q models. The cluster distributions and the fluctuations behave in a manner reminiscent of critical phenomena and percolation. The findings of long-lived metastable states and extremely slow relaxational behavior in the transition region are shown to be caused by the presence of intermediate lipid chain conformational states which kinetically stabilize the cluster distribution and the effective phase coexistence. This has as its macroscopic consequence that the first-order transition apperas as a continuous transition, as invariably observed in all experiments on uncharged lecithin bilayer membranes. The results also suggest an explanation of the non-horizontal isotherms of lipid monolayers. Possible implications of lipid bilayer softening and enhanced passive permeability for the functioning of biological membranes are discussed.Abbreviations PC phosphatidvlcholine - DMPC dimyristoyl PC - DPPC dipalmitoyl PC - ac alternating current - DSC differential scanning calorimetry - T _m lipid gel-to-fluid phase transition temperature - TEMPO 2,2,6,6-tetramethylpiperidine-N-oxyl Supported by the Danish Natural Science Research Council and A/S De Danske Spritfabrikkers JubilæumslegatSupported in part by the NSERC of Canada and Le FCAC du Quebec     

Kinetic analysis of bacterial reaction centers reconstituted in lipid bilayers

(1) Reaction center-lipid complexes were extracted into octane solutions. Different methods for generating an assymetric membrane distribution of reaction centers are discussed, which allow the measurement of electrical signals upon illumination. (2) The dichroism of the chromophoric groups in the reaction centers was investigated in planar lipid bilayers and the angle β between each transition moment and the normal to the membrane could be determined to be β(757 nm) = 29.5 ± 1.2, β(801 nm) = 34 ± 1.0 and β(860 nm) = 41.3 ± 0.9°. (3) The kinetics of the reaction centers from Rhodopseudomonas sphaeroides were analysed by electrical measurements and the relevant rate constants could be determined. In addition, the interaction between reaction centers and the intramembrane, ubiquinone-containing pool was investigated and described in a kinetic model. (4) The interaction between the electron-donating ferrocytochromes exhibited two distinguishable sources, a fast accessible, membrane-bound pool, which is limited by diffusion, and a pool consisting of an aqueous solution of ferrocytochrome c, which is accessible with a slower rate constant.     

Characterization of the solubilization of lipid bilayers by surfactants

This communication addresses the state of aggregation of lipid-detergent mixed dispersions. Analysis of recently published data suggest that for any given detergent-lipid mixture the most important factor in determining the type of aggregates (mixed vesicles or mixed micelles) and the size of the aggregate is the detergent to lipid molar ratio in these aggregates, herein denoted the effective ratio, R_e. For mixed bilayers this effective ratio has been previously shown to be a function of the lipid and detergent concentrations and of an equilibrium partition coefficient, K, which describes the distribution of the detergent between the bilayers and the aqueous phase. We show that, similar to mixed bilayers, the size of mixed micelles is also a function of the effective ratio, but for these dispersions the distribution of detergent between the mixed micelles and the aqueous medium obeys a much higher partition coefficient. In practical terms, the detergent concentration in the mixed micelles is equal to the difference between the total detergent concentration and the critical micelle concentration (cmc). Thus, the effective ratio is equal to this difference divided by the lipid concentration. Transformation of mixed bilayers to mixed micelles, commonly denoted solubilization, occurs when the surfactant to lipid effective ratio reaches a critical value. Experimental evaluation of this critical ratio can be based on the linear dependence of detergent concentration, required for solubilization, on the lipid concentration. According to the ‘equilibrium partition model’, the dependence of the ‘solubilizing detergent concentration’ on the lipid concentration intersects with the lipid axis at −1/K, while the slope of this dependence is the critical effective ratio. On the other hand, assuming that when solubilization occurs the detergent concentration in the aqueous phase is approximately equal to the critical micelle concentration, implies that the above dependence intersects with the detergent axis at the critical micelle concentration, while its slope, again, is equal to the critical effective ratio. Analysis of existing data suggests that within experimental error both these distinctively different approaches are valid, indicating that the critical effective ratio at which solubilization occurs is approximately equal to the product of the critical micelle concentration and the distribution coefficient K. Since the nature of detergent affects K and the critical micelle concentration in opposite directions, the critical (‘solubilizing’) effective ratio depends upon the nature of detergent less than any of these two factors.     

Toad bladder amiloride-sensitive channels reconstituted into planar lipid bilayers

In the present study we used established methods to obtain apical membrane vesicles from the toad urinary bladder and incorporated these membrane fragments to solvent-free planar lipid bilayer membranes. This resulted in the appearance of a macroscopic conductance highly sensitive to the diuretic amiloride added to the cis side. The blockage is voltage dependent and well described by a model which assumes that the drug binds to sites in the channel lumen. This binding site is localized at about 15% of the electric field across the membrane. The apparent inhibition constant (K(0)) is equal to 0.98 microM. Ca2+, in the micromolar range on the cis side, is a potent blocker of this conductance. The effect of the divalent has a complex voltage dependence and is modulated by pH. At the unitary level we have found two distinct amiloride-blockable channels with conductances of 160 pS (more frequent) and 120 pS. In the absence of the drug the mean open time is around 0.5 sec for both channels and is not dependent on voltage. The channels are cation selective (PNa/PCl = 15) and poorly discriminate between Na+ and K+ (PNa/PK = 2). Amiloride decreases the lifetime in the open state of both channels and also the conductance of the 160-pS channel.     

Reconstitution of membrane proteins. Spontaneous association of integral membrane proteins with preformed unilamellar lipid bilayers

We have developed a simple method for reconstituting pure, integral membrane proteins into phospholipid-protein vesicles. The method does not depend on use of detergents or sonication. It has been used successfully with three different types of integral membrane proteins: UDPglucuronosyltransferase (EC 2.4.1.17) from pig liver microsomes, cytochrome oxidase (EC 1.9.3.1) from pig heart, and bacteriorhodopsin from Halobacterium halobium. The method depends on preparing unilamellar vesicles of dimyristoylphosphatidylcholine (DMPC) that contain a small amount of myristate as fusogen. Under conditions that the vesicles of DMPC have the property of fusing, all of the above proteins incorporated into bilayers. Two events appear to be involved in forming the phospholipid-protein complexes. The first is a rapid insertion of all proteins into a small percentage of total vesicles. The second is slower but continued fusion of the remaining phospholipid-protein vesicles, or proteoliposomes, with small unilamellar vesicles of DMPC. This latter process was inhibited by conditions under which vesicles of DMPC themselves would not fuse. On the basis of proton pumping by bacteriorhodopsin and negative staining, the vesicles were unilamellar and large. The data suggest that insertion of the above integral membrane proteins into vesicles occurred independently of fusion between vesicles.     

A new method for membrane reconstitution: Fusion of protein-containing vesicles with planar bilayer membranes below lipid phase transition temperature

A new method of membrane reconstitution was developed by fusion of channel protein containing vesicles with planar bilayer membranes. The fusion process only occurred below and near the phase transition temperature of the lipid used. We obtained the following results: 1. Our system is solvent-free and vesicles do not come into contact with the air/water interface. This obviates a possible denaturation of hydrophobic proteins. 2. Channel forming proteins and protein complexes, respectively, are active in a frozen lipid matrix. 3. We detected an unknown channel in cilia fragments. 4. Purified acetylcholine receptors form fluctuating channels in a membrane consisting of a pure synthetic lecithin (two component system).     

Nanomechanics of lipid bilayers by force spectroscopy with AFM: A perspective

Lipid bilayers determine the architecture of cell membranes and regulate a myriad of distinct processes that are highly dependent on the lateral organization of the phospholipid molecules that compose the membrane. Indeed, the mechanochemical properties of the membrane are strongly correlated with the function of several membrane proteins, which demand a very specific, highly localized physicochemical environment to perform their function. Several mesoscopic techniques have been used in the past to investigate the mechanical properties of lipid membranes. However, they were restricted to the study of the ensemble properties of giant bilayers. Force spectroscopy with AFM has emerged as a powerful technique able to provide valuable insights into the nanomechanical properties of supported lipid membranes at the nanometer/nanonewton scale in a wide variety of systems. In particular, these measurements have allowed direct measurement of the molecular interactions arising between neighboring phospholipid molecules and between the lipid molecules and the surrounding solvent environment. The goal of this review is to illustrate how these novel experiments have provided a new vista on membrane mechanics in a confined area within the nanometer realm, where most of the specific molecular interactions take place. Here we report in detail the main discoveries achieved by force spectroscopy with AFM on supported lipid bilayers, and we also discuss on the exciting future perspectives offered by this growing research field.     

The resealing process of lipid bilayers after reversible electrical breakdown

The resealing process of lipid bilayer membranes after reversible electrical breakdown was investigated using two voltage pulses switched on together. Electrical breakdown of the membranes was induced with a voltage pulse of high intensity and short duration. The time course of the change in membrane conductance after the application of the high (short) voltage pulse was measured with a longer voltage pulse of low amplitude. The decrease in membrane conductance during the resealing process could be fitted to a single exponential curve with a time constant of 10-2 μs in the temperature range between 2 and 20°C. The activation energy for this exponential decay process was found to be about 50 kJ/mol, which might indicate a diffusion process. Above 25°C the resealing process is controlled by two exponential processes.The data obtained for the time course of the resealing process can be explained in terms of pore formation in the membranes in response to the high electrical field strength. A radius of about 4 nm is calculated for the initial pore size. From the assumed exponential change of the pore area with progressive resealing time a diffusion constant of 10^?8 cm²/s for lateral lipid diffusion can be estimated.     

Capturing the nanoscale complexity of cellular membranes in supported lipid bilayers

The lateral mobility of cell membranes plays an important role in cell signaling, governing the rate at which embedded proteins can interact with other biomolecules. The past two decades have seen a dramatic transformation in understanding of this environment, as the mechanisms and potential implications of nanoscale structure of these systems has become accessible to theoretical and experimental investigation. In particular, emerging micro- and nano-scale fabrication techniques have made possible the direct manipulation of model membranes at the scales relevant to these biological processes. This review focuses on recent advances in nanopatterning of supported lipid bilayers, capturing the impact of membrane nanostructure on molecular diffusion and providing a powerful platform for further investigation of the role of this spatial complexity on cell signaling.     

The diffusion-solubility of oxygen in lipid bilayers

The product, D_oα, of the oxygen diffusion coefficient, D_o, and the oxygen solubility, α, is determined in phosphatidylcholine bilayers at temperatures above the lipid phase transitions from ESR spin-exchange measurements. The resulting values of D_oα are in good agreement with those obtained from fluorescence-quenching experiments. The use of fatty acid spin labels makes it possible to measure D_oα as a function of the coordinate perpendicular to the bilayer surface. The results indicate that D_oα is a strong function of this coordinate; it is greatest in the bilayer center and least near the bilayer head groups.     

Pure lipid vesicles can induce channel-like conductances in planar bilayers

Summary Typical channel-like current fluctuations were observed in planar lipid bilayers following brief exposure to large concentrations of lipid vesiclesdevoid of protein. Vesicles, formed by sonication of pure lipids suspended in 150mm salt solutions, were ejected 0.5 mm from a planar bilayer with a pipette. Over the next several minutes the bilayer conductance changed in ways usually considered to be indicative of reconstituted protein channels including step conductance changes (both up and down), flickering, ion selectivity, and inactivation. This observation demonstrates the need for caution in interpreting conductance changes which occur following ejection of channel-containing vesicles near a membrane.     

Implicit membrane treatment of buried charged groups: Application to peptide translocation across lipid bilayers

The energetic cost of burying charged groups in the hydrophobic core of lipid bilayers has been controversial, with simulations giving higher estimates than certain experiments. Implicit membrane approaches are usually deemed too simplistic for this problem. Here we challenge this view. The free energy of transfer of amino acid side chains from water to the membrane center predicted by IMM1 is reasonably close to all-atom free energy calculations. The shape of the free energy profile, however, for the charged side chains needs to be modified to reflect the all-atom simulation findings (IMM1-LF). Membrane thinning is treated by combining simulations at different membrane widths with an estimate of membrane deformation free energy from elasticity theory. This approach is first tested on the voltage sensor and the isolated S4 helix of potassium channels. The voltage sensor is stably inserted in a transmembrane orientation for both the original and the modified model. The transmembrane orientation of the isolated S4 helix is unstable in the original model, but a stable local minimum in IMM1-LF, slightly higher in energy than the interfacial orientation. Peptide translocation is addressed by mapping the effective energy of the peptide as a function of vertical position and tilt angle, which allows identification of minimum energy pathways and transition states. The barriers computed for the S4 helix and other experimentally studied peptides are low enough for an observable rate. Thus, computational results and experimental studies on the membrane burial of peptide charged groups appear to be consistent. This article is part of a Special Issue entitled: Interfacially Active Peptides and Proteins. Guest Editors: William C. Wimley and Kalina Hristova.     

The modulation of thermal properties of vinblastine by cholesterol in membrane bilayers

It has been shown that the partitioning of vinblastine in 1,2-dipalmitoyl-sn-glycero-3-phosphatidylcholine (DPPC) single and multiple bilayer dispersions induces partial interdigitation of the lipid alkyl chains. Similar behavior has been observed for abietic and ursodeoxycholic acids and may well be generalized for the partitioning of bulky amphoteric molecules, which tend to localize in the vicinity of the polar heads. For the present study, differential scanning calorimetry (DSC) has been employed to investigate the role of lipid molecular characteristics such as the alkyl chain length and the polarity of the head-group, as well as the impact of cholesterol upon vinblastine-induced interdigitation. It is found that vinblastine does not induce interdigitation in lipids with either shorter or longer alkyl chains than DPPC, or having head-groups of different polarity. In addition, it is shown that the presence of cholesterol in the lipid bilayer tends to modulate the phase behavior of the lipid/vinblastine bilayer system. Preliminary studies show that such properties directly affect the encapsulation efficiency and the pharmacokinetics of liposomes.     

Interaction of fatty acids with lipid bilayers

We present a method by which it is possible to describe the binding of fatty acids to phospholipid bilayers. Binding constants for oleic acid and a number of fatty acids used as spectroscopic probes are deduced from electrophoresis measurements. There is a large shift in $\text{p}\text{K}$ value for the fatty acids on binding to the phospholipid bilayers, consistent with stronger binding of the uncharged form of the fatty acid. For dansylundecanoic acid, fluorescence titrations are consistent with the binding constants derived from the electrophoresis experiments. For 12-(9-anthroyloxy)stearic acid, fluorescence and electrophoresis data are inconsistent, and we attribute this to quenching of fluorescence at high molar ratios of 12-anthroylstearic acid to phospholipid in the bilayer.