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).     

Application of SVM to predict membrane protein types

As a continuous effort to develop automated methods for predicting membrane protein types that was initiated by Chou and Elrod (PROTEINS: Structure, Function, and Genetics, 1999, 34, 137-153), the support vector machine (SVM) is introduced. Results obtained through re-substitution, jackknife, and independent data set tests, respectively, have indicated that the SVM approach is quite a promising one, suggesting that the covariant discriminant algorithm (Chou and Elrod, Protein Eng. 12 (1999) 107) and SVM, if effectively complemented with each other, will become a powerful tool for predicting membrane protein types and the other protein attributes as well.     

Effects of membrane composition and lipid structure on the photopolymerization of lipid diacetylenes in bilayer membranes

Molecules analogous to biological and synthetic lipids have been prepared with conjugated diacetylene moieties in the long alkyl chain. These lipid diacetylenes form bilayer structures when suspended in aqueous buffers. Ultraviolet light (254 nm) exposure initiates the polymerization of the diacetylenes in the lipid bilayer to give a fully conjugated, highly colored product. The reaction is topotactic, and its efficiency depends on the correct alignment of the monomeric units. Thus, the lipid diacetylenes are photopolymerizable if the hydrocarbon chains are in a regular lattice found at temperatures below the lipid transition temperature; polymerization is inhibited above this transition. The photopolymerization of a diacetylenic glycerophosphocholine in lipid bilayer membranes was observed in two-component mixtures with a nonpolymerizable lipid, either dioleoylphosphatidylcholine or distearoylphosphatidylcholine. The photochemical and thermochemical characteristics suggest that the diacetylenic glycerophosphocholine exists largely in separate domains in the mixed bilayers. Lipid diacetylenes analogous to a dialkyldimethylammonium salt and to a dialkyl phosphate have a plane of symmetry, which suggests that both chains penetrate equally into the bilayer. The photopolymerization of these symmetrical synthetic species is more than 10³-times more efficient than that of the diacetylenic glycerophosphocholine. These differences are interpretable in terms of the expected conformational preference of the lipid molecules.     

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.     

On the transport noise of hydrophobic ions in lipid bilayer membranes

The effect of transit time on the electrical transport noise of a closed one-barrier model at equilibrium as proposed by Kolb and Läuger [6] is studied using the master-equation approach. A transit time is the time for an ion to cross the energy barrier (membrane interior) when the energy of the ion reaches the barrier height. Both the time correlation function and the noise power spectrum are obtained as functions of the transit time of the ions. Possible effects of transit time on the time correlation function of transport of dipicrylamine ions in lipid bilayers as reported by Bruner and Hall [13] and on the noise power spectrum as reported by Kolb and Läuger [6] are discussed.     

A synthetic peptide forms voltage-gated porin-like ion channels in lipid bilayer membranes

Design of simple protein structures represents the essential first step toward novel macromolecules and understanding the basic principles of protein folding. Our work focuses on the ion channel formation and structure of peptides having a repeated pattern of glycine residues. Investigation of the ion channel properties of a glycine repeat peptide, VSLGLSIGFSVGVSIGWSFGRSRG revealed the formation of porin-like high conductance, multimeric, non-selective voltage-gated channels in phospholipid bilayer membranes. ATR-IR and CD spectroscopic studies showed an anti-parallel beta sheet structure in membranes. The formation of porin-like ion channels by a beta sheet peptide suggests spontaneous assembly into a beta barrel structure through oligomerization as in pore forming bacterial toxins. The present work is the first example of a short synthetic peptide mimicking the pore characteristics of a complex beta barrel protein and demonstrates that smaller peptides are capable of mimicking the complex functional properties of natural ion channels. This will have implications in understanding the folding of beta sheet proteins in membranes, the mechanism of two state voltage gating, and the role of glycine residues in beta barrel proteins.     

Pore forming activity of the major outer membrane protein of Rhodobacter capsulatus in lipid bilayer membranes

Porin of the outer membrane of Rhodobacter capsulatus St. Louis (ATCC 23782) was isolated and reconstituted into lipid bilayer membranes. The porin was obtained either by the sodium dodecyl sulfate treatment of cell envelopes (SDS-porin) or by saline extraction of whole cells (NaCl-porin). Nanomolar concentrations of both porin preparations resulted in a strong conductance increase of the lipid bilayer membranes by many orders of magnitude. At small protein concentrations the conductance increased in a stepwise fashion, the average single channel conductance being about 0.35 nS in 0.1 M KCl for SDS-porin and NaCl-porin as well. The single channel conductance was a linear function of the specific conductance of the aqueous phase. The results were consistent with the assumption that the porin formed large water-filled transmembrane channels in the membrane. From the average value of the single channel conductance in 0.1 M KCl an effective channel diameter of about 1.5 nm was estimated for both types of porins.Abbreviations EDTA ethylenediamine tetraacetic acid - SDS sodium dodecyl sulfate     

A simple theory of peptide interactions on a membrane surface: excluded volume and entropic order

A simple theory of the interactions of peptides bound onto a lipid membrane is developed, modeling the peptides as rods on a surface. At low peptide surface-concentration, excluded volume dominates the peptide-peptide interactions and the orientation of the peptides is random, resulting in an isotropic configuration. However, at high peptide density on the membrane, the peptides become orientationally ordered, resulting in an anisotropic configuration. This effect is entirely entropic in origin, and simply reflects the fact that peptides can be exchanged more easily on the surface if they are equally aligned, resulting in a larger number of possible configurations. In three dimensions, this phenomenon corresponds to the well-known isotropic-nematic phase transition. In two dimensions, the problem is not as well understood. The theoretical treatment presented here yields a simple, manageable expression which can be compared with experimental data. Two-dimensional ordering results in an increase in the apparent binding constant of peptides to membranes at high concentration of peptides relative to what is expected from the effect of excluded volume alone. The possible implications of side-by-side alignment for several biological processes, such as peptide translocation across membranes and plaque formation in Alzheimer's disease, are discussed.     

E(z), a depth-dependent potential for assessing the energies of insertion of amino acid side-chains into membranes: derivation and applications to determining the orientation of transmembrane and interfacial helices

We have developed an empirical residue-based potential (E(z) potential) for protein insertion in lipid membranes. Propensities for occurrence as a function of depth in the bilayer were calculated for the individual amino acid types from their distribution in known structures of helical membrane proteins. The propensities were then fit to continuous curves and converted to a potential using a reverse-Boltzman relationship. The E(z) potential demonstrated a good correlation with experimental data such as amino acid transfer free energy scales (water to membrane center and water to interface), and it incorporates transmembrane helices of varying composition in the membrane with trends similar to those obtained with translocon-mediated insertion experiments. The potential has a variety of applications in the analysis of natural membrane proteins as well as in the design of new ones. It can help in calculating the propensity of single helices to insert in the bilayer and estimate their tilt angle with respect to the bilayer normal. It can be utilized to discriminate amphiphilic helices that assume a parallel orientation at the membrane interface, such as those of membrane-active peptides. In membrane protein design applications, the potential allows an environment-dependent selection of amino acid identities.     

Galectin-3 interacts with membrane lipids and penetrates the lipid bilayer

The precise mechanism by which galectin-3 and other cytosolic proteins that lack signal peptides are secreted is yet to be elucidated. In the present analyses, we determined that galectin-3, a beta-galactoside binding protein, can interact directly with membrane lipids in solid phase binding assays. More interestingly, we determined by spectrophotometric methods that it can spontaneously penetrate the lipid bilayer of liposomes in either direction. These findings suggest that galectin-3 on its own has the capacity to traverse the lipid bilayer. Whereas the situation is rather simplified in liposomes, the interaction of galectin-3 with the plasma membrane may involve cholesterol-rich membrane domains where galectin-3 can be concentrated and form multimers or interact covalently with other proteins.     

Assembly of transmembrane helices of simple polytopic membrane proteins from sequence conservation patterns

The transmembrane (TM) domains of most membrane proteins consist of helix bundles. The seemingly simple task of TM helix bundle assembly has turned out to be extremely difficult. This is true even for simple TM helix bundle proteins, i.e., those that have the simple form of compact TM helix bundles. Herein, we present a computational method that is capable of generating native-like structural models for simple TM helix bundle proteins having modest numbers of TM helices based on sequence conservation patterns. Thus, the only requirement for our method is the presence of more than 30 homologous sequences for an accurate extraction of sequence conservation patterns. The prediction method first computes a number of representative well-packed conformations for each pair of contacting TM helices, and then a library of tertiary folds is generated by overlaying overlapping TM helices of the representative conformations. This library is scored using sequence conservation patterns, and a subsequent clustering analysis yields five final models. Assuming that neighboring TM helices in the sequence contact each other (but not that TM helices A and G contact each other), the method produced structural models of Calpha atom root-mean-square deviation (CA RMSD) of 3-5 A from corresponding crystal structures for bacteriorhodopsin, halorhodopsin, sensory rhodopsin II, and rhodopsin. In blind predictions, this type of contact knowledge is not available. Mimicking this, predictions were made for the rotor of the V-type Na(+)-adenosine triphosphatase without such knowledge. The CA RMSD between the best model and its crystal structure is only 3.4 A, and its contact accuracy reaches 55%. Furthermore, the model correctly identifies the binding pocket for sodium ion. These results demonstrate that the method can be readily applied to ab initio structure prediction of simple TM helix bundle proteins having modest numbers of TM helices.     

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.     

Differentiation between transmembrane helices and peripheral helices by the deconvolution of circular dichroism spectra of membrane proteins.

The interpretation of the circular dichroism (CD) spectra of proteins to date requires additional secondary structural information of the proteins to be analyzed, such as X-ray or NMR data. Therefore, these methods are inappropriate for a CD database whose secondary structures are unknown, as in the case of the membrane proteins. The convex constraint analysis algorithm (Perczel, A., Hollósi, M., Tusnády, G., & Fasman, G. D., 1991, Protein Eng. 4, 669-679), on the other hand, operates only on a collection of spectral data to extract the common spectral components with their spectral weights. The linear combinations of these derived "pure" CD curves can reconstruct the original data set with great accuracy. For a membrane protein data set, the five-component spectra so obtained from the deconvolution consisted of two different types of alpha helices (the alpha helix in the soluble domain and the alpha T helix, for the transmembrane alpha helix), a beta-pleated sheet, a class C-like spectrum related to beta turns, and a spectrum correlated with the unordered conformation. The deconvoluted CD spectrum for the alpha T helix was characterized by a positive red-shifted band in the range 195-200 nm (+95,000 deg cm2 dmol-1), with the intensity of the negative band at 208 nm being slightly less negative than that of the 222-nm band (-50,000 and -60,000 deg cm2 dmol-1, respectively) in comparison with the regular alpha helix, with a positive band at 190 nm and two negative bands at 208 and 222 nm with magnitudes of +70,000, -30,000, and -30,000 deg cm2 dmol-1, respectively.     

A knowledge-based scale for amino acid membrane propensity

In this article, a membrane-propensity scale for amino acids is derived using only two ingredients: (i) a set of transmembrane helices segments from membrane protein crystal structures and (ii) the request that each component of the set has a free energy lower than that of a typical soluble protein sequence of the same length. Although the most widely used hydropathy scales satisfy this request, we use an optimization procedure that allows for extraction of an optimal scale, which correlates equally well with those scales. We show that, if the choice of the sequence database is accurate, significant knowledge-based scales, which are robust with respect to changes in the learning set, can be easily derived. The obtained scales can be used for transmembrane helices prediction. The predictive power of one of these scales is tested on membrane proteins, soluble proteins, and signal peptides databases, finding that its performances is comparable with those of the hydropathy scales.     

A highly unusual palindromic transmembrane helical hairpin formed by SARS coronavirus E protein

The agent responsible for the recent severe acute respiratory syndrome (SARS) outbreak is a previously unidentified coronavirus. While there is a wealth of epidemiological studies, little if any molecular characterization of SARS coronavirus (SCoV) proteins has been carried out. Here we describe the molecular characterization of SCoV E protein, a critical component of the virus responsible for virion envelope morphogenesis. We conclusively show that SCoV E protein contains an unusually short, palindromic transmembrane helical hairpin around a previously unidentified pseudo-center of symmetry, a structural feature which seems to be unique to SCoV. The hairpin deforms lipid bilayers by way of increasing their curvature, providing for the first time a molecular explanation of E protein's pivotal role in viral budding. The molecular understanding of this critical component of SCoV may represent the beginning of a concerted effort aimed at inhibiting its function, and consequently, viral infectivity.     

Fission of the bilayer lipid tube

The fusion of two black lipid membranes results in the formation of peculiar bilayer lipid tubes (‘cylindrical’) membranes (Neher, E. (1974) Biochim. Biophys. Acta 373, 328–336 and Melikyan, G.B., Abidor, L.G., Chernomordik, L.V. and Chailakhyan, L.M. (1983) Biochim. Biophys. Acta 730, 395–398). The mechanical stability of such tubes has been investigated experimentally and theoretically. With increasing hydrostatic pressure on the outside of the tube the radius of its middle part decreases. After this radius has reached a critical value, which constitutes 0.55 of the radius of the tube base, there occurs a collapse of the tube and its disintegration into two planar bilayers (fission). Expressions are obtained which relate the transmembrane difference of the hydrostatic pressure, causing the collapse, to the geometrical characteristics of the tube (its length and the radius of its base) and to the tension of the lipid bilayer. A method for measuring the membrane tension is proposed on the basis of the phenomenon considered.     

Organization and dynamics of NBD-labeled lipids in lipid bilayer analyzed by FRET using the small membrane fluorescent probe AHBA as donor

Fluorescent probes are employed to investigate natural and model membranes. It is important to know probe location and extent of perturbations they cause into the lipid bilayer. Förster Resonance Energy Transfer (FRET) is a useful tool to investigate phenomena involving plasma membranes, and reports in literature used relatively large fluorophores like 1,6-diphenylhexatriene, located at the center of the hydrophobic region, 4-aminophthalimide-based molecules located at lipid/water interfaces and BODIPY-labeled phosphatidylcholine. In this work we explored FRET process in 1,2-dimyristoyl-L-α-GPC large unilamellar vesicles, in gel and fluid phase, using as donor the very small group o-Abz bound to hexadecyl chain (2-amino-N-hexadecyl-benzamide - AHBA) and 7-nitro-2-1,3-benzoxadiazol-4-yl (NBD) labeled lipids as acceptor. From the intensity decay of donor in presence of acceptors, the FRET efficiency was calculated, and used to fit the model proposed by Fung and Stryer to that efficiency. Using lipid bilayer structural data, the procedure allowed the determination of Förster distance for each donor-acceptor pair in vesicles, without imposing any value for the orientational factor κ². From distance distributions between o-Abz in AHBA and NBD in lipid bilayer obtained using the program CONTIN, we obtained donor-acceptor populations having different separation distances. The populations reflect the occurrence of FRET involving probes in the same or in opposite leaflet. A dynamic picture emerged showing how relative position of the probes is dependent on the structural thermal phase of the DMPC bilayer. The results emphasize the need of careful analysis in order to understand processes involving fluorescent probes in model membranes.     

Regulation of channel function due to physical energetic coupling with a lipid bilayer

Regulation of membrane protein functions due to hydrophobic coupling with a lipid bilayer has been investigated. An energy formula describing interactions between lipid bilayer and integral ion channels with different structures, which is based on the screened Coulomb interaction approximation, has been developed. Here the interaction energy is represented as being due to charge-based interactions between channel and lipid bilayer. The hydrophobic bilayer thickness channel length mismatch is found to induce channel destabilization exponentially while negative lipid curvature linearly. Experimental parameters related to channel dynamics are consistent with theoretical predictions. To measure comparable energy parameters directly in the system and to elucidate the mechanism at an atomistic level we performed molecular dynamics (MD) simulations of the ion channel forming peptide–lipid complexes. MD simulations indicate that peptides and lipids experience electrostatic and van der Waals interactions for short period of time when found within each other’s proximity. The energies from these two interactions are found to be similar to the energies derived theoretically using the screened Coulomb and the van der Waals interactions between peptides (in ion channel) and lipids (in lipid bilayer) due to mainly their charge properties. The results of in silico MD studies taken together with experimental observable parameters and theoretical energetic predictions suggest that the peptides induce ion channels inside lipid membranes due to peptide–lipid physical interactions. This study provides a new insight helping better understand of the underlying mechanisms of membrane protein functions in cell membrane leading to important biological implications.     

Oxygen permeability of the lipid bilayer membrane made of calf lens lipids

The oxygen permeability coefficient across the membrane made of the total lipid extract from the plasma membrane of calf lens was estimated from the profile of the oxygen transport parameter (local oxygen diffusion-concentration product) and compared with those estimated for membranes made of an equimolar 1-palmitoyl-2-oleoylphosphatidylcholine/cholesterol (POPC/Chol) mixture and of pure POPC. Profiles of the oxygen transport parameter were obtained by observing the collision of molecular oxygen with nitroxide radical spin labels placed at different depths in the membrane using the saturation-recovery EPR technique and were published by us earlier (J. Widomska, M. Raguz, J. Dillon, E. R. Gaillard, W. K. Subczynski, Biochim. Biophys. Acta. 1768 (2007) 1454-1465). At 35 °C, the estimated oxygen permeability coefficients were 51.3, 49.7, and 157.4 cm/s for lens lipid, POPC/Chol, and POPC membranes, respectively (compared with 53.3 cm/s for a water layer with the same thickness as a membrane). Membrane permeability significantly decreases at lower temperatures. In the lens lipid membrane, resistance to the oxygen transport is located in and near the polar headgroup region of the membrane to the depth of the ninth carbon, which is approximately where the steroid-ring structure of cholesterol reaches into the membrane. In the central region of the membrane, oxygen transport is enhanced, significantly exceeding that in bulk water. It is concluded that the high level of cholesterol in lens lipids is responsible for these unique membrane properties.     

Properties of the large ion-permeable pores formed from protein F of Pseudomonas aeruginosa in lipid bilayer membranes

The incorporation of porin protein F from the outer membrane of Pseudomonas aeruginosa into artificial lipid bilayers results in an increase of the membrane conductance by many orders of magnitude. The membrane conductance is caused by the formation of large ion-permeable channels with a single-channel conductance in the order of 5 nS for 1 M alkali chlorides. The conductance has an ohmic current vs. voltage relationship. Further information on the structure of the pore formed by protein F was obtained by determining the single-channel conductance for various species differing in charge and size, and from zero-current potential measurements. The channel was found to be permeable for large organic ions (Tris⁺, N(C₂H₅)₄⁺, Hepes^?) and a channel diameter of 2.2 nm could be estimated from the conductance data (pore length of 7.5 nm). At neutral pH the pore is about two times more permeable for cations than for anions, possibly caused by negative charges in the pore. The consistent observation of large water filled pores formed by porin protein F in model membrane systems is discussed in the light of the known low permeability of the Ps. aeruginosa outer membrane towards antibiotics. It is suggested that this results from a relatively low proportion of open functional porin protein F pores in vivo.