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.     

Extracellular loops 1 and 3 and their associated transmembrane regions of the calcitonin receptor-like receptor are needed for CGRP receptor function

The first and third extracellular loops (ECL) of G protein-coupled receptors (GPCRs) have been implicated in ligand binding and receptor function. This study describes the results of an alanine/leucine scan of ECLs 1 and 3 and loop-associated transmembrane (TM) domains of the secretin-like GPCR calcitonin receptor-like receptor which associates with receptor activity modifying protein 1 to form the CGRP receptor. Leu195Ala, Val198Ala and Ala199Leu at the top of TM2 all reduced αCGRP-mediated cAMP production and internalization; Leu195Ala and Ala199Leu also reduced αCGRP binding. These residues form a hydrophobic cluster within an area defined as the “minor groove” of rhodopsin-like GPCRs. Within ECL1, Ala203Leu and Ala206Leu influenced the ability of αCGRP to stimulate adenylate cyclase. In TM3, His219Ala, Leu220Ala and Leu222Ala have influences on αCGRP binding and cAMP production; they are likely to indirectly influence the binding site for αCGRP as well as having an involvement in signal transduction. On the exofacial surfaces of TMs 6 and 7, a number of residues were identified that reduced cell surface receptor expression, most noticeably Leu351Ala and Glu357Ala in TM6. The residues may contribute to the RAMP1 binding interface. Ile360Ala impaired αCGRP-mediated cAMP production. Ile360 is predicted to be located close to ECL2 and may facilitate receptor activation. Identification of several crucial functional loci gives further insight into the activation mechanism of this complex receptor system and may aid rational drug design.     

TIM23-mediated insertion of transmembrane α-helices into the mitochondrial inner membrane

While overall hydrophobicity is generally recognized as the main characteristic of transmembrane (TM) α-helices, the only membrane system for which there are detailed quantitative data on how different amino acids contribute to the overall efficiency of membrane insertion is the endoplasmic reticulum (ER) of eukaryotic cells. Here, we provide comparable data for TIM23-mediated membrane protein insertion into the inner mitochondrial membrane of yeast cells. We find that hydrophobicity and the location of polar and aromatic residues are strong determinants of membrane insertion. These results parallel what has been found previously for the ER. However, we see striking differences between the effects elicited by charged residues flanking the TM segments when comparing the mitochondrial inner membrane and the ER, pointing to an unanticipated difference between the two insertion systems.     

Helix packing moments reveal diversity and conservation in membrane protein structure

Helical membrane proteins are more tightly packed and the packing interactions are more diverse than those found in helical soluble proteins. Based on a linear correlation between amino acid packing values and interhelical propensity, we propose the concept of a helix packing moment to predict the orientation of helices in helical membrane proteins and membrane protein complexes. We show that the helix packing moment correlates with the helix interfaces of helix dimers of single pass membrane proteins of known structure. Helix packing moments are also shown to help identify the packing interfaces in membrane proteins with multiple transmembrane helices, where a single helix can have multiple contact surfaces. Analyses are described on class A G protein-coupled receptors (GPCRs) with seven transmembrane helices. We show that the helix packing moments are conserved across the class A family of GPCRs and correspond to key structural contacts in rhodopsin. These contacts are distinct from the highly conserved signature motifs of GPCRs and have not previously been recognized. The specific amino acid types involved in these contacts, however, are not necessarily conserved between subfamilies of GPCRs, indicating that the same protein architecture can be supported by a diverse set of interactions. In GPCRs, as well as membrane channels and transporters, amino acid residues with small side-chains (Gly, Ala, Ser, Cys) allow tight helix packing by mediating strong van der Waals interactions between helices. Closely packed helices, in turn, facilitate interhelical hydrogen bonding of both weakly polar (Ser, Thr, Cys) and strongly polar (Asn, Gln, Glu, Asp, His, Arg, Lys) amino acid residues. We propose the use of the helix packing moment as a complementary tool to the helical hydrophobic moment in the analysis of transmembrane sequences.     

Elucidation of the role of peptide linker in calcium-sensing receptor activation process

Family 3 G-protein-coupled receptors (GPCRs), which includes metabotropic glutamate receptors (mGluRs), sweet and "umami" taste receptors (T1Rs), and the extracellular calcium-sensing receptor (CaR), represent a distinct group among the superfamily of GPCRs characterized by large amino-terminal extracellular ligand-binding domains (ECD) with homology to bacterial periplasmic amino acid-binding proteins that are responsible for signal detection and receptor activation through as yet unresolved mechanism(s) via the seven-transmembrane helical domain (7TMD) common to all GPCRs. To address the mechanism(s) by which ligand-induced conformational changes are conveyed from the ECD to the 7TMD for G-protein activation, we altered the length and composition of a 14-amino acid linker segment common to all family 3 GPCRs except GABA(B) receptor, in the CaR by insertion, deletion, and site-directed mutagenesis of specific highly conserved residues. Small alterations in the length and composition of the linker impaired cell surface expression and abrogated signaling of the chimeric receptors. The exchange of nine amino acids within the linker of CaR with the homologous sequence of mGluR1, however, preserved receptor function. Ala substitution for the four highly conserved residues within this amino acid sequence identified a Leu at position 606 of the CaR critical for cell surface expression and signaling. Substitution of Leu(606) for Ala resulted in impaired cell surface expression. However, Ile and Val substitutions displayed strong activating phenotypes. Disruption of the linker by insertion of nine amino acids of a random-coiled structure uncoupled the ECD from regulating the 7TMD. These data are consistent with a model of receptor activation in which the peptide linker, and particularly Leu(606), provides a critical interaction for the CaR signal transmission, a finding likely to be relevant for all family 3 GPCRs containing this conserved motif.     

Conformational preferences of non-polar amino acid residues: an additional factor in amyloid formation

Amyloid consists of β-sheet polymers and is associated with disease and with functional assemblies. Amyloid-forming proteins differ widely in native structures and sequences. We describe here how conformational preferences of non-polar amino acid residues can affect amyloid formation. The most non-polar residues promote either β-strands (Val, Ile, Phe, and Cys, VIFC) or α-helices (Leu, Ala, and Met, LAM), while the most polar residues promote only α-helices. For 12 proteins associated with disease, the localizations of the amyloid core regions are known. Eleven of these contain segments that are biased for VIFC, but essentially lack segments that are biased for LAM. For the amyloid β-peptide associated with Alzheimer’s disease and an amyloidogenic fragment of the prion protein, observed effects of mutations support that VIFC bias favors formation of β-sheet aggregates and amyloid, while LAM bias prevents it. VIFC and LAM profiles combine information on secondary structure propensities and polarity, and add a simple criterion to the prediction of amyloidogenic regions.     

Empirical lipid propensities of amino acid residues in multispan alpha helical membrane proteins

Characterizing the interactions between amino acid residues and lipid molecules is important for understanding the assembly of transmembrane helices and for studying membrane protein folding. In this study we develop TMLIP (TransMembrane helix-LIPid), an empirically derived propensity of individual residue types to face lipid membrane based on statistical analysis of high-resolution structures of membrane proteins. Lipid accessibilities of amino acid residues within the transmembrane (TM) region of 29 structures of helical membrane proteins are studied with a spherical probe of radius of 1.9 A. Our results show that there are characteristic preferences for residues to face the headgroup region and the hydrocarbon core region of lipid membrane. Amino acid residues Lys, Arg, Trp, Phe, and Leu are often found exposed at the headgroup regions of the membrane, where they have high propensity to face phospholipid headgroups and glycerol backbones. In the hydrocarbon core region, the strongest preference for interacting with lipids is observed for Ile, Leu, Phe and Val. Small and polar amino acid residues are usually buried inside helical bundles and are strongly lipophobic. There is a strong correlation between various hydrophobicity scales and the propensity of a given residue to face the lipids in the hydrocarbon region of the bilayer. Our data suggest a possibly significant contribution of the lipophobic effect to the folding of membrane proteins. This study shows that membrane proteins have exceedingly apolar exteriors rather than highly polar interiors. Prediction of lipid-facing surfaces of boundary helices using TMLIP1 results in a 54% accuracy, which is significantly better than random (25% accuracy). We also compare performance of TMLIP with another lipid propensity scale, kPROT, and with several hydrophobicity scales using hydrophobic moment analysis.     

Multiple distant amino acid residues present in the serpentine region of the follitropin receptor modulate the rate of agonist-induced internalization

The amino acid sequences of the human (h) and rat (r) follitropin receptors (FSHR) are approximately 89% identical, but the half-time of internalization of agonist mediated by the rFSHR is approximately 3 times faster than that of the hFSHR. Chimeras of the hFSHR and the rFSHR showed that this difference in rate is dictated mostly by the serpentine domain. Further analysis identified six residues, two non-contiguous residues in the transmembrane helix 4 (Leu/Thr in the rFSHR and Met/Ile in the hFSHR), three non-contiguous residues in the third intracellular loop (Thr/Thr/Lys in the rFSHR and Ile/Asn/Arg in the hFSHR), and one in transmembrane helix 7 (Tyr in the rFSHR and His in the hFSHR) that are fully responsible for the difference in the rates of internalization of the hFSHR and the rFSHR.     

Comprehensive analysis of the numbers,lengths and amino acid compositions of transmembrane helices in prokaryotic,eukaryotic and viral integral membrane proteins of high-resolution structure

We report a comprehensive analysis of the numbers, lengths and amino acid compositions of transmembrane helices in 235 high-resolution structures of integral membrane proteins. The properties of 1551 transmembrane helices in the structures were compared with those obtained by analysis of the same amino acid sequences using topology prediction tools. Explanations for the 81 (5.2%) missing or additional transmembrane helices in the prediction results were identified. Main reasons for missing transmembrane helices were mis-identification of N-terminal signal peptides, breaks in α-helix conformation or charged residues in the middle of transmembrane helices and transmembrane helices with unusual amino acid composition. The main reason for additional transmembrane helices was mis-identification of amphipathic helices, extramembrane helices or hairpin re-entrant loops. Transmembrane helix length had an overall median of 24 residues and an average of 24.9 ± 7.0 residues and the most common length was 23 residues. The overall content of residues in transmembrane helices as a percentage of the full proteins had a median of 56.8% and an average of 55.7 ± 16.0%. Amino acid composition was analysed for the full proteins, transmembrane helices and extramembrane regions. Individual proteins or types of proteins with transmembrane helices containing extremes in contents of individual amino acids or combinations of amino acids with similar physicochemical properties were identified and linked to structure and/or function. In addition to overall median and average values, all results were analysed for proteins originating from different types of organism (prokaryotic, eukaryotic, viral) and for subgroups of receptors, channels, transporters and others.     

Identification of amino acids involved in protein structural uniqueness: implication for de novo protein design

Structural uniqueness is characteristic of native proteins and is essential to express their biological functions. The major factors that bring about the uniqueness are specific interactions between hydrophobic residues and their unique packing in the protein core. To find the origin of the uniqueness in their amino acid sequences, we analyzed the distribution of the side chain rotational isomers (rotamers) of hydrophobic amino acids in protein tertiary structures and derived deltaS(contact), the conformational-entropy changes of side chains by residue-residue contacts in each secondary structure. The deltaS(contact) values indicate distinct tendencies of the residue pairs to restrict side chain conformation by inter-residue contacts. Of the hydrophobic residues in alpha-helices, aliphatic residues (Leu, Val, Ile) strongly restrict the side chain conformations of each other. In beta-sheets, Met is most strongly restricted by contact with Ile, whereas Leu, Val and Ile are less affected by other residues in contact than those in alpha-helices. In designed and native protein variants, deltaS(contact) was found to correlate with the folding-unfolding cooperativity. Thus, it can be used as a specificity parameter for designing artificial proteins with a unique structure.     

Expression and function of the gonadotropin-releasing hormone receptor are dependent on a conserved apolar amino acid in the third intracellular loop

The coupling of agonist-activated heptahelical receptors to their cognate G proteins is often dependent on the amino-terminal region of the third intracellular loop. Like many G protein-coupled receptors, the gonadotropin-releasing hormone (GnRH) receptor contains an apolar amino acid in this region at a constant distance from conserved Pro and Tyr/Asn residues in the fifth transmembrane domain (TM V). An analysis of the role of this conserved residue (Leu(237)) in GnRH receptor function revealed that the binding affinities of the L237I and L237V mutant receptors were unchanged, but their abilities to mediate GnRH-induced inositol phosphate signaling, G protein coupling, and agonist-induced internalization were significantly impaired. Receptor expression at the cell surface was reduced by replacement of Leu(237) with Val, and abolished by replacement with Ala, Arg, or Asp residues. These results are consistent with molecular modeling of the TM V and VI regions of the GnRH receptor, which predicts that Leu(237) is caged by several apolar amino acids (Ile(233), Ile(234), and Val(240) in TM V, and Leu(262), Leu(265), and Val(269) in TM VI) to form a tight hydrophobic cluster. These findings indicate that the conserved apolar residue (Leu(237)) in the third intracellular loop is an important determinant of GnRH receptor expression and activation, and possibly that of other G protein-coupled receptors.     

Binary patterning of polar and nonpolar amino acids in the sequences and structures of native proteins

Protein sequences can be represented as binary patterns of polar (○) and nonpolar (?) amino acids. These binary sequence patterns are categorized into two classes: Class A patterns match the structural repeat of an idealized amphiphilic α-helix (3.6 residues per turn), and class B patterns match the structural repeat of an idealized amphiphilic β-strand (2 residues per turn). The difference between these two classes of sequence patterns has led to a strategy for de novo protein design based on binary patterning of polar and nonpolar amino acids. Here we ask whether similar binary patterning is incorporated in the sequences and structures of natural proteins. Analysis of the Protein Data Bank demonstrates the following. (1) Class A sequence patterns occur considerably more frequently in the sequences of natural proteins than would be expected at random, but class B patterns occur less often than expected. (2) Each pattern is found predominantly in the secondary structure expected from the binary strategy for protein design. Thus, class A patterns are found more frequently in α-helices than in β-strands, and class B patterns are found more frequently in β-strands than in α-helices. (3) Among the α-helices of natural proteins, the most commonly used binary patterns are indeed the class A patterns. (4) Among all β-strands in the database, the most commonly used binary patterns are not the expected class B patterns. (5) However, for solvent-exposed β-strands, the correlation is striking: All β-strands in the database that contain the class B patterns are exposed to solvent. (6) The bias of class A patterns for α-structure over β-structure and the bias of class B patterns for β-structure over α-structure are significant, not merely when compared to other binary patterns of polar (○) and nonpolar (?) amino acids, but also when compared to the full range of sequences in the database. The implications for the design of novel proteins are discussed.     

Packing and hydrophobicity effects on protein folding and stability: effects of beta-branched amino acids, valine and isoleucine, on the formation and stability of two-stranded alpha-helical coiled coils/leucine zippers.

The aim of this study was to examine the differences between hydrophobicity and packing effects in specifying the three-dimensional structure and stability of proteins when mutating hydrophobes in the hydrophobic core. In DNA-binding proteins (leucine zippers), Leu residues are conserved at positions "d," and beta-branched amino acids, Ile and Val, often occur at positions "a" in the hydrophobic core. In order to discern what effect this selective distribution of hydrophobes has on the formation and stability of two-stranded alpha-helical coiled coils/leucine zippers, three Val or three Ile residues were simultaneously substituted for Leu at either positions "a" (9, 16, and 23) or "d" (12, 19, and 26) in both chains of a model coiled coil. The stability of the resulting coiled coils was monitored by CD in the presence of Gdn.HCl. The results of the mutations of Ile to Val at either positions "a" or "d" in the reduced or oxidized coiled coils showed a significant hydrophobic effect with the additional methylene group in Ile stabilizing the coiled coil (delta delta G values range from 0.45 to 0.88 kcal/mol/mutation). The results of mutations of Leu to Ile or Val at positions "a" in the reduced or oxidized coiled coils showed a significant packing effect in stabilizing the coiled coil (delta delta G values range from 0.59 to 1.03 kcal/mol/mutation). Our results also indicate the subtle control hydrophobic packing can have not only on protein stability but on the conformation adopted by the amphipathic alpha-helices. These structural findings correlate with the observation that in DNA-binding proteins, the conserved Leu residues at positions "d" are generally less tolerant of amino acid substitutions than the hydrophobic residues at positions "a."     

Structural and mechanistic investigations on Salmonella typhimurium acetate kinase (AckA): identification of a putative ligand binding pocket at the dimeric interface

Background

A large number of studies have been carried out to obtain amino acid propensities for ??-helices and ??-sheets. The obtained propensities for ??-helices are consistent with each other, and the pair-wise correlation coefficient is frequently high. On the other hand, the ??-sheet propensities obtained by several studies differed significantly, indicating that the context significantly affects ??-sheet propensity.

Results

We calculated amino acid propensities for ??-helices and ??-sheets for 39 and 24 protein folds, respectively, and addressed whether they correlate with the fold. The propensities were also calculated for exposed and buried sites, respectively. Results showed that ??-helix propensities do not differ significantly by fold, but ??-sheet propensities are diverse and depend on the fold. The propensities calculated for exposed sites and buried sites are similar for ??-helix, but such is not the case for the ??-sheet propensities. We also found some fold dependence on amino acid frequency in ??-strands. Folds with a high Ser, Thr and Asn content at exposed sites in ??-strands tend to have a low Leu, Ile, Glu, Lys and Arg content (correlation coefficient = ?0.90) and to have flat ??-sheets. At buried sites in ??-strands, the content of Tyr, Trp, Gln and Ser correlates negatively with the content of Val, Ile and Leu (correlation coefficient = ?0.93). "All-??" proteins tend to have a higher content of Tyr, Trp, Gln and Ser, whereas "??/??" proteins tend to have a higher content of Val, Ile and Leu.

Conclusions

The ??-helix propensities are similar for all folds and for exposed and buried residues. However, ??-sheet propensities calculated for exposed residues differ from those for buried residues, indicating that the exposed-residue fraction is one of the major factors governing amino acid composition in ??-strands. Furthermore, the correlations we detected suggest that amino acid composition is related to folding properties such as the twist of a ??-strand or association between two ?? sheets.     

Follicle-stimulating hormone interacts with exoloop 3 of the receptor

The human follicle-stimulating hormone (FSH) receptor consists of two distinct domains of approximately 330 amino acids, the N-terminal extracellular exodomain and membrane-associated endodomain including three exoloops and seven transmembrane helices. The exodomain binds the hormone with high affinity, and the resulting hormone/exodomain complex modulates the endodomain where receptor activation occurs. It has been an enigma whether the hormone interacts with the endodomain. In a step to address the question, exoloop 3 of (580)KVPLITVSKAK(590) was examined by Ala scan, multiple substitution, assays for hormone binding, cAMP and inositol phosphate (IP) induction, and photoaffinity labeling. We present the evidence for the interaction of FSH and exoloop 3. A peptide mimic of exoloop 3 specifically and saturably photoaffinity-labels FSH alpha but not FSH beta. This is in contrast to photoaffinity labeling of FSH beta by the peptide mimic of the N-terminal region of the receptor. Leu(583) and Ile(584) are crucial for the interaction of FSH and exoloop 3. Substitutions of these two residues enhanced the hormone binding affinity. This is due to the loss of the original side chains but not the introduction of new side chains. The Leu(583) and Ile(584) side chains appear to project in opposite directions. Ile(584) appears to be so specific and to require flexibility and stereo specificity so that no other amino acids can fit into its place. Leu(583) is less specific. The improvement in hormone binding by substitutions was offset by the severe impairment of signal generation of cAMP and/or inositol phosphate. For example, the Phe or Tyr substitution of Leu(583) improved the hormone binding and cAMP induction but impaired IP induction. On the other hand, the substitutions for Ile(584) and Lys(590) abolished the cAMP and IP induction. Our results open a logical question whether Leu(583), Ile(584), and Lys(590) interact with the exodomain and/or the hormone. The answers will provide new insights into the mechanisms of hormone binding and signal generation.     

Functional role of the extracellular N-terminal domain of neuropeptide Y subfamily receptors in membrane integration and agonist-stimulated internalization

The N terminus is the most variable element in G protein-coupled receptors (GPCRs), ranging from seven residues up to approximately 5900 residues. For family B and C GPCRs it is described that at least part of the ligand binding site is located within the N terminus. Here we investigated the role of the N terminus in the neuropeptide Y receptor family, which belongs to the class A of GPCRs. We cloned differentially truncated Y receptor mutants, in which the N terminus was partially or completely deleted. We found, that eight amino acids are sufficient for full ligand binding and signal transduction activity. Interestingly, we could show that no specific amino acids but rather the extension of the first transmembrane helix by any residues is sufficient for receptor activity but also for membrane integration in case of the hY(1) and the hY(4) receptors. In contrast, the complete deletion of the N terminus in the hY(2) receptors resulted in a mutant that is fully integrated in the membrane but does not bind the ligand very well and internalizes much slower compared to the wild type receptor. Interestingly, also these effects could be reverted by any N-terminal extension. Accordingly, the most important function of the N termini seems to be the stabilization of the first transmembrane helix to ensure the correct receptor structure, which obviously is essential for ligand binding, integration into the cell membrane and receptor internalization.     

Structural and functional roles of small group-conserved amino acids present on helix-H7 in the β(2)-adrenergic receptor

Sequence analysis of the class A G protein-coupled receptors (GPCRs) reveals that most of the highly conserved sites are located in the transmembrane helices. A second level of conservation exists involving those residues that are conserved as a group characterized by small and/or weakly polar side chains (Ala, Gly, Ser, Cys, Thr). These positions can have group conservation levels of up to 99% across the class A GPCRs and have been implicated in mediating helix-helix interactions in membrane proteins. We have previously shown that mutation of group-conserved residues present on transmembrane helices H2-H4 in the β(2)-adrenergic receptor (β(2)-AR) can influence both receptor expression and function. We now target the group-conserved sites, Gly315(7.42) and Ser319(7.46), on H7 for structure-function analysis. Replacing Ser319(7.46) with smaller amino acids (Ala or Gly) did not influence the ability of the mutant receptors to bind to the antagonist dihydroalprenolol (DHA) but resulted in ~15-20% agonist-independent activity. Replacement of Ser319(7.46) with the larger amino acid leucine lowered the expression of the S319L mutant and its ability to bind DHA. Both the G315A and G315S mutants also exhibited agonist-independent signaling, while the G315L mutant did not show specific binding to DHA. These data indicate that Gly315(7.42) and Ser319(7.46) are stabilizing β(2)-AR in an inactive conformation. We discuss our results in the context of van der Waals interactions of Gly315(7.42) with Trp286(6.48) and hydrogen bonding interactions of Ser319(7.46) with amino acids on H1-H2-H7 and with structural water.     

Hydrophobicity and prediction of the secondary structure of membrane proteins and peptides.

Reliability of the hydropathy method to predict the formation of membrane-spanning alpha-helices by integral membrane proteins and peptides whose structure is known from X-ray crystallography is analysed. It is shown that Kyte-Doolittle hydropathy plots do not predict accurately 22 transmembrane alpha-helices in the reaction centres (RC) of the photosynthetic bacteria Rhodopseudomonas viridis and Rhodobacter sphaeroides (R-26). The accuracy of prediction for these proteins was improved using an optimised Kyte-Doolittle hydrophobicity scale. However, this hydrophobicity scale did not improve the predictions for the alphabeta-peptides of the B800-850 (LH2) complexes of the photosynthetic bacteria Rhodopseudomonas acidophila and Rhodospirillum molischianum, which were excluded from the optimisation procedure. The best and worst predictions of membrane-spanning alpha-helices for the RC proteins and LH2 peptides, respectively, were obtained with a propensity scale (PRC) calculated from the amino acid sequences and X-ray data for the RC proteins. A propensity scale (PLH) obtained using the amino acid sequences and X-ray data for the alphabeta-peptides of the LH2 complexes did not give an acceptable prediction of the transmembrane segments in the LH2 peptides; moreover, it markedly contradicted the PRC scale. Amino acids have been concluded to have no significant preference to localisation in transmembrane segments. Therefore, the predictive ability of the hydropathy methodology appears to be limited: the number of transmembrane segments can be correctly calculated for the best case only, and the lengths and positions of membrane-spanning alpha-helices in a protein amino acid sequence can not be predicted exactly.