Implicit solvent model studies of the interactions of the influenza hemagglutinin fusion peptide with lipid bilayers

The "fusion peptide," a segment of approximately 20 residues of the influenza hemagglutinin (HA), is necessary and sufficient for HA-induced membrane fusion. We used mean-field calculations of the free energy of peptide-membrane association (DeltaG(tot)) to deduce the most probable orientation of the fusion peptide in the membrane. The main contributions to DeltaG(tot) are probably from the electrostatic (DeltaG(el)) and nonpolar (DeltaG(np)) components of the solvation free energy; these were calculated using continuum solvent models. The peptide was described in atomic detail and was modeled as an alpha-helix based on spectroscopic data. The membrane's hydrocarbon region was described as a structureless slab of nonpolar medium embedded in water. All the helix-membrane configurations, which were lower in DeltaG(tot) than the isolated helix in the aqueous phase, were in the same (wide) basin in configurational space. In each, the helix was horizontally adsorbed at the water-bilayer interface with its principal axis parallel to the membrane plane, its hydrophobic face dissolved in the bilayer, and its polar face in the water. The associated DeltaG(tot) value was approximately -8 to -10 kcal/mol (depending on the rotameric state of one of the phenylalanine residues). In contrast, the DeltaG(tot) values associated with experimentally observed oblique orientations were found to be near zero, suggesting they are marginally stable at best. The theoretical model did not take into account the interactions of the polar headgroups with the peptide and peptide-induced membrane deformation effects. Either or both may overcompensate for the DeltaG(tot) difference between the horizontal and oblique orientations.     

Amino acid side-chain partition energies and distribution of residues in soluble proteins

Energies required to transfer amino acid side chains from water to less polar environments were calculated from results of several studies and compared with several statistical analyses of residue distributions in soluble proteins. An analysis that divides proteins into layers parallel with their surfaces is more informative than those that simply classify residues as exposed or buried. Most residues appear to be distributed as a function of the distance from the protein-water interface in a manner consistent with partition energies calculated from partitioning of amino acids between water and octanol phases and from solubilities of amino acids in water, ethanol, and methanol. Lys, Arg, Tyr, and Trp residues tend to concentrate near the water-protein interface where their apolar side-chain components are more buried than their polar side-chain components. Residue distributions calculated in this manner do not correlate well with side-chain solvation energies calculated from vapor pressures of side-chain analogs over a water phase. Results of statistical studies that classify residues as exposed to solvent or buried inside the protein interior appear to depend on the method used to classify residues. Data from some of these studies correlate better with solvation energies, but other data correlate better with partition energies. Most other statistical methods that have been used to evaluate effects of water on residue distributions yield results that correlate better with partition energies than with solvation energies.     

Quantification and characterization of P-glycoprotein-substrate interactions

It is generally accepted that P-glycoprotein binds its substrates in the lipid phase of the membrane. Quantification and characterization of the lipid-transporter binding step are, however, still a matter of debate. We therefore selected 15 structurally diverse drugs and measured the binding constants from water to the activating (inhibitory) binding region of P-glycoprotein, K(tw(1)) (K(tw(2))), as well as the lipid-water partition coefficients, K(lw). The former were obtained by measuring the concentrations of half-maximum activation (inhibition), K(1) (K(2)), in living NIH-MDR-G185 mouse embryo fibroblasts using a Cytosensor microphysiometer, and the latter were derived from surface activity measurements. This allowed determination of the membrane concentration of drugs at half-maximum P-glycoprotein activation (C(b(1)) = (0.02 to 67) mmol/L lipid), which is much higher than the corresponding aqueous concentration (K(1) = (0.02 to 376) microM). Moreover we determined the free energy of drug binding from water to the activating binding region of the transporter (DeltaG degrees (tw(1)) = (-30 to -54) kJ/mol), the free energy of drug partitioning into the lipid membrane (DeltaG degrees (lw) = (-23 to -34) kJ/mol), and, as the difference of the two, the free energy of drug binding from the lipid membrane to the activating binding region of the transporter (DeltaG degrees (tl(1)) = (-7 to -27) kJ/mol). For the compounds tested DeltaG degrees (tl(1)) was less negative than DeltaG degrees (lw) but varied more strongly. The free energies of substrate binding to the transporter within the lipid phase, DeltaG degrees (tl(1)), are consistent with a modular binding concept, where the energetically most efficient binding module comprises two hydrogen bond acceptor groups.     

Partitioning of ABA into bilayers of Di-saturated phosphatidylcholines as measured by DSC

Using differential scanning calorimetry, we have investigated partitioning of the plant hormone abscisic acid into a homologous series of di-saturated phosphatidylcholines increasing in chain length from C(14) to C(19). Partition coefficients calculated from the shift in T(m) range from 1280 for DiC(14)PC to 480 for DiC(19)PC. The free energy of transfer is chain-length independent with a value of DeltaG = -17.4 kJ/mol and an enthalpic contribution of DeltaH = -22.6 kJ/mol. The low net entropic contribution of -TDeltaS = -5.2 J/mol agrees with the concept of the bilayer effect, but differs from that of the entropy-driven classic hydrophobic effect valid for partitioning between bulk solvents. Preferential location of the hormone in the outer region of the membrane is indicated by characteristic changes in the transition profiles and by comparison with partitioning into organic solvents whose dielectric constants model the interior and exterior regions of the bilayer. Differences in partitioning and surface pKa between the biologically active ct-ABA and the inactive tt-isomer are discussed for biological relevance.     

Protein chemistry at membrane interfaces: non-additivity of electrostatic and hydrophobic interactions

Non-specific binding of proteins and peptides to charged membrane interfaces depends upon the combined contributions of hydrophobic (DeltaG(HPhi)) and electrostatic (DeltaG(ES)) free energies. If these are simply additive, then the observed free energy of binding (DeltaG(obs)) will be given by DeltaG(obs)=DeltaG(HPhi)+DeltaG(ES), where DeltaG(HPhi)=-sigma(NP)A(NP) and DeltaG(ES)=zFphi. In these expressions, A(NP) is the non-polar accessible area, sigma(NP) the non-polar solvation parameter, z the formal peptide valence, F the Faraday constant, and phi the membrane surface potential. But several lines of evidence suggest that hydrophobic and electrostatic binding free energies of proteins at membrane interfaces, such as those associated with cell signaling, are not simply additive. In order to explore this issue systematically, we have determined the interfacial partitioning free energies of variants of indolicidin, a cationic proline-rich antimicrobial peptide. The synthesized variants of the 13 residue peptide covered a wide range of hydrophobic free energies, which allowed us to examine the effect of hydrophobicity on electrostatic binding to membranes formed from mixtures of neutral and anionic lipids. Although DeltaG(obs) was always a linear function of DeltaG(HPhi), the slope depended upon anionic lipid content: the slope was 1.0 for pure, zwitterionic phosphocholine bilayers and 0.3 for pure phosphoglycerol membranes. DeltaG(obs) also varied linearly with surface potential, but the slope was smaller than the expected value, zF. As observed by others, this suggests an effective peptide valence z(eff) that is smaller than the formal valence z. Because of our systematic approach, we were able to establish a useful rule-of-thumb: z(eff) is reduced relative to z by about 20 % for each 3 kcal mol(-1) (1 kcal=4.184 kJ) favorable increase in DeltaG(HPhi). For neutral phosphocholine interfaces, we found that DeltaG(obs) could be predicted with remarkable accuracy using the Wimley-White experiment-based interfacial hydrophobicity scale.     

Drug-membrane interactions studied in phospholipid monolayers adsorbed on nonporous alkylated microspheres

Characterization of interactions with phospholipids is an integral part of the in vitro profiling of drug candidates because of the roles the interactions play in tissue accumulation and passive diffusion. Currently used test systems may inadequately emulate the bilayer core solvation properties (immobilized artificial membranes [IAM]), suffer from potentially slow transport of some chemicals (liposomes in free or immobilized forms), and require a tedious separation (if used for free liposomes). Here the authors introduce a well-defined system overcoming these drawbacks: nonporous octadecylsilica particles coated with a self-assembled phospholipid monolayer. The coating mimics the structure of the headgroup region, as well as the thickness and properties of the hydrocarbon core, more closely than IAM. The monolayer has a similar transition temperature pattern as the corresponding bilayer. The particles can be separated by filtration or a mild centrifugation. The partitioning equilibria of 81 tested chemicals were dissected into the headgroup and core contributions, the latter using the alkane/water partition coefficients. The deconvolution allowed a successful prediction of the bilayer/water partition coefficients with the standard deviation of 0.26 log units. The plate-friendly assay is suitable for high-throughput profiling of drug candidates without sacrificing the quality of analysis or details of the drug-phospholipid interactions.     

Thermodynamics of sodium dodecyl sulfate partitioning into lipid membranes

The partition equilibria of sodium dodecyl sulfate (SDS) and lithium dodecyl sulfate between water and bilayer membranes were investigated with isothermal titration calorimetry and spectroscopic methods (light scattering, (31)P-nuclear magnetic resonance) in the temperature range of 28 degrees C to 56 degrees C. The partitioning of the dodecyl sulfate anion (DS(-)) into the bilayer membrane is energetically favored by an exothermic partition enthalpy of Delta H(O)(D) = -6.0 kcal/mol at 28 degrees C. This is in contrast to nonionic detergents where Delta H(O)(D) is usually positive. The partition enthalpy decreases linearly with increasing temperature and the molar heat capacity is Delta C(O)(P) = -50 +/- 3 cal mol(-1) K(-1). The partition isotherm is nonlinear if the bound detergent is plotted versus the free detergent concentration in bulk solution. This is caused by the electrostatic repulsion between the DS(-) ions inserted into the membrane and those free in solution near the membrane surface. The surface concentration of DS(-) immediately above the plane of binding was hence calculated with the Gouy-Chapman theory, and a strictly linear relationship was obtained between the surface concentration and the extent of DS(-) partitioning. The surface partition constant K describes the chemical equilibrium in the absence of electrostatic effects. For the SDS-membrane equilibrium K was found to be 1.2 x 10(4) M(-1) to 6 x 10(4) M(-1) for the various systems and conditions investigated, very similar to data available for nonionic detergents of the same chain length. The membrane-micelle phase diagram was also studied. Complete membrane solubilization requires a ratio of 2.2 mol SDS bound per mole of total lipid at 56 degrees C. The corresponding equilibrium concentration of SDS free in solution is C (sat)(D,F) approximately 1.7 mM and is slightly below the critical micelles concentration (CMC) = 2.1 mM (at 56 degrees C and 0.11 M buffer). Membrane saturation occurs at approximately 0.3 mol SDS per mol lipid and the equilibrium SDS concentration is C (sat)(D,F)approximately equal 2.2 mM +/- 0.6 mM. SDS translocation across the bilayer is slow at ambient temperature but increases at high temperatures.     

High-performance liquid chromatographic chiral separation of beta2-homoamino acids

Reversed-phase high-performance liquid chromatographic methods were developed for the separation of enantiomers of eleven unnatural beta(2)-homoamino acids on chiral stationary phases containing macrocyclic glycopeptides (teicoplanin-containing Chirobiotic T and T2) or the macrocyclic peptide teicoplanin aglycone (Chirobiotic TAG) as chiral selectors. The effects of the organic modifier, the mobile phase composition, temperature, and the structures of the analytes on the separations were investigated. Separations were carried out at constant mobile phase compositions in temperature range 7-45 degrees C and the changes in enthalpy, Delta(DeltaH(o)), entropy, Delta(DeltaS(o)), and free energy, Delta(DeltaG(o)), were calculated. The -Delta(DeltaG(o)) values obtained on the three columns indicated that Chirobiotic TAG, without sugar units, may promote the interactions of the enantiomers of beta(2)-homoamino acids with branched alkyl or aryl side-chains, whereas for beta(2)-homoamino acids with alkyl side-chains Chirobiotic T and T2 seem to be more favorable. The elution sequence was determined in some cases and was observed to be R < S.     

Transport methods for probing the barrier domain of lipid bilayer membranes.

Two experimental techniques have been utilized to explore the barrier properties of lecithin/decane bilayer membranes with the aim of determining the contributions of various domains within the bilayer to the overall barrier. The thickness of lecithin/decane bilayers was systematically varied by modulating the chemical potential of decane in the annulus surrounding the bilayer using different mole fractions of squalene in decane. The dependence of permeability of a model permeant (acetamide) on the thickness of the solvent-filled region of the bilayer was assessed in these bilayers to determine the contribution of this region to the overall barrier. The flux of acetamide was found to vary linearly with bilayer area with Pm = (2.9 +/- 0.3) x 10(-4) cm s-1, after correcting for diffusion through unstirred water layers. The ratio between the overall membrane permeability coefficient and that calculated for diffusion through the hydrocarbon core in membranes having maximum thickness was 0.24, suggesting that the solvent domain contributes only slightly to the overall barrier properties. Consistent with these results, the permeability of acetamide was found to be independent of bilayer thickness. The relative contributions of the bilayer interface and ordered hydrocarbon regions to the transport barrier may be evaluated qualitatively by exploring the effective chemical nature of the barrier microenvironment. This may be probed by comparing functional group contributions to transport with those obtained for partitioning between water and various model bulk solvents ranging in polarity or hydrogen-bonding potential. A novel approach is described for obtaining group contributions to transport using ionizable permeants and pH adjustment. Using this approach, bilayer permeability coefficients of p-toluic acid and p-hydroxymethyl benzoic acid were determined to be 1.1 +/- 0.2 cm s-1 and (1.6 +/- 0.4) x 10(-3) cm s-1, respectively. From these values, the -OH group contribution to bilayer transport [delta(delta G0-OH)] was found to be 3.9 kcal/mol. This result suggests that the barrier region of the bilayer does not resemble the hydrogen-bonding environment found in octanol, but is somewhat less selective (more polar) than a hydrocarbon solvent.     

Hydrophobic ion hydration and the magnitude of the dipole potential

The magnitude of the dipole potential of lipid membranes is often estimated from the difference in conductance between the hydrophobic ions, tetraphenylborate, and tetraphenylarsonium or tetraphenylphosphonium. The calculation is based on the tetraphenylarsonium-tetraphenylborate hypothesis that the magnitude of the hydration energies of the anions and cations are equal (i.e., charge independent), so that their different rates of transport across the membrane are solely due to differential interactions with the membrane phase. Here we investigate the validity of this assumption by quantum mechanical calculations of the hydration energies. Tetraphenylborate (Delta G(hydr) = -168 kJ mol(-1)) was found to have a significantly stronger interaction with water than either tetraphenylarsonium (Delta G(hydr) = -145 kJ mol(-1)) or tetraphenylphosphonium (Delta G(hydr) = -157 kJ mol(-1)). Taking these differences into account, literature conductance data were recalculated to yield values of the dipole potential 57 to 119 mV more positive in the membrane interior than previous estimates. This may partly account for the discrepancy of at least 100 mV generally observed between dipole potential values calculated from lipid monolayers and those determined on bilayers.     

An evaluation of Poisson-Boltzmann electrostatic free energy calculations through comparison with experimental mutagenesis data

For systems involving highly and oppositely charged proteins, electrostatic forces dominate association and contribute to biomolecular complex stability. Using experimental or theoretical alanine-scanning mutagenesis, it is possible to elucidate the contribution of individual ionizable amino acids to protein association. We evaluated our electrostatic free energy calculations by comparing calculated and experimental data for alanine mutants of five protein complexes. We calculated Poisson-Boltzmann electrostatic free energies based on a thermodynamic cycle, which incorporates association in a reference (Coulombic) and solvated (solution) state, as well as solvation effects. We observe that Coulombic and solvation free energy values correlate with experimental data in highly and oppositely charged systems, but not in systems comprised of similarly charged proteins. We also observe that correlation between solution and experimental free energies is dependent on dielectric coefficient selection for the protein interior. Free energy correlations improve as protein dielectric coefficient increases, suggesting that the protein interior experiences moderate dielectric screening, despite being shielded from solvent. We propose that higher dielectric coefficients may be necessary to more accurately predict protein-protein association. Additionally, our data suggest that Coulombic potential calculations alone may be sufficient to predict relative binding of protein mutants.     

Mutational and structural-based analyses of the osmolyte effect on protein stability

It is known that several naturally occurring substances known as osmolytes increase the conformational stability of proteins. Bolen and co-worker proposed the osmophobic theory, which asserts the osmolyte effect occurs because of an unfavorable interaction of osmolytes mainly with the protein backbone, based on the results on the transfer Gibbs energy of amino acids (Deltag) [Bolen and Baskakov (2001) J. Mol. Biol. 310, 955-963]. In this paper, we report the effect of sarcosine on the conformational stability (DeltaG) of RNase Sa (96 residues and one disulfide bond) and four mutant proteins. The thermal denaturation curves for RNase Sa in sarcosine fitted a two-state model on nonlinear least-squares analysis. All the RNase Sa proteins were stabilized by sarcosine. For example, the increase in stability of the wild-type protein in 4 M sarcosine due to the osmolyte effect (Delta(o)DeltaG) is 3.2 kcal/mol. Mutational analysis of the osmolyte effect indicated that the changed Delta(o)DeltaG values upon mutation (Delta(m)Delta(o)DeltaG), as estimated from the Deltag values, are similar to the experimental values. Structural-based analysis of the osmolyte effect was also performed using model denatured structures: (a) a fully extended model (single chain) with no disulfide bond, (b) two-part, unfolded models (two chains) with a disulfide bond constructed through molecular dynamic (MD) simulation, and (c) a two-part, folded model (two chains). The two-part, unfolded models were expected to be more suitable as denatured structures. The Delta(o)DeltaG values calculated using the two-part, unfolded models were more consistent with experimental values than those calculated using the fully extended and two-part, folded models. This suggests that MD simulation is useful for testing denatured structures. These results indicate that the osmophobic theory can explain the osmolyte effect on protein stability.     

A kinetic folding intermediate probed by native state hydrogen exchange

Stopped-flow fluorescence studies on the N-terminal domain of rat CD2 (CD2.d1) have demonstrated that folding from the fully denatured state (U) proceeds via the transient accumulation of an apparent intermediate (I) in a so-called burst phase that precedes the rate-limiting transition leading to the native state (N). A previous pH-dependent equilibrium hydrogen exchange (HX) study identified a subset of amides in CD2.d1 which, under EX2 conditions, exchange from N with free energies greater than or equal to the free energy difference between the N and I states calculated from the stopped-flow data. Under EX1 conditions the rates of HX for these amides tend towards an asymptote that matches the global unfolding rate calculated from the stopped-flow data, suggesting that exchange for these amides requires traversing the N-to-I transition state barrier. Exchange for these amides presumably occurs from exchange-competent forms comprising the kinetic burst phase therefore. To explore this idea further, native state HX (NHX) data have been collected for CD2.d1 under EX2 conditions using denaturant concentrations which span either side of the denaturant concentration where, according to the stopped-flow data, the apparent U and I states are iso-energetic. The data fit to a two-component, sub-global (sg)/global (g) NHX mechanism, yielding Delta G and m value parameters (where the m value is a measure of hydrocarbon solvation). Regression analysis demonstrates that the (m(sg), Delta G(sg)) and (m(g), Delta G(g)) values calculated for this subset of amides correspond with those describing the kinetic burst phase transition. This result confirms the ability of the NHX technique to explore the structural and energetic properties of kinetic folding intermediates.     

Molecular regions controlling the activity of CNG channels

The alpha subunits of CNG channels of retinal photoreceptors (rod) and olfactory neurons (olf) are proteins that consist of a cytoplasmic NH(2) terminus, a transmembrane core region (including the segments S1-S6), and a cytoplasmic COOH terminus. The COOH terminus contains a cyclic nucleotide monophosphate binding domain NBD) that is linked by the C-linker (CL) to the core region. The binding of cyclic nucleotides to the NBD promotes channel opening by an allosteric mechanism. We examined why the sensitivity to cGMP is 22 times higher in olf than in rod by constructing chimeric channels and determining the [cGMP] causing half maximum channel activity (EC(50)). The characteristic difference in the EC(50) value between rod and olf was introduced by the NH(2) terminus and the core-CL region, whereas the NBD showed a paradoxical effect. The difference of the free energy difference Delta(DeltaG) was determined for each of these three regions with all possible combinations of the other two regions. For rod regions with respect to corresponding olf regions, the open channel conformation was destabilized by the NH(2) terminus (Delta(DeltaG) = -1.0 to -2.0 RT) and the core-CL region (Delta(DeltaG) = -2.0 to -2.9 RT), whereas it was stabilized by the NBD (Delta(DeltaG) = 0.3 to 1.1 RT). The NH(2) terminus deletion mutants of rod and olf differed by Delta(DeltaG) of only 0.9 RT, whereas the wild-type channels differed by the much larger value of 3.1 RT. The results show that in rod and olf, the NH(2) terminus, the core-CL region, and the NBD differ by characteristic Delta(DeltaG) values that do not depend on the specific composition of the other two regions and that the NH(2) terminus generates the main portion of Delta(DeltaG) between the wild-type channels.     

Partitioning of amphiphiles between coexisting ordered and disordered phases in two-phase lipid bilayer membranes

The partition coefficients (K(P)) of a series of single-chain and double-chain fluorescent amphiphiles, between solid ordered (P(beta') and L(beta)) and liquid disordered (L(alpha) of the type l(d)) lipid phases coexisting in the same lipid bilayer, was studied using steady-state fluorescence emission anisotropy. The single-chain amphiphiles were N-(7-nitrobenzoxa-2, 3-diazol-4-yl)-alkylamines, and the double-chain amphiphiles were N-(7-nitrobenzoxa-2, 3-diazol-4-yl)-phosphatidylethanolamines with chain lengths of 12-18 carbon atoms. Saturated 18-carbon alkyl/acyl chain compounds were also compared with Delta(9)-cis unsaturated chains of the same chain length. The fluorescence anisotropy of the probes was examined in lipid bilayers (multilamellar vesicles) prepared from an equimolar mixture of dilauroylphosphatidylcholine and distearoylphosphatidylcholine and studied as a function of temperature through the entire temperature range of coexistence of ordered gel phases and a disordered fluid phase in this system. The unsaturated chain amphiphiles partitioned exclusively into the fluid phase whenever this phase was present, as did the saturated chain amphiphiles with the shortest chains (C(12:0)), while K(P) ranges between 1 and 2, in favor of the L(beta) solid phase, for the amphiphiles with long saturated (C(18:0)) alkyl/acyl chains, with intermediate behavior for the intermediate chain lengths. All probes appeared to be totally excluded from P(beta') solid (gel) phases. The technique was also used to determine partitioning of some of the probes between coexisting liquid ordered (cholesterol-containing) (l(o)) and liquid disordered (l(d)) L(alpha) phases. In this case the ratio of signal amplitude to noise allowed us to obtain a qualitative, but not quantitative, measure of the phase partitioning of the probes. We conclude that the partitioning behavior of the probes examined between coexisting l(o) and l(d) phases is qualitatively similar to that observed between solid ordered and liquid disordered phases.     

Prediction of the aqueous solvation free energy of organic compounds by using autocorrelation of molecular electrostatic potential surface properties combined with response surface analysis

Several quantitative structure-property relationship (QSPR) approaches have been explored for the prediction of aqueous solubility or aqueous solvation free energies, DeltaG(sol), as crucial parameter affecting the pharmacokinetic profile and toxicity of chemical compounds. It is mostly accepted that aqueous solvation free energies can be expressed quantitatively in terms of properties of the molecular surface electrostatic potentials of the solutes. In the present study we have introduced autocorrelation molecular electrostatic potential (autoMEP) vectors in combination with nonlinear response surface analysis (RSA) as alternative 3D-QSPR strategy to evaluate the aqueous solvation free energy of organic compounds. A robust QSPR model (r(cv)=0.93) has been obtained by using a collection of 248 organic chemicals. An external test set based on 23 molecules confirmed the good predictivity of the autoMEP/RSA model suggesting its further applicability in the in silico prediction of water solubility of large organic compound libraries.     

Melittin-induced bilayer leakage depends on lipid material properties: evidence for toroidal pores

The membrane-lytic peptide melittin has previously been shown to form pores in lipid bilayers that have been described in terms of two different structural models. In the "barrel stave" model the bilayer remains more or less flat, with the peptides penetrating across the bilayer hydrocarbon region and aggregating to form a pore, whereas in the "toroidal pore" melittin induces defects in the bilayer such that the bilayer bends sharply inward to form a pore lined by both peptides and lipid headgroups. Here we test these models by measuring both the free energy of melittin transfer (DeltaG degrees ) and melittin-induced leakage as a function of bilayer elastic (material) properties that determine the energetics of bilayer bending, including the area compressibility modulus (K(a)), bilayer bending modulus (k(c)), and monolayer spontaneous curvature (R(o)). The addition of cholesterol to phosphatidylcholine (PC) bilayers, which increases K(a) and k(c), decreases both DeltaG degrees and the melittin-induced vesicle leakage. In contrast, the addition to PC bilayers of molecules with either positive R(o), such as lysoPC, or negative R(o), such as dioleoylglycerol, has little effect on DeltaG degrees , but produces large changes in melittin-induced leakage, from 86% for 8:2 PC/lysoPC to 18% for 8:2 PC/dioleoylglycerol. We observe linear relationships between melittin-induced leakage and both K(a) and 1/R(o)(2). However, in contrast to what would be expected for a barrel stave model, there is no correlation between observed leakage and bilayer hydrocarbon thickness. All of these results demonstrate the importance of bilayer material properties on melittin-induced leakage and indicate that the melittin-induced pores are defects in the bilayer lined in part by lipid molecules.     

The Barrier Domain for Solute Permeation Varies With Lipid Bilayer Phase Structure

The chemical selectivities of the transport barriers in lipid bilayers varying in composition and phase structure (gel-phase DPPC and DHPC bilayers and liquid-crystalline DPPC/CHOL/50:50 mol% bilayers) have been investigated by determining functional group contributions to transport of a series of α-substituted p-toluic acid analogs obtained in vesicle efflux experiments. Linear free energy relationships are established between the free energies of transfer for this series of compounds from water to the barrier domain and corresponding values for their transfer from water into six model bulk solvents (hexadecane, hexadecene, decadiene, chlorobutane, butyl ether, and octanol) determined in partitioning experiments to compare the barrier microenvironment to that in these model solvents. The barrier microenvironment in all bilayers studied is substantially more hydrophobic than octanol, thus establishing the location of the barrier beyond the hydrated headgroup interfacial region, as the interface is expected to be more hydrophilic than octanol. The chemical nature of the barrier domain microenvironment varies with bilayer phase structure. The barrier regions in non-interdigitated DPPC and interdigitated DHPC gel-phase bilayers exhibit some degree of hydrogen-bond acceptor capacity as may occur if these domains lie in the vicinity of the ester/ether linkages between the headgroups and the acyl chains. Intercalation of 50 mol% cholesterol into DPPC bilayers, which induces a phase transition to a liquid-crystalline phase, substantially increases the apparent barrier domain hydrophobicity relative to gel-phase bilayers to a nonhydrogen bonding, hydrocarbonlike environment resembling hexadecene. This result, combined with similar observations in liquid-crystalline egg-PC bilayers (J. Pharm. Sci. (1994), 83:1511–1518), supports the notion that the transition from the gel-phase to liquid-crystalline phase shifts the barrier domain further into the bilayer interior (i.e., deeper within the ordered chain region). Received: 16 September 1997/Revised: 14 May 1998     

Bilayer interfacial properties modulate the binding of amphipathic peptides

The free energy of transfer (DeltaG degrees ) from water to lipid bilayers was measured for two amphipathic peptides, the presequence of the mitochondrial peptide rhodanese (MPR) and melittin. Experiments were designed to determine the effects on peptide partitioning of the addition of lipids that produce structural modifications to the bilayer/water interface. In particular, the addition of cholesterol or the cholesterol analog 6-ketocholestanol increases the bilayer area compressibility modulus, indicating that these molecules modify lipid-lipid interactions in the plane of the bilayer. The addition of 6-ketocholestanol or lipids with attached polyethylene glycol chains (PEG-lipids) modify the effective thickness of the interfacial region; 6-ketocholestanol increases the width of hydrophilic headgroup region in the direction of the acyl chains whereas the protruding PEG chains of PEG-lipids increase the structural width of the headgroup region into the surrounding aqueous phase. The incorporation of PEG-lipids with PEG molecular weights of 2000 or 5000 had no appreciable effect on peptide partitioning that could not be accounted for by the presence of surface charge. However, for both MPR and melittin DeltaG degrees decreased linearly with increasing bilayer compressibility modulus, demonstrating the importance of bilayer mechanical properties in the binding of amphipathic peptides.