waveTM: wavelet-based transmembrane segment prediction

waveTM is a web tool for the prediction of transmembrane segments in alpha-helical membrane proteins. Prediction is performed by a dynamic programming algorithm on wavelet-denoised 'hydropathy' signals. Users submit a protein sequence and receive interactively the results. Topology prediction can also be obtained in conjunction with the algorithm OrienTM. A web server that implements the waveTM algorithm is freely available at http://bioinformatics.biol.uoa.gr/waveTM.     

STAM: simple transmembrane alignment method

MOTIVATION: The database of transmembrane protein (TMP) structures is still very small. At the same time, more and more TMP sequences are being determined. Molecular modeling is an interim answer that may bridge the gap between the two databases. The first step in homology modeling is to achieve a good alignment between the target sequences and the template structure. However, since most algorithms to obtain the alignments were constructed with data derived from globular proteins, they perform poorly when applied to TMPs. In our application, we automate the alignment procedure and design it specifically for TMP. We first identify segments likely to form transmembrane alpha-helices. We then apply different sets of criteria for transmembrane and non-transmembrane segments. For example, the penalty for insertion/deletions in the transmembrane segments is much higher than that of a penalty in the loop region. Different substitution matrices are used since the frequencies of occurrence of the various amino acids differ for transmembrane segments and water-soluble domains. RESULTS: This program leads to better models since it does not treat the protein as a single entity with the same properties, but accounts for the different physical properties of the various segments. STAM is the first multisequence alignment program that is directly targeted at transmembrane proteins. AVAILABILITY: Source code and installation package are available on request from the authors. Web access is currently implemented.     

On filtering false positive transmembrane protein predictions

While helical transmembrane (TM) region prediction tools achieve high (>90%) success rates for real integral membrane proteins, they produce a considerable number of false positive hits in sequences of known nontransmembrane queries. We propose a modification of the dense alignment surface (DAS) method that achieves a substantial decrease in the false positive error rate. Essentially, a sequence that includes possible transmembrane regions is compared in a second step with TM segments in a sequence library of documented transmembrane proteins. If the performance of the query sequence against the library of documented TM segment-containing sequences in this test is lower than an empirical threshold, it is classified as a non-transmembrane protein. The probability of false positive prediction for trusted TM region hits is expressed in terms of E-values. The modified DAS method, the DAS-TMfilter algorithm, has an unchanged high sensitivity for TM segments ( approximately 95% detected in a learning set of 128 documented transmembrane proteins). At the same time, the selectivity measured over a non-redundant set of 526 soluble proteins with known 3D structure is approximately 99%, mainly because a large number of falsely predicted single membrane-pass proteins are eliminated by the DAS-TMfilter algorithm.     

Genome-wide analysis of integral membrane proteins from eubacterial,archaean, and eukaryotic organisms.

We have carried out detailed statistical analyses of integral membrane proteins of the helix-bundle class from eubacterial, archaean, and eukaryotic organisms for which genome-wide sequence data are available. Twenty to 30% of all ORFs are predicted to encode membrane proteins, with the larger genomes containing a higher fraction than the smaller ones. Although there is a general tendency that proteins with a smaller number of transmembrane segments are more prevalent than those with many, uni-cellular organisms appear to prefer proteins with 6 and 12 transmembrane segments, whereas Caenorhabditis elegans and Homo sapiens have a slight preference for proteins with seven transmembrane segments. In all organisms, there is a tendency that membrane proteins either have many transmembrane segments with short connecting loops or few transmembrane segments with large extra-membraneous domains. Membrane proteins from all organisms studied, except possibly the archaeon Methanococcus jannaschii, follow the so-called "positive-inside" rule; i.e., they tend to have a higher frequency of positively charged residues in cytoplasmic than in extra-cytoplasmic segments.     

Analysis of dihedral angle preferences for alanine and glycine residues in alpha and beta transmembrane regions

For the past 50?years, the Ramachandran map has been used effectively to study the protein structure and folding. However, though extensive analysis has been done on dihedral angle preferences of residues in globular proteins, related studies and reports of membrane proteins are limited. It is of interest to explore the conformational preferences of residues in transmembrane regions of membrane proteins which are involved in several important and diverse biological processes. Hence, in the present work, a systematic comparative computational analysis has been made on dihedral angle preferences of alanine and glycine in alpha and beta transmembrane regions (the two major classes of transmembrane proteins) with the aid of the Ramachandran map. Further, the conformational preferences of residues in transmembrane regions were compared with the non-transmembrane regions. We have extracted cation-pi interacting residues present in transmembrane regions and explored the dihedral angle preferences. From our observations, we reveal the higher percentage of occurrences of glycine in alpha and beta transmembrane regions than other hydrophobic residues. Further, we noted a clear shift in ψ-angle preferences of glycine residues from negative bins in alpha transmembrane regions to positive bins in beta transmembrane regions. Also, cation-pi interacting residues in beta transmembrane regions avoid preferring ψ-angles in the range of ?59° to ?30°. In this article, we insist that the studies on preferences of dihedral angles in transmembrane regions, thorough understanding of structure and folding of transmembrane proteins, can lead to modeling of novel transmembrane regions towards designing membrane proteins.     

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.     

The yeast mitochondrial transport proteins: new sequences and consensus residues, lack of direct relation between consensus residues and transmembrane helices, expression patterns of the transport protein genes, and protein-protein interactions with other proteins

Mitochondrial transport proteins (MTP) typically are homodimeric with a 30-kDa subunit with six transmembrane helices. The subunit possesses a sequence motif highly similar to Pro X Asp/Glu X X Lys/Arg X Arg within each of its three similar 10-kDa segments. Four (YNL083W, YFR045W, YPR021C, YDR470C) of the 35 yeast (S. cerevisiae) MTP genes were resequenced since the masses of their proteins deviate significantly from the typical 30 kDa. We now find these four proteins to have 545, 285, 902, and 502 residues, respectively. Together with only four other MTPs, the sequences of YPR021C and YDR470C show substitutions of some of the five residues that are absolutely conserved among the 12 MTPs with identified transport function and 17 other MTPs. We do now find these five consensus residues also in the new sequences of YNL083W and YFR045W. Additional analyses of the 35 yeast MTPs show that the location of transmembrane helix sequences do not correlate with the general consensus residues of the MTP family; protein segments connecting the six transmembrane helices and facing the intermembrane space are not uniformly short (about 20 residues) or long (about 40 residues) when facing the matrix; most MTPs have at least one transmembrane helix for which the sum of the negative hydropathy values of all residues yields a very small negative value, suggesting a membrane location bordering polar faces of other transmembrane helices or a non-transmembrane location. The extra residues of the three large MTPs are hydrophilic and at the N-terminal. The 200-residue N-terminal segment of YNL083W has four putative Ca2+-binding sites. The 500-residue N-terminal segment of YPR021C shows sequence similarity to enzymes of nucleic acid metabolism. cDNA microarray data show that YNL083W is expressed solely during sporulation, while the expressions of YFR045W, YPR021C, and YDR470C are induced by various stress situations. These results also show that the 35 MTP genes are expressed under a rather diverse set of metabolic conditions that may help identify the function of the proteins. Interestingly, yeast two-hybrid screens, that will also be useful in identifying the function of MTPs, indicate that MIR1, AAC3, YOR100C, and YPR011C do interact with non-MTPs.     

An improved hidden Markov model for transmembrane protein detection and topology prediction and its applications to complete genomes

MOTIVATION: Knowledge of the transmembrane helical topology can help identify binding sites and infer functions for membrane proteins. However, because membrane proteins are hard to solubilize and purify, only a very small amount of membrane proteins have structure and topology experimentally determined. This has motivated various computational methods for predicting the topology of membrane proteins. RESULTS: We present an improved hidden Markov model, TMMOD, for the identification and topology prediction of transmembrane proteins. Our model uses TMHMM as a prototype, but differs from TMHMM by the architecture of the submodels for loops on both sides of the membrane and also by the model training procedure. In cross-validation experiments using a set of 83 transmembrane proteins with known topology, TMMOD outperformed TMHMM and other existing methods, with an accuracy of 89% for both topology and locations. In another experiment using a separate set of 160 transmembrane proteins, TMMOD had 84% for topology and 89% for locations. When utilized for identifying transmembrane proteins from non-transmembrane proteins, particularly signal peptides, TMMOD has consistently fewer false positives than TMHMM does. Application of TMMOD to a collection of complete genomes shows that the number of predicted membrane proteins accounts for approximately 20-30% of all genes in those genomes, and that the topology where both the N- and C-termini are in the cytoplasm is dominant in these organisms except for Caenorhabditis elegans. AVAILABILITY: http://liao.cis.udel.edu/website/servers/TMMOD/     

Hiding behind hydrophobicity. Transmembrane segments in mass spectrometry

Proteomics of membrane proteins is essential for the understanding of cellular function. However, mass spectrometric analysis of membrane proteomes has been less successful than the proteomic determination of soluble proteins. To elucidate the mystery of transmembrane proteins in mass spectrometry, we present a detailed statistical analysis of experimental data derived from chloroplast membranes. This approach was further accomplished by the analysis of the Arabidopsis thaliana proteome after in silico digestion. We demonstrate that both the length and the hydrophobicity of the proteolytic fragments containing transmembrane segments are major determinants for detection by mass spectrometry. Based on a comparative analysis, we discuss possibilities to overcome the problem and provide possible protocols to shift the hydrophobicity of transmembrane segment-containing peptides to facilitate their detection.     

Identification and characterization of hepsin/-TM, a non-transmembrane hepsin isoform

Type II transmembrane serine proteases (TTSPs), including hepsin, are a new class of cell surface catalytic enzymes. In the present study, a non-transmembrane isoform of hepsin, named hepsin/-TM that originates from alternative splicing, was identified. Unlike the transmembrane hepsin isoform, this non-transmembrane isoform was distributed within the cytoplasm. Real-time PCR experiments revealed that while hepsin was expressed in all tested human tissues, hepsin/-TM was restricted in kidney, brain and lung tissues. Significantly, hepsin/-TM was not expressed in liver where hepsin was originally identified. However, hepsin/-TM was highly expressed in brain where hepsin was expressed at a relatively lower level. Moreover, these two isoforms showed different expression patterns in a number of colon adenocarcinoma cell lines. In addition, in contrast to hepsin, expression of hepsin/-TM in vivo does not exert any apparent inhibitory effect on mammalian cell growth.     

A heptad motif of leucine residues found in membrane proteins can drive self-assembly of artificial transmembrane segments

Specific interactions between alpha-helical transmembrane segments are important for folding and/or oligomerization of membrane proteins. Previously, we have shown that most transmembrane helix-helix interfaces of a set of crystallized membrane proteins are structurally equivalent to soluble leucine zipper interaction domains. To establish a simplified model of these membrane-spanning leucine zippers, we studied the homophilic interactions of artificial transmembrane segments using different experimental approaches. Importantly, an oligoleucine, but not an oligoalanine, se- quence efficiently self-assembled in membranes as well as in detergent solution. Self-assembly was maintained when a leucine zipper type of heptad motif consisting of leucine residues was grafted onto an alanine host sequence. Analysis of point mutants or of a random sequence confirmed that the heptad motif of leucines mediates self-recognition of our artificial transmembrane segments. Further, a data base search identified degenerate versions of this leucine motif within transmembrane segments of a variety of functionally different proteins. For several of these natural transmembrane segments, self-interaction was experimentally verified. These results support various lines of previously reported evidence where these transmembrane segments were implicated in the oligomeric assembly of the corresponding proteins.     

Dislocation by the m-AAA Protease Increases the Threshold Hydrophobicity for Retention of Transmembrane Helices in the Inner Membrane of Yeast Mitochondria

Sorting of mitochondrial inner membrane proteins is a complex process in which translocons and proteases function in a concerted way. Many inner membrane proteins insert into the membrane via the TIM23 translocon, and some are then further acted upon by the mitochondrial m-AAA protease, a molecular motor capable of dislocating proteins from the inner membrane. This raises the possibility that the threshold hydrophobicity for the retention of transmembrane segments in the inner membrane is different depending on whether they belong to membrane proteins that are m-AAA protease substrates or not. Here, using model transmembrane segments engineered into m-AAA protease-dependent proteins, we show that the threshold hydrophobicity for membrane retention measured in yeast cells in the absence of a functional m-AAA protease is markedly lower than that measured in its presence. Whether a given hydrophobic segment in a mitochondrial inner membrane protein will ultimately form a transmembrane helix may therefore depend on whether or not it will be exposed to the pulling force exerted by the m-AAA protease during biogenesis.     

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     

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.     

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.     

Getting greasy: how transmembrane polypeptide segments integrate into the lipid bilayer

Many integral membrane proteins use the same translocation machinery for membrane insertion as secretory proteins use to get across the membrane. This requires that transmembrane segments can be discriminated from other parts of the protein during membrane translocation, and further requires that the transmembrane segments can be moved laterally out of the translocation channel into the surrounding lipid. The molecular basis for this remarkable intramembraneous sorting event is a major focus of current studies of membrane protein biogenesis.     

Role of the Vam3p transmembrane segment in homodimerization and SNARE complex formation

Intracellular membrane fusion in eukaryotic cells is mediated by SNARE (soluble N-ethylmaleimide sensitive factor (NSF) attachment protein receptor) proteins and is known to involve assembly of cognate subunits to heterooligomeric complexes. For synaptic SNAREs, it has previously been shown that the transmembrane segments drive homotypic and support heterotypic interactions. Here, we demonstrate that a significant fraction of the yeast vacuolar SNARE Vam3p is a homodimer in detergent extracts of vacuolar membranes. This homodimer exists in parallel to the heterooligomeric SNARE complex. A Vam3p homodimer also formed from the isolated recombinant protein. Interestingly, homodimerization depended on the transmembrane segment. In contrast, formation of the quaternary SNARE complex from recombinant Vam3p, Nyv1p, Vti1p, and Vam7p subunits did not depend on the transmembrane segment of Vam3p nor on the transmembrane segments of its partner proteins. We conclude that Vam3p homodimerization, but not quaternary SNARE complex formation, is promoted by TMS-TMS interaction. As the transmembrane segments of Vam3p and other SNARE homologues were previously shown to be critical for membrane fusion downstream of membrane apposition, our results may shed light on the functional significance of SNARE TMS-TMS interactions.     

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.     

Membrane-protein integration and the role of the translocation channel

Most eukaryotic membrane proteins are integrated into the lipid bilayer during their synthesis at the endoplasmic reticulum (ER). Their integration occurs with the help of a protein-conducting channel formed by the heterotrimeric Sec61 membrane-protein complex. The crystal structure of an archaeal homolog of the complex suggests mechanisms that enable the channel to open across the membrane and to release laterally hydrophobic transmembrane segments of nascent membrane proteins into lipid. Many aspects of membrane-protein integration remain controversial and poorly understood, but new structural data provide testable hypotheses. We propose a model of how the channel recognizes transmembrane segments, orients them properly with respect to the plane of the membrane and releases them into lipid. We also discuss how the channel would prevent small molecules from crossing the lipid bilayer while it is integrating proteins.