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.     

Functional estimation of loop-helix boundaries in the lactose permease of Escherichia coli by single amino acid deletion analysis

Mutants with single amino acid deletions in the loops of lactose permease retain activity, while mutants with single deletions in transmembrane helices are inactive, and the loop--helix boundaries of helices IV, V, VII, VIII, and IX have been approximated functionally by the systematic deletion of single residues [Wolin, C. D., and Kaback, H. R. (1999) Biochemistry 38, 8590-8597]. The experimental approach is applied here to the remainder of the permease. Periplasmic and cytoplasmic loop-helix boundaries for helices I, II, X, XI, and XII and the cytoplasmic boundary of helix III are in reasonable agreement with structural predictions. In contrast, the periplasmic end of helix III appears to be five to eight residues further into the transmembrane domain than predicted. Taken together with the previous findings, the analysis estimates that 11 of the 12 transmembrane helices have an average length of 21 residues. Surprisingly, deletion analysis of loop V/VI, helix VI, and loop VI/VII does not yield an activity profile typical of the rest of the protein, as individual deletion of only three residues in this region abolishes activity. Thus, transmembrane domain VI which is probably on the periphery of the 12-helix bundle may make few functionally important contacts.     

Proportion of membrane proteins in proteomes of 15 single-cell organisms analyzed by the SOSUI prediction system

A software system, SOSUI, was previously developed for discriminating between soluble and membrane proteins and predicting transmembrane regions (Hirokawa et al., Bioinformatics, 14 (1998) 378-379). The performance of the system was 99% for the discrimination between two types of proteins and 96% for the prediction of transmembrane helices. When all of the amino acid sequences from 15 single-cell organisms were analyzed by SOSUI, the proportion of predicted polytopic membrane proteins showed an almost constant value of 15-20%, irrespective of the total genome size. However, single-cell organisms appeared to be categorized in terms of the preference of the number of transmembrane segments: species with small genomes were characterized by a significant peak at a helix number of approximately six or seven; species with large genomes showed a peak at 10 or 11 helices; and species with intermediate genome sizes showed a monotonous decrease of the population of membrane proteins against the number of transmembrane helices.     

Stabilization of secondary structure elements by specific combinations of hydrophilic and hydrophobic amino acid residues is more important for proteins encoded by GC-poor genes

Stabilization of secondary structure elements by specific combinations of hydrophobic and hydrophilic amino acids has been studied by the way of analysis of pentapeptide fragments from twelve partial bacterial proteomes. PDB files describing structures of proteins from species with extremely high and low genomic GC-content, as well as with average G + C were included in the study. Amino acid residues in 78,009 pentapeptides from alpha helices, beta strands and coil regions were classified into hydrophobic and hydrophilic ones. The common propensity scale for 32 possible combinations of hydrophobic and hydrophilic amino acid residues in pentapeptide has been created: specific pentapeptides for helix, sheet and coil were described. The usage of pentapeptides preferably forming alpha helices is decreasing in alpha helices of partial bacterial proteomes with the increase of the average genomic GC-content in first and second codon positions. The usage of pentapeptides preferably forming beta strands is increasing in coil regions and in helices of partial bacterial proteomes with the growth of the average genomic GC-content in first and second codon positions. Due to these circumstances the probability of coil-sheet and helix-sheet transitions should be increased in proteins encoded by GC-rich genes making them prone to form amyloid in certain conditions. Possible causes of the described fact that importance of alpha helix and coil stabilization by specific combinations of hydrophobic and hydrophilic amino acids is growing with the decrease of genomic GC-content have been discussed.     

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.     

Structures of the intradiskal loops and amino terminus of the G-protein receptor, rhodopsin.

The intradiskal surface of the transmembrane protein, rhodopsin, consists of the amino terminal domain and three loops connecting six of the seven transmembrane helices. This surface corresponds to the extracellular surface of other G-protein receptors. Peptides that represent each of the extramembraneous domains on this surface (three loops and the amino terminus) were synthesized. These peptides also included residues which, based on a hydrophobic plot, could be expected to be part of the transmembrane helix. The structure of each of these peptides in solution was then determined using two-dimensional 1H nuclear magnetic resonance. All peptide domains showed ordered structures in solution. The structures of each of the peptides from intradiskal loops of rhodopsin exhibited a turn in the central region of the peptide. The ends of the peptides show an unwinding of the transmembrane helices to form this turn. The amino terminal domain peptide exhibited alpha-helical regions with breaks and bends at proline residues. This region forms a compact domain. Together, the structures for the loop and amino terminus domains indicate that the intradiskal surface of rhodopsin is ordered. These data further suggest a structural motif for short loops in transmembrane proteins. The ordered structures of these loops, in the absence of the transmembrane helices, indicate that the primary sequences of these loops are sufficient to code for the turn.     

Defining the structural basis for assembly of a transmembrane cytochrome

To define the structural basis for cofactor binding to membrane proteins, we introduce a manageable model system, which allows us, for the first time, to study the influence of individual transmembrane helices and of single amino acid residues on the assembly of a transmembrane cytochrome. In vivo as well as in vitro analyses indicate central roles of single amino acid residues for either interaction of the transmembrane helices or for binding of the cofactor. The results clearly show that interaction of the PsbF transmembrane helix is independent from binding of the heme cofactor. On the other hand, binding of the cofactor highly depends on helix-helix interactions. By site-directed mutagenesis critical amino acid residues were identified, which are involved in the assembly of a functional transmembrane cytochrome. Especially, a highly conserved glycine residue is critical for interaction of the transmembrane helices and assembly of the cytochrome. Based on the two-stage-model of alpha-helical membrane protein folding, the presented results clearly indicate a third stage of membrane protein folding, in which a cofactor binds to a pre-assembled transmembrane protein.     

Proline-induced distortions of transmembrane helices

Proline residues in the transmembrane (TM) alpha-helices of integral membrane proteins have long been suspected to play a key role for helix packing and signal transduction by inducing regions of helix distortion and/or dynamic flexibility (hinges). In this study we try to characterise the effect of proline on the geometric properties of TM alpha-helices. We have examined 199 transmembrane alpha-helices from polytopic membrane proteins of known structure. After examining the location of proline residues within the amino acid sequences of TM helices, we estimated the helix axes either side of a hinge and hence identified a hinge residue. This enabled us to calculate helix kink and swivel angles. The results of this analysis show that proline residues occur with a significant concentration in the centre of sequences of TM alpha-helices. In this location, they may induce formation of molecular hinges, located on average about four residues N-terminal to the proline residue. A superposition of proline-containing TM helices structures shows that the distortion induced is anisotropic and favours certain relative orientations (defined by helix kink and swivel angles) of the two helix segments.     

Folding amphipathic helices into membranes: amphiphilicity trumps hydrophobicity

High amphiphilicity is a hallmark of interfacial helices in membrane proteins and membrane-active peptides, such as toxins and antimicrobial peptides. Although there is general agreement that amphiphilicity is important for membrane-interface binding, an unanswered question is its importance relative to simple hydrophobicity-driven partitioning. We have examined this fundamental question using measurements of the interfacial partitioning of a family of 17-residue amidated-acetylated peptides into both neutral and anionic lipid vesicles. Composed only of Ala, Leu, and Gln residues, the amino acid sequences of the peptides were varied to change peptide amphiphilicity without changing total hydrophobicity. We found that peptide helicity in water and interface increased linearly with hydrophobic moment, as did the favorable peptide partitioning free energy. This observation provides simple tools for designing amphipathic helical peptides. Finally, our results show that helical amphiphilicity is far more important for interfacial binding than simple hydrophobicity.     

Thermostability of membrane protein helix-helix interaction elucidated by statistical analysis

A prerequisite for the survival of (micro)organisms at high temperatures is an adaptation of protein stability to extreme environmental conditions. In contrast to soluble proteins, where many factors have already been identified, the mechanisms by which the thermostability of membrane proteins is enhanced are almost unknown. The hydrophobic membrane environment constrains possible stabilizing factors for transmembrane domains, so that a difference might be expected between soluble and membrane proteins. Here we present sequence analysis of predicted transmembrane helices of the genomes from eight thermophilic and 12 mesophilic organisms. A comparison of the amino acid compositions indicates that more polar residues can be found in the transmembrane helices of thermophilic organisms. Particularly, the amino acids aspartic acid and glutamic acid replace the corresponding amides. Cysteine residues are found to be significantly decreased by about 70% in thermophilic membrane domains suggesting a non-specific function of most cysteine residues in transmembrane domains of mesophilic organisms. By a pair-motif analysis of the two sets of transmembrane helices, we found that the small residues glycine and serine contribute more to transmembrane helix-helix interactions in thermophilic organisms. This may result in a tighter packing of the helices allowing more hydrogen bond formation.     

Computational differentiation of N-terminal signal peptides and transmembrane helices

Signal peptides and transmembrane helices both contain a stretch of hydrophobic amino acids. This common feature makes it difficult for signal peptide and transmembrane helix predictors to correctly assign identity to stretches of hydrophobic residues near the N-terminal methionine of a protein sequence. The inability to reliably distinguish between N-terminal transmembrane helix and signal peptide is an error with serious consequences for the prediction of protein secretory status or transmembrane topology. In this study, we report a new method for differentiating protein N-terminal signal peptides and transmembrane helices. Based on the sequence features extracted from hydrophobic regions (amino acid frequency, hydrophobicity, and the start position), we set up discriminant functions and examined them on non-redundant datasets with jackknife tests. This method can incorporate other signal peptide prediction methods and achieve higher prediction accuracy. For Gram-negative bacterial proteins, 95.7% of N-terminal signal peptides and transmembrane helices can be correctly predicted (coefficient 0.90). Given a sensitivity of 90%, transmembrane helices can be identified from signal peptides with a precision of 99% (coefficient 0.92). For eukaryotic proteins, 94.2% of N-terminal signal peptides and transmembrane helices can be correctly predicted with coefficient 0.83. Given a sensitivity of 90%, transmembrane helices can be identified from signal peptides with a precision of 87% (coefficient 0.85). The method can be used to complement current transmembrane protein prediction and signal peptide prediction methods to improve their prediction accuracies.     

A single amino acid substitution in rhodopsin (lysine 248----leucine) prevents activation of transducin

In structure-function studies on bovine rhodopsin by in vitro site-specific mutagenesis, we have prepared three mutants in the cytoplasmic loop between the putative transmembrane helices E and F. In each mutant, charged amino acid residues were replaced by neutral residues: mutant 1, Glu239----Gln; mutant 2, Lys248----Leu; and mutant 3, Glu247----Gln, Lys248----Leu, and Glu249----Gln. The mutant rhodopsin genes were expressed in monkey kidney (COS-1) cells. After the addition of 11-cis-retinal to the cells, the rhodopsin mutants were purified by immunoaffinity adsorption. Each mutant gave a wild-type rhodopsin visible absorption spectrum. The mutants were assayed for their ability to stimulate the GTPase activity of transducin in a light-dependent manner. While mutants 1 and 3 showed wild-type activity, mutant 2 (Lys248----Leu) was inactive.     

Structural classification and prediction of reentrant regions in alpha-helical transmembrane proteins: application to complete genomes

Alongside the well-studied membrane spanning helices, alpha-helical transmembrane (TM) proteins contain several functionally and structurally important types of substructures. Here, existing 3D structures of transmembrane proteins have been used to define and study the concept of reentrant regions, i.e. membrane penetrating regions that enter and exit the membrane on the same side. We find that these regions can be divided into three distinct categories based on secondary structure motifs, namely long regions with a helix-coil-helix motif, regions of medium length with the structure helix-coil or coil-helix and regions of short to medium length consisting entirely of irregular secondary structure. The residues situated in reentrant regions are significantly smaller on average compared to other regions and reentrant regions can be detected in the inter-transmembrane loops with an accuracy of approximately 70% based on their amino acid composition. Using TOP-MOD, a novel method for predicting reentrant regions, we have scanned the genomes of Escherichia coli, Saccharomyces cerevisiae and Homo sapiens. The results suggest that more than 10% of transmembrane proteins contain reentrant regions and that the occurrence of reentrant regions increases linearly with the number of transmembrane regions. Reentrant regions seem to be most commonly found in channel proteins and least commonly in signal receptors.     

The sensor protein KdpD inserts into the Escherichia coli membrane independent of the Sec translocase and YidC.

KdpD is a sensor kinase protein in the inner membrane of Escherichia coli containing four transmembrane regions. The periplasmic loops connecting the transmembrane regions are intriguingly short and protease mapping allowed us to only follow the translocation of the second periplasmic loop. The results show that neither the Sec translocase nor the YidC protein are required for membrane insertion of the second loop of KdpD. To study the translocation of the first periplasmic loop a short HA epitope tag was genetically introduced into this region. The results show that also the first loop was translocated independently of YidC and the Sec translocase. We conclude that KdpD resembles a new class of membrane proteins that insert into the membrane without enzymatic assistance by the known translocases. When the second periplasmic loop was extended by an epitope tag to 27 amino acid residues, the membrane insertion of this loop of KdpD depended on SecE and YidC. To test whether the two periplasmic regions are translocated independently of each other, the KdpD protein was split between helix 2 and 3 into two approximately equal-sized fragments. Both constructed fragments, which contained KdpD-N (residues 1-448 of KdpD) and the KdpD-C (residues 444-894 of KdpD), readily inserted into the membrane. Similar to the epitope-tagged KdpD protein, only KdpD-C depended on the presence of the Sec translocase and YidC. This confirms that the four transmembrane helices of KdpD are inserted pairwise, each translocation event involving two transmembrane helices and a periplasmic loop.     

Interactions between charged amino acid residues within transmembrane helices in the sulfate transporter SHST1

The aim of this study was to identify charged amino acid residues important for activity of the sulfate transporter SHST1. We mutated 10 charged amino acids in or near proposed transmembrane helices and expressed the resulting mutants in a sulfate transport-deficient yeast strain. Mutations affecting four residues resulted in a complete loss of sulfate transport; these residues were D107 and D122 in helix 1 and R354 and E366 in helix 8. All other mutants showed some reduction in transport activity. The E366Q mutant was unusual in that expression of the mutant protein was toxic to yeast cells. The R354Q mutant showed reduced trafficking to the plasma membrane, indicating that the protein was misfolded. However, transporter function (to a low level) and wild-type trafficking could be recovered by combining the R354Q mutation with either the E175Q or E270Q mutations. This suggested that R354 interacts with both E175 and E270. The triple mutant E175Q/E270Q/R354Q retained only marginal sulfate transport activity but was trafficked at wild-type levels, suggesting that a charge network between these three residues may be involved in the transport pathway, rather than in folding. D107 was also found to be essential for the ion transport pathway and may form a charge pair with R154, both of which are highly conserved. The information obtained on interactions between charged residues provides the first evidence for the possible spatial arrangement of transmembrane helices within any member of this transporter family. This information is used to develop a model for SHST1 tertiary structure.     

Genetic mapping of the human C5a receptor. Identification of transmembrane amino acids critical for receptor function

Many hormones and sensory stimuli signal through a superfamily of seven transmembrane-spanning receptors to activate heterotrimeric G proteins. How the seven transmembrane segments of the receptors (a molecular architecture of bundled alpha-helices conserved from yeast to man) work as "on/off" switches remains unknown. Previously, we used random saturation mutagenesis coupled with a genetic selection in yeast to determine the relative importance of amino acids in four of the seven transmembrane segments of the human C5a receptor (Baranski, T. J., Herzmark, P., Lichtarge, O., Gerber, B. O., Trueheart, J., Meng, E. C., Iiri, T., Sheikh, S. P., and Bourne, H. R. (1999) J. Biol. Chem. 274, 15757-15765). In this study, we evaluate helices I, II, and IV, thereby furnishing a complete mutational map of the seven transmembrane helices of the human C5a receptor. Our analysis identified 19 amino acid positions resistant to non-conservative substitutions. When combined with the 25 essential residues previously identified in helices III and V-VII, they delineate two distinct components of the receptor switch: a ligand-binding surface at or near the extracellular surface of the helix bundle and a core cluster in the cytoplasmic half of the bundle. In addition, we found critical amino acids in the first and second helices that are predicted to face the lipid membrane. These residues form an extended surface that might mediate interactions with lipids and other membrane proteins or function as an oligomerization domain with other receptors.     

Random coil structures in bacterial proteins. Relationships of their amino acid compositions to flanking structures and corresponding genic base compositions

In this study we classified regions of random coil into four types: coil between alpha helix and beta strand, coil between beta strand and alpha helix, coil between two alpha helices and coil between two beta strands. This classification may be considered as natural. We used 610 3D structures of proteins collected from the Protein Data Bank from bacteria with low, average and high genomic GC-content. Relatively short regions of coil are not random: certain amino acid residues are more or less frequent in each of the types of coil. Namely, hydrophobic amino acids with branched side chains (Ile, Val and Leu) are rare in coil between two beta strands, unlike some acrophilic amino acids (Asp, Asn and Gly). In contrast, coil between two alpha helices is enriched by Leu. Regions of coil between alpha helix and beta strand are enriched by positively charged amino acids (Arg and Lys), while the usage of residues with side chains possessing hydroxyl group (Ser and Thr) is low in them, in contrast to the regions of coil between beta strand and alpha helix. Regions of coil between beta strand and alpha helix are significantly enriched by Cys residues. The response to the symmetric mutational pressure (AT-pressure or GC-pressure) is also quite different for four types of coil. The most conserved regions of coil are “connecting bridges” between beta strand and alpha helix, since their amino acid content shows less strong dependence on GC-content of genes than amino acid contents of other three types of coil. Possible causes and consequences of the described differences in amino acid content distribution between different types of random coil have been discussed.     

High resolution NMR analysis of the seven transmembrane domains of a heptahelical receptor in organic-aqueous medium

The NMR properties of seven peptides representing the transmembrane domains of the alpha-factor receptor from Saccharomyces cerevisiae were examined in trifluoroethanol/water (4:1) at 10 to 55 degrees C. The parameters extracted indicated all peptides were helical in this membrane mimetic solvent. Using chemical shift indices as the criterion, helicity varied from 64 to 83%. The helical residues in the peptides corresponded to the region predicted to cross the hydrocarbon interior of the bilayer. A study of a truncated 25-residue peptide corresponding to domain 2 gave evidence that the helix extended all the way to the N-terminus of this peptide, indicating that sequence and not chain end effects are very important in helix termination for our model peptides. Both nuclear Overhauser effect spectroscopy (NOESY) connectivities and chemical shift indices revealed significant perturbations around prolyl residues in the helices formed by transmembrane domains 6 and 7. Molecular models of the transmembrane domains indicate that helices for domains 6 and 7 are severely kinked at these prolyl residues. The helix perturbation around proline 258 in transmembrane domain 6 correlates with mutations that cause phenotypic changes in this receptor.     

A study of the membrane-water interface region of membrane proteins

The most conspicuous structural characteristic of the alpha-helical membrane proteins is their long transmembrane alpha-helices. However, other structural elements, as yet largely ignored in statistical studies of membrane protein structure, are found in those parts of the protein that are located in the membrane-water interface region. Here, we show that this region is enriched in irregular structure and in interfacial helices running roughly parallel with the membrane surface, while beta-strands are extremely rare. The average amino acid composition is different between the interfacial helices, the parts of the transmembrane helices located in the interface region, and the irregular structures. In this region, hydrophobic and aromatic residues tend to point toward the membrane and charged/polar residues tend to point away from the membrane. The interface region thus imposes different constraints on protein structure than do the central hydrocarbon core of the membrane and the surrounding aqueous phase.     

TM-MOTIF: an alignment viewer to annotate predicted transmembrane helices and conserved motifs in aligned set of sequences

Multiple sequence alignments become biologically meaningful only if conserved and functionally important residues and secondary structural elements preserved can be identified at equivalent positions. This is particularly important for transmembrane proteins like G-protein coupled receptors (GPCRs) with seven transmembrane helices. TM-MOTIF is a software package and an effective alignment viewer to identify and display conserved motifs and amino acid substitutions (AAS) at each position of the aligned set of homologous sequences of GPCRs. The key feature of the package is to display the predicted membrane topology for seven transmembrane helices in seven colours (VIBGYOR colouring scheme) and to map the identified motifs on its respective helices /loop regions. It is an interactive package which provides options to the user to submit query or pre-aligned set of GPCR sequences to align with a reference sequence, like rhodopsin, whose structure has been solved experimentally. It also provides the possibility to identify the nearest homologue from the available inbuilt GPCR or Olfactory Receptor cluster dataset whose association is already known for its receptor type. AVAILABILITY: The database is available for free at mini@ncbs.res.in.