Effect of sequence hydrophobicity and bilayer width upon the minimum length required for the formation of transmembrane helices in membranes

The minimum hydrophobic length necessary to form a transmembrane (TM) helix in membranes was investigated using model membrane-inserted hydrophobic helices. The fluorescence of a Trp at the center of the sequence and its sensitivity to quenching were used to ascertain helix position within the membrane. Peptides with hydrophobic cores composed of poly(Leu) were compared to sequences containing a poly 1:1 Leu:Ala core (which have a hydrophobicity typical of natural TM helices). Studies varying bilayer width revealed that the poly(Leu) core peptides predominately formed a TM state when the bilayer width exceeded hydrophobic sequence length by (i.e. when negative mismatch was) up to ∼ 11-12 Å (e.g. the case of a 11-12 residue hydrophobic sequence in bilayers with a biologically relevant width, i.e. dioleoylphosphatidylcholine (DOPC) bilayers), while poly(LeuAla) core peptides formed predominantly TM state with negative mismatch of up to 9 Å (a 13 residue hydrophobic sequence in DOPC bilayers). This indicates that minimum length necessary to form a predominating amount of a TM state (minimum TM length) is only modestly hydrophobicity-dependent for the sequences studied here, and a formula that defines the minimum TM length as a function of hydrophobicity for moderately-to-highly hydrophobic sequences was derived. The minimum length able to form a stable TM helix for alternating LeuAla sequences, and that for sequences with a Leu block followed by an Ala block, was similar, suggesting that a hydrophobicity gradient along the sequence may not be an important factor in TM stability. TM stability was also similar for sequences flanked by different charged ionizable residues (Lys, His, Asp). However, ionizable flanking residues destabilized the TM configuration much more when charged than when uncharged. The ability of short hydrophobic sequences to form TM helices in membranes in the presence of substantial negative mismatch implies that lipid bilayers have a considerable ability to adjust to negative mismatch, and that short TM helices may be more common than generally believed. Factors that modulate the ability of bilayers to adjust to mismatch may strongly affect the configuration of short hydrophobic helices.     

Molecular modeling of the amphipathic helices of the plasma apolipoproteins.

In this paper we propose a classification of the amphipathic helical repeats occurring in the plasma apolipoprotein sequences. It is based upon the calculation of the molecular hydrophobicity potential around the helical segments. The repeats were identified using a new autocorrelation matrix, based upon similarities of hydrophobic and hydrophilic properties of the amino acid residues within the apolipoprotein sequences. The helices were constructed by molecular modeling, the molecular hydrophobicity potential was calculated, and isopotential contour lines drawn around the helices yielded a three-dimensional visualization of the hydrophobicity potential. Two classes of apolipoproteins could be differentiated by comparing the hydrophobic angles obtained by projection of the isopotential contour lines on a plane perpendicular to the long axis of the helix. The isopotential contour lines around apo AI, AIV, and E are more hydrophilic than hydrophobic, whereas they are of similar intensity for apo AII, CI, and CIII. In both cases discoidal lipid-protein complexes are generated, with the amphipathic helices around the edge of the lipid core. The long axis of the helices is oriented parallel to the phospholipid acyl chains and the hydrophilic side of the helix toward the aqueous phase. As a result of the differences in hydrophobicity potential, the contact between the hydrophobic side of the helices and the phospholipid acyl chains is larger for apo AII, CI, and CIII than for the other apolipoproteins. This might account for the greater stability of the discoidal complexes generated between phospholipids and these apoproteins.     

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.     

Membrane Partitioning of the Pore-Forming Domain of Colicin A. Role of the Hydrophobic Helical Hairpin

The colicins are bacteriocins that target Escherichia coli and kill bacterial cells through different mechanisms. Colicin A forms ion channels in the inner membranes of nonimmune bacteria. This activity resides exclusively in its C-terminal fragment (residues 387–592). The soluble free form of this domain is a 10 α-helix bundle. The hydrophobic helical hairpin, H8–H9, is buried inside the structure and shielded by eight amphipathic surface helices. The interaction of the C-terminal colicin A domain and several chimeric variants with lipidic vesicles was examined here by isothermal titration calorimetry. In the mutant constructions, natural sequences of the hydrophobic helices H8 and H9 were either removed or substituted by polyalanine or polyleucine. All the constructions fully associated with DOPG liposomes including the mutant that lacked helices H8 and H9, indicating that amphipathic rather than hydrophobic helices were the major determinants of the exothermic binding reactions. Alanine is not specially favored in the lipid-bound form; the chimeric construct with polyalanine produced lower enthalpy gain. On the other hand, the large negative heat capacities associated with partitioning, a characteristic feature of the hydrophobic effect, were found to be dependent on the sequence hydrophobicity of helices H8 and H9.     

Membrane Partitioning of the Pore-Forming Domain of Colicin A. Role of the Hydrophobic Helical Hairpin

The colicins are bacteriocins that target Escherichia coli and kill bacterial cells through different mechanisms. Colicin A forms ion channels in the inner membranes of nonimmune bacteria. This activity resides exclusively in its C-terminal fragment (residues 387–592). The soluble free form of this domain is a 10 α-helix bundle. The hydrophobic helical hairpin, H8–H9, is buried inside the structure and shielded by eight amphipathic surface helices. The interaction of the C-terminal colicin A domain and several chimeric variants with lipidic vesicles was examined here by isothermal titration calorimetry. In the mutant constructions, natural sequences of the hydrophobic helices H8 and H9 were either removed or substituted by polyalanine or polyleucine. All the constructions fully associated with DOPG liposomes including the mutant that lacked helices H8 and H9, indicating that amphipathic rather than hydrophobic helices were the major determinants of the exothermic binding reactions. Alanine is not specially favored in the lipid-bound form; the chimeric construct with polyalanine produced lower enthalpy gain. On the other hand, the large negative heat capacities associated with partitioning, a characteristic feature of the hydrophobic effect, were found to be dependent on the sequence hydrophobicity of helices H8 and H9.     

The transmembrane helices of the L, M, and N subunits of Complex I from E. coli can be assigned on the basis of conservation and hydrophobic moment analysis

An assignment of the transmembrane helices of subunits L, M, and N of the Escherichia coli Complex I has been made from the helices as determined in a recent crystal structure [Efromov et al., Nature (2010) 465, 441-446]. The amino acid sequences of the three subunits were evaluated for hydrophobicity, and hydrophobic moments, to identify the helices that are likely to be in contact with membrane lipids. Using 29 closely related species, a similar analysis of average conservation, and conservation moments was performed. In each subunit, transmembrane helices 9 and 12 are predicted to form the discontinuous helices, which are likely to play a key role in function.     

Frequencies of amino acid strings in globular protein sequences indicate suppression of blocks of consecutive hydrophobic residues

Patterns of hydrophobic and hydrophilic residues play a major role in protein folding and function. Long, predominantly hydrophobic strings of 20-22 amino acids each are associated with transmembrane helices and have been used to identify such sequences. Much less attention has been paid to hydrophobic sequences within globular proteins. In prior work on computer simulations of the competition between on-pathway folding and off-pathway aggregate formation, we found that long sequences of consecutive hydrophobic residues promoted aggregation within the model, even controlling for overall hydrophobic content. We report here on an analysis of the frequencies of different lengths of contiguous blocks of hydrophobic residues in a database of amino acid sequences of proteins of known structure. Sequences of three or more consecutive hydrophobic residues are found to be significantly less common in actual globular proteins than would be predicted if residues were selected independently. The result may reflect selection against long blocks of hydrophobic residues within globular proteins relative to what would be expected if residue hydrophobicities were independent of those of nearby residues in the sequence.     

Prediction of protein helices with a derivative of the strip-of-helix hydrophobicity algorithm

The strip-of-helix hydrophobicity algorithm was devised to identify protein sequences which, when coiled as alpha or 3(10) helices, had one axial, hydrophobic strip and otherwise variably hydrophilic residues. The strip-of-helix hydrophobicity algorithm also ranked such sequences according to an index, the mean hydrophobicity of amino acids in the axial strip. This algorithm well predicted T cell-presented fragments of antigenic proteins. A derivative of this algorithm (the structural helices algorithm (SHA] was tested for the prediction of helices in crystallographically defined proteins. For the SHA, eight amino acid sequences, 2 cycles plus one amino acid in an alpha helix, with strip-of-helix hydrophobicity indices greater than 2.5, were selected with overlapping segments joined. These selections were terminated according to simple "capping rules," which took into account the roles of N-terminal Asn or Pro and C-terminal Gly in the stability of helices. In analyses of 35 crystallographically defined proteins with known alpha and 3(10) helices, the predictions with the SHA overlapped (had overlap indices x greater than or equal to 0.5) with 34% of known helices, touched (had overlap indices 0.5 greater than x greater than 0) or overlapped with 66% of known helices, or were neighboring (came within 6 residues) or touched or overlapped with 82% of known helices. At each level of judging the quality of prediction, the SHA was usually less sensitive (correct predictions/total number of known helices) and more efficient (correct predictions/total number of predictions) than the Chou-Fasman and Garnier-Robson methods. It was simpler in design and calculation. The chemical mechanisms underlying these algorithms appear to apply both to protein folding and to selection of T cell-presented antigenic sequences.     

Topography of helices 5-7 in membrane-inserted diphtheria toxin T domain: identification and insertion boundaries of two hydrophobic sequences that do not form a stable transmembrane hairpin

The T domain of diphtheria toxin undergoes a low pH-induced conformational change that allows it to penetrate cell membranes. T domain hydrophobic helices 8 and 9 can adopt two conformations, one close to the membrane surface (P state) and a second in which they apparently form a transmembrane hairpin (TM state). We have now studied T domain helices 5-7, a second cluster of hydrophobic helices, using Cys-scanning mutagenesis. After fluorescently labeling a series of Cys residues, penetration into a non-polar environment, accessibility to externally added antibodies, and relative depth in the bilayer were monitored. It was found that helices 5-7 insert shallowly in the P state and deeply in the TM state. Thus, the conformational changes in helices 5-7 are both similar and somehow linked to those in helices 8 and 9. The boundaries of deeply inserting sequences were also identified. One deeply inserted segment was found to span residues 270 to 290, which overlaps helix 5, and a second spanned residues 300 to 320, which includes most of helix 6 and all of helix 7. This indicates that helices 6 and 7 form a continuous hydrophobic segment despite their separation by a Pro-containing kink. Additionally, it is found that in the TM state some residues in the hydrophilic loop between helices 5 and 6 become more highly exposed than they are in the P state. Their exposure to external solution in the TM state indicates that helices 5-7 do not form a stable transmembrane hairpin. However, helix 5 and/or helices 6 plus 7 could form transmembrane structures that are in equilibrium with non-transmembrane states, or be kinetically prevented from forming a transmembrane structure. How helices 5-7 might influence the mechanism by which the T domain aids translocation of the diphtheria toxin A chain across membranes is discussed.     

Analyzing topography of membrane-inserted diphtheria toxin T domain using BODIPY-streptavidin: at low pH, helices 8 and 9 form a transmembrane hairpin but helices 5-7 form stable nonclassical inserted segments on the cis side of the bilayer

Low pH-induced membrane insertion by diphtheria toxin T domain is crucial for A chain translocation into the cytoplasm. To define the membrane topography of the T domain, the exposure of biotinylated Cys residues to the cis and trans bilayer surfaces was examined using model membrane vesicles containing a deeply inserted T domain. To do this, the reactivity of biotin with external and vesicle-entrapped BODIPY-labeled streptavidin was measured. The T domain was found to insert with roughly 70-80% of the molecules in the physiologically relevant orientation. In this orientation, residue 349, located in the loop between hydrophobic helices 8 and 9, was exposed to the trans side of the bilayer, while other solution-exposed residues along the hydrophobic helices 5-9 region of the T domain located near the cis surface. A protocol developed to detect the movement of residues back and forth across the membranes demonstrated that T domain sequences did not rapidly equilibrate between the cis and the trans sides of the bilayer. Binding streptavidin to biotinylated residues prior to membrane insertion only inhibited T domain pore formation for residues in the loop between helices 8 and 9. Pore formation experiments used an approach avoiding interference from transient membrane defects/leakage that may occur upon the initial insertion of protein. Combined, these results indicate that at low pH hydrophobic helices 8 and 9 form a transmembrane hairpin, while hydrophobic helices 5-7 form a nonclassical deeply inserted nontransmembraneous state. We propose that this represents a novel pre-translocation state that is distinct from a previously defined post-translocation state.     

Differentiation of lipid-associating helices by use of three-dimensional molecular hydrophobicity potential calculations.

Several types of lipid-associating helices exist: transmembrane helices such as in receptor proteins, pore-forming helices in ion channel proteins, fusion-inducing peptides in viral proteins, and amphipathic helices such as in plasma apolipoproteins. In order to propose a classification of these helices according to their molecular properties, we introduce the concept of molecular hydrophobicity potential for such helical segments. The calculation of this parameter for alpha-helices enables the visualization of the hydrophobic and hydrophilic envelopes around the peptide and their three-dimensional representation by molecular graphics. We have used this parameter to differentiate between pore-forming helices with a hydrophobic envelope larger than the hydrophilic component, membrane-spanning helices surrounded almost entirely by an hydrophobic envelope, fusiogenic peptides with an hydrophobicity gradient both around the helix and along the axis, and finally, amphipathic helices with a predominantly hydrophilic envelope. The structure of the lipid-protein complexes is determined by a number of different interactions: the hydrophobic interaction of the apolar faces of the helices with lipids, the polar interaction of the hydrophilic sides of different helices with each other, and the interaction of hydrophilic residues with the aqueous solvent. The relative magnitude of the hydrophobic and hydrophilic envelopes accounts for the differences in the structure of the lipid-protein complexes. Purely hydrophobic interactions stabilize transmembrane helical segments, while hydrophobic interactions with the lipid phase and with each other are involved in the stabilization of the pore-forming helices. In contrast, both hydrophobic interactions with the lipids and hydrophilic interactions with the aqueous phase contribute to the arrangement of amphipathic helices around the edges of the discoidal lipid-apoprotein complexes.     

Structural clues in the sequences of the aquaporins

The large number of sequences available for the aquaporin family represents a valuable source of information to incorporate into three-dimensional structure determination. Phylogenetic analysis was used to define type sequences to avoid extreme over-representation of some subfamilies, and as a measure of the quality of multiple sequence alignment. Inspection of the sequence alignment suggested eight conserved segments that define the core architecture of six transmembrane helices and two functional loops, B and E, projecting into the plane of the membrane. The sum of the core segments and the minimum lengths of the interlinking loops constitute the 208 residues necessary to satisfy the aquaporin architecture. Analysis of hydrophobic and conservation periodicity and of correlated mutations across the alignment indicated the likely assignment and orientation of the helices in the bilayer. This assignment is examined with respect to the structure of the erythrocyte aquaporin 1 determined by electron crystallography. The aquaporin 1 tetramer is described as three rings of helices, each ring with a different exposure to the lipid environment. The sequence analysis clearly suggests that two helices are exposed along their whole lengths, two helices are exposed only at their N termini, and two helices are not exposed to lipid. It is further proposed that, besides loops B and E, the highly conserved motifs on helices 1 and 4, ExxxTxxF/L, could line the water channel.     

Sequence context and modified hydrophobic moment plots help identify 'horizontal' surface helices in transmembrane protein structure prediction

Transmembrane proteins make up at least one-fifth of the genome of most organisms and are critical components of key pathways for cell survival and interactions with the environment. The function of helices found at the membrane surface in transmembrane proteins has not been greatly explored, but it is likely that they play an ancillary role to membrane spanning helices and are analogous to the surface active helices of peripheral membrane proteins, being involved in: lipid association, membrane perturbation, transmembrane signal transduction and regulation, and transmembrane helical bundle formation. Due to the difficulties in obtaining high-resolution structural data for this class of proteins, structure-from-sequence predictive methods continue to be developed as a means to obtain structural models for these largely intractable systems. A simple but effective variant of the hydrophobic moment analysis of amino acid sequences is described here as part of a protocol for distinguishing helical sequences that are parallel to or 'horizontal' at the membrane bilayer/aqueous phase interface from helices that are membrane-embedded or located in extra-membranous domains. This protocol when tested on transmembrane spanning protein amino acid sequences not used in its development, was found to be 84-91% accurate when the results were compared to the partition locations in the corresponding structures determined by X-ray crystallography, and 72% accurate in determining which helices lie horizontal or near horizontal at the lipid interface.     

The Positive Inside Rule Is Stronger When Followed by a Transmembrane Helix

The translocon recognizes transmembrane helices with sufficient level of hydrophobicity and inserts them into the membrane. However, sometimes less hydrophobic helices are also recognized. Positive inside rule, orientational preferences of and specific interactions with neighboring helices have been shown to aid in the recognition of these helices, at least in artificial systems. To better understand how the translocon inserts marginally hydrophobic helices, we studied three naturally occurring marginally hydrophobic helices, which were previously shown to require the subsequent helix for efficient translocon recognition. We find no evidence for specific interactions when we scan all residues in the subsequent helices. Instead, we identify arginines located at the N-terminal part of the subsequent helices that are crucial for the recognition of the marginally hydrophobic transmembrane helices, indicating that the positive inside rule is important. However, in two of the constructs, these arginines do not aid in the recognition without the rest of the subsequent helix; that is, the positive inside rule alone is not sufficient. Instead, the improved recognition of marginally hydrophobic helices can here be explained as follows: the positive inside rule provides an orientational preference of the subsequent helix, which in turn allows the marginally hydrophobic helix to be inserted; that is, the effect of the positive inside rule is stronger if positively charged residues are followed by a transmembrane helix. Such a mechanism obviously cannot aid C-terminal helices, and consequently, we find that the terminal helices in multi-spanning membrane proteins are more hydrophobic than internal helices.     

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.     

Frequencies of hydrophobic and hydrophilic runs and alternations in proteins of known structure

Patterns of alternation of hydrophobic and polar residues are a profound aspect of amino acid sequences, but a feature not easily interpreted for soluble proteins. Here we report statistics of hydrophobicity patterns in proteins of known structure in a current protein database as compared with results from earlier, more limited structure sets. Previous studies indicated that long hydrophobic runs, common in membrane proteins, are underrepresented in soluble proteins. Long runs of hydrophobic residues remain significantly underrepresented in soluble proteins, with none longer than 16 residues observed. These long runs most commonly occur as buried alpha helices, with extended hydrophobic strands less common. Avoiding aggregation of partially folded intermediates during intracellular folding remains a viable explanation for the rarity of long hydrophobic runs in soluble proteins. Comparison between database editions reveals robustness of statistics on aqueous proteins despite an approximately twofold increase in nonredundant sequences. The expanded database does now allow us to explain several deviations of hydrophobicity statistics from models of random sequence in terms of requirements of specific secondary structure elements. Comparison to prior membrane-bound protein sequences, however, shows significant qualitative changes, with the average hydrophobicity and frequency of long runs of hydrophobic residues noticeably increasing between the database editions. These results suggest that the aqueous proteins of solved structure may represent an essentially complete sample of the universe of aqueous sequences, while the membrane proteins of known structure are not yet representative of the universe of membrane-associated proteins, even by relatively simple measures of hydrophobic patterns.     

Residues in the longitudinal, hydrophobic strip-of-helix relate to terminations and crossings of alpha-helices.

An alpha-helix terminates when the virtual extension of its most hydrophobic, longitudinal strip containing Leu, Ile, Val, Phe, and Met lacks those residues. In each of 247 helices a template was fitted to maximize the mean hydrophobicity of positions forming a longitudinal strip-of-helix. The template was then extended into sequences beyond the ends of the helices. Leu, Ile, Val, Phe, and Met occurred in positions in the longitudinal strip-of-helix at an increased frequency (p less than 0.001), but in the first and second positions beyond either end of each true helix, they occurred at the same frequency as for their empirical distribution over all the proteins. Excesses of Asp and Glu were found in the N-terminal loop, and of Arg, His, and Lys in specific positions about the C terminus of helices. The longitudinal hydrophobic strip, the smallest amino acid in that strip, and charged amino acids in that strip, related to rotational and longitudinal orientation of alpha-helices in 15 proteins. Adjacent helices generally crossed through their longitudinal hydrophobic strips. They usually crossed through the smallest residue in the strip. Charged residues, when they occurred in the strips, were excluded from the crossing regions.     

A Calpha model for the transmembrane alpha helices of gap junction intercellular channels

Gap junction channels connect the cytoplasms of apposed cells via an intercellular conduit formed by the end-to-end docking of two hexameric hemichannels called connexons. We used electron cryomicroscopy to derive a three-dimensional density map at 5.7 angstroms in-plane and 19.8 angstroms vertical resolution, allowing us to identify the positions and tilt angles for the 24 alpha helices within each hemichannel. The four hydrophobic segments in connexin sequences were assigned to the alpha helices in the map based on biochemical and phylogenetic data. Analyses of evolutionary conservation and compensatory mutations in connexin evolution identified the packing interfaces between the helices. The final model, which specifies the coordinates of Calpha atoms in the transmembrane domain, provides a structural basis for understanding the different physiological effects of almost 30 mutations and polymorphisms in terms of structural deformations at the interfaces between helices, revealing an intimate connection between molecular structure and disease.     

Critical functional role of the COOH-terminal ends of longitudinal hydrophobic strips in alpha-helices of T4 lysozyme.

The sensitivity of bacteriophage T4 lysozyme function to amino acid substitutions at defined positions in and around the longitudinal, hydrophobic strips of 9 alpha-helices was assessed after systematic replacement of each residue in the protein with a series of 13 amino acids. The hydrophobic strips were defined by identifying the longitudinal sectors in the helices with the highest mean residue hydrophobicities. Sensitivity to mutation (the percentage of replacements leading to loss of function) was calculated for each residue in the following positions: whole protein, helices, hydrophobic strips, other positions within the helices, and various positions within the hydrophobic strips as well as their extensions beyond the helices. Substitutions at positions in the hydrophobic strips led more frequently to loss of function than substitutions in the protein as a whole. One subset, the COOH-terminal hydrophobic strip residues, is apparently critical; substitutions of these residues (but not of their NH2-terminal counterparts) led at least as frequently to loss of function as substitutions of solvent-inaccessible residues, and nearly as frequently as substitutions of the most highly conserved residues.