Cumulative effects of amino acid substitutions and hydrophobic mismatch upon the transmembrane stability and conformation of hydrophobic alpha-helices

The effects of amino acid substitutions upon the behavior of poly(Leu)-rich alpha-helices inserted into model membrane vesicles were investigated. One or two consecutive Leu residues in the hydrophobic core of the helix were substituted with A, F, G, S, D, K, H, P, GG, SS, PG, PP, KK, or DD residues. A Trp placed at the center of the sequence allowed assessment of peptide behavior via fluorescence emission lambda(max) and dual quenching analysis of Trp depth [Caputo, G. A., and London, E. (2003) Biochemistry 42, 3265-3274]. In vesicles composed of dioleoylphosphatidylcholine (DOPC), all of the peptides with single substitutions adopted a transmembrane (TM) state. Experiments were also performed in thicker bilayers composed of dierucoylphosphatidylcholine (DEuPC). In DEuPC vesicles TM states were destabilized by mismatch between helix length and bilayer thickness. Nevertheless, in DEuPC vesicles TM states were still prevalent for peptides with single substitutions, although less so for peptides with P, K, H, or D substitutions. In contrast to single substitutions, certain consecutive double substitutions strongly interfered with formation of TM states. In both DOPC and DEuPC vesicles DD and KK substitutions abolished the normal TM state, but GG and SS substitutions had little effect. In even wider bilayers, a SS substitution reduced the formation of a TM state. A peptide with a PP substitution maintained the TM state in DOPC vesicles, but in DEuPC vesicles the level of formation of the TM state was significantly reduced. Upon disruption of normal TM insertion peptides moved close to the bilayer surface, with the exception of the KK-substituted peptide in DOPC vesicles, which formed a truncated TM segment. These studies begin to provide a detailed relationship between sequence and the stability of TM insertion and show that the influence of insertion-destabilizing residues upon hydrophobic helices can be strongly modulated by properties such as mismatch. For certain helix-forming hydrophobic sequences, sensitivity to lipid structure may be sufficient to induce large conformational changes in vivo.     

Position and ionization state of Asp in the core of membrane-inserted alpha helices control both the equilibrium between transmembrane and nontransmembrane helix topography and transmembrane helix positioning

The behavior of model-membrane-inserted polyLeu-rich peptides containing Asp residues located at various positions in their hydrophobic core was investigated. The topography of the bilayer-inserted alpha helices formed by these peptides was evaluated by measuring the emission lambda(max) and quenching the fluorescence of a Trp at the center of the peptide sequence. When Asp residues were protonated (at low pH), peptides that were incorporated into vesicles composed of dioleoylphosphatidylcholine (DOPC) adopted a topography in which the polyLeu sequence predominantly formed a normal transmembrane (TM) helix. When Asp residues were ionized (at neutral or high pH), topography was altered in a manner that would allow the charged Asp residues to reside near the bilayer surface. In DOPC vesicles, most peptides repositioned so that the longest segment of consecutive hydrophobic residues (12 residue minimum) formed a truncated/shifted TM structure. However, peptides with one or two charged Asp residues close to the center of the hydrophobic sequence and thus lacking even a 12-residue continuous hydrophobic segment, formed a helical non-TM state locating near the bilayer surface. At low pH, incorporation of the peptides into thicker bilayers composed of dierucoylphosphatidylcholine (DEuPC) resulted in the formation of a mixture of the normal TM state and the non-TM helical state located near the bilayer surface. In DEuPC vesicles at high pH, the non-TM state tended to predominate. How Asp-ionization-dependent shifts in helix topography may regulate the function of membrane proteins exposed to environments with differing pH in vivo (e.g., endosomes) is discussed.     

Revisiting Hydrophobic Mismatch with Free Energy Simulation Studies of Transmembrane Helix Tilt and Rotation

Protein-lipid interaction and bilayer regulation of membrane protein functions are largely controlled by the hydrophobic match between the transmembrane (TM) domain of membrane proteins and the surrounding lipid bilayer. To systematically characterize responses of a TM helix and lipid adaptations to a hydrophobic mismatch, we have performed a total of 5.8-μs umbrella sampling simulations and calculated the potentials of mean force (PMFs) as a function of TM helix tilt angle under various mismatch conditions. Single-pass TM peptides called WALPn (n = 16, 19, 23, and 27) were used in two lipid bilayers with different hydrophobic thicknesses to consider hydrophobic mismatch caused by either the TM length or the bilayer thickness. In addition, different flanking residues, such as alanine, lysine, and arginine, instead of tryptophan in WALP23 were used to examine their influence. The PMFs, their decomposition, and trajectory analysis demonstrate that 1), tilting of a single-pass TM helix is the major response to a hydrophobic mismatch; 2), TM helix tilting up to ∼10° is inherent due to the intrinsic entropic contribution arising from helix precession around the membrane normal even under a negative mismatch; 3), the favorable helix-lipid interaction provides additional driving forces for TM helix tilting under a positive mismatch; 4), the minimum-PMF tilt angle is generally located where there is the hydrophobic match and little lipid perturbation; 5), TM helix rotation is dependent on the specific helix-lipid interaction; and 6), anchoring residues at the hydrophilic/hydrophobic interface can be an important determinant of TM helix orientation.     

Effect of lipid composition on the topography of membrane-associated hydrophobic helices: stabilization of transmembrane topography by anionic lipids

To investigate the effect of lipid structure upon the membrane topography of hydrophobic helices, the behavior of hydrophobic peptides was studied in model membrane vesicles. To define topography, fluorescence and fluorescence quenching methods were used to determine the location of a Trp at the center of the hydrophobic sequence. For peptides with cationic residues flanking the hydrophobic sequence, the stability of the transmembrane (TM) configuration (relative to a membrane-bound non-TM state) increased as a function of lipid composition on the order: 1:1 (mol:mol) 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC):1-palmitoyl-2-oleoyl-sn-glycero-3-phosphoethanolamine ∼ 6:4 POPC:cholesterol < POPC ∼ dioleoylphosphatidylcholine (DOPC) < 1,2-dioleoyl-sn-glycero-3-[phospho-rac-(1-glycerol)] sodium salt (DOPG) ≤ 1,2-dioleoyl-sn-glycero-3-[phospho-l-serine] sodium salt (DOPS), indicating that the anionic lipids DOPG and DOPS most strongly stabilized the TM configuration. TM stabilization was near maximal at 20-30 mol% anionic lipid, which are physiologically relevant values. TM stabilization by anionic lipid was observed for hydrophobic sequences with a diverse set of sequences (including polyAla), diverse lengths (from 12 to 22 residues), and various cationic flanking residues (H, R, or K), but not when the flanking residues were uncharged. TM stabilization by anionic lipid was also dependent on the number of cationic residues flanking the hydrophobic sequence, but was still significant with only one cationic residue flanking each end of the peptide. These observations are consistent with TM-stabilizing effects being electrostatic in origin. However, Trp located more deeply in DOPS vesicles relative to DOPG vesicles, and peptides in DOPS vesicles showed increased helix formation relative to DOPG and all other lipid compositions. These observations fit a model in which DOPS anchors flanking residues near the membrane surface more strongly than does DOPG and/or increases the stability of the TM state to a greater degree than DOPG. We conclude that anionic lipids can have significant and headgroup structure-specific effects upon membrane protein topography.     

Tilt angle of a trans-membrane helix is determined by hydrophobic mismatch

In order to investigate the compensation mechanism of a trans-membrane helix in response to hydrophobic mismatch, the tilt and rotation angles of the trans-membrane helix of Vpu aligned in lipid bilayers of various thickness were determined using orientation-dependent frequencies obtained from solid-state NMR experiments of aligned samples. A tilt angle of 18 degrees was observed in 18:1-O-PC/DOPG (9:1) lipid bilayers, which have a hydrophobic thickness that approximately matches the hydrophobic length of the trans-membrane helix of Vpu. Upon decreasing the hydrophobic thickness of lipid bilayers, no significant change in rotation angle was observed. However, the tilt angle increased systematically with increasing positive mismatch to 27 degrees in 14:0-O-PC/DMPG (9:1), 35 degrees in 12:0-O-PC/DLPG (9:1), and 51 degrees in 10:0 PC/10:0 PG (9:1) lipid bilayers, indicating that the change in tilt angle of the trans-membrane helix is a principal compensation mechanism for hydrophobic mismatch. In addition, the distinctive kink in the middle of the helix observed in 18:1 bilayers disappears in thinner bilayers. Although the opposite of what might be expected, this finding suggests that a helix kink may also be a part of the hydrophobic matching mechanism for trans-membrane 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.     

Hydrophobic Mismatch Drives the Interaction of E5 with the Transmembrane Segment of PDGF Receptor

The oncogenic E5 protein from bovine papillomavirus is a short (44 amino acids long) integral membrane protein that forms homodimers. It activates platelet-derived growth factor receptor (PDGFR) β in a ligand-independent manner by transmembrane helix-helix interactions. The nature of this recognition event remains elusive, as numerous mutations are tolerated in the E5 transmembrane segment, with the exception of one hydrogen-bonding residue. Here, we examined the conformation, stability, and alignment of the E5 protein in fluid lipid membranes of substantially varying bilayer thickness, in both the absence and presence of the PDGFR transmembrane segment. Quantitative synchrotron radiation circular dichroism analysis revealed a very long transmembrane helix for E5 of ∼26 amino acids. Oriented circular dichroism and solid-state ¹⁵N-NMR showed that the alignment and stability of this unusually long segment depend critically on the membrane thickness. When reconstituted alone in exceptionally thick DNPC lipid bilayers, the E5 helix was found to be inserted almost upright. In moderately thick bilayers (DErPC and DEiPC), it started to tilt and became slightly deformed, and finally it became aggregated in conventional DOPC, POPC, and DMPC membranes due to hydrophobic mismatch. On the other hand, when E5 was co-reconstituted with the transmembrane segment of PDGFR, it was able to tolerate even the most pronounced mismatch and was stabilized by binding to the receptor, which has the same hydrophobic length. As E5 is known to activate PDGFR within the thin membranes of the Golgi compartment, we suggest that the intrinsic hydrophobic mismatch of these two interaction partners drives them together. They seem to recognize each other by forming a closely packed bundle of mutually aligned transmembrane helices, which is further stabilized by a specific pair of hydrogen-bonding residues.     

Polar/Ionizable Residues in Transmembrane Segments: Effects on Helix-Helix Packing

The vast majority of membrane proteins are anchored to biological membranes through hydrophobic α-helices. Sequence analysis of high-resolution membrane protein structures show that ionizable amino acid residues are present in transmembrane (TM) helices, often with a functional and/or structural role. Here, using as scaffold the hydrophobic TM domain of the model membrane protein glycophorin A (GpA), we address the consequences of replacing specific residues by ionizable amino acids on TM helix insertion and packing, both in detergent micelles and in biological membranes. Our findings demonstrate that ionizable residues are stably inserted in hydrophobic environments, and tolerated in the dimerization process when oriented toward the lipid face, emphasizing the complexity of protein-lipid interactions in biological membranes.     

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.     

Hydrophobic organization of α-helix membrane bundle in bacteriorhodopsin

The hydrophobic organization of the intramembraneα-helical bundle in bacteriorhodopsin (BRh) was assessed based on a new approach to characterization of spatial hydrophobic properties of transmembrane (TM)α-helical peptides. The method employs two independent techniques: Monte Carlo simulations of nonpolar solvent around TM peptides and analysis of molecular hydrophobicity potential on their surfaces. The results obtained by the two methods agree with each other and permit precise hydrophobicity mapping of TM peptides. Superimposition of such data on the experimentally derived spatial model of the membrane moiety together with 2D maps of hydrophobic hydrophilic contacts provide considerable insight into the hydrophobic organization of BRh. The helix bundle is stabilized to a large extent by hydrophobic interactions between helices—neighbors in the sequence of BRh, by long-range interactions in helix pairs C-E, C-F, and C-G, and by nonpolar contracts between retinal and helices C, D, E, F. Unlike globular proteins, no polar contacts between residues distantly separated in the sequence of BRh were found in the bundle. One of the most striking results of this study is the finding that the hydrophobic organization of BRh is significantly different from those in bacterial photoreaction centers. Thus, TMα-helices in BRh expose their most nonpolar sides to the bilayer as well as to the neighboring helices and to the interior of the bundle. Some of them contact lipids with their relatively hydrophilic surfaces. No correlation was found between disposition of the most hydrophobic and the most variable sides of the TM helices.     

Control of the transmembrane orientation and interhelical interactions within membranes by hydrophobic helix length

We examined the effect of the length of the hydrophobic core of Lys-flanked poly(Leu) peptides on their behavior when inserted into model membranes. Peptide structure and membrane location were assessed by the fluorescence emission lambdamax of a Trp residue in the center of the peptide sequence, the quenching of Trp fluorescence by nitroxide-labeled lipids (parallax analysis), and circular dichroism. Peptides in which the hydrophobic core varied in length from 11 to 23 residues were found to be largely alpha-helical when inserted into the bilayer. In dioleoylphosphatidylcholine (diC18:1PC) bilayers, a peptide with a 19-residue hydrophobic core exhibited highly blue-shifted fluorescence, an indication of Trp location in a nonpolar environment, and quenching localized the Trp to the bilayer center, an indication of transmembrane structure. A peptide with an 11-residue hydrophobic core exhibited emission that was red-shifted, suggesting a more polar Trp environment, and quenching showed the Trp was significantly displaced from the bilayer center, indicating that this peptide formed a nontransmembranous structure. A peptide with a 23-residue hydrophobic core gave somewhat red-shifted fluorescence, but quenching demonstrated the Trp was still close to the bilayer center, consistent with a transmembrane structure. Analogous behavior was observed when the behavior of individual peptides was examined in model membranes with various bilayer widths. Other experiments demonstrated that in diC18:1PC bilayers the dilution of the membrane concentration of the peptide with a 23-residue hydrophobic core resulted in a blue shift of fluorescence, suggesting the red-shifted fluorescence at higher peptide concentrations was due to helix oligomerization. The intermolecular self-quenching of rhodamine observed when the peptide was rhodamine-labeled, and the concentration dependence of self-quenching, supported this conclusion. These studies indicate that the mismatch between helix length and bilayer width can control membrane location, orientation, and helix-helix interactions, and thus may mismatch control both membrane protein folding and the interactions between membrane proteins.     

Highly restricted distributions of hydrophobic and charged amino acids in longitudinal quadrants of alpha-helices

Helix formation in folding proteins is stabilized by binding of recurrent hydrophobic side chains in one longitudinal quadrant against the locally most hydrophobic region of the protein. To test this hypothesis, we fitted sequences of 247 alpha-helices of 55 proteins to the circular (infinite) template (symbol; see text) to maximize the strip-of-helix hydrophobicity index (the mean hydrophobicity of residues in (symbol; see text) positions). These template-predicted configurations closely matched crystallographic structures in 87% of four- or five-turn helices compared. We determined the longitudinal quadrant distributions of amino acids in the template-fitted, sheet projections of alpha-helices with respect to the best longitudinal, hydrophobic strip on each helix and to the N and C termini, interiors, and entire helices. Amino acids Leu, Ile, Val, and Phe were concentrated in one longitudinal quadrant (p less than 0.001). Lys, Arg, Asp, and Glu were not in the quadrant of Leu, Ile, Val, and Phe (p less than 0.001). Significant quadrant distributions for other amino acids and for termini of the helices were also found.     

Local hydrophobicity stabilizes secondary structures in proteins

The probability of occurrence of helix and β-sheet residues in 47 globular proteins was determined as a function of local hydrophobicity, which was defined by the sum of the Nozaki-Tanford transfer free energies at two nearest-neighbors on both sides of the amino acid sequence. In general, hydrophilic amino acids favor neither helix nor β-sheet formations when neighbor residues are also hydrophilic but favor helix formation at higher local hydrophobicity. On the other hand, some hydrophobic amino acids such as Met, Leu, and Ile favor helix formation when neighbor residues are hydrophilic. None of the hydrophobic amino acids favor β-sheet formation with hydrophilic neighbors, but most of them strongly favor β-sheet formation at high local hydrophobicity. When the average of 20 amino acids is taken, both helix and β-sheet residue probabilities are higher at higher local hydrophobicity, although the increase is steeper for β-sheets. Therefore, β-sheet formation is more influenced by local hydrophobicity than helix formation. Generally, helices are nearer the surface and tend to have hydrophilic and hydrophobic faces at opposite sides. The tendency of alternating regions of hydrophilic and hydrophobic residues in a helical sequence was revealed by calculating the correlation of the Nozaki-Tanford values. Such amphipathic helices may be important in protein–protein and protein–lipid interactions and in forming hydrophilic channels in the membrane. The choice of 30 nonhomologous proteins as the data set did not alter the above results.     

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.     

Studies of the minimum hydrophobicity of alpha-helical peptides required to maintain a stable transmembrane association with phospholipid bilayer membranes

The effects of the hydrophobicity and the distribution of hydrophobic residues on the surfaces of some designed alpha-helical transmembrane peptides (acetyl-K2-L(m)-A(n)-K2-amide, where m + n = 24) on their solution behavior and interactions with phospholipids were examined. We find that although these peptides exhibit strong alpha-helix forming propensities in water, membrane-mimetic media, and lipid model membranes, the stability of the helices decreases as the Leu content decreases. Also, their binding to reversed phase high-performance liquid chromatography columns is largely determined by their hydrophobicity and generally decreases with decreases in the Leu/Ala ratio. However, the retention of these peptides by such columns is also affected by the distribution of hydrophobic residues on their helical surfaces, being further enhanced when peptide helical hydrophobic moments are increased by clustering hydrophobic residues on one side of the helix. This clustering of hydrophobic residues also increases peptide propensity for self-aggregation in aqueous media and enhances partitioning of the peptide into lipid bilayer membranes. We also find that the peptides LA3LA2 [acetyl-K2-(LAAALAA)3LAA-K2-amide] and particularly LA6 [acetyl-K2-(LAAAAAA)3LAA-K2-amide] associate less strongly with and perturb the thermotropic phase behavior of phosphatidylcholine bilayers much less than peptides with higher L/A ratios. These results are consistent with free energies calculated for the partitioning of these peptides between water and phospholipid bilayers, which suggest that LA3LA2 has an equal tendency to partition into water and into the hydrophobic core of phospholipid model membranes, whereas LA6 should strongly prefer the aqueous phase. We conclude that for alpha-helical peptides of this type, Leu/Ala ratios of greater than 7/17 are required for stable transmembrane associations with phospholipid bilayers.     

Impact of histidine residues on the transmembrane helices of viroporins

Abstract

The role of histidine in channel-forming transmembrane (TM) helices was investigated by comparing the TM helices from Virus protein ‘u' (Vpu) and the M2 proton channel. Both proteins are members of the viroporin family of small membrane proteins that exhibit ion channel activity, and have a single TM helix that is capable of forming oligomers. The TM helices from both proteins have a conserved tryptophan towards the C-terminus. Previously, alanine 18 of Vpu was mutated to histidine in order to artificially introduce the same HXXXW motif that is central to the proton channel activity of M2. Interestingly, the mutated Vpu TM resulted in an increase in helix tilt angle of 11° in lipid bilayers compared to the wild-type Vpu TM. Here, we find the reverse, when histidine 37 of the HXXXW motif in M2 was mutated to alanine, it decreased the helix tilt by 10° from that of wild-type M2. The tilt change is independent of both the helix length and the presence of tryptophan. In addition, compared to wild-type M2, the H37A mutant displayed lowered sensitivity to proton concentration. We also found that the solvent accessibility of histidine-containing M2 is greater than without histidine. This suggests that the TM helix may increase the solvent exposure by changing its tilt angle in order to accommodate a polar/charged residue within the hydrophobic membrane region. The comparative results of M2, Vpu and their mutants demonstrated the significance of histidine in a transmembrane helix and the remarkable plasticity of the function and structure of ion channels stemming from changes at a single amino acid site.     

Structure stability of lytic peptides during their interactions with lipid bilayers.

In this work, molecular dynamics simulations were used to examine the consequences of a variety of analogs of cecropin A on lipid bilayers. Analog sequences were constructed by replacing either the N- or C-terminal helix with the other helix in native or reverse sequence order, by making palindromic peptides based on both the N- and C-terminal helices, and by deleting the hinge region. The structure of the peptides was monitored throughout the simulation. The hinge region appeared not to assist in maintaining helical structure but help in motion flexibility. In general, the N-terminal helix of peptides was less stable than the C-terminal one during the interaction with anionic lipid bilayers. Sequences with hydrophobic helices tended to regain helical structure after an initial loss while sequences with amphipathic helices were less able to do this. The results suggests that hydrophobic design peptides have a high structural stability in an anionic membrane and are the candidates for experimental investigation.     

The nature of the hydrophobic binding of small peptides at the bilayer interface: implications for the insertion of transbilayer helices

One method of obtaining useful information about the physical chemistry of peptide/bilayer interactions is to relate thermodynamic parameters of the interactions to structural parameters obtained by diffraction methods. We report here the results of the application of this approach to interactions of hydrophobic tripeptides of the form Ala-X-Ala-O-tert-butyl with lipid bilayers. The thermodynamic constants (delta Gt, delta Ht, and delta St) for the transfer of the tripeptides from water into DMPC vesicles were determined for X = Leu, Phe, and Trp and found to be consistent with those expected for hydrophobic interactions above the phase transition of DMPC. Combining these results with the earlier ones of Jacobs and White [(1986) Biochemistry 25, 2605-2612], the favorable free energies of transfer with different amino acids in the -X- position increase in the order Gly less than Ala less than Leu less than Phe less than Trp in agreement with the Nozaki and Tanford [(1971) J. Biol. Chem. 246, 2211-2217] hydrophobicity scale. Determination of the location of Ala-[2H5]Trp-Ala-O-tert-butyl in oriented DOPC bilayers by neutron diffraction shows that the most hydrophobic peptide of the series is confined to the bilayer headgroup/water region. Refinement of the diffraction measurements shows that only 13% of the tryptophan is associated with the hydrocarbon core. The distribution of the water tends to mirror that of the peptide. Unlike peptide-free bilayers, 5% of the water penetrates the hydrocarbon, which is about 100-fold greater than expected. A quantitative thermodynamic analysis of the interfacial binding of the peptides suggests that (1) the hydrophobic interactions are 60-70% complete upon binding at the bilayer interface, (2) the interface is likely to play an important role in helix formation and insertion, (3) the hydrogen bond status of amino acid side chains is crucial to insertion, and (4) an a priori lack of knowledge of the status of such bonds could limit the precision of hydrophobicity plots. We introduce an interfacial hydrophobicity scale, IFH(h), with a variable hydrogen bond parameter (h) that permits one to consider explicitly hydrogen bonding in transbilayer helix searches.     

Strong Correlation Between Statistical Transmembrane Tendency and Experimental Hydrophobicity Scales for Identification of Transmembrane Helices

Direct physical chemistry measurements of the hydrophobicity of amino acids or their derivatives have often been used to estimate the propensity of amino acids to participate in transmembrane helices. In this short note, it is found that there is a very high degree of correlation (r = 0.944–0.965) between an average physical chemistry hydrophobicity scale (an average of scales derived, e.g., from the solubility of amino acid derivatives in organic solvents versus water or their binding to hydrophobic particles) and the statistically based transmembrane tendency scale (derived from the relative abundance of residues in known transmembrane and soluble protein sequences (Zhao and London, Protein Sci 15:1987–2001, 2006)). This correlation indicates that, other than hydrophobicity, amino acid properties/interactions that promote or inhibit transmembrane helix formation in a specific membrane protein largely cancel out when averaged over all transmembrane sequences. In other words, other than hydrophobicity, there are no properties of a specific amino acid residue within a hydrophobic segment that have a strong systematic effect upon transmembrane helix formation independent of the remainder of the sequence in that hydrophobic segment. However, proline is an exception to this rule.