- Research article
- Open Access
Transmembrane protein topology prediction using support vector machines
© Nugent and Jones; licensee BioMed Central Ltd. 2009
- Received: 26 September 2008
- Accepted: 26 May 2009
- Published: 26 May 2009
Alpha-helical transmembrane (TM) proteins are involved in a wide range of important biological processes such as cell signaling, transport of membrane-impermeable molecules, cell-cell communication, cell recognition and cell adhesion. Many are also prime drug targets, and it has been estimated that more than half of all drugs currently on the market target membrane proteins. However, due to the experimental difficulties involved in obtaining high quality crystals, this class of protein is severely under-represented in structural databases. In the absence of structural data, sequence-based prediction methods allow TM protein topology to be investigated.
We present a support vector machine-based (SVM) TM protein topology predictor that integrates both signal peptide and re-entrant helix prediction, benchmarked with full cross-validation on a novel data set of 131 sequences with known crystal structures. The method achieves topology prediction accuracy of 89%, while signal peptides and re-entrant helices are predicted with 93% and 44% accuracy respectively. An additional SVM trained to discriminate between globular and TM proteins detected zero false positives, with a low false negative rate of 0.4%. We present the results of applying these tools to a number of complete genomes. Source code, data sets and a web server are freely available from http://bioinf.cs.ucl.ac.uk/psipred/.
The high accuracy of TM topology prediction which includes detection of both signal peptides and re-entrant helices, combined with the ability to effectively discriminate between TM and globular proteins, make this method ideally suited to whole genome annotation of alpha-helical transmembrane proteins.
- Signal Peptide
- Unlabelled Data
- Matthews Correlation Coefficient
- Signal Peptide Prediction
- Topology Prediction
Alpha-helical transmembrane (TM) proteins constitute roughly 30% of a typical genome and are involved in a wide variety of important biological processes including cell signaling, transport of membrane-impermeable molecules and cell recognition. However, due to the experimental difficulties involved in obtaining high quality crystals, this class of protein is severely under represented in structural databases, making up only 1.7% of known structures in the PDB [1, 2]. Given the biological and pharmacological importance of TM proteins, an understanding of their topology – the total number of TM helices, their boundaries and in/out orientation relative to the membrane – is essential for structural and functional analysis, and directing further experimental work. In the absence of structural data, sequence-based prediction methods allow TM protein topology to be investigated.
Early prediction methods, based on the physicochemical principle of a sliding window of hydrophobicity combined with the 'positive-inside' rule , have been superseded by machine learning approaches which prevail due to their statistical formulation. These include Hidden Markov models (HMMs), neural networks (NNs) and more recently, support vector machines (SVMs). While NNs and HMMs are capable of producing multiple outputs, SVMs are predominantly binary classifiers therefore multiple SVMs must be employed to classify the numerous residue preferences before being combined into a probabilistic framework. Like NNs and HMMs, SVMs are capable of learning complex relationships among the amino acids within a given window with which they are trained, particularly when provided with evolutionary information. However, they are considered more resilient to the problem of over-training compared to other machine learning methods.
One problem faced by modern topology predictors is the discrimination between TM helices and other features composed largely of hydrophobic residues. These include targeting motifs such as signal peptides and signal anchors, amphipathic helices, and re-entrant helices – membrane penetrating helices that enter and exit the membrane on the same side, common in many ion channel families. The high similarity between such features and the hydrophobic profile of a TM helix frequently leads to crosstalk between the different types of predictions. Should these elements be predicted as TM helices, the ensuing topology prediction is likely to be corrupted. Some prediction methods, such as SignalP  and TargetP , are effective in identifying signal peptides, and may be used as a pre-filter prior to analysis using a TM topology predictor. Phobius  uses a HMM to successfully address the problem of signal peptides in TM protein topology prediction, while PolyPhobius  further increases accuracy by including homology information. Other methods such as TMLOOP , TOP-MOD  and OCTOPUS  have attempted to identify re-entrant regions, the latter two in combination with a TM topology predictor, but there is significant room for improvement.
TM topology predictors also exist that are able to use experimentally derived information in order to guide topology prediction. With reliable experimental data, prediction accuracy is likely to benefit substantially. Methods include HMM-TM , HMMTOP  and TMHMMfix . Tools such as SOSUI  and PRED-CLASS  are designed to discriminate between globular and TM proteins, while others such as PRED-TMBB  specialise in the discrimination and prediction of beta-barrel TM proteins.
A key element when constructing any prediction method is the use of a high quality data set for both training and validation purposes. Previously described TM data sets such as the Möller set  have contained relatively few sequences with structures available (33), but substantially more with TM region annotation based on varying types of biochemical characterisation (166). A number of experimental methods, including glycosylation analysis, insertion tags, antibody studies and fusion protein constructs, allow the topological location of a region to be identified. However, such studies are often conflicting [18, 19] and also risk upsetting the natural topology by altering the protein sequence. As a result, orientation and helix boundary errors in databases are not infrequent and add an element of noise. While such noise is often well tolerated by machine learning methods, the problem is more significant in smaller data sets. With a view to unifying, updating and verifying existing resources while minimising topological error, a new TM data set has been created containing protein sequences with topologies derived solely from crystal structures.
We thus present a new TM topology predictor trained and benchmarked with full cross-validation on a novel data set of 131 sequences with crystal structures. The method uses evolutionary information and four SVMs, combining the outputs using a dynamic programming algorithm, to return a list of predicted topologies ranked by overall likelihood, and incorporates signal peptide and re-entrant helix prediction. Overall, the method predicted the correct topology and location of TM helices for 89% of the test set, a significant improvement on our previous NN-based method MEMSAT3 . An additional SVM has been trained to discriminate between TM and globular proteins with zero false positives and a low false negative rate of 0.4%, making this method highly suitable for whole genome analysis.
Support vector machine performance
SVM per residue accuracy.
TM Helix/¬TM Helix
Inside Loop/Outside Loop
Re-entrant Helix/¬Re-entrant Helix
Signal Peptide/¬Signal Peptide
TM Protein/Globular Protein
The inside loop/outside loop SVM was the only SVM to perform optimally using a polynomial kernel, which justifies our use of multiple SVMs to classify each of the residue preferences rather than a single multiclass ranking SVM. The highest MCC value we could achieve using a radial basis function (RBF) kernel for this SVM was 0.35, significantly lower than the value of 0.63 achieved using a third-order polynomial kernel, therefore demonstrating that no single kernel function is capable of optimally separating all the data classes and suggesting the structure of loop data is strongly favoured by this kernel.
Detection of re-entrant helices remains challenging compared to other regions, with lack of training data a significant issue. Despite the addition of 64 proteins to the training set, all were homologous to one of the original 11 re-entrant helix-containing proteins and were therefore removed from the respective training files. Their contribution was therefore reflected by a low false positive rate of 0.008, but a low true positive rate of 0.478 owing to the lack of positive training examples.
In contrast, the signal peptide/¬signal peptide and TM protein/globular protein SVM performance was close to that of the TM helix/¬TM helix SVM, aided by sufficient quantities of training data. While largely driven by hydrophobicity, the signal peptide/¬signal peptide SVM must accurately discriminate between signal peptides, which contain a 7–15 residue long hydrophobic helix, and an equally hydrophobic but slightly longer TM helix. For all signal peptide-containing proteins, the residues ranked highest by the SVM appear to be close to the C-terminal end of the signal peptide region, suggesting the SVM is efficiently detecting the polar and uncharged 3–8 amino acid residue long C-region and the neutral residues that lie adjacent to the cleavage point . Similarly, the TM protein/globular protein SVM must discriminate between hydrophobic residues that compose TM helices and those that form the core of globular proteins, a challenge reflected by the difference in MCC compared to the TM helix/¬TM helix SVM.
Overall topology prediction accuracy
Benchmark results for the SVM-based method ('MEMSAT-SVM') against a selection of leading topology predictors
Correct helix count
Correct helix locations
Correct SP topology
Correct RE topology
NN + HMM
NN + HMM
To assess overall topology prediction accuracy, correct prediction of 3 components were required: the N-terminal location, number of TM helices and TM helix locations, based on an overlap of at least 5 residues with boundary definitions. Correct signal peptide and re-entrant helix predictions were not required for a correct overall topology prediction, though failure to predict these features was likely to result in an incorrect topology. Based on this definition, MEMSAT-SVM correctly predicts topology in 89% (116 out of 131) of cases, a 10% improvement on OCTOPUS which predicted 79% (103) of cases correctly (column 10). Using a more stringent criterion of a 10-residue helix overlap, the margin increases to 11% (MEMSAT-SVM 87%, OCTOPUS 76%), suggesting good segment end point prediction. In terms of the 3 individual components, MEMSAT-SVM is consistently better than all other methods (columns 3–5), and in particular performs well at predicting the correct number of TM helices (95% accuracy). MEMSAT-SVM also had a balanced number of over- and under predictions (columns 6–7) which is favourable to avoid bias towards either type of prediction, and suggests good sensitivity while avoiding over predicting helices. Since this work was completed, an extension to the OCTOPUS method which incorporates signal peptide prediction, SPOCTOPUS , has been released. This method achieved 87% accuracy on the test set, largely addressing the poor performance of OCTOPUS on sequences containing signal peptides (column 8).
Signal peptide and re-entrant helix prediction
MEMSAT-SVM correctly predicts the topology of 93% (13 out of 14) of proteins which contain signal peptides, a substantial improvement on the limited signal peptide prediction capability of our previous method MEMSAT3. In all 13 cases, signal peptides were also predicted. This accuracy is matched by PHOBIUS, the only other method that is specially trained to identify signal peptides in TM proteins. Amongst proteins that did not contain signal peptides, no false positive signal peptides were predicted. Proteins containing re-entrant helices proved much harder to predict, with only 64% (7 out of 11) correctly predicted. This is matched by MEMSAT3 and HMMTOP, though is slightly lower than the 73% (8) accuracy achieved by OCTOPUS. However, this additional correct prediction could well be attributed to the overlap between the test and training sets, as, in the absence of cross-validation, MEMSAT-SVM is able to predict 82% (9) topologies correctly. In terms of predicting re-entrant helices, MEMSAT-SVM identifies 44% (8 out of 18) with 2 false positive predictions, which compares favourably with OCTOPUS results of 22% (4) with 4 false positives. Since the numbers of proteins containing re-entrant helices and signal peptides are relatively small (14 and 11 respectively), care should be taken when interpreting these results as a relatively large percentage difference in performance may only reflect the correct prediction of one additional sequence.
MEMSAT-SVM incorrectly predicts topologies in 15 cases. Four of these correspond to proteins containing re-entrant helices that are erroneously predicted as TM helices – ABC transporter BtuCD, Proton Glutamate Symport protein, Aquaporin Z and Clc chloride channel (PDB: 1l7vB, 1xfhA, 2abmH and 2feeB) – accounting for the majority of over predicted TM helices. The remaining over prediction is due to a highly hydrophobic N-terminal region within a chain from Cytochrome bc1 (PDB: 1sqxD).
In seven cases, incorrect topologies are a result of under predicted TM helices – Photosystem I (chains A, L and K), Steryl-sulfatase, Light-Harvesting Complex II, Particulate Methane Monooxygenase and Sodium/proton antiporter 1 (PDB: 1jb0A, 1jb0L, 1jb0K, 1p49A, 1vcrA, 1yewB and 1zcdB). These under predictions fall into two categories; weakly predicted helices (PDB: 1jb0A, 1jb0L, 1jb0K, 1p49A and 1vcrA) or prediction of one helix rather than two shorter ones (PDB: 1yewB and 1zcdB). Of the weakly predicted helix errors, sequence analysis indicates low hydrophobicity for many of these helices, often due to a large fraction of charged residues. Such helices are therefore extremely difficult to predict and suggest a novel membrane insertion mechanism. Other helices appear sufficiently hydrophobic to be detected; errors are possibly the results of PSI-BLAST alignment which reduce their detectability.
The remaining three incorrect predictions are all single TM helix proteins that are inverted – Photosystem I, Cytochrome bc1 and Cytochrome b6f (PDB: 1jb0I, 1p84I and 1q90N). In all three cases, the confidence of the prediction compared to the correct topology (measured by the difference between the two scores) is extremely small. With no clear signal to differentiate between either orientation, interplay with other chains from the same protein may influence the final conformation.
Prediction accuracy using the Möller and TOPDB data sets
Prediction performance using the Möller and TOPDB data sets
We then tested performance using the TOPDB  data set, a comprehensive collection of TM protein containing experimentally derived topology information (Table 3). It currently contains records for 1452 alpha-helical TM proteins. Using this data set, MEMSAT-SVM achieved 67% accuracy, MEMSAT3 66%, OCTOPUS 64% and PHOBIUS 62%. The data set also contains 317 sequences containing signal peptides. Of these, MEMSAT-SVM correctly predicted the topologies for 77% of cases. This value was lower than that of PHOBIUS which achieved 85% accuracy. However, the MEMSAT-SVM false positive rate for signal peptide prediction is 7%, half the PHOBIUS value of 14%. These results show that on this data set, MEMSAT-SVM signal peptide performance is below that of PHOBIUS, though MEMSAT-SVM overall prediction accuracy is 5% higher due to the relatively poor performance of PHOBIUS on sequences that do not contain signal peptides (a substantially larger fraction) – 54% accuracy compared to 63% for MEMSAT-SVM. These results should again be treated with caution as they were not cross-validated.
These results are clearly lower than those attained using the crystal structure-based data set, and we believe this is likely due to errors in TOPDB. We analysed sequences from the original, uncorrected Möller set that at the time did not have crystal structures. 55 of these sequences now have a homologous PDB structure (E-value < 0.001), and of these only 38 (69%) of the original Möller topologies are correct based on current OPM definitions (taking into account only the N-terminal location and TM helix count). There is no reason to believe that the error rate in other data sets such as TOPDB, composed predominantly of sequences whose topologies were determined by biochemical means, should be significantly different. Perfect prediction methods are therefore unlikely to be able to achieve results higher than this, while older methods trained on erroneous topologies have the potential to achieve higher scores but may in reality be poorer predictors, a fact likely to be highlighted when tested against a crystal structure-based set.
Discriminating between globular and transmembrane proteins
Results for TM/globular protein discrimination rates.
False positive rate
False negative rate
NN + HMM
Application to a number of complete genomes
The fraction of proteins predicted as transmembrane, and to contain re-entrant helices and signal peptides, in a number of complete genomes.
Fraction of genome predicted as TM proteins
Fraction of TM proteins predicted to contain re-entrant helices
Fraction of TM proteins predicted to contain signal peptides
We then carried out full topology prediction on sequences predicted to be TM proteins and analysed these for the presence of re-entrant helices and signal peptides. In most species, re-entrant helices were detected in at least 2% of TM proteins, with more than 3% detected in Takifugu rubripes which can be explained by the extensive channel heterogeneity discussed above. However, given the low true positive rate of 44%, this figure is likely to be an underestimate. A positive predictive value (PPV) of 0.8 suggests a value in the range 3–4.5% is more realistic. This range is close to one previous estimate of 5%  but below another of 10% , although the latter was based on a broader definition of re-entrant regions that did not necessarily contain helical secondary structure.
In this paper we have implemented a novel SVM-based TM protein topology predictor, an area previously dominated by HMM and NN-based machine learning approaches, and have shown that it outperforms a selection of the best performing prediction methods when fully cross-validated on a novel high resolution data set of 131 protein sequences. This data set includes proteins containing both re-entrant helices and signal peptides, features that this method is also able to predict. The method has also been benchmarked on the Möller data set, which contains a higher fraction of eukaryotic sequences, improving on the best current methods. And we have achieved extremely low false positive and false negative rates for TM/globular protein discrimination. Using these tools, we have estimated the fraction of TM proteins, re-entrant helices and signal peptides in a number of complete genomes. Overall, our results suggest that MEMSAT-SVM is ideally suited to whole genome annotation of alpha-helical TM proteins.
Assembling a novel data set of transmembrane proteins
The novel data set was based solely on crystal structure data. Additional information was collected from MPTOPO , OPM , PDB_TM , SWISS-PROT  and from the literature. SWISS-PROT files were parsed for entries containing the keyword 'TRANSMEM' in feature table (FT) lines. N-terminal data was also extracted using keyword 'TOPO_DOM' where available. To avoid partial sequences being included, any entries containing keywords such as 'FRAGMENT' were excluded. Sequences were then scanned against the PDB in order to identify entries for which the TM region had complete structural coverage. Alignments occasionally highlighted chain breaks. In these cases, the sequence was excluded unless a visual inspection ensured the topology could not be cast in doubt by the break. This left a redundant data set containing 944 sequences which was then homology reduced at the 40% sequence identity level. A number of sequences were then removed. These included colicins (e.g. PDB: 1col) and bee venom (PDB: 2mlt) which are bilayer disrupting and thus are not native integral membrane proteins, sequences labelled as 'secreted protein', and sequences where the N-terminal location or topology could not be verified.
OPM was then used to define TM helix boundaries, or in the absence of an OPM entry, PDB_TM was used. OPM uses a theoretical multi-feature approach to position proteins in a membrane which has been shown to be in good agreement with experimental studies of 24 TM proteins. In some cases where a visual inspection appeared to indicate an incorrect placement of the membrane, PDB_TM helix boundary definitions were used instead. For example, OPM lists Mechanosensitive channel protein MscS (PDB: 2oau) as having two TM helices, neither of which fully cross the membrane, whereas the PDB_TM definition of 3 TM helices appears more plausible.
PDB_TM was also used to annotate proteins containing re-entrant helices. A re-entrant helix was defined as a helix-containing region that enters and exits the membrane on the same side, penetrating at least 6Å but not more than 6Å from the opposite membrane face. Re-entrant regions that did not contain a helix formed by at least three contiguous amino acid residues were excluded. Sequences containing signal peptides were then labelled according to SWISS-PROT annotations.
Data set composition.
Number in set
Single-spanning TM segment
Multiple-spanning TM segments
Contains re-entrant helix
Contains signal peptide
Support vector machine training
As Support Vector Machines (SVMs) are binary classifiers , we chose to combine multiple SVMs to classify each of the residue preferences found in TM proteins. Although multiclass ranking SVMs do exist, they are generally considered unreliable since in many cases no single mathematical function exists to separate all classes of data from one another . We therefore trained four SVMs to classify TM helix/¬TM helix, inside loop/outside loop, re-entrant helix/¬re-entrant helix and signal peptide/¬signal peptide. Residue labelling was performed according to our data set definitions.
x is the raw score to be normalised. μ and σ are the mean and standard deviation PSI-BLAST score for each of the 20 amino acid, generated using profiles for all 131 sequences.
In order to accentuate the contribution of re-entrant helices for which data is particularly sparse, the sequences of 64 proteins, all homologous to the 11 re-entrant helix-containing sequences in our initial data set, were also used to train the TM helix/¬TM helix and re-entrant helix/¬re-entrant helix SVMs. Helix, loop and re-entrant helix boundaries were determined by PDB_TM definitions.
We also attempted to train the TM helix/¬TM helix SVM using unlabelled data via transduction. In transduction, the learning task is to assign labels to unlabelled data as accurately as possible . SVMs can perform transduction by finding the hyperplane that maximises the margin relative to both the labelled and unlabelled data, in order to improve the generalisation performance. We selected sequences from SWISS-PROT identified by the MEMSAT3 TM/globular protein discriminator as TM proteins. Sequences with greater than 40% sequence identity to sequences in the labelled data set were removed, as were those with signal peptides predicted by SignalP. Of those remaining, 135 sequences were used as unlabelled training data.
For training the signal peptide/¬signal peptide SVM, we included data from the Phobius training set which contains 2654 well annotated examples of TM and globular proteins, with and without signal peptides. This was supplemented by a search of SWISS-PROT for sequences labelled with the keyword 'SIGNAL' (but excluding entries labelled 'POTENTIAL' or 'BY SIMILARITY') to add to the signal peptide set, and sequences without keyword 'SIGNAL' to add to the non-signal peptide set. The combined set was then homology reduced at the 40% sequence identity level, leaving 3205 (1222 with signal peptides, 1983 without signal peptides) sequences for which PSI-BLAST profiles were then generated as outlined above.
Stringent cross validation was performed using a jack knife test (leave-one-out cross validation) for the TM helix/¬TM helix, inside loop/outside loop and re-entrant helix/¬re-entrant helix SVMs. In training, the target sequence, along with any other sequences with greater than 25% sequence identity, were excluded. For the signal peptide/¬signal peptide SVM we used 10-fold cross validation, again excluding sequences from the training set with greater than 25% sequence identity to any sequence in the test set. For training and classification, SVM-Light  was used. The performance of several kernels was investigated in combination with a comprehensive grid search of SVM parameters. To determine optimal windows sizes, the data set was split randomly into two and the highest scoring window which ranked equally in each split was selected, therefore demonstrating consistency between data sets and reducing the risk of overfitting. We used the MCC to optimise these values which is a more robust measure than using recall or precision alone .
To calculate a list of topologies ranked by overall likelihood, the TM helix/¬TM helix, inside loop/outside loop and signal peptide/¬signal peptide raw SVM outputs were combined in a modified version of the dynamic programming algorithm used in the original MEMSAT method . The algorithm was simplified by treating TM helices as discrete units, rather than separating them into inside, outside and middle components, though a signal peptide state was added. Loop regions between predicted TM helices were scanned for re-entrant helices using the re-entrant helix/¬re-entrant helix raw SVM output and a simple scoring function. For evaluating signal peptide preference, residues with positive signal peptide scores up to position 40 in a target sequence were added to the outside loop score and subtracted from the inside loops score where positive, in order to direct prediction towards a non-cytoplasmic amino-terminus. The value was also scaled by a factor of 10 and subtracted from the TM helix SVM score to prevent TM helix prediction. Residues were therefore predicted to lie in one of five different topological regions: inside loop, outside loop, TM helix, re-entrant helix and signal peptide.
To evaluate performance, four metrics were used. Firstly, correct location of the amino terminus; secondly, correct number of TM helices; thirdly, correct number and location of TM helices (based on an overlap of at least five residues with the helix boundaries in our data set) and fourthly, correct overall topology. For comparison, we also evaluated a number of other leading topology predictors. For this method and MEMSAT3, the appropriate cross-validated training data was used in assessing performance. Where equivalent data was unavailable for the other methods, performance is likely to be overestimated as it is likely that there is significant overlap between test and training sets. We also assessed performance of the method against proteins containing signal peptides and re-entrant helices.
We also trained an additional SVM to discriminate between TM and globular proteins, to be used as a pre-filter prior to TM topology prediction. For SVM training, we used the data set of 131 TM proteins and 416 globular proteins from non-redundant PDB chains as used by MEMSAT3. To accurately compare with MEMSAT3 we used exactly the same test set consisting of 184 TM proteins from the Möller data set and a separate set of 2269 non-redundant globular protein chains, giving a total of 2453 test cases. PSI-BLAST profiles were generated for all sequences and 10-fold cross validation was used to assess performance, again removing sequences from the training fold with greater than 25% sequences identity to any sequence in the test fold.
For whole genome analysis, ten genomes – nine eukaryotic and one prokaryotic – were downloaded from the Ensembl  and NCBI  websites. Protein sequences were extracted and PSI-BLAST profiles were generated using the SWISS-PROT database. The TM/globular predictor was used to identify TM proteins, which were then subject to full topology prediction.
MEMSAT-SVM is available as downloadable source code and as a web server from the URL below and is free for non-commercial use. All data sets are also available, and cross-validation SVM model files are available on request. The software has been tested on a Linux operating system. In order to compile and run, the gcc compiler, Perl interpreter, and NCBI tools are required.
Part of this work was supported by the BioSapiens project, which is funded by the European Commission within its FP6 Programme, under the thematic area "Life sciences, genomics and biotechnology for health", contract number LSHG-CT-2003-503265. Additional funding was provided by the Biotechnology and Biological Sciences Research Council.
- Berman HM, Henrick K, Nakamura H: Announcing the worldwide Protein Data Bank. Nature Structural Biology 2003, 10(12):980. 10.1038/nsb1203-980View ArticlePubMedGoogle Scholar
- White S: The progress of membrane protein structure determination. Protein Sci 2004, 13: 1948–1949. 10.1110/ps.04712004PubMed CentralView ArticlePubMedGoogle Scholar
- von Heijne G: Membrane Protein Structure Prediction, Hydrophobicity Analysis and the Positive-inside Rule. Mol Biol 1992, 225(2):487–494. 10.1016/0022-2836(92)90934-CView ArticleGoogle Scholar
- Bendtsen JD, Nielsen H, von Heijne G, Brunak S: Improved prediction of signal peptides: SignalP 3.0. J Mol Biol 2004, 340: 783–795. 10.1016/j.jmb.2004.05.028View ArticlePubMedGoogle Scholar
- Emanuelsson O, Brunak S, von Heijne G, Nielsen H: Locating proteins in the cell using TargetP, SignalP and related tools. Nat Protoc 2007, 2(4):953–971. 10.1038/nprot.2007.131View ArticlePubMedGoogle Scholar
- Käll L, Krogh A, Sonnhammer EL: A combined transmembrane topology and signal peptide prediction method. J Mol Biol 2004, 338(5):1027–1036. 10.1016/j.jmb.2004.03.016View ArticlePubMedGoogle Scholar
- Käll L, Krogh A, Sonnhammer EL: An HMM posterior decoder for sequence feature prediction that includes homology information. Bioinformatics 2005, 21: i251–257. 10.1093/bioinformatics/bti1014View ArticlePubMedGoogle Scholar
- Lasso G, Antoniw JF, Mullins JG: A combinatorial pattern discovery approach for the prediction of membrane dipping (re-entrant) loops. Bioinformatics 2006, 22(14):290–297. 10.1093/bioinformatics/btl209View ArticleGoogle Scholar
- Viklund H, Granseth E, Elofsson A: Structural classification and prediction of reentrant regions in alpha-helical transmembrane proteins: application to co mplete genomes. J Mol Biol 2006, 361(3):591–603. 10.1016/j.jmb.2006.06.037View ArticlePubMedGoogle Scholar
- Viklund H, Elofsson A: OCTOPUS: Improving topology prediction by two-track ANN-based preference scores and an extended topological grammar. Bioinformatics 2008, 24(15):1662–1668. 10.1093/bioinformatics/btn221View ArticlePubMedGoogle Scholar
- Bagos PG, Liakopoulos TD, Hamodrakas SJ: Algorithms for incorporating prior topological information in HMMs: application to transmembrane proteins. BMC Bioinformatics 2006, 7: 189. 10.1186/1471-2105-7-189PubMed CentralView ArticlePubMedGoogle Scholar
- Tusnady GE, Simon I: The HMMTOP transmembrane topology prediction server. Bioinformatics 2001, 17(9):849–850. 10.1093/bioinformatics/17.9.849View ArticlePubMedGoogle Scholar
- Melén K, Krogh A, von Heijne G: Reliability measures for membrane protein topology prediction algorithms. J Mol Biol 2003, 327(3):735–744. 10.1016/S0022-2836(03)00182-7View ArticlePubMedGoogle Scholar
- Hirokawa T, Boon-Chieng S, Mitaku S: SOSUI: classification and secondary structure prediction system for membrane proteins. Bioinformatics 1998, 14: 378–379. 10.1093/bioinformatics/14.4.378View ArticlePubMedGoogle Scholar
- Pasquier C, Promponas VJ, Hamodrakas SJ: PRED-CLASS: cascading neural networks for generalized protein classification and genome-wide applications. Proteins 2001, 44(3):361–369. 10.1002/prot.1101View ArticlePubMedGoogle Scholar
- Bagos PG, Liakopoulos TD, Spyropoulos IC, Hamodrakas SJ: PRED-TMBB: a web server for predicting the topology of beta-barrel outer membrane proteins. Nucleic Acids Res 2004, 32: W400–404. 10.1093/nar/gkh417PubMed CentralView ArticlePubMedGoogle Scholar
- Möller S, Kriventseva EV, Apweiler R: A collection of well characterised integral membrane proteins. Bioinformatics 2000, 16(12):1159–1160. 10.1093/bioinformatics/16.12.1159View ArticlePubMedGoogle Scholar
- Kyttälä A, Ihrke G, Vesa J, Schell MJ, Luzio JP: Two motifs target Batten disease protein CLN3 to lysosomes in transfected nonneuronal and neuronal cells. Mol Biol Cell 2004, 15(3):1313–1323. 10.1091/mbc.E03-02-0120PubMed CentralView ArticlePubMedGoogle Scholar
- Mao Q, Foster BJ, Xia H, Davidson BL: Membrane topology of CLN3, the protein underlying Batten disease. FEBS Lett 2003, 541(1–3):40–46. 10.1016/S0014-5793(03)00284-9View ArticlePubMedGoogle Scholar
- Jones DT: Improving the accuracy of transmembrane protein topology prediction using evolutionary information. Bioinformatics 2007, 23(5):538–544. 10.1093/bioinformatics/btl677View ArticlePubMedGoogle Scholar
- von Heijne G: Patterns of amino acids near signal-sequence cleavage sites. Eur J Biochem 1983, 133: 17. 10.1111/j.1432-1033.1983.tb07424.xView ArticlePubMedGoogle Scholar
- Viklund H, Bernsel A, Skwark M, Elofsson A: SPOCTOPUS: A combined predictor of signal peptides and membrane protein topology. Bioinformatics 2008, 24: 2928–2929. 10.1093/bioinformatics/btn550View ArticlePubMedGoogle Scholar
- Granseth E, Seppälä S, Rapp M, Daley DO, von Heijne G: Membrane protein structural biology – how far can the bugs take us? Mol Membr Biol 2007, 24(5–6):329–332. 10.1080/09687680701413882View ArticlePubMedGoogle Scholar
- Tusnády GE, Kalmár L, Simon I: TOPDB: topology data bank of transmembrane proteins. Nucleic Acids Res 2008, D234–239.Google Scholar
- Perlman D, Halvorson HO: A putative signal peptidase recognition site and sequence in eukaryotic and prokaryotic signal peptides. J Mol Biol 1983, 167: 391. 10.1016/S0022-2836(83)80341-6View ArticlePubMedGoogle Scholar
- Viklund H, Elofsson A: Best alpha-helical transmembrane protein topology predictions are achieved using hidden Markov models and evolutionary information. Protein Sci 2004, 13(7):1908–1917. 10.1110/ps.04625404PubMed CentralView ArticlePubMedGoogle Scholar
- Rost B, Fariselli P, Casadio R: Topology prediction for helical transmembrane proteins at 86% accuracy. Protein Sci 1996, 5(8):1704–1718. 10.1002/pro.5560050824PubMed CentralView ArticlePubMedGoogle Scholar
- Martelli PL, Fariselli P, Casadio R: An ENSEMBLE machine learning approach for the prediction of all-alpha membrane proteins. Bioinformatics 2003, 19(1):i205–211. 10.1093/bioinformatics/btg1027View ArticlePubMedGoogle Scholar
- Amico M, Finelli M, Rossi I, Zauli A, Elofsson A, Viklund H, von Heijne G, Jones D, Krogh A, Fariselli P, Luigi Martelli P, Casadio R: PONGO: a web server for multiple predictions of all-alpha transmembrane proteins. Nucleic Acids Res 2006, 34: W169–172. 10.1093/nar/gkl208PubMed CentralView ArticlePubMedGoogle Scholar
- Lo A, Chiu HS, Sung TY, Lyu PC, Hsu WL: Enhanced membrane protein topology prediction using a hierarchical classification method and a new scoring function. Journal of Proteome Research 2008, 7(2):487–496. 10.1021/pr0702058View ArticlePubMedGoogle Scholar
- Chen CP, Kernytsky A, Rost B: Transmembrane helix predictions revisited. Protein Sci 2002, 11(12):2774–2791. 10.1110/ps.0214502PubMed CentralView ArticlePubMedGoogle Scholar
- Wallin E, von Heijne G: Genome-wide analysis of integral membrane proteins from eubacterial, archaean, and eukaryotic organisms. Protein Sci 1998, 7(4):1029–1038.PubMed CentralView ArticlePubMedGoogle Scholar
- Wong E, Yu WP, Yap WH, Venkatesh B, Soong TW: Comparative genomics of the human and Fugu voltage-gated calcium channel alpha1-subunit gene family reveals greater diversity in Fugu. Gene 2006, 366(1):117–127. 10.1016/j.gene.2005.08.022View ArticlePubMedGoogle Scholar
- Larsson TA, Olsson F, Sundström G, Brenner S, Venkatesh B, Larhammar D: Pufferfish and zebrafish have five distinct NPY receptor subtypes, but have lost appetite receptors Y1 and Y5. Ann N Y Acad Sci 2005, 1040: 375–377. 10.1196/annals.1327.066View ArticlePubMedGoogle Scholar
- 35. Liu Y, Engelman DM, Gerstein M: Genomic analysis of membrane protein families: abundance and conserved motifs. Genome Biol 2002.,19;3(10)Google Scholar
- Liu Y, Gerstein M, Engelman DM: Transmembrane protein domains rarely use covalent domain recombination as an evolutionary mechanism. Proc Natl Acad Sci USA 2004, 101(10):3495–3497. 10.1073/pnas.0307330101PubMed CentralView ArticlePubMedGoogle Scholar
- Petty HR: Molecular Biology of Membranes: Structure and Function. New York: Springer; 1993.View ArticleGoogle Scholar
- Cuthbertson JM, Doyle DA, Sansom MS: Transmembrane helix prediction: a comparative evaluation and analysis. Protein Eng Des Sel 2005, 18(6):295–308. 10.1093/protein/gzi032View ArticlePubMedGoogle Scholar
- Gafvelin G, Sakaguchi M, Andersson H, von Heijne G: Topological rules for membrane protein assembly in eukaryotic cells. J Biol Chem 1997, 272(10):6119–6127. 10.1074/jbc.272.10.6119View ArticlePubMedGoogle Scholar
- Krogh A, Larsson B, von Heijne G, Sonnhammer EL: Predicting transmembrane protein topology with a hidden Markov model: application to complete genomes. J Mol Biol 2001, 305(3):567–580. 10.1006/jmbi.2000.4315View ArticlePubMedGoogle Scholar
- Jayasinghe S, Hristova K, White SH: MPtopo: A database of membrane protein topology. Protein Sci 2001, 10(2):455–458. 10.1110/ps.43501PubMed CentralView ArticlePubMedGoogle Scholar
- Lomize MA, Lomize AL, Pogozheva ID, Mosberg HI: OPM: Orientations of Proteins in Membranes database. Bioinformatics 2006, 22: 623–625. 10.1093/bioinformatics/btk023View ArticlePubMedGoogle Scholar
- Tusnády GE, Dosztányi Z, Simon I: PDB_TM: selection and membrane localization of transmembrane proteins in the protein data bank. Nucleic Acids Res 2005, 33: D275–278. 10.1093/nar/gki002PubMed CentralView ArticlePubMedGoogle Scholar
- Boeckmann B, Bairoch A, Apweiler R, Blatter M-C, Estreicher A, Gasteiger E, Martin MJ, Michoud K, O'Donovan C, Phan I, Pilbout S, Schneider M: The SWISS-PROT protein knowledgebase and its supplement TrEMBL in 2003. Nucleic Acids Res 2003, 31: 365–370. 10.1093/nar/gkg095PubMed CentralView ArticlePubMedGoogle Scholar
- Vapnik V: Statistical Learning Theory. New York: Wiley-Interscience; 1998.Google Scholar
- Abe S: Analysis of Multiclass Support Vector Machines. International Conference on Computational Intelligence for Modelling Control and Automation 2003, 385–396.Google Scholar
- Altschul SF, Madden TL, Schäffer AA, Zhang J, Zhang Z, Miller W, Lipman DJ: Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res 1997, 25(17):3389–3402. 10.1093/nar/25.17.3389PubMed CentralView ArticlePubMedGoogle Scholar
- Hedman M, Deloof H, von Heijne G, Elofsson A: Improved detection of homologous membrane proteins by inclusion of information from topology predictions. Protein Sci 2002, 11(3):652–658. 10.1110/ps.39402PubMed CentralView ArticlePubMedGoogle Scholar
- Chen Y, Wang G, Dong S: Learning with progressive transductive support vector machine. Pattern Recognition Letters 2003, 24(12):1845–1855. 10.1016/S0167-8655(03)00008-4View ArticleGoogle Scholar
- Joachims T: Making large-scale SVM learning practical. In Advances in Kernel Methods-Support Vector Learning. Edited by: Schölkopf B, Burges C, Smola A. Cambridge: MIT; 1998.Google Scholar
- Matthews BW: Comparison of the predicted and observed secondary structure of T4 phage lysozyme. Biochim Biophys Acta 1975, 405(2):442–451.View ArticlePubMedGoogle Scholar
- Jones DT, Taylor WR, Thornton JM: A model recognition approach to the prediction of all-helical membrane protein structure and topology. Biochemistry 1994, 33: 3038–3048. 10.1021/bi00176a037View ArticlePubMedGoogle Scholar
- Flicek P, Aken BL, Beal K, Ballester B, Caccamo M, Chen Y, Clarke L, Coates G, Cunningham F, Cutts T, Down T, Dyer SC, Eyre T, Fitzgerald S, Fernandez-Banet J, Gräf S, Haider S, Hammond M, Holland R, Howe KL, Howe K, Johnson N, Jenkinson A, Kähäri A, Keefe D, Kokocinski F, Kulesha E, Lawson D, Longden I, Megy K, Meidl P, Overduin B, Parker A, Pritchard B, Prlic A, Rice S, Rios D, Schuster M, Sealy I, Slater G, Smedley D, Spudich G, Trevanion S, Vilella AJ, Vogel J, White S, Wood M, Birney E, Cox T, Curwen V, Durbin R, Fernandez-Suarez XM, Herrero J, Hubbard TJ, Kasprzyk A, Proctor G, Smith J, Ureta-Vidal A, Searle S: Ensembl 2008. Nucleic Acids Res 2008, 36: D707–714. 10.1093/nar/gkm988PubMed CentralView ArticlePubMedGoogle Scholar
- Benson DA, Karsch-Mizrachi I, Lipman DJ, Ostell J, Wheeler DL: GenBank. Nucleic Acids Res 2008, 36: D25–30. 10.1093/nar/gkm929PubMed CentralView ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.