Experimental evidence for hydrophobic matching and membrane-mediated interactions in lipid bilayers containing gramicidin.

Hydrophobic matching, in which transmembrane proteins cause the surrounding lipid bilayer to adjust its hydrocarbon thickness to match the length of the hydrophobic surface of the protein, is a commonly accepted idea in membrane biophysics. To test this idea, gramicidin (gD) was embedded in 1, 2-dilauroyl-sn-glycero-3-phosphocholine (DLPC) and 1, 2-myristoyl-sn-glycero-3-phosphocholine (DMPC) bilayers at the peptide/lipid molar ratio of 1:10. Circular dichroism (CD) was measured to ensure that the gramicidin was in the beta6.3 helix form. The bilayer thickness (the phosphate-to-phosphate distance, or PtP) was measured by x-ray lamellar diffraction. In the Lalpha phase near full hydration, PtP is 30.8 A for pure DLPC, 32.1 A for the DLPC/gD mixture, 35.3 A for pure DMPC, and 32.7 A for the DMPC/gD mixture. Gramicidin apparently stretches DLPC and thins DMPC toward a common thickness as expected by hydrophobic matching. Concurrently, gramicidin-gramicidin correlations were measured by x-ray in-plane scattering. In the fluid phase, the gramicidin-gramicidin nearest-neighbor separation is 26.8 A in DLPC, but shortens to 23.3 A in DMPC. These experiments confirm the conjecture that when proteins are embedded in a membrane, hydrophobic matching creates a strain field in the lipid bilayer that in turn gives rise to a membrane-mediated attractive potential between proteins.     

Influence of hydrophobic mismatch on structures and dynamics of gramicidin a and lipid bilayers

Gramicidin A (gA) is a 15-amino-acid antibiotic peptide with an alternating L-D sequence, which forms (dimeric) bilayer-spanning, monovalent cation channels in biological membranes and synthetic bilayers. We performed molecular dynamics simulations of gA dimers and monomers in all-atom, explicit dilauroylphosphatidylcholine (DLPC), dimyristoylphosphatidylcholine (DMPC), dioleoylphosphatidylcholine (DOPC), and 1-palmitoyl-2-oleoyl-phosphatidylcholine (POPC) bilayers. The variation in acyl chain length among these different phospholipids provides a way to alter gA-bilayer interactions by varying the bilayer hydrophobic thickness, and to determine the influence of hydrophobic mismatch on the structure and dynamics of both gA channels (and monomeric subunits) and the host bilayers. The simulations show that the channel structure varied little with changes in hydrophobic mismatch, and that the lipid bilayer adapts to the bilayer-spanning channel to minimize the exposure of hydrophobic residues. The bilayer thickness, however, did not vary monotonically as a function of radial distance from the channel. In all simulations, there was an initial decrease in thickness within 4–5 Å from the channel, which was followed by an increase in DOPC and POPC or a further decrease in DLPC and DMPC bilayers. The bilayer thickness varied little in the monomer simulations—except one of three independent simulations for DMPC and all three DLPC simulations, where the bilayer thinned to allow a single subunit to form a bilayer-spanning water-permeable pore. The radial dependence of local lipid area and bilayer compressibility is also nonmonotonic in the first shell around gA dimers due to gA-phospholipid interactions and the hydrophobic mismatch. Order parameters, acyl chain dynamics, and diffusion constants also differ between the lipids in the first shell and the bulk. The lipid behaviors in the first shell around gA dimers are more complex than predicted from a simple mismatch model, which has implications for understanding the energetics of membrane protein-lipid interactions.     

Molecular dynamics simulations of model trans-membrane peptides in lipid bilayers: a systematic investigation of hydrophobic mismatch

Hydrophobic mismatch, which is the difference between the hydrophobic length of trans-membrane segments of a protein and the hydrophobic width of the surrounding lipid bilayer, is known to play a role in membrane protein function. We have performed molecular dynamics simulations of trans-membrane KALP peptides (sequence: GKK(LA)nLKKA) in phospholipid bilayers to investigate hydrophobic mismatch alleviation mechanisms. By varying systematically the length of the peptide (KALP15, KALP19, KALP23, KALP27, and KALP31) and the lipid hydrophobic length (DLPC, DMPC, and DPPC), a wide range of mismatch conditions were studied. Simulations of durations of 50-200 ns show that under positive mismatch, the system alleviates the mismatch predominantly by tilting the peptide and to a smaller extent by increased lipid ordering in the immediate vicinity of the peptide. Under negative mismatch, alleviation takes place by a combination of local bilayer bending and the snorkeling of the lysine residues of the peptide. Simulations performed at a higher peptide/lipid molar ratio (1:25) reveal slower dynamics of both the peptide and lipid relative to those at a lower peptide/lipid ratio (1:128). The lysine residues have favorable interactions with specific oxygen atoms of the phospholipid headgroups, indicating the preferred localization of these residues at the lipid/water interface.     

Protein-induced bilayer perturbations: Lipid ordering and hydrophobic coupling

The host lipid bilayer is increasingly being recognized as an important non-specific regulator of membrane protein function. Despite considerable progress the interplay between hydrophobic coupling and lipid ordering is still elusive. We use electron spin resonance (ESR) to study the interaction between the model protein gramicidin and lipid bilayers of varying thickness. The free energy of the interaction is up to −6 kJ/mol; thus not strongly favored over lipid-lipid interactions. Incorporation of gramicidin results in increased order parameters with increased protein concentration and hydrophobic mismatch. Our findings also show that at high protein:lipid ratios the lipids are motionally restricted but not completely immobilized. Both exchange on and off rate values for the lipid ↔ gramicidin interaction are lowest at optimal hydrophobic matching. Hydrophobic mismatch of few Å results in up to 10-fold increased exchange rates as compared to the ‘optimal’ match situation pointing to the regulatory role of hydrophobic coupling in lipid-protein interactions.     

Modulation of gramicidin channel conformation and organization by hydrophobic mismatch in saturated phosphatidylcholine bilayers

The matching of hydrophobic lengths of integral membrane proteins and the surrounding lipid bilayer is an important factor that influences both structure and function of integral membrane proteins. The ion channel gramicidin is known to be uniquely sensitive to membrane properties such as bilayer thickness and membrane mechanical properties. The functionally important carboxy terminal tryptophan residues of gramicidin display conformation-dependent fluorescence which can be used to monitor gramicidin conformations in membranes [S.S. Rawat, D.A. Kelkar, A. Chattopadhyay, Monitoring gramicidin conformations in membranes: a fluorescence approach, Biophys. J. 87 (2004) 831-843]. We have examined the effect of hydrophobic mismatch on the conformation and organization of gramicidin in saturated phosphatidylcholine bilayers of varying thickness utilizing the intrinsic conformation-dependent tryptophan fluorescence. Our results utilizing steady state and time-resolved fluorescence spectroscopic approaches, in combination with circular dichroism spectroscopy, show that gramicidin remains predominantly in the channel conformation and gramicidin tryptophans are at the membrane interfacial region over a range of mismatch conditions. Interestingly, gramicidin conformation shifts toward non-channel conformations in extremely thick gel phase membranes although it is not excluded from the membrane. In addition, experiments utilizing self quenching of tryptophan fluorescence indicate peptide aggregation in thicker gel phase membranes.     

Modulation of gramicidin channel conformation and organization by hydrophobic mismatch in saturated phosphatidylcholine bilayers

The matching of hydrophobic lengths of integral membrane proteins and the surrounding lipid bilayer is an important factor that influences both structure and function of integral membrane proteins. The ion channel gramicidin is known to be uniquely sensitive to membrane properties such as bilayer thickness and membrane mechanical properties. The functionally important carboxy terminal tryptophan residues of gramicidin display conformation-dependent fluorescence which can be used to monitor gramicidin conformations in membranes [S.S. Rawat, D.A. Kelkar, A. Chattopadhyay, Monitoring gramicidin conformations in membranes: a fluorescence approach, Biophys. J. 87 (2004) 831-843]. We have examined the effect of hydrophobic mismatch on the conformation and organization of gramicidin in saturated phosphatidylcholine bilayers of varying thickness utilizing the intrinsic conformation-dependent tryptophan fluorescence. Our results utilizing steady state and time-resolved fluorescence spectroscopic approaches, in combination with circular dichroism spectroscopy, show that gramicidin remains predominantly in the channel conformation and gramicidin tryptophans are at the membrane interfacial region over a range of mismatch conditions. Interestingly, gramicidin conformation shifts toward non-channel conformations in extremely thick gel phase membranes although it is not excluded from the membrane. In addition, experiments utilizing self quenching of tryptophan fluorescence indicate peptide aggregation in thicker gel phase membranes.     

Solvent history dependence of gramicidin A conformations in hydrated lipid bilayers.

Reconstituted lipid bilayers of dimyristoylphosphatidylcholine (DMPC) and gramicidin A' have been prepared by cosolubilizing gramicidin and DMPC in one of three organic solvent systems followed by vacuum drying and hydration. The conformational state of gramicidin as characterized by 23Na NMR, circular dichroism, and solid state 15N NMR is dependent upon the cosolubilizing solvent system. In particular, two conformational states are described; a state in which Na+ has minimal interactions with the polypeptide, referred to as a nonchannel state, and a state in which Na+ interacts very strongly with the polypeptide, referred to as the channel state. Both of these conformations are intimately associated with the hydrophobic core of the lipid bilayer. Furthermore, both of these states are stable in the bilayer at neutral pH and at a temperature above the bilayer phase transition temperature. These results with gramicidin suggest that the conformation of membrane proteins may be dictated by the conformation before membrane insertion and may be dependent upon the mechanism by which the insertion is accomplished.     

A one-dimensional continuum elastic model for membrane-embedded gramicidin dimer dissociation

Membrane elastic properties, which are subject to alteration by compounds such as cholesterol, lipid metabolites and other amphiphiles, as well as pharmaceuticals, can have important effects on membrane proteins. A useful tool for measuring some of these effects is the gramicidin A channels, which are formed by transmembrane dimerization of non-conducting subunits that reside in each bilayer leaflet. The length of the conducting channels is less than the bilayer thickness, meaning that channel formation is associated with a local bilayer deformation. Electrophysiological studies have shown that the dimer becomes increasingly destabilized as the hydrophobic mismatch between the channel and the host bilayer increases. That is, the bilayer imposes a disjoining force on the channel, which grows larger with increasing hydrophobic mismatch. The energetic analysis of the channel-bilayer coupling is usually pursued assuming that each subunit, as well as the subunit-subunit interface, is rigid. Here we relax the latter assumption and explore how the bilayer junction responds to changes in this disjoining force using a simple one-dimensional energetic model, which reproduces key features of the bilayer regulation of gramicidin channel lifetimes.     

Interfacial polar interactions affect gramicidin channel kinetics

Critical to biological processes such as secretion and transport, protein-lipid interactions within the membrane and at the membrane-water interface still raise many questions. Here we examine the role of lipid headgroups in these interactions by using gramicidin A (gA) channels in planar bilayers as a probe. We show that although headgroup demethylation from phosphatidylcholine (DOPC) to phosphatidylethanolamine decreases the lifetime of gA channels by an order of magnitude in accordance with the currently accepted hydrophobic mismatch mechanism, our findings with diether-DOPC suggest the importance of the headgroup-peptide interactions. According to our x-ray diffraction measurements, this lipid has the same hydrophobic thickness as DOPC but increases gA lifetime by a factor of 2. Thus we demonstrate that peptide-headgroup interactions may dominate over the effect of hydrophobic mismatch in regulating protein function.     

Membrane topology of a 14-mer model amphipathic peptide: a solid-state NMR spectroscopy study

We have investigated the interaction between a synthetic amphipathic 14-mer peptide and model membranes by solid-state NMR. The 14-mer peptide is composed of leucines and phenylalanines modified by the addition of crown ethers and forms a helical amphipathic structure in solution and bound to lipid membranes. To shed light on its membrane topology, 31P, 2H, 15N solid-state NMR experiments have been performed on the 14-mer peptide in interaction with mechanically oriented bilayers of dilauroylphosphatidylcholine (DLPC), dimyristoylphosphatidylcholine (DMPC), and dipalmitoylphosphatidylcholine (DPPC). The 31P, 2H, and 15N NMR results indicate that the 14-mer peptide remains at the surface of the DLPC, DMPC, and DPPC bilayers stacked between glass plates and perturbs the lipid orientation relative to the magnetic field direction. Its membrane topology is similar in DLPC and DMPC bilayers, whereas the peptide seems to be more deeply inserted in DPPC bilayers, as revealed by the greater orientational and motional disorder of the DPPC lipid headgroup and acyl chains. 15N{31P} rotational echo double resonance experiments have also been used to measure the intermolecular dipole-dipole interaction between the 14-mer peptide and the phospholipid headgroup of DMPC multilamellar vesicles, and the results indicate that the 14-mer peptide is in contact with the polar region of the DMPC lipids. On the basis of these studies, the mechanism of membrane perturbation of the 14-mer peptide is associated to the induction of a positive curvature strain induced by the peptide lying on the bilayer surface and seems to be independent of the bilayer hydrophobic thickness.     

Regulation of sodium channel function by bilayer elasticity: the importance of hydrophobic coupling. Effects of Micelle-forming amphiphiles and cholesterol

Membrane proteins are regulated by the lipid bilayer composition. Specific lipid-protein interactions rarely are involved, which suggests that the regulation is due to changes in some general bilayer property (or properties). The hydrophobic coupling between a membrane-spanning protein and the surrounding bilayer means that protein conformational changes may be associated with a reversible, local bilayer deformation. Lipid bilayers are elastic bodies, and the energetic cost of the bilayer deformation contributes to the total energetic cost of the protein conformational change. The energetics and kinetics of the protein conformational changes therefore will be regulated by the bilayer elasticity, which is determined by the lipid composition. This hydrophobic coupling mechanism has been studied extensively in gramicidin channels, where the channel-bilayer hydrophobic interactions link a "conformational" change (the monomer<-->dimer transition) to an elastic bilayer deformation. Gramicidin channels thus are regulated by the lipid bilayer elastic properties (thickness, monolayer equilibrium curvature, and compression and bending moduli). To investigate whether this hydrophobic coupling mechanism could be a general mechanism regulating membrane protein function, we examined whether voltage-dependent skeletal-muscle sodium channels, expressed in HEK293 cells, are regulated by bilayer elasticity, as monitored using gramicidin A (gA) channels. Nonphysiological amphiphiles (beta-octyl-glucoside, Genapol X-100, Triton X-100, and reduced Triton X-100) that make lipid bilayers less "stiff", as measured using gA channels, shift the voltage dependence of sodium channel inactivation toward more hyperpolarized potentials. At low amphiphile concentration, the magnitude of the shift is linearly correlated to the change in gA channel lifetime. Cholesterol-depletion, which also reduces bilayer stiffness, causes a similar shift in sodium channel inactivation. These results provide strong support for the notion that bilayer-protein hydrophobic coupling allows the bilayer elastic properties to regulate membrane protein function.     

Orientation and lipid-peptide interactions of gramicidin A in lipid membranes: polarized attenuated total reflection infrared spectroscopy and spin-label electron spin resonance

Gramicidin A was incorporated at a peptide/lipid ratio of 1:10 mol/mol in aligned bilayers of dimyristoyl phosphatidylcholine (DMPC), phosphatidylserine (DMPS), phosphatidylglycerol (DMPG), and phosphatidylethanolamine (DMPE), from trifluoroethanol. Orientations of the peptide and lipid chains were determined by polarized attenuated total reflection infrared spectroscopy. Lipid-peptide interactions with gramicidin A in DMPC bilayers were studied with different spin-labeled lipid species by using electron spin resonance spectroscopy. In DMPC membranes, the orientation of the lipid chains is comparable to that in the absence of peptide, in both gel and fluid phases. In gel-phase DMPC, the effective tilt of the peptide exceeds that of the lipid chains, but in the fluid phase both are similar. For gramicidin A in DMPS, DMPG, and DMPE, the degree of orientation of the peptide and lipid chains is less than in DMPC. In the fluid phase of DMPS, DMPG, and DMPE, gramicidin A is also less well oriented than are the lipid chains. In DMPE especially, gramicidin A is largely disordered. In DMPC membranes, three to four lipids per monomer experience direct motional restriction on interaction with gramicidin A. This is approximately half the number of lipids expected to contact the intramembranous perimeter of the gramicidin A monomer. A selectivity for certain negatively charged lipids is found in the interaction with gramicidin A in DMPC. These results are discussed in terms of the integration of gramicidin A channels in lipid bilayers, and of the interactions of lipids with integral membrane proteins.     

Gramicidin structure and disposition in highly curved membranes

With a view to deciphering aspects of the mechanism of membrane protein crystallization in lipidic mesophases (in meso crystallization), an examination of the structure and disposition of the pore-forming peptide, gramicidin, in the lipidic cubic phase was undertaken. At its simplest, the cubic phase consists of lipid and water in the form of a molecular 'sponge.' The lipid exists as a continuous, highly curved bilayer that divides the aqueous component into two interpenetrating but non-contacting channels. In this study, we show that gramicidin reconstitutes into the lipid bilayer of the cubic phase and that it adopts the channel, or helical dimer, conformation therein. Fluorescence quenching with brominated lipid was used to establish the bilayer location of the peptide. Electronic absorption and emission spectroscopies corroborated this finding. Peptide conformation in the cubic phase membrane was determined by circular dichroism. The identity and microstructure of the mesophases, and their capacity to accommodate gramicidin and other system components (sodium dodecyl sulfate, trifluoroethanol), was established by small-angle X-ray diffraction. Beyond a limiting concentration, gramicidin destabilized the cubic phase in favor of the inverted hexagonal phase. While gramicidin remained bilayer bound as membrane thickness changed, its conformation responded to the degree of bilayer mismatch with the hydrophobic surface of the peptide. These findings support the hypothesis that reconstitution into the lipid bilayer is an integral part of the in meso crystallization process as applied to membrane proteins. They also suggest ways for improving the process of membrane protein crystallogenesis.     

X-ray scattering with momentum transfer in the plane of membrane. Application to gramicidin organization.

We demonstrate a technique for measuring x-ray (or neutron) scattering with the momentum transfer confined in the plane of membrane, for the purpose of studying lateral organization of proteins and peptides in membrane. Unlike freeze-fracture electron microscopy or atomic force microscopy which requires the membrane to be frozen or fixed, in-plane x-ray scattering can be performed with the membrane maintained in the liquid crystalline state. As an example, the controversial question of whether gramicidin forms aggregates in membrane was investigated. We used dilauroylphosphatidylcholine (DLPC) bilayers containing gramicidin in the molar ratio of 10:1. Very clear scattering curves reflecting gramicidin channel-channel correlation were obtained, even for the sample containing no heavy atoms. Thallium ions bound to gramicidin channels merely increase the magnitude of the scattering curve. Analysis of the data shows that the channels were randomly distributed in the membrane, similar to a computer simulation of freely moving disks in a plane. We suggest that oriented proteins may provide substantial x-ray contrast against the lipid background without requiring heavy-atom labeling. This should open up many possible new experiments.     

Molecular sorting of lipids by bacteriorhodopsin in dilauroylphosphatidylcholine/distearoylphosphatidylcholine lipid bilayers.

A combined experimental and theoretical study is performed on binary dilauroylphosphatidylcholine/distearoylphosphatidylcholine (DLPC/DSPC) lipid bilayer membranes incorporating bacteriorhodopsin (BR). The system is designed to investigate the possibility that BR, via a hydrophobic matching principle related to the difference in lipid bilayer hydrophobic thickness and protein hydrophobic length, can perform molecular sorting of the lipids at the lipid-protein interface, leading to lipid specificity/selectivity that is controlled solely by physical factors. The study takes advantage of the strongly nonideal mixing behavior of the DLPC/DSPC mixture and the fact that the average lipid acyl-chain length is strongly dependent on temperature, particularly in the main phase transition region. The experiments are based on fluorescence energy transfer techniques using specifically designed lipid analogs that can probe the lipid-protein interface. The theoretical calculations exploit a microscopic molecular interaction model that embodies the hydrophobic matching as a key parameter. At low temperatures, in the gel-gel coexistence region, experimental and theoretical data consistently indicate that BR is associated with the short-chain lipid DLPC. At moderate temperatures, in the fluid-gel coexistence region, BR remains in the fluid phase, which is mainly composed of short-chain lipid DLPC, but is enriched at the interface between the fluid and gel domains. At high temperatures, in the fluid phase, BR stays in the mixed lipid phase, and the theoretical data suggest a preference of the protein for the long-chain DSPC molecules at the expense of the short-chain DLPC molecules. The combined results of the experiments and the calculations provide evidence that a molecular sorting principle is active because of hydrophobic matching and that BR exhibits physical lipid selectivity. The results are discussed in the general context of membrane organization and compartmentalization and in terms of nanometer-scale lipid-domain formation.     

Hydrophobic matching controls the tilt and stability of the dimeric platelet-derived growth factor receptor (PDGFR) β transmembrane segment

The platelet-derived growth factor receptor β is a member of the cell surface receptor tyrosine kinase family and dimerizes upon activation. We determined the structure of the transmembrane segment in dodecylphosphocholine micelles by liquid-state NMR and found that it forms a stable left-handed helical dimer. Solid-state NMR and oriented circular dichroism were used to measure the tilt angle of the helical segments in macroscopically aligned model membranes with different acyl chain lengths. Both methods showed that decreasing bilayer thickness (DEPC-POPC-DMPC) led to an increase in the helix tilt angle from 10° to 30° with respect to the bilayer normal. At the same time, reconstitution of the comparatively long hydrophobic segment became less effective, eventually resulting in complete protein aggregation in the short-chain lipid DLPC. Unrestrained molecular dynamics simulations of the dimer were carried out in explicit lipid bilayers (DEPC, POPC, DMPC, sphingomyelin), confirming the observed dependence of the helix tilt angle on bilayer thickness. Notably, molecular dynamics revealed that the left-handed dimer gets tilted en bloc, whereas conformational transitions to alternative (e.g. right-handed dimeric) states were not supported. The experimental data along with the simulation results demonstrate a pronounced interplay between the platelet-directed growth factor receptor β transmembrane segment and the bilayer thickness. The effect of hydrophobic mismatch might play a key role in the redistribution and activation of the receptor within different lipid microdomains of the plasma membrane in vivo.     

Thermodynamics of mechanosensitivity

Mechanosensitivity of ion channels is conventionally interpreted as being driven by a change of their in-plane cross-sectional area A(msc). This, however, does not include any factors relating to membrane stiffness, thickness, spontaneous curvature or changes in channel shape, length or stiffness. Because the open probability of a channel is sensitive to all these factors, we constructed a general thermodynamic formalism. These equations provide the basis for the analysis of the behaviour of mechanosensitive channels in lipids of different geometric and chemical properties such as the hydrophobic mismatch at the boundary between the protein and lipid or the effects of changes in the bilayer intrinsic curvature caused by the adsorption of amphipaths. This model predicts that the midpoint gamma(1/2) and the slope(1/2) of the gating curve are generally not independent. Using this relationship, we have predicted the line tension at the channel/lipid border of MscL as approximately 10 pN, and found it to be much less than the line tension of aqueous pores in pure lipid membranes. The MscL channel appears quite well matched to its lipid environment. Using gramicidin as a model system, we have explained its observed conversion from stretch-activated to stretch-inactivated gating as a function of bilayer thickness and composition.     

Comparative Molecular Dynamics Simulation Studies of Protegrin-1 Monomer and Dimer in Two Different Lipid Bilayers

Antimicrobial peptides interact specifically with the membrane of a pathogen and kill the pathogen by releasing its cellular contents. Protegrin-1 (PG-1), a β-hairpin antimicrobial peptide, is known to exist as a transmembrane monomer in a 1,2-dilauroylphosphatidylcholine (DLPC) bilayer and shows concentration-dependent oligomerization in a 1-palmitoyl-2-oleoylphosphatidylcholine (POPC) bilayer. To examine its structure, dynamics, orientation, and interaction in membranes, we performed comparative molecular dynamics simulations of PG-1 monomer and dimer in DLPC and POPC bilayers for a total of 840 ns. The PG-1 monomer exhibits larger tilting in DLPC than in POPC due to a hydrophobic mismatch. PG-1 tilting is dependent on its rotation angle. The specific orientation of PG-1 in membranes is governed by the interactions of its aromatic residues with lipid headgroups. The calculated ¹⁵N and ¹³CO chemical shifts of Val¹⁶ in DLPC reveal that there are different sets of tilt and rotation angles that satisfy the experimental values reasonably, suggesting that more experiments are needed to determine its orientation. The dimer simulations show that the dimer interface is better preserved in POPC than in DLPC because POPC's greater hydrophobic thickness causes reduced flexibility of the C-terminal strands. Both monomer and dimer simulations show membrane thinning around PG-1, largely due to arginine-lipid interactions.     

Stability of an ion channel in lipid bilayers: implicit solvent model calculations with gramicidin

Gramicidin is a helical peptide, 15 residues in length, which dimerizes to form ion-conducting channels in lipid bilayers. Here we report calculations of its free energy of transfer from the aqueous phase into bilayers of different widths. The electrostatic and nonpolar contributions to the desolvation free energy were calculated using implicit solvent models, in which gramicidin was described in atomic detail and the hydrocarbon region of the membrane was described as a slab of hydrophobic medium embedded in water. The free energy penalties from the lipid perturbation and membrane deformation effects, and the entropy loss associated with gramicidin immobilization in the bilayer, were estimated from a statistical thermodynamic model of the bilayer. The calculations were carried out using two classes of experimentally observed conformations: a head-to-head dimer of two single-stranded (SS) beta-helices and a double-stranded (DS) intertwined double helix. The calculations showed that gramicidin is likely to partition into the bilayer in all of these conformations. However, the SS conformation was found to be significantly more stable than the DS in the bilayer, in agreement with most of the experimental data. We tested numerous transmembrane and surface orientations of gramicidin in bilayers of various widths. Our calculations indicate that the most favorable orientation is transmembrane, which is indeed to be expected from a channel-forming peptide. The calculations demonstrate that gramicidin insertion into the membrane is likely to involve a significant deformation of the bilayer to match the hydrophobic width of the peptide (22 A), again in good agreement with experimental data. Interestingly, deformation of the bilayer was induced by all of the gramicidin conformations.     

A correlation between lipid domain shape and binary phospholipid mixture composition in free standing bilayers: A two-photon fluorescence microscopy study

Giant unilamellar vesicles (GUVs) composed of different phospholipid binary mixtures were studied at different temperatures, by a method combining the sectioning capability of the two-photon excitation fluorescence microscope and the partition and spectral properties of 6-dodecanoyl-2-dimethylamino-naphthalene (Laurdan) and Lissamine rhodamine B 1,2-dihexadecanoyl-sn-glycero-3-phosphoethanolamine (N-Rh-DPPE). We analyzed and compared fluorescence images of GUVs composed of 1,2-dilauroyl-sn-glycero-3-phosphocholine/1, 2-dipalmitoyl-sn-glycero-3-phosphocholine (DLPC/DPPC), 1, 2-dilauroyl-sn-glycero-3-phosphocholine/1, 2-distearoyl-sn-glycero-3-phosphocholine (DLPC/DSPC), 1, 2-dilauroyl-sn-glycero-3-phosphocholine/1, 2-diarachidoyl-sn-glycero-3-phosphocholine (DLPC/DAPC), 1, 2-dimyristoyl-sn-glycero-3-phosphocholine/1, 2-distearoyl-sn-glycero-3-phosphocholine (DMPC/DSPC) (1:1 mol/mol in all cases), and 1,2-dimyristoyl-sn-glycero-3-phosphoethanolamine/1, 2-dimyristoyl-sn-glycero-3-phosphocholine (DMPE/DMPC) (7:3 mol/mol) at temperatures corresponding to the fluid phase and the fluid-solid phase coexistence. In addition, we studied the solid-solid temperature regime for the DMPC/DSPC and DMPE/DMPC mixtures. From the Laurdan intensity images the generalized polarization function (GP) was calculated at different temperatures to characterize the phase state of the lipid domains. We found a homogeneous fluorescence distribution in the GUV images at temperatures corresponding to the fluid region for all of the lipid mixtures. At temperatures corresponding to phase coexistence we observed concurrent fluid and solid domains in the GUVs independent of the lipid mixture. In all cases the lipid solid domains expanded and migrated around the vesicle surface as we decreased the temperature. The migration of the solid domains decreased dramatically at temperatures close to the solid-fluid-->solid phase transition. For the DLPC-containing mixtures, the solid domains showed line, quasicircular, and dendritic shapes as the difference in the hydrophobic chain length between the components of the binary mixture increases. In addition, for the saturated PC-containing mixtures, we found a linear relationship between the GP values for the fluid and solid domains and the difference between the hydrophobic chain length of the binary mixture components. Specifically, at the phase coexistence temperature region the difference in the GP values, associated with the fluid and solid domains, increases as the difference in the chain length of the binary mixture component increases. This last finding suggests that in the solid-phase domains, the local concentration of the low melting temperature phospholipid component increases as the hydrophobic mismatch decreases. At the phase coexistence temperature regime and based on the Laurdan GP data, we observe that when the hydrophobic mismatch is 8 (DLPC/DAPC), the concentration of the low melting temperature phospholipid component in the solid domains is negligible. This last observation extends to the saturated PE/PC mixtures at the phase coexistence temperature range. For the DMPC/DSPC we found that the nonfluorescent solid regions gradually disappear in the solid temperature regime of the phase diagram, suggesting lipid miscibility. This last result is in contrast with that found for DMPE/DMPC mixtures, where the solid domains remain on the GUV surface at temperatures corresponding to that of the solid region. In all cases the solid domains span the inner and outer leaflets of the membrane, suggesting a strong coupling between the inner and outer monolayers of the lipid membrane. This last finding extends previous observations of GUVs composed of DPPE/DPPC and DLPC/DPPC mixtures (, Biophys. J. 78:290-305).