Comparing function and structure between entire proteomes

More than 30 organisms have been sequenced entirely. Here, we applied a variety of simple bioinformatics tools to analyze 29 proteomes for representatives from all three kingdoms: eukaryotes, prokaryotes, and archaebacteria. We confirmed that eukaryotes have relatively more long proteins than prokaryotes and archaes, and that the overall amino acid composition is similar among the three. We predicted that approximately 15%-30% of all proteins contained transmembrane helices. We could not find a correlation between the content of membrane proteins and the complexity of the organism. In particular, we did not find significantly higher percentages of helical membrane proteins in eukaryotes than in prokaryotes or archae. However, we found more proteins with seven transmembrane helices in eukaryotes and more with six and 12 transmembrane helices in prokaryotes. We found twice as many coiled-coil proteins in eukaryotes (10%) as in prokaryotes and archaes (4%-5%), and we predicted approximately 15%-25% of all proteins to be secreted by most eukaryotes and prokaryotes. Every tenth protein had no known homolog in current databases, and 30%-40% of the proteins fell into structural families with >100 members. A classification by cellular function verified that eukaryotes have a higher proportion of proteins for communication with the environment. Finally, we found at least one homolog of experimentally known structure for approximately 20%-45% of all proteins; the regions with structural homology covered 20%-30% of all residues. These numbers may or may not suggest that there are 1200-2600 folds in the universe of protein structures. All predictions are available at http://cubic.bioc.columbia.edu/genomes.     

Large Tilts in Transmembrane Helices Can Be Induced during Tertiary Structure Formation

While early structural models of helix-bundle integral membrane proteins posited that the transmembrane α-helices [transmembrane helices (TMHs)] were orientated more or less perpendicular to the membrane plane, there is now ample evidence from high-resolution structures that many TMHs have significant tilt angles relative to the membrane. Here, we address the question whether the tilt is an intrinsic property of the TMH in question or if it is imparted on the TMH during folding of the protein. Using a glycosylation mapping technique, we show that four highly tilted helices found in multi-spanning membrane proteins all have much shorter membrane-embedded segments when inserted by themselves into the membrane than seen in the high-resolution structures. This suggests that tilting can be induced by tertiary packing interactions within the protein, subsequent to the initial membrane-insertion step.     

A predictor of transmembrane α-helix domains of proteins based on neural networks

Back-propagation, feed-forward neural networks are used to predict a-helical transmembrane segments of proteins. The networks are trained on the few membrane proteins whose transmembrane -helix domains are known to atomic or nearly atomic resolution. When testing is performed with a jackknife procedure on the proteins of the training set, the fraction of total correct assignments is as high as 0.87, with an average length for the transmembrane segments of 20 residues. The method correctly fails to predict any transmembrane domain for porin, whose transmembrane segments are -sheets. When tested on globular proteins, lower and upper limits of 1.6 and 3.5% for a total of 26826 residues are determined for the mispredicted cases, indicating that the predictor is highly specific for -helical domains of membrane proteins. The predictor is also tested on 37 membrane proteins whose transmembrane topology is partially known. The overall accuracy is 0.90, two percentage points higher than that obtained with statistical methods. The reliability of the prediction is 100% for 60% of the total 18242 predicted residues of membrane proteins. Our results show that the local directional information automatically extracted by the neural networks during the training phase plays a key role in determining the accuracy of the prediction. Correspondence to: R. Casadio     

Analysis of the membrane organization of an Escherichia coli protein translocator, HlyB, a member of a large family of prokaryote and eukaryote surface transport proteins

Haemolysin B (HlyB) is essential for secretion of the 107 x 10(3) Mr haemolysin A protein from Escherichia coli and is a member of a family of highly conserved, apparently ATP-dependent surface proteins in many organisms. We have shown in this study that both HlyB and HlyD fractionate primarily with the cytoplasmic membrane of E. coli and are accessible to proteases after removal of the outer membrane. We have measured experimentally the topological organization of HlyB within the membrane by construction of fusions to beta-lactamase as a reporter. The predicted folding of HlyB, with a minimum of six transmembrane segments, does not always coincide with regions of highest average hydrophobicity. This suggests that HlyB may have a novel organization within the bilayer. From our data and comparative sequence analysis, we have been able to predict very similar topological models for the other members of the HlyB family.     

Single-spanning transmembrane domains in cell growth and cell-cell interactions: More than meets the eye?

As a whole, integral membrane proteins represent about one third of sequenced genomes, and more than 50% of currently available drugs target membrane proteins, often cell surface receptors. Some membrane protein classes, with a defined number of transmembrane (TM) helices, are receiving much attention because of their great functional and pharmacological importance, such as G protein-coupled receptors possessing 7 TM segments. Although they represent roughly half of all membrane proteins, bitopic proteins (with only 1 TM helix) have so far been less well characterized. Though they include many essential families of receptors, such as adhesion molecules and receptor tyrosine kinases, many of which are excellent targets for biopharmaceuticals (peptides, antibodies, et al.). A growing body of evidence suggests a major role for interactions between TM domains of these receptors in signaling, through homo and heteromeric associations, conformational changes, assembly of signaling platforms, etc. Significantly, mutations within single domains are frequent in human disease, such as cancer or developmental disorders. This review attempts to give an overview of current knowledge about these interactions, from structural data to therapeutic perspectives, focusing on bitopic proteins involved in cell signaling.Key words: bitopic membrane proteins, transmembrane domains, transmembrane signaling, helix-helix interactions, receptors     

Estimating intrinsic structural preferences of de novo emerging random-sequence proteins: Is aggregation the main bottleneck?

Present-day proteins are believed to have evolved features to reduce the risk of aggregation. However, proteins can emerge de novo by translation of non-coding DNA segments. In this study we assess the aggregation, disorder and transmembrane propensity of protein sequences generated by translating random nucleotide sequences of varying GC-content. Potential de novo random-sequence proteins translated from regions with GC content between 40% and 60% do not show stronger aggregation propensity than existing ones and exhibit similar tendency to be disordered. We suggest that de novo emerging proteins do not mean an unavoidable aggregation threat to evolving organisms.     

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.     

NMR structures of polytopic integral membrane proteins

Membrane proteins represent up to 30% of the proteins in all organisms, they are involved in many biological processes and are the molecular targets for around 50% of validated drugs. Despite this, membrane proteins represent less than 1% of all high-resolution protein structures due to various challenges associated with applying the main biophysical techniques used for protein structure determination. Recent years have seen an explosion in the number of high-resolution structures of membrane proteins determined by NMR spectroscopy, especially for those with multiple transmembrane-spanning segments. This is a review of the structures of polytopic integral membrane proteins determined by NMR spectroscopy up to the end of the year 2010, which includes both β-barrel and α-helical proteins from a number of different organisms and with a range in types of function. It also considers the challenges associated with performing structural studies by NMR spectroscopy on membrane proteins and how some of these have been overcome, along with its exciting potential for contributing new knowledge about the molecular mechanisms of membrane proteins, their roles in human disease, and for assisting drug design.     

Role of medium--and long-range interactions in discriminating globular and membrane proteins.

The analysis of inter-residue interactions in protein structures provides considerable insight to understand their folding and stability. We have previously analyzed the role of medium- and long-range interactions in the folding of globular proteins. In this work, we study the distinct role of such interactions in the three-dimensional structures of membrane proteins. We observed a higher number of long-range contacts in the termini of transmembrane helical (TMH) segments, implying their role in the stabilization of helix-helix interactions. The transmembrane strand (TMS) proteins are having appreciably higher long-range contacts than that in all-beta class of globular proteins, indicating closer packing of the strands in TMS proteins. The residues in membrane spanning segments of TMH proteins have 1.3 times higher medium-range contacts than long-range contacts whereas that of TMS proteins have 14 times higher long-range contacts than medium-range contacts. Residue-wise analysis indicates that in TMH proteins, the residues Cys, Glu, Gly, Pro, Gln, Ser and Tyr have higher long-range contacts than medium-range contacts in contrast with all-alpha class of globular proteins. The charged residue pairs have higher medium-range contacts in all-alpha proteins, whereas hydrophobic residue pairs are dominant in TMH proteins. The information on the preference of residue pairs to form medium-range contacts has been successfully used to discriminate the TMH proteins from all-alpha proteins. The statistical significance of the results obtained from the present study has been verified using randomized structures of TMH and TMS protein templates.     

Protein Disorder and Short Conserved Motifs in Disordered Regions Are Enriched near the Cytoplasmic Side of Single-Pass Transmembrane Proteins

Intracellular juxtamembrane regions of transmembrane proteins play pivotal roles in cell signalling, mediated by protein-protein interactions. Disordered protein regions, and short conserved motifs within them, are emerging as key determinants of many such interactions. Here, we investigated whether disorder and conserved motifs are enriched in the juxtamembrane area of human single-pass transmembrane proteins. Conserved motifs were defined as short disordered regions that were much more conserved than the adjacent disordered residues. Human single-pass proteins had higher mean disorder in their cytoplasmic segments than their extracellular parts. Some, but not all, of this effect reflected the shorter length of the cytoplasmic tail. A peak of cytoplasmic disorder was seen at around 30 residues from the membrane. We noted a significant increase in the incidence of conserved motifs within the disordered regions at the same location, even after correcting for the extent of disorder. We conclude that elevated disorder within the cytoplasmic tail of many transmembrane proteins is likely to be associated with enrichment for signalling interactions mediated by conserved short motifs.     

Recognition of alpha-helix transmembrane domains with an amphipathy scale generated by molecular dynamics using only the primary sequence of proteins

We recently developed an amphipathy scale, elaborated from molecular dynamics data that can be used for the identification of hydrophobic or hydrophilic regions in proteins. This amphipathy scale reflects side chain/water molecule interaction energies. We have now used this amphipathy scale to find candidates for transmembrane segments, by examining a large sample of membrane proteins with alpha-helix segments. The candidates were selected based on an amphipathy coefficient value range and the minimum number of residues in a segment. We compared our results with the transmembrane segments previously identified in the PDB_TM database by the TMDET algorithm. We expected that the hydrophobic segments would be identified using only the primary structures of the proteins and the amphipathy scale. However, some of these hydrophobic segments may pertain to hydrophobic pockets not included in transmembrane regions. We found that our amphipathy scale could identify alpha-helix transmembrane regions with a probability of success of 76% when all segments were included and 90% when all membrane proteins were included.     

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.     

Size comparisons among integral membrane transport protein homologues in bacteria, Archaea, and Eucarya

Integral membrane proteins from over 20 ubiquitous families of channels, secondary carriers, and primary active transporters were analyzed for average size differences between homologues from the three domains of life: Bacteria, Archaea, and Eucarya. The results showed that while eucaryotic homologues are consistently larger than their bacterial counterparts, archaeal homologues are significantly smaller. These size differences proved to be due primarily to variations in the sizes of hydrophilic domains localized to the N termini, the C termini, or specific loops between transmembrane alpha-helical spanners, depending on the family. Within the Eucarya domain, plant homologues proved to be substantially smaller than their animal and fungal counterparts. By contrast, extracytoplasmic receptors of ABC-type uptake systems in Archaea proved to be larger on average than those of their bacterial homologues, while cytoplasmic enzymes from different organisms exhibited little or no significant size differences. These observations presumably reflect evolutionary pressure and molecular mechanisms that must have been operative since these groups of organisms diverged from each other.     

NMR structures of polytopic integral membrane proteins

Abstract

Membrane proteins represent up to 30% of the proteins in all organisms, they are involved in many biological processes and are the molecular targets for around 50% of validated drugs. Despite this, membrane proteins represent less than 1% of all high-resolution protein structures due to various challenges associated with applying the main biophysical techniques used for protein structure determination. Recent years have seen an explosion in the number of high-resolution structures of membrane proteins determined by NMR spectroscopy, especially for those with multiple transmembrane-spanning segments. This is a review of the structures of polytopic integral membrane proteins determined by NMR spectroscopy up to the end of the year 2010, which includes both β-barrel and α-helical proteins from a number of different organisms and with a range in types of function. It also considers the challenges associated with performing structural studies by NMR spectroscopy on membrane proteins and how some of these have been overcome, along with its exciting potential for contributing new knowledge about the molecular mechanisms of membrane proteins, their roles in human disease, and for assisting drug design.     

Plant P-type ATPases

P-type ATPases are a superfamily of membrane proteins involved in many physiological processes that are fundamental for all living organisms. Using ATP, they can transport a variety of ions and other substances across all types of cell membranes against a concentration electrochemical gradient. P-type ATPases form a phosphorylated intermediate and are sensitive to vanadate. Based on evolutionary relations and sequence homology, P-type ATPases are divided into five major families. All P-type ATPases share a simple structure and mechanism, but also possess domains characteristic for each family, which are crucial for substrate specificity. These proteins usually have a single subunit with eight to twelve transmembrane segments, a large central cytoplasmic domain with the conservative ATP binding site along with N and C termini exposed to the cytoplasm. Because of variety of proteins that belong to P-type ATPase superfamily, in this review the comparison of functional and structure properties of plant cells P-type ATPases is presented, as well as their important role in adaptation to environmental stress.     

Membrane-integration Characteristics of Two ABC Transporters, CFTR and P-glycoprotein

To what extent do corresponding transmembrane helices in related integral membrane proteins have different membrane-insertion characteristics? Here, we compare, side-by-side, the membrane insertion characteristics of the 12 transmembrane helices in the adenosine triphosphate-binding cassette (ABC) transporters, P-glycoprotein (P-gp) and the cystic fibrosis transmembrane conductance regulator (CFTR). Our results show that 10 of the 12 CFTR transmembrane segments can insert independently into the ER membrane. In contrast, only three of the P-gp transmembrane segments are independently stable in the membrane, while the majority depend on the presence of neighboring loops and/or transmembrane segments for efficient insertion. Membrane-insertion characteristics can thus vary widely between related proteins.     

Sequence similarity as a predictor of the transmembrane topology of membrane-intrinsic subunits of bacterial respiratory chain enzymes

Integral membrane proteins usually have a predominantly alpha-helical secondary structure in which transmembrane segments are connected by membrane-extrinsic loops. Although a number of membrane protein structures have been reported in recent years, in most cases transmembrane topologies are initially predicted using a variety of theoretical techniques, including hydropathy analyses and the "positive inside" rule. We have explored the use of plots of the distribution of sequence similarity within families of membrane proteins comprising homeomorphic domains as a new method for the prediction/verification of the orientation of transmembrane topology models within certain families of multimeric respiratory chain enzymes. Within such proteins, analyses of sequence similarity can: i) identify heme and/or quinol binding sites; ii) identify potential electron-transfer conduits to/from prosthetic groups; and iii) locate regions defining potential subunit-subunit interactions. We mined emerging bioinformatic data for sequences of 11 families of membrane-intrinsic proteins that are part of multimeric respiratory chain complexes that also have membrane-extrinsic subunits. The sequences of each family were then aligned and the resultant alignments converted into a graphical format recording an empirical measure of the sequence similarity plotted versus residue position. In each case, this plot was compared to the predicted transmembrane topology. With one exception, there is a strong correlation between the existence     

An artificial transmembrane segment directs SecA, SecB, and electrochemical potential-dependent translocation of a long amino-terminal tail

Many integral membrane proteins contain an amino-terminal segment, often referred to as an N-tail, that is translocated across a membrane. In many cases, translocation of the N-tail is initiated by a cleavable, amino-terminal signal peptide. For N-tail proteins lacking a signal peptide, translocation is initiated by a transmembrane segment that is carboxyl to the translocated segment. The mechanism of membrane translocation of these segments, although poorly understood, has been reported to be independent of the protein secretion machinery. In contrast, here we describe alkaline phosphatase mutants containing artificial transmembrane segments that demonstrate that translocation of a long N-tail across the membrane is dependent upon SecA, SecB, and the electrochemical potential in the absence of a signal peptide. The corresponding mutants containing signal peptides also use the secretion machinery but are less sensitive to inhibition of its components. We present evidence that inhibition of SecA by sodium azide is incomplete even at high concentrations of inhibitor, which suggests why SecA-dependent translocation may not have been detected in other systems. Furthermore, by varying the charge around the transmembrane segment, we find that in the absence of a signal peptide, the orientation of the membrane-bound alkaline phosphatase is dictated by the positive inside rule. However, the presence of a signal peptide is an overriding factor in membrane orientation and renders all mutants in an Nout-Cin orientation.     

Characteristics affecting expression and solubilization of yeast membrane proteins

Biochemical and structural analysis of membrane proteins often critically depends on the ability to overexpress and solubilize them. To identify properties of eukaryotic membrane proteins that may be predictive of successful overexpression, we analyzed expression levels of the genomic complement of over 1000 predicted membrane proteins in a recently completed Saccharomyces cerevisiae protein expression library. We detected statistically significant positive and negative correlations between high membrane protein expression and protein properties such as size, overall hydrophobicity, number of transmembrane helices, and amino acid composition of transmembrane segments. Although expression levels of membrane and soluble proteins exhibited similar negative correlations with overall hydrophobicity, high-level membrane protein expression was positively correlated with the hydrophobicity of predicted transmembrane segments. To further characterize yeast membrane proteins as potential targets for structure determination, we tested the solubility of 122 of the highest expressed yeast membrane proteins in six commonly used detergents. Almost all the proteins tested could be solubilized using a small number of detergents. Solubility in some detergents depended on protein size, number of transmembrane segments, and hydrophobicity of predicted transmembrane segments. These results suggest that bioinformatic approaches may be capable of identifying membrane proteins that are most amenable to overexpression and detergent solubilization for structural and biochemical analyses. Bioinformatic approaches could also be used in the redesign of proteins that are not intrinsically well-adapted to such studies.     

Stability of loops in the structure of lactose permease

Structural analysis of peptide fragments has provided useful information on the secondary structure of integral membrane proteins built from a helical bundle (up to seven transmembrane segments). Comparison of those results to recent X-ray crystallographic results showed agreement between the structures of the fragments and the structures of the intact proteins. Lactose permease of Escherichia coli (lac Y) offers an opportunity to test that hypothesis on a substantially larger integral membrane protein. Lac Y contains a bundle of 12 transmembrane segments connected by 11 loops. Eleven segments, each corresponding to one of the loops in this protein, were studied. Five of these segments form defined structures in solution as determined by multidimensional nuclear magnetic resonance. Four peptides form turns, and one peptide reveals the end of one of the transmembrane helices. These results suggest that some loops in helical bundles are stabilized by short-range interactions, particularly in smaller bundles, and such intrinsically stable loops may contribute to protein stability and influence the pathway of folding. Greater conformational flexibility may be found in large integral membrane proteins.