Simulation studies of protein-induced bilayer deformations, and lipid-induced protein tilting, on a mesoscopic model for lipid bilayers with embedded proteins

Biological membranes are complex and highly cooperative structures. To relate biomembrane structure to their biological function it is often necessary to consider simpler systems. Lipid bilayers composed of one or two lipid species, and with embedded proteins, provide a model system for biological membranes. Here we present a mesoscopic model for lipid bilayers with embedded proteins, which we have studied with the help of the dissipative particle dynamics simulation technique. Because hydrophobic matching is believed to be one of the main physical mechanisms regulating lipid-protein interactions in membranes, we considered proteins of different hydrophobic length (as well as different sizes). We studied the cooperative behavior of the lipid-protein system at mesoscopic time- and lengthscales. In particular, we correlated in a systematic way the protein-induced bilayer perturbation, and the lipid-induced protein tilt, with the hydrophobic mismatch (positive and negative) between the protein hydrophobic length and the pure lipid bilayer hydrophobic thickness. The protein-induced bilayer perturbation was quantified in terms of a coherence length, xi(P), of the lipid bilayer hydrophobic thickness profile around the protein. The dependence on temperature of xi(P), and the protein tilt-angle, were studied above the main-transition temperature of the pure system, i.e., in the fluid phase. We found that xi(P) depends on mismatch, i.e., the higher the mismatch is, the longer xi(P) becomes, at least for positive values of mismatch; a dependence on the protein size appears as well. In the case of large model proteins experiencing extreme mismatch conditions, in the region next to the so-called lipid annulus, there appears an undershooting (or overshooting) region where the bilayer hydrophobic thickness is locally lower (or higher) than in the unperturbed bilayer, depending on whether the protein hydrophobic length is longer (or shorter) than the pure lipid bilayer hydrophobic thickness. Proteins may tilt when embedded in a too-thin bilayer. Our simulation data suggest that, when the embedded protein has a small size, the main mechanism to compensate for a large hydrophobic mismatch is the tilt, whereas large proteins react to negative mismatch by causing an increase of the hydrophobic thickness of the nearby bilayer. Furthermore, for the case of small, peptidelike proteins, we found the same type of functional dependence of the protein tilt-angle on mismatch, as was recently detected by fluorescence spectroscopy measurements.     

Perturbation of the hydrophobic core of lipid bilayers by the human antimicrobial peptide LL-37

LL-37 is a cationic, amphipathic alpha-helical antimicrobial peptide found in humans that kills cells by disrupting the cell membrane. To disrupt membranes, antimicrobial peptides such as LL-37 must alter the hydrophobic core of the bilayer. Differential scanning calorimetry and deuterium ((2)H) NMR experiments on acyl chain perdeuterated lipids demonstrate that LL-37 inserts into the hydrophobic region of the bilayer and alters the chain packing and cooperativity. The results show that hydrophobic interactions between LL-37 and the hydrophobic acyl chains are as important for the ability of this peptide to disrupt lipid bilayers as its electrostatic interactions with the polar headgroups. The (2)H NMR data are consistent with the previously determined surface orientation of LL-37 (Henzler Wildman, K. A., et al. (2003) Biochemistry 42, 6545) with an estimated 5-6 A depth of penetration of the hydrophobic face of the amphipathic helix into the hydrophobic interior of the bilayer. LL-37 also alters the material properties of lipid bilayers, including the area per lipid, hydrophobic thickness, and coefficient of thermal expansion in a manner that varies with lipid type and temperature. Comparison of the effect of LL-37 on 1-palmitoyl-2-oleoyl-phosphatidylcholine (POPC-d(31)) and 1,2-dimyristoyl-phosphatidylcholine (DMPC-d(54)) at different temperatures demonstrates the importance of bilayer order in determining the type and extent of disordering and disruption of the hydrophobic core by LL-37. One possible explanation, which accounts for both the (2)H NMR data presented here and the known surface orientation of LL-37 under identical conditions, is that bilayer order influences the depth of insertion of LL-37 into the hydrophobic/hydrophilic interface of the bilayer, altering the balance of electrostatic and hydrophobic interactions between the peptide and the lipids.     

Combined influence of cholesterol and synthetic amphiphillic peptides upon bilayer thickness in model membranes.

Deuterium (2H) NMR was used to study bilayer hydrophobic thickness and mechanical properties when cholesterol and/or synthetic amphiphillic polypeptides were added to deuterated POPC lipid bilayer membranes in the liquid-crystalline (fluid) phase. Smoothed acyl chain orientational order profiles were used to calculate bilayer hydrophobic thickness. Addition of 30 mol% cholesterol to POPC at 25 degrees C increased the bilayer thickness from 2.58 to 2.99 nm. The peptides were chosen to span the bilayers with more or less mismatch between the hydrophobic peptide length and membrane hydrophobic thickness. The average thickness of the pure lipid bilayers was significantly perturbed upon addition of peptide only in cases of large mismatch, being increased (decreased) when the peptide hydrophobic length was greater (less) than that of the pure bilayer, consistent with the "mattress" model of protein lipid interactions (Mouritsen, O.G., and M. Bloom. 1984. Biophys. J. 46:141-153). The experimental results were also used to examine the combined influence of the polypeptides and cholesterol on the orientational order profile and thickness expansivity of the membranes. A detailed model for the spatial distribution of POPC and cholesterol molecules in the bilayers was proposed to reconcile the general features of these measurements with micromechanical measurements of area expansivity in closely related systems. Experiments to test the model were proposed.     

Interaction of a peptide model of a hydrophobic transmembrane alpha-helical segment of a membrane protein with phosphatidylcholine bilayers: differential scanning calorimetric and FTIR spectroscopic studies.

High-sensitivity differential scanning calorimetry (DSC) and Fourier transform infrared (FTIR) spectroscopy were used to study the interaction of a synthetic model hydrophobic peptide, Lys2-Gly-Leu24-Lys2-Ala-amide, and members of the homologous series of n-saturated diacylphosphatidylcholines. In the low range of peptide mole fractions, the DSC thermograms exhibited by the lipid/peptide mixtures are resolvable into two components. One of these components is fairly narrow, highly cooperative, and exhibits properties which are similar to but not identical with those of the pure lipid. In addition, the fractional contribution of this component to the total enthalpy change, the peak transition temperature, and cooperativity decrease with an increase in peptide concentration, more or less independently of acyl chain length. The other component is very broad and predominates in the high range of peptide concentration. These two components have been assigned to the chain-melting phase transitions of populations of bulk lipid and peptide-associated lipid, respectively. Moreover, when the mean hydrophobic thickness of the PC bilayer is less than the peptide hydrophobic length, the peptide-associated lipid melts at higher temperatures than does the bulk lipid and vice versa. In addition, the chain-melting enthalpy of the broad endotherm does not decrease to zero even at high peptide concentrations, suggesting that this peptide reduces but do not abolish the cooperative gel/liquid-crystalline phase transition of the lipids with which it is in contact. Our DSC results indicate that the width of the phase transition observed at high peptide concentration is inversely but discontinuously related to hydrocarbon chain length and that gel phase immiscibility occurs when the hydrophobic thickness of the bilayer greatly exceeds the hydrophobic length of the peptide. The FTIR spectroscopic data indicate that the peptide forms a very stable alpha-helix under all of our experimental conditions but that small distortions of its alpha-helical conformation are induced in response to any mismatch between peptide hydrophobic length and bilayer hydrophobic thickness. These results also indicate that the peptide alters the conformational disposition of the acyl chains in contact with it and that the resultant conformational changes in the lipid hydrocarbon chains tend to minimize the extent of mismatch of peptide hydrophobic length and bilayer hydrophobic thickness.     

Haemolysin-derived synthetic peptides with pore-forming and haemolytic activity

Escherichia coli haemolysin (Hlya) is a pore-forming protein which belongs to the family of 'Repeat-toxins' (RTX) (Lo et al., 1987; Lally et al., 1989; Kraig et al., 1990). A model for the pore-forming structure of HlyA has been proposed (Ludwig et al., 1991) which consists of eight transmembrane segments all present in this hydrophobic region of HlyA. We report here that two synthetic peptides of 10 and 8 amino acids in length (Pep1 and Pep2, respectively), which are derived from transmembrane segment V, are able to form pores in an artificial lipid bilayer. In addition, Pep1 exhibits strong haemolytic activity when tested on human red blood cells (HRBCs). The haemolytic activity of Pep1 and of E. coli haemolysin is completely inhibited by antibodies raised against Pep1.     

Effect of variations in the structure of a polyleucine-based alpha-helical transmembrane peptide on its interaction with phosphatidylcholine bilayers

High-sensitivity differential scanning calorimetry (DSC) and Fourier transform infrared (FTIR) spectroscopy were used to study the interaction of an alpha-helical transmembrane peptide, acetyl-Lys2-Leu24-Lys2-amide (L24), and odd-chain members of the homologous series of n-saturated diacylphosphatidylcholines. An analogue of L24, in which the lysine residues were all replaced by 2,3-diaminopropionic acid, and another, in which a leucine residue at each end of the polyLeu sequence was replaced by a tryptophan, were also studied. At low peptide concentrations, the DSC thermograms exhibited by these lipid/peptide mixtures are resolvable into two components. One of these components is fairly narrow, highly cooperative, and exhibits properties which are similar to but not identical with those of the pure lipid. In addition, the transition temperature and cooperativity of this component, and its fractional contribution to the total enthalpy change, decrease with an increase in peptide concentration, more or less independently of phospholipid acyl chain length. The other component is very broad and predominates at high peptide concentrations. These two components have been assigned to the chain-melting phase transitions of populations of peptide-poor and peptide-enriched lipid domains, respectively. Moreover, when the mean hydrophobic thickness of the PC bilayer is less than the peptide hydrophobic length, the peptide-associated lipid melts at higher temperatures than does the bulk lipid and vice versa. In addition, the chain-melting enthalpy of the broad endotherm does not decrease to zero even at high peptide concentrations, suggesting that these peptides reduce somewhat but do not abolish the cooperative gel/liquid-crystalline phase transition of the lipids with which it is in contact. Our DSC results indicate that the width of the broad phase transition observed at high peptide concentration is inversely but discontinuously related to hydrocarbon chain length. Our FTIR spectroscopic data indicate that these peptides form a very stable alpha-helix under all of our experimental conditions but that small distortions of their alpha-helical conformation are induced in response to mismatch between peptide hydrophobic length and gel-state bilayer hydrophobic thickness. We also present evidence that these distortions are localized to the N- and C-terminal regions of these peptides. Interestingly, replacing the terminal Lys residues of L24 by 2,3-diaminopropionic acid residues actually attenuates the hydrophobic mismatch effects of the peptide on the thermotropic phase behavior of the host PC bilayer, in contrast to the predictions of the snorkel hypothesis. We rationalize this attenuated hydrophobic mismatch effect by postulating that the 2,3-diaminopropionic acid residues are too short to engage in significant electrostatic and hydrogen-bonding interactions with the polar headgroups of the host phospholipid bilayer, even in the absence of any hydrophobic mismatch between incorporated peptide and the bilayer. Similarly, the reduced hydrophobic mismatch effect also observed when the two terminal Leu residues of L24 are replaced by Trp residues is rationalized by considering the lower energetic cost of exposing the Trp as opposed to the Leu residues to the aqueous phase in thin PC bilayers and the higher cost of inserting the Trp as opposed to the Leu residues into the hydrophobic cores of thick PC bilayers.     

Influence of transmembrane peptides on bilayers of phosphatidylcholines with different acyl chain lengths studied by solid-state NMR

The molecular orientation in a lipid membrane of the peptide fragment VEYAGIALFFVAAVLTLWSMLQYLSAAR (phosphatidylglycerophosphate synthase (Pgs) peptide E) of an integral membrane protein, Pgs, in Escherichia coli has been investigated by solid-state 15N nuclear magnetic resonance (NMR) on macroscopically aligned lipid bilayers. The secondary structure of the peptide in lipid vesicles was determined by circular dichroism spectroscopy. Furthermore, the phase behaviour of the Pgs peptide E/dierucoylphosphatidylcholine (DEruPC)/water system was determined by (2)H, (31)P and 15N solid-state NMR spectroscopy. The phase behaviour obtained was then compared to that of the Pgs peptide E solubilised in dioleoylphosphatidylcholine and water that was previously studied by Morein et al. [Biophys. J. 73 (1997) 3078-3088]. This was aimed to answer the question whether a difference in the length of the hydrophobic part of this peptide and the hydrophobic thickness of the lipid bilayer (hydrophobic mismatch) will affect the phase behaviour. The peptide mostly has a transmembrane orientation and is in an alpha-helical conformation. An isotropic phase is formed in DEruPC with high peptide content (peptide/lipid molar ratio (p/l) > or =1:15) and high water content (> or =50%, w/w) at 35 degrees C. At 55 and 65 degrees C an isotropic phase is induced at high water content (> or =50%, w/w) at all peptide contents studied (no isotropic phase forms in the lipid/water system under the conditions in this study). At high peptide contents (p/l> or =1:15) an isotropic phase forms at 20 and 40% (w/w) of water at 55 and 65 degrees C. A comparison of the phase behaviour of the two homologous lipid systems reveals striking similarities, although the thicknesses of the two lipid bilayers differ by 7 A. This suggests that the rationalisation of the phase behaviour in terms of the hydrophobic mismatch is not applicable to these systems. The C-terminus of Pgs peptide E is amphiphilic and a considerable part of the peptide is situated outside the hydrophobic part of the bilayer, a property of the peptide that to a large extent will affect the lipid/peptide phase behaviour.     

Solid-state NMR investigation of the depth of insertion of protegrin-1 in lipid bilayers using paramagnetic Mn2+

The depth of insertion of an antimicrobial peptide, protegrin-1 (PG-1), in lipid bilayers is investigated using solid-state NMR. Paramagnetic Mn(2+) ions bind to the surface of lipid bilayers and induce distance-dependent dipolar relaxation of nuclear spins. By comparing the signal dephasing of the peptide with that of the lipids, whose segmental depths of insertion are known, we determined the depths of several residues of PG-1 in 1,2 dilauryl-sn-glycero-3-phosphotidylcholine (DLPC) bilayers. We found that residues G2 at the N-terminus and F12 at the beta-turn of the peptide reside near the membrane surface, whereas L5 and V16 are embedded in the acyl chain region. The depths increase in the order of G2 < F12 < L5 < V16. These intensity-dephasing results are confirmed by direct measurement of the paramagnetically enhanced (13)C transverse relaxation rates. The relative depths indicate that PG-1 is tilted from the bilayer normal, which is consistent with independent solid-state NMR measurements of PG-1 orientation in the same lipids (Yamaguchi et al., 2001). They also indicate that PG-1 is fully immersed in the lipid bilayer. However, a quantitative mismatch between the bilayer thickness and PG-1 length suggests a local thinning of the DLPC bilayer by 8-10 A. The depth sensitivity of this Mn(2+) dephasing technique is tunable with the Mn(2+) concentration to focus on different regions of the lipid bilayer.     

Transmembrane helices of membrane proteins may flex to satisfy hydrophobic mismatch

A novel mechanism for membrane modulation of transmembrane protein structure, and consequently function, is suggested in which mismatch between the hydrophobic surface of the protein and the hydrophobic interior of the lipid bilayer induces a flexing or bending of a transmembrane segment of the protein. Studies on model hydrophobic transmembrane peptides predict that helices tilt to submerge the hydrophobic surface within the lipid bilayer to satisfy the hydrophobic effect if the helix length exceeds the bilayer width. The hydrophobic surface of transmembrane helix 1 (TM1) of lactose permease, LacY, is accessible to the bilayer, and too long to be accommodated in the hydrophobic portion of a typical lipid bilayer if oriented perpendicular to the membrane surface. Hence, nuclear magnetic resonance (NMR) data and molecular dynamics simulations show that TM1 from LacY may flex as well as tilt to satisfy the hydrophobic mismatch with the bilayer. In an analogous study of the hydrophobic mismatch of TM7 of bovine rhodopsin, similar flexing of the transmembrane segment near the conserved NPxxY sequence is observed. As a control, NMR data on TM5 of lacY, which is much shorter than TM1, show that TM5 is likely to tilt, but not flex, consistent with the close match between the extent of hydrophobic surface of the peptide and the hydrophobic thickness of the bilayer. These data suggest mechanisms by which the lipid bilayer in which the protein is embedded modulates conformation, and thus function, of integral membrane proteins through interactions with the hydrophobic transmembrane helices.     

Transmembrane helices of membrane proteins may flex to satisfy hydrophobic mismatch

A novel mechanism for membrane modulation of transmembrane protein structure, and consequently function, is suggested in which mismatch between the hydrophobic surface of the protein and the hydrophobic interior of the lipid bilayer induces a flexing or bending of a transmembrane segment of the protein. Studies on model hydrophobic transmembrane peptides predict that helices tilt to submerge the hydrophobic surface within the lipid bilayer to satisfy the hydrophobic effect if the helix length exceeds the bilayer width. The hydrophobic surface of transmembrane helix 1 (TM1) of lactose permease, LacY, is accessible to the bilayer, and too long to be accommodated in the hydrophobic portion of a typical lipid bilayer if oriented perpendicular to the membrane surface. Hence, nuclear magnetic resonance (NMR) data and molecular dynamics simulations show that TM1 from LacY may flex as well as tilt to satisfy the hydrophobic mismatch with the bilayer. In an analogous study of the hydrophobic mismatch of TM7 of bovine rhodopsin, similar flexing of the transmembrane segment near the conserved NPxxY sequence is observed. As a control, NMR data on TM5 of lacY, which is much shorter than TM1, show that TM5 is likely to tilt, but not flex, consistent with the close match between the extent of hydrophobic surface of the peptide and the hydrophobic thickness of the bilayer. These data suggest mechanisms by which the lipid bilayer in which the protein is embedded modulates conformation, and thus function, of integral membrane proteins through interactions with the hydrophobic transmembrane helices.     

Stability of KcsA tetramer depends on membrane lateral pressure

The potassium channel KcsA forms an extremely stable tetramer. Despite this high stability, it has been shown that the membrane-mimicking solvent 2,2,2-trifluoroethanol (TFE) can induce tetramer dissociation [Valiyaveetil, F. I., et al. (2002) Biochemistry 41, 10771-7, and Demmers, J. A. A., et al. (2003) FEBS Lett. 541, 69-77]. Here we have studied the effect of TFE on the structure and oligomeric state of the KcsA tetramer, reconstituted in different lipid systems. It was found that TFE changes the secondary and tertiary structure of KcsA and that it can dissociate the KcsA tetramer in all systems used. The tetramer is stabilized by a lipid bilayer as compared to detergent micelles. The extent of stabilization was found to depend on the nature of the lipids: a strong stabilizing effect of the nonbilayer lipid phosphatidylethanolamine (PE) was observed, but no effect of the charged phoshosphatidylglycerol (PG) as compared to phosphatidylcholine (PC) was found. To understand how lipids stabilize KcsA against TFE-induced tetramer dissociation, we also studied the effect of TFE on the bilayer organization in the various lipid systems, using (31)P and (2)H NMR. The observed lipid dependency was similar as was found for tetramer stabilization: PE increased the bilayer stability as compared to PC, while PG behaved similar to PC. Furthermore, it was found that TFE has a large effect on the acyl chain ordering. The results indicate that TFE inserts primarily in the membrane interface. We suggest that the lipid bilayer stabilizes the KcsA tetramer by the lateral pressure in the acyl chain region and that this stabilizing effect increases when a nonbilayer lipid like PE is present.     

Tilt angles of transmembrane model peptides in oriented and non-oriented lipid bilayers as determined by 2H solid-state NMR

Solid-state NMR methods employing (2)H NMR and geometric analysis of labeled alanines (GALA) were used to study the structure and orientation of the transmembrane alpha-helical peptide acetyl-GWW(LA)(8)LWWA-amide (WALP23) in phosphatidylcholine (PC) bilayers of varying thickness. In all lipids the peptide was found to adopt a transmembrane alpha-helical conformation. A small tilt angle of 4.5 degrees was observed in di-18:1-PC, which has a hydrophobic bilayer thickness that approximately matches the hydrophobic length of the peptide. This tilt angle increased slightly but systematically with increasing positive mismatch to 8.2 degrees in di-C12:0-PC, the shortest lipid used. This small increase in tilt angle is insufficient to significantly change the effective hydrophobic length of the peptide and thereby to compensate for the increasing hydrophobic mismatch, suggesting that tilt of these peptides in a lipid bilayer is energetically unfavorable. The tilt and also the orientation around the peptide axis were found to be very similar to the values previously reported for a shorter WALP19 peptide (GWW(LA)(6)LWWA). As also observed in this previous study, the peptide rotates rapidly around the bilayer normal, but not around its helix axis. Here we show that these properties allow application of the GALA method not only to macroscopically aligned samples but also to randomly oriented samples, which has important practical advantages. A minimum of four labeled alanine residues in the hydrophobic transmembrane sequence was found to be required to obtain accurate tilt values using the GALA method.     

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.     

Structure, topology, and tilt of cell-signaling peptides containing nuclear localization sequences in membrane bilayers determined by solid-state NMR and molecular dynamics simulation studies

Cell-signaling peptides have been extensively used to transport functional molecules across the plasma membrane into living cells. These peptides consist of a hydrophobic sequence and a cationic nuclear localization sequence (NLS). It has been assumed that the hydrophobic region penetrates the hydrophobic lipid bilayer and delivers the NLS inside the cell. To better understand the transport mechanism of these peptides, in this study, we investigated the structure, orientation, tilt of the peptide relative to the bilayer normal, and the membrane interaction of two cell-signaling peptides, SA and SKP. Results from CD and solid-state NMR experiments combined with molecular dynamics simulations suggest that the hydrophobic region is helical and has a transmembrane orientation with the helical axis tilted away from the bilayer normal. The influence of the hydrophobic mismatch, between the hydrophobic length of the peptide and the hydrophobic thickness of the bilayer, on the tilt angle of the peptides was investigated using thicker POPC and thinner DMPC bilayers. NMR experiments showed that the hydrophobic domain of each peptide has a tilt angle of 15 +/- 3 degrees in POPC, whereas in DMPC, 25 +/- 3 degree and 30 +/- 3 degree tilts were observed for SA and SKP peptides, respectively. These results are in good agreement with molecular dynamics simulations, which predict a tilt angle of 13.3 degrees (SA in POPC), 16.4 degrees (SKP in POPC), 22.3 degrees (SA in DMPC), and 31.7 degrees (SKP in DMPC). These results and simulations on the hydrophobic fragment of SA or SKP suggest that the tilt of helices increases with a decrease in bilayer thickness without changing the phase, order, and structure of the lipid bilayers.     

Peptide Nanopores and Lipid Bilayers: Interactions by Coarse-Grained Molecular-Dynamics Simulations

A set of 49 protein nanopore-lipid bilayer systems was explored by means of coarse-grained molecular-dynamics simulations to study the interactions between nanopores and the lipid bilayers in which they are embedded. The seven nanopore species investigated represent the two main structural classes of membrane proteins (α-helical and β-barrel), and the seven different bilayer systems range in thickness from ∼28 to ∼43 Å. The study focuses on the local effects of hydrophobic mismatch between the nanopore and the lipid bilayer. The effects of nanopore insertion on lipid bilayer thickness, the dependence between hydrophobic thickness and the observed nanopore tilt angle, and the local distribution of lipid types around a nanopore in mixed-lipid bilayers are all analyzed. Different behavior for nanopores of similar hydrophobic length but different geometry is observed. The local lipid bilayer perturbation caused by the inserted nanopores suggests possible mechanisms for both lipid bilayer-induced protein sorting and protein-induced lipid sorting. A correlation between smaller lipid bilayer thickness (larger hydrophobic mismatch) and larger nanopore tilt angle is observed and, in the case of larger hydrophobic mismatches, the simulated tilt angle distribution seems to broaden. Furthermore, both nanopore size and key residue types (e.g., tryptophan) seem to influence the level of protein tilt, emphasizing the reciprocal nature of nanopore-lipid bilayer interactions.     

Solid-state NMR investigations of peptide-lipid interaction and orientation of a beta-sheet antimicrobial peptide,protegrin

Protegrin-1 (PG-1) is a broad-spectrum beta-sheet antimicrobial peptide found in porcine leukocytes. The mechanism of action and the orientation of PG-1 in lipid bilayers are here investigated using (2)H, (31)P, (13)C, and (15)N solid-state NMR spectroscopy. (2)H spectra of mechanically aligned and chain-perdeuterated 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphatidylcholine (POPC) bilayers indicate that PG-1 at high concentrations destroys the orientational order of the aligned lamellar bilayer. The conformation of the lipid headgroups in the unoriented region is significantly altered, as seen from the (31)P spectra of POPC and the (2)H spectra of headgroup-deuterated 1,2-dipalmitoyl-sn-glycero-3-phosphatidylcholine. These observations indicate that PG-1 disrupts microbial membranes by breaking the extended bilayer into smaller disks, where a significant fraction of lipids is located in the edges of the disks with a distribution of orientations. These edges allow the lipid bilayer to bend back on itself as in toroidal pores. Interestingly, this loss of bilayer orientation occurs only in long-chain lipids such as POPC and not in shorter chain lipids such as 1,2-dilauroyl-sn-glycero-3-phosphatidylcholine (DLPC). To understand the mode of binding of PG-1 to the lipid bilayer, we determined the orientation of PG-1 in DLPC bilayers. The (13)CO and (15)N chemical shifts of Val-16 labeled PG-1 indicate that the beta-strand axis is tilted by 55 degrees +/- 5 degrees from the bilayer normal while the normal of the beta-sheet plane is 48 degrees +/- 5 degrees from the bilayer normal. This orientation favors interaction of the hydrophobic backbone of the peptide with the hydrophobic core of the bilayer and positions the cationic Arg side chains to interact with the anionic phosphate groups. This is the first time that the orientation of a disulfide-stabilized beta-sheet membrane peptide has been determined by solid-state NMR.     

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.     

Energetics of inclusion-induced bilayer deformations.

The material properties of lipid bilayers can affect membrane protein function whenever conformational changes in the membrane-spanning proteins perturb the structure of the surrounding bilayer. This coupling between the protein and the bilayer arises from hydrophobic interactions between the protein and the bilayer. We analyze the free energy cost associated with a hydrophobic mismatch, i.e., a difference between the length of the protein's hydrophobic exterior surface and the average thickness of the bilayer's hydrophobic core, using a (liquid-crystal) elastic model of bilayer deformations. The free energy of the deformation is described as the sum of three contributions: compression-expansion, splay-distortion, and surface tension. When evaluating the interdependence among the energy components, one modulus renormalizes the other: e.g., a change in the compression-expansion modulus affects not only the compression-expansion energy but also the splay-distortion energy. The surface tension contribution always is negligible in thin solvent-free bilayers. When evaluating the energy per unit distance (away from the inclusion), the splay-distortion component dominates close to the bilayer/inclusion boundary, whereas the compression-expansion component is more prominent further away from the boundary. Despite this complexity, the bilayer deformation energy in many cases can be described by a linear spring formalism. The results show that, for a protein embedded in a membrane with an initial hydrophobic mismatch of only 1 A, an increase in hydrophobic mismatch to 1.3 A can increase the Boltzmann factor (the equilibrium distribution for protein conformation) 10-fold due to the elastic properties of the bilayer.     

Structure and dynamics of the lipid modifications of a transmembrane α-helical peptide determined by ²H solid-state NMR spectroscopy

The fusion of biological membranes is mediated by integral membrane proteins with α-helical transmembrane segments. Additionally, those proteins are often modified by the covalent attachment of hydrocarbon chains. Previously, a series of de novo designed α-helical peptides with mixed Leu/Val sequences was presented, mimicking fusiogenically active transmembrane segments in model membranes (Hofmann et al., Proc. Natl. Acad. Sci. USA 101 (2004) 14776-14781). From this series, we have investigated the peptide LV16 (KKKW LVLV LVLV LVLV LVLV KKK), which was synthesized featuring either a free N-terminus or a saturated N-acylation of 2, 8, 12, or 16 carbons. We used 2H and 31P NMR spectroscopy to investigate the structure and dynamics of those peptide lipid modifications in POPC and DLPC bilayers and compared them to the hydrocarbon chains of the surrounding membrane. Except for the C2 chain, all peptide acyl chains were found to insert well into the membrane. This can be explained by the high local lipid concentrations the N-terminal lipid chains experience. Further, the insertion of these peptides did not influence the membrane structure and dynamics as seen from the 2H and 31P NMR data. In spite of the fact that the longer acyl chains insert into the membrane, they do not adapt their lengths to the thickness of the bilayer. Even the C16 lipid chain on the peptide, which could match the length of the POPC palmitoyl chain, exhibited lower order parameters in the upper chain, which get closer and finally reach similar values in the lower chain region. 2H NMR square law plots reveal motions of slightly larger amplitudes for the peptide lipid chains compared to the surrounding phospholipids. In spite of the significantly different chain lengths of the acylations, the fraction of gauche defects in the inserted chains is constant.     

Theoretical analysis of hydrophobic matching and membrane-mediated interactions in lipid bilayers containing gramicidin.

We present a quantitative analysis of the effects of hydrophobic matching and membrane-mediated protein-protein interactions exhibited by gramicidin embedded in dimyristoylphosphatidylcholine (DMPC) and dilauroylphosphatidylcholine (DLPC) bilayers (Harroun et al., 1999. Biophys. J. 76:937-945). Incorporating gramicidin, at 1:10 peptide/lipid molar ratio, decreases the phosphate-to-phosphate (PtP) peak separation in the DMPC bilayer from 35.3 A without gramicidin to 32.7 A. In contrast, the same molar ratio of gramicidin in DLPC increases the PtP from 30.8 A to 32.1 A. Concurrently, x-ray in-plane scattering showed that the most probable nearest-neighbor separation between gramicidin channels was 26.8 A in DLPC, but reduced to 23.3 A in DMPC. In this paper we review the idea of hydrophobic matching in which the lipid bilayer deforms to match the hydrophobic surface of the embedded proteins. We use a simple elasticity theory, including thickness compression, tension, and splay terms to describe the membrane deformation. The energy of membrane deformation is compared with the energy cost of hydrophobic mismatch. We discuss the boundary conditions between a gramicidin channel and the lipid bilayer. We used a numerical method to solve the problem of membrane deformation profile in the presence of a high density of gramicidin channels and ran computer simulations of 81 gramicidin channels to find the equilibrium distributions of the channels in the plane of the bilayer. The simulations contain four parameters: bilayer thickness compressibility 1/B, bilayer bending rigidity Kc, the channel-bilayer mismatch Do, and the slope of the interface at the lipid-protein boundary s. B, Kc, and Do were experimentally measured; the only free parameter is s. The value of s is determined by the requirement that the theory produces the experimental values of bilayer thinning by gramicidin and the shift in the peak position of the in-plane scattering due to membrane-mediated channel-channel interactions. We show that both hydrophobic matching and membrane-mediated interactions can be understood by the simple elasticity theory.