Solution NMR resonance assignment strategies for β‐barrel membrane proteins

Membrane proteins in detergent micelles are large and dynamic complexes that present challenges for solution NMR investigations such as spectral overlap and line broadening. In this study, multiple methods are introduced to facilitate resonance assignment of β‐barrel membrane proteins using Opa₆₀ from Neisseria gonorrhoeae as a model system. Opa₆₀ is an eight‐stranded β‐barrel with long extracellular loops (～63% of the protein) that engage host receptors and induce engulfment of the bacterium. The NMR spectra of Opa₆₀ in detergent micelles exhibits significant spectral overlap and resonances corresponding to the loop regions had variable line widths, which interfered with a complete assignment of the protein. To assign the β‐barrel residues, trypsin cleavage was used to remove much of the extracellular loops while preserving the detergent solubilized β‐barrel. The removal of the loop resonances significantly improved the assignment of the Opa₆₀ β‐barrel region (97% of the resonances corresponding to the β‐barrel and periplasmic turns were assigned). For the loop resonance assignments, two strategies were implemented; modulating temperature and synthetic peptides. Lowering the temperature broadened many peaks beyond detection and simplified the spectra to only the most dynamic regions of the loops facilitating 27 loop resonances to be assigned. To further assign functionally important and unstructured regions of the extracellular loops, a synthetic 20 amino acid peptide was synthesized and had nearly complete spectral overlap with the full‐length protein allowing 17 loop resonances to be assigned. Collectively, these strategies are effective tools that may accelerate solution NMR structure determination of β‐barrel membrane proteins.     

Preformed structural elements feature in partner recognition by intrinsically unstructured proteins

Intrinsically unstructured proteins (IUPs) are devoid of extensive structural order but often display signs of local and limited residual structure. To explain their effective functioning, we reasoned that such residual structure can be crucial in their interactions with their structured partner(s) in a way that preformed structural elements presage their final conformational state. To check this assumption, a database of 24 IUPs with known 3D structures in the bound state has been assembled and the distribution of secondary structure elements and backbone torsion angles have been analysed. The high proportion of residues in coil conformation and with phi, psi angles in the disallowed regions of the Ramachandran map compared to the reference set of globular proteins shows that IUPs are not fully ordered even in their bound form. To probe the effect of partner proteins on IUP folding, inherent conformational preferences of IUP sequences have been assessed by secondary structure predictions using the GOR, ALB and PROF algorithms. The accuracy of predicting secondary structure elements of IUPs is similar to that of their partner proteins and is significantly higher than the corresponding values for random sequences. We propose that strong conformational preferences mark regions in IUPs (mostly helices), which correspond to their final structural state, while regions with weak conformational preferences represent flexible linkers between them. In our interpretation, preformed elements could serve as initial contact points, the binding of which facilitates the reeling of the flexible regions onto the template. This finding implies that IUPs draw a functional advantage from preformed structural elements, as they enable their facile, kinetically and energetically less demanding, interaction with their physiological partner.     

On residues in the disallowed region of the Ramachandran map.

An analysis of the occurrence of nonglycyl residues in conformations disallowed in the Ramachandran plot is presented. Ser, Asn, Thr, and Cys have the highest propensities to exhibit such conformations, and the branched aliphatic residues the lowest. Residues cluster in five regions and there are some trends in the types of residues and their side-chain conformations (chi(1)) occupying these. Majority of the residues are found at the edge of helices and strands and in short loops, and are involved in different types of weak, stabilizing interactions. A structural motif has been identified where a residue in disallowed conformation occurs as the first residue of a short 3(10)-helix. On the basis of the types of neighboring residues, the location in the three-dimensional structure and accessibility, there are similarities with the occurrence of cis peptide bonds in protein structures.     

The Effect of Loops on the Structural Organization of α-Helical Membrane Proteins

Loops connecting the transmembrane (TM) α-helices in membrane proteins are expected to affect the structural organization of the thereby connected helices and the helical bundles as a whole. This effect, which has been largely ignored previously, is studied here by analyzing the x-ray structures of 41 α-helical membrane proteins. First we define the loop flexibility ratio, R, and find that 53% of the loops are stretched, where a stretched loop constrains the distance between the two connected helices. The significance of this constraining effect is supported by experiments carried out with bacteriorhodopsin and rhodopsin, in which cutting or eliminating their (predominately stretched) loops has led to a decrease in protein stability, and for rhodopsin, in most cases, also to the destruction of the structure. We show that for nonstretched loops in the extramembranous regions, the fraction of hydrophobic residues is comparable to that for soluble proteins; furthermore (as is also the case for soluble proteins), the hydrophobic residues in these regions are preferentially buried. This is expected to lead to the compact structural organization of the loops, which is transferred to the TM helices, causing them to assemble. We argue that a soluble protein complexed with a membrane protein similarly promotes compactness; other properties of such complexes are also studied. We calculate complementary attractive interactions between helices, including hydrogen bonds and van der Waals interactions of sequential motifs, such as GXXXG. The relative and combined effects of all these factors on the association of the TM helices are discussed and protein structures with only a few of these factors are analyzed. Our study emphasizes the need for classifying membrane proteins into groups according to structural organization. This classification should be considered when procedures for structural analysis or prediction are developed and applied. Detailed analysis of each structure is provided at http://flan.blm.cs.cmu.edu/memloop/     

The dual role of a loop with low loop contact distance in folding and domain swapping

Alpha helices, beta strands, and loops are the basic building blocks of protein structure. The folding kinetics of alpha helices and beta strands have been investigated extensively. However, little is known about the formation of loops. Experimental studies show that for some proteins, the formation of a single loop is the rate-determining step for folding, whereas for others, a loop (or turn) can misfold to serve as the hinge loop region for domain-swapped species. Computer simulations of an all-atom model of fragment B of Staphylococcal protein A found that the formation of a single loop initiates the dominant folding pathway. On the other hand, the stability analysis of intermediates suggests that the same loop is a likely candidate to serve as a hinge loop for domain swapping. To interpret the simulation result, we developed a simple structural parameter: the loop contact distance (LCD), or the sequence distance of contacting residues between a loop and the rest of the protein. The parameter is applied to a number of other proteins, including SH3 domains and prion protein. The results suggest that a locally interacting loop (low LCD) can either promote folding or serve as the hinge region for domain swapping. Thus, there is an intimate connection between folding and domain swapping, a possible cause of misfolding and aggregation.     

Dynamic alpha-helices: conformations that do not conform

Structural transitions are important for the stability and function of proteins, but these phenomena are poorly understood. An extensive analysis of Protein Data Bank entries reveals 103 regions in proteins with a tendency to transform from helical to nonhelical conformation and vice versa. We find that these dynamic helices, unlike other helices, are depleted in hydrophobic residues. Furthermore, the dynamic helices have higher surface accessibility and conformational mobility (P-value = 3.35e-07) than the rigid helices. Contact analyses show that these transitions result from protein-ligand, protein-nucleic acid, and crystal-contacts. The immediate structural environment differs quantitatively (P-value = 0.003) as well as qualitatively in the two alternate conformations. Often, dynamic helix experiences more contacts in its helical conformation than in the nonhelical counterpart (P-value = 0.001). There is differential preference for the type of short contacts observed in two conformational states. We also demonstrate that the regions in protein that can undergo such large conformational transitions can be predicted with a reasonable accuracy using logistic regression model of supervised learning. Our findings have implications in understanding the molecular basis of structural transitions that are coupled with binding and are important for the function and stability of the protein. Based on our observations, we propose that several functionally relevant regions on the protein surface can switch over their conformation from coil to helix and vice-versa, to regulate the recognition and binding of their partner and hence these may work as "molecular switches" in the proteins to regulate certain biological process. Our results supports the idea that protein structure-function paradigm should transform from static to a highly dynamic one.     

Structure and dynamics of a membrane protein in micelles from three solution NMR experiments

Three solution NMR experiments on a uniformly ¹⁵N labeled membrane protein in micelles provide sufficient information to describe the structure, topology, and dynamics of its helices, as well as additional information that characterizes the principal features of residues in terminal and inter-helical loop regions. The backbone amide resonances are assigned with an HMQC-NOESY experiment and the backbone dynamics are characterized by a ¹H-¹⁵N heteronuclear NOE experiment, which clearly distinguishes between the structured helical residues and the more mobile residues in the terminal and interhelical loop regions of the protein. The structure and topology of the helices are described by Dipolar waves and PISA wheels derived from experimental measurements of residual dipolar couplings (RDCs) and residual chemical shift anisotropies (RCSAs). The results show that the membrane-bound form of Pf1 coat protein has a 20-residue trans-membrane hydrophobic helix with an orientation that differs by about 90° from that of an 8-residue amphipathic helix. This combination of three-experiments that yields Dipolar waves and PISA wheels has the potential to contribute to high-throughput structural characterizations of membrane proteins.     

NMR structure and binding studies confirm that PA4608 from Pseudomonas aeruginosa is a PilZ domain and a c-di-GMP binding protein

PA4608 is a 125 residue protein from Pseudomonas aeruginosa with a recent identification as a PilZ domain and putative bis-(3'-5')-cyclic dimeric guanosine monophosphate (c-di-GMP) adaptor protein that plays a role in bacterial second-messenger regulated processes. The nuclear magnetic resonance (NMR) structure of PA4608 has been determined and c-di-GMP binding has been confirmed by NMR titration studies. The monomeric core structure of PA4608 contains a six-stranded anti-parallel beta barrel flanked by three helices. Conserved surface residues among PA4608 homologs suggest the c-di-GMP binding site is at one end of the barrel and includes residues in the helices as well as in the unstructured N-terminus. Chemical shift changes in PA4608 resonances upon titration with c-di-GMP confirm binding. This evidence supports the hypothesis that proteins containing PilZ domains are the long-sought c-di-GMP adaptor proteins.     

Modeling of protein loops by simulated annealing.

A method is presented to model loops of protein to be used in homology modeling of proteins. This method employs the ESAP program of Higo et al. (Higo, J., Collura, V., & Garnier, J., 1992, Biopolymers 32, 33-43) and is based on a fast Monte Carlo simulation and a simulated annealing algorithm. The method is tested on different loops or peptide segments from immunoglobulin, bovine pancreatic trypsin inhibitor, and bovine trypsin. The predicted structure is obtained from the ensemble average of the coordinates of the Monte Carlo simulation at 300 K, which exhibits the lowest internal energy. The starting conformation of the loop prior to modeling is chosen to be completely extended, and a closing harmonic potential is applied to N, CA, C, and O atoms of the terminal residues. A rigid geometry potential of Robson and Platt (1986, J. Mol. Biol. 188, 259-281) with a united atom representation is used. This we demonstrate to yield a loop structure with good hydrogen bonding and torsion angles in the allowed regions of the Ramachandran map. The average accuracy of the modeling evaluated on the eight modeled loops is 1 A root mean square deviation (rmsd) for the backbone atoms and 2.3 A rmsd for all heavy atoms.     

Long-range effects of tag sequence on marginally stabilized structure in HIV-1 p24 capsid protein monitored using NMR

N-terminal domain of HIV-1 p24 capsid protein is a globular fold composed of seven helices and two β-strands with a flexible structure including the α4–5 loop and both N- and C-terminal ends. However, the protein shows a high tendency (48%) for an intrinsically disordered structure based on the PONDR VL-XT prediction from the primary sequence. To assess the possibility of marginally stabilized structure under physiological conditions, the N-terminal domain of p24 was destabilized by the addition of an artificial flexible tag to either N- or C-terminal ends, and it was analyzed using T₁, T₂, hetero-nuclear NOE, and amide-proton exchange experiments. When the C-terminal tag (12 residues) was attached, the regions of the α3–4 loop and helix 6 as well as the α4–5 loop attained the flexible structures. Furthermore, in the protein containing the N-terminal tag (27 residues), helix 4 in addition to the above-mentioned area including α3–4 and α4–5 loops as well as helix 6 exhibited highly disordered structures. Thus, the long-range effects of the existence of tag sequence was observed in the stepwise manner of the appearance of disordered structures (step 1: α4–5 loop, step 2: α3–4 loop and helix 6, and step 3: helix 4). Furthermore, the disordered regions in tagged proteins were consistent with the PONDR VL-XT disordered prediction. The dynamic structure located in the middle part (α3–4 loop to helix 6) of the protein shown in this study may be related to the assembly of the viral particle.     

Functional Analysis of the Autographa californica Multiple Nucleopolyhedrovirus GP64 Terminal Fusion Loops and Interactions with Membranes

The Autographa californica multiple nucleopolyhedrovirus (AcMNPV) glycoprotein GP64 is the major envelope protein of the budded virus (BV). GP64 is a class III fusion protein that mediates BV attachment to the cell surface and low-pH-triggered membrane fusion between the BV envelope and the endosome membrane during entry. Class III fusion proteins contain terminal looped structures that are believed to interact with membranes. To examine the functions of 3 loops found at the apex of the GP64 postfusion structure, we generated 2-alanine substitutions that scanned the two so-called fusion loops (loop 1 and loop 2) plus an adjacent loop structure (loop 3) that is closely attached to loop 2 and is also found at the apex of the GP64 postfusion structure. We identified essential residues from Y75 to T86 (loop 1) and N149 to H156 (loop 2) that are required for fusion activity, but no essential residues in loop 3. Further analysis revealed that critical fusion loop residues fall within two groups that are associated with either membrane merger (hemifusion) or fusion pore expansion. We next examined the interactions of soluble GP64 proteins and BV with membranes composed of various phospholipids. BV interacted directly with small unilamellar vesicles (SUVs) comprised of phospholipids phosphatidylcholine and phosphatidic acid (PC/PA) or phosphatidylcholine and phosphatidylserine (PC/PS) under neutral and acidic pH. We also examined the interactions of soluble GP64 constructs containing substitutions of the most hydrophobic residues within each of the two fusion loops. We found that a 2-residue substitution in either single loop (loop 1 [positions 81 and 82] or loop 2 [positions 153 and 154]) was not sufficient to substantially reduce the GP64-liposome interaction, but the same substitutions in both fusion loops severely reduced the GP64-liposome association at neutral pH. These results suggest that critical hydrophobic residues in both fusion loops may be involved in the interaction of GP64 with host cellular membranes and direct GP64-membrane interactions may represent a receptor-binding step prior to a low-pH-triggered conformational change.     

Calcium- and membrane-induced changes in the structure and dynamics of three helical hairpins in annexin B12

Annexins are a family of soluble proteins that can undergo reversible Ca(2+)-dependent interaction with the interfacial region of phospholipid membranes. The helical hairpins on the convex face of the crystal structure of soluble annexins are proposed to mediate binding to membranes, but the mechanism is not defined. For this study, we used a site-directed spin labeling (SDSL) experimental approach to investigate Ca(2+) and membrane-induced structural and dynamic changes that occurred in the helical hairpins encompassing three of the four D and E helices of annexin B12. Electron paramagnetic resonance (EPR) parameters were analyzed for the soluble and Ca(2+)-dependent membrane-bound states of the following nitroxide scans of annexin B12: a continuous 24-residue scan of the D and E helices in the third repeat (residues 219-242) and short scans encompassing the D-E loop regions of the first repeat (residues 68-74) and the fourth repeat (300-305). EPR mobility and accessibility parameters of most sites were similar when the protein was in solution or in the membrane-bound state, and both sets of data were consistent with the crystal structure of the protein. However, membrane-induced changes in mobility and accessibility were observed in all three loop regions, with the most dramatic changes noted at sites corresponding to the highly conserved serine and glycine residues in the loops. EPR accessibility parameters clearly established that nitroxide side chains placed at these sites made direct contact with the bilayer. EPR mobility parameters showed that these sites were very mobile in solution, but immobilized on the EPR time scale in the membrane-bound state. Since the headgroup regions of bilayer phospholipids are relatively mobile in the absence of annexins, Ca(2+)-dependent binding of annexin B12 appears to form a complex in which the mobility of the D-E loop region of the protein and the headgroup region of the phospholipid are highly constrained. Possible biological consequences of annexin-induced restriction of membrane mobility are discussed.     

The NMR structure of a stable and compact all-beta-sheet variant of intestinal fatty acid-binding protein

Intestinal fatty acid-binding protein (I-FABP) has a clam-shaped structure that may serve as a scaffold for the design of artificial enzymes and drug carriers. In an attempt to optimize the scaffold for increased access to the interior-binding cavity, several helix-less variants of I-FABP have been engineered. The solution-state NMR structure of the first generation helix-less variant, known as Delta17-SG, revealed a larger-than-expected and structurally ill-defined loop flanking the deletion site. We hypothesized that the presence of this loop, on balance, was energetically unfavorable for the stability of the protein. The structure exhibited no favorable pairwise or nonpolar interactions in the loop that could offset the loss of configurational entropy associated with the folding of this region of the protein. As an attempt to generate a more stable protein, we engineered a second-generation helix-less variant of I-FABP (Delta27-GG) by deleting 27 contiguous residues of the wild-type protein and replacing them with a G-G linker. The deletion site of this variant (D9 through N35) includes the 10 residues spanning the unstructured loop of Delta17-SG. Chemical denaturation experiments using steady-state fluorescence spectroscopy showed that the second-generation helix-less variant is energetically more stable than Delta17-SG. The three-dimensional structure of apo-Delta27-GG was solved using triple-resonance NMR spectroscopy along with the structure calculation and refinement protocols contained in the program package ARIA/CNS. In spite of the deletion of 27 residues, the structure assumes a compact all-beta-sheet fold with no unstructured loops and open access to the interior cavity.     

Conformational analysis and clustering of short and medium size loops connecting regular secondary structures: a database for modeling and prediction.

Loops are regions of nonrepetitive conformation connecting regular secondary structures. We identified 2,024 loops of one to eight residues in length, with acceptable main-chain bond lengths and peptide bond angles, from a database of 223 protein and protein-domain structures. Each loop is characterized by its sequence, main-chain conformation, and relative disposition of its bounding secondary structures as described by the separation between the tips of their axes and the angle between them. Loops, grouped according to their length and type of their bounding secondary structures, were superposed and clustered into 161 conformational classes, corresponding to 63% of all loops. Of these, 109 (51% of the loops) were populated by at least four nonhomologous loops or four loops sharing a low sequence identity. Another 52 classes, including 12% of the loops, were populated by at least three loops of low sequence similarity from three or fewer nonhomologous groups. Loop class suprafamilies resulting from variations in the termini of secondary structures are discussed in this article. Most previously described loop conformations were found among the classes. New classes included a 2:4 type IV hairpin, a helix-capping loop, and a loop that mediates dinucleotide-binding. The relative disposition of bounding secondary structures varies among loop classes, with some classes such as beta-hairpins being very restrictive. For each class, sequence preferences as key residues were identified; those most frequently at these conserved positions than in proteins were Gly, Asp, Pro, Phe, and Cys. Most of these residues are involved in stabilizing loop conformation, often through a positive phi conformation or secondary structure capping. Identification of helix-capping residues and beta-breakers among the highly conserved positions supported our decision to group loops according to their bounding secondary structures. Several of the identified loop classes were associated with specific functions, and all of the member loops had the same function; key residues were conserved for this purpose, as is the case for the parvalbumin-like calcium-binding loops. A significant number, but not all, of the member loops of other loop classes had the same function, as is the case for the helix-turn-helix DNA-binding loops. This article provides a systematic and coherent conformational classification of loops, covering a broad range of lengths and all four combinations of bounding secondary structure types, and supplies a useful basis for modelling of loop conformations where the bounding secondary structures are known or reliably predicted.     

Clustering of protein structural fragments reveals modular building block approach of nature

Structures of peptide fragments drawn from a protein can potentially occupy a vast conformational continuum. We co-ordinatize this conformational space with the help of geometric invariants and demonstrate that the peptide conformations of the currently available protein structures are heavily biased in favor of a finite number of conformational types or structural building blocks. This is achieved by representing a peptides' backbone structure with geometric invariants and then clustering peptides based on closeness of the geometric invariants. This results in 12,903 clusters, of which 2207 are made up of peptides drawn from functionally and/or structurally related proteins. These are termed "functional" clusters and provide clues about potential functional sites. The rest of the clusters, including the largest few, are made up of peptides drawn from unrelated proteins and are termed "structural" clusters. The largest clusters are of regular secondary structures such as helices and beta strands as well as of beta hairpins. Several categories of helices and strands are discovered based on geometric differences. In addition to the known classes of loops, we discover several new classes, which will be useful in protein structure modeling. Our algorithm does not require assignment of secondary structure and, therefore, overcomes the limitations in loop classification due to ambiguity in secondary structure assignment at loop boundaries.     

Ribosome binding of a single copy of the SecY complex: implications for protein translocation

The SecY complex associates with the ribosome to form a protein translocation channel in the bacterial plasma membrane. We have used cryo-electron microscopy and quantitative mass spectrometry to show that a nontranslating E. coli ribosome binds to a single SecY complex. The crystal structure of an archaeal SecY complex was then docked into the electron density maps. In the resulting model, two cytoplasmic loops of SecY extend into the exit tunnel near proteins L23, L29, and L24. The loop between transmembrane helices 8 and 9 interacts with helices H59 and H50 in the large subunit RNA, while the 6/7 loop interacts with H7. We also show that point mutations of basic residues within either loop abolish ribosome binding. We suggest that SecY binds to this primary site on the ribosome and subsequently captures and translocates the nascent chain.     

A novel method for protein-protein interaction site prediction using phylogenetic substitution models

Protein-protein binding events mediate many critical biological functions in the cell. Typically, functionally important sites in proteins can be well identified by considering sequence conservation. However, protein-protein interaction sites exhibit higher sequence variation than other functional regions, such as catalytic sites of enzymes. Consequently, the mutational behavior leading to weak sequence conservation poses significant challenges to the protein-protein interaction site prediction. Here, we present a phylogenetic framework to capture critical sequence variations that favor the selection of residues essential for protein-protein binding. Through the comprehensive analysis of diverse protein families, we show that protein binding interfaces exhibit distinct amino acid substitution as compared with other surface residues. On the basis of this analysis, we have developed a novel method, BindML, which utilizes the substitution models to predict protein-protein binding sites of protein with unknown interacting partners. BindML estimates the likelihood that a phylogenetic tree of a local surface region in a query protein structure follows the substitution patterns of protein binding interface and nonbinding surfaces. BindML is shown to perform well compared to alternative methods for protein binding interface prediction. The methodology developed in this study is very versatile in the sense that it can be generally applied for predicting other types of functional sites, such as DNA, RNA, and membrane binding sites in proteins.     

Bulged residues promote the progression of a loop-loop interaction to a stable and inhibitory antisense-target RNA complex

In several groups of bacterial plasmids, antisense RNAs regulate copy number through inhibition of replication initiator protein synthesis. These RNAs are characterized by a long hairpin structure interrupted by several unpaired residues or bulged loops. In plasmid R1, the inhibitory complex between the antisense RNA (CopA) and its target mRNA (CopT) is characterized by a four-way junction structure and a side-by-side helical alignment. This topology facilitates the formation of a stabilizer intermolecular helix between distal regions of both RNAs, essential for in vivo control. The bulged residues in CopA/CopT were shown to be required for high in vitro binding rate and in vivo activity. This study addresses the question of why removal of bulged nucleotides blocks stable complex formation. Structure mapping, modification interference, and molecular modeling of bulged-less mutant CopA–CopT complexes suggests that, subsequent to loop–loop contact, helix propagation is prevented. Instead, a fully base paired loop–loop interaction is formed, inducing a continuous stacking of three helices. Consequently, the stabilizer helix cannot be formed, and stable complex formation is blocked. In contrast to the four-way junction topology, the loop–loop interaction alone failed to prevent ribosome binding at its loading site and, thus, inhibition of RepA translation was alleviated.     

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.