An automatic method for identifying surface proteins in bacteria: SLEP
© Giombini et al; licensee BioMed Central Ltd. 2010
Received: 13 September 2009
Accepted: 20 January 2010
Published: 20 January 2010
Bacterial infections represent a global health challenge. The identification of novel antibacterial targets for both therapy and vaccination is needed on a constant basis because resistance continues to spread worldwide at an alarming rate. Even infections that were once easy to treat are becoming difficult or, in some cases, impossible to cure. Ideal targets for both therapy and vaccination are bacterial proteins exposed on the surface of the organism, which are often involved in host-pathogen interaction. Their identification can greatly benefit from technologies such as bioinformatics, proteomics and DNA microarrays.
Here we describe a pipeline named SLEP (Surface Localization Extracellular Proteins), based on an automated optimal combination and sequence of usage of reliable available tools for the computational identification of the surfome, i.e. of the subset of proteins exposed on the surface of a bacterial cell.
The tool not only simplifies the usage of these methods, but it also improves the results by selecting the specifying order and combination of the instruments. The tool is freely available at http://www.caspur.it/slep.
Bacterial diseases are among the major causes of mortality and morbidity in humans. Antibiotics are the first line of defence against bacteria, however more and more bacteria are antibiotic resistant and the phenomenon is spreading at an alarming rate [1, 2]. Many diseases are becoming increasingly difficult to fight. There are several examples of microbial infections that are becoming resistant to all existing therapies and for which a vaccination strategy is deemed to be appropriate, such as gonorrhea, tuberculosis, pneumonia, septicaemia and childhood ear infections [3–7].
Among the proteins encoded by bacteria, secreted and surface proteins are particularly important in bacterial pathogenesis. The former can be involved in host cell toxicity and lead to more or less subtle alterations of the host cell for the benefit of the pathogen. Bacterial surface proteins play a fundamental role in the interaction with the cell environment [8–12]. They can be involved in adhesion and invasion of the host cells as well as in defending against host responses. Because of this, surface proteins are potential drug targets . Moreover, surface proteins are likely to interact with the host immune system and are ideal candidates for vaccine development [14–16].
Surface proteins include integral or transmembrane proteins that span the membrane and have a hydrophilic cytosolic domain, which interacts with internal molecules, a hydrophobic membrane-spanning domain that anchors it within the cell membrane, and a hydrophilic extracellular domain that interacts with external molecules. Lipid anchored proteins are instead covalently-bound to one or more lipid molecules. Other membrane proteins are peripheral, i.e. they are attached to integral membrane proteins, or associated with regions of the lipid bilayer.
Gram-positive bacteria possess a thick cell wall containing many layers of peptidoglycan and teichoic acids. In contrast, Gram-negative bacteria have a relatively thin cell wall consisting of a few layers of peptidoglycan surrounded by a second lipid membrane containing lipopolysaccharides and lipoproteins. This is reflected in their membrane protein composition. Cell wall proteins are found in Gram+ bacteria while β-barrel membrane proteins are only found in the outer membranes of Gram- organisms, in mitochondria and chloroplast .
Despite the biological relevance of bacterial surface proteins, their characterization is still incomplete. There are two main routes to identify surface proteins. In one approach, membrane and cell wall fractions are separated from the cytoplasmic fraction and then proteins are identified by two-dimensional (2D)-electrophoresis or 2D-chromatography coupled to mass spectrometry [see for example [18–23]]. The other possibility is to take advantage of bioinformatics and attempt their prediction on the basis of one of the many specifically developed algorithms.
There is a plethora of available tools for predicting the membrane localization and topology of a protein and the presence of specific localization signals in its sequence, but not every method is equally accurate and, especially, an end user is not always well informed about novel developments in the field. The order in which these tools are used might also make a difference, as we will show here. Furthermore, each of them tends to use different input formats and not always self explanatory output formats.
The aim of the work described here is to bring these tools in a coordinated and easy-to use form to the bench scientists who, on one side, should not need to be familiar with the ins and outs of each and every tool, but, on the other, should be given sufficient information to assess the reliability of the methods.
SLEP and all the related tools have been implemented locally on a linux SLES 10 server.
If the user inputs a genome, putative genes need to be identified and translated into their amino acid sequence. This is achieved using Glimmer, a gene finding program based on Interpolated Markov Models (IMMs) [24–26]. The accuracy of gene identification by Glimmer depends upon the length and the GC-content the genome and is reported at http://www.cbcb.umd.edu/software/glimmer/.
The translated gene products, or the input proteins (if the user selected to start with a known proteome) are analysed for the presence of transmembrane regions using TMHMM [27, 28] and prodiv-HMM , two independent Hidden Markov Model-based prediction methods. It is a known problem in the field that signal peptides might often be mispredicted as transmembrane helices and vice versa. To alleviate this problem, we only assign the "membrane protein" tag to proteins for which more than three transmembrane helix are predicted by at least one method. As described later, proteins for which no signal peptide is identified are re-submitted to the transmembrane prediction tools.
Comparison between SLEP and LipoP
The next step consists in running pSORTb [31, 32] on the remaining set of proteins. pSORTb is used for recognizing cell wall proteins (in Gram+ bacteria) and outer membrane proteins (in Gram- bacteria) as well as exported proteins. The remaining proteins are reanalysed by TMHMM and prodiv-HMM in order to identify proteins with a single membrane spanning helix. As mentioned before, we remove clearly detectable membrane proteins before attempting the prediction of the presence of signal peptides. Only if no signal peptide has been identified in the sequence, we look for single membrane spanning helix.
Where TP, TN, FP and FN are the number of True Positive, True Negative, False Positive and False Negative results, respectively.
Results and discussion
We tested the accuracy of the procedure using the manually curated database SwissProt. This dataset, named SP, contained 18,510 protein sequences of known localization (as reported in the SUBCELLULAR LOCATION field), roughly equally populated by protein from Gram- and Gram+ bacteria (9,946 and 9,564, respectively). For Gram+ bacteria we used Enterococcus (EN, 228 proteins), Listeria (LI, 749 proteins), Staphylococcus (SP, 3981 proteins), Streptococcus (ST, 179 proteins) and a pool of Bacilli Gram+ organisms (B+, 4427 proteins). For Gram-, the datasets contained data from E. coli (EC, 3891 proteins), Legionella (LE, 421 proteins), Pseudomonas (PS, 3369 proteins), Salmonella (SA, 2019 proteins) and a pool of Bacilli Gram- organisms (B-, 246 proteins).
SLEP overall accuracy
Gram+ Total number of proteins: 9,564
Gram-. Total number of proteins: 9,946
Comparison between SLEP and other available tools
Cell wall/Outer membrane
Bioinformatics tools are extremely useful for the bench scientists and most of them are mature enough to be considered part of a toolbox that should be readily and easily accessible to all.
The appropriate usage of the tools is however essential. This is far from being trivial: one of the most cogent problems in bioinformatics is that way too often obsolete tools remain available and are used by experimentalists who are unaware of more recent developments. Users are confronted with too many available tools, not all properly benchmarked and updated and this can result in a waste of time and effort. The problem is even more relevant when the methods need to be used as start points of a set of experiments where an incorrect selection/usage of the methods can seriously affect the end results.
The initial selection of the set of transcripts/proteins from a pathogen to be used as targets for the development of vaccines and/or inhibitor screening is one such case and yet no comprehensive easy-to-use system was available so far. Perhaps the most complete resource available is Augur  which includes a precompiled list of protein localizations and other useful features, but does not allow users to supply their own genome/proteome or set of proteins as is the case in SLEP and is limited to Gram negative bacteria.
We have described here an automatic procedure designed to achieve an accurate prediction of bacterial protein localization via an appropriate sequence of usage of the available methods that is, at the same time, extremely easy to use.
SLEP uses a combination of state of the art methods that have been shown to be the most accurate available [29, 30, 32]. The specific order of usage of these programs has been designed to reduce the chance of misclassification by each of the tools.
The system relieves the bench scientists from the burden of selecting the most accurate programs for the task at hand. SLEP will be continuously updated to reflect novel developments and plans to be the one-stop shop for the analysis of bacterial protein localization that is perhaps the most important aspect of therapeutic target selection.
Availability and requirements
This work is supported by the FIRB projects ITALBIONET, Proteomica and by the Italian Ministry of Labour, Health and Social Policies contract no.onc_ord 25/07.
- Cohen R: The need for prudent use of antibiotics and routine use of vaccines. Clin Microbiol Infect 2009, 15(Suppl 3):21–23. 10.1111/j.1469-0691.2009.02727.xView ArticlePubMedGoogle Scholar
- Salisbury D, Ramsay M, Noakes K, (eds): Immunisation against infectious disease. 3rd edition. London: The Stationery Office; 2006.
- Balestrazzi P, Banchini G, Menozzi MG, Frassi C, De Angelis GL, Giovannelli G: [Yersinia enterocolitica infection in thalassemia. Report of one case (author's transl)]. Pediatr Med Chir 1981, 3(4):269–271.PubMedGoogle Scholar
- Obaro SK, Monteil MA, Henderson DC: The pneumococcal problem. Bmj 1996, 312(7045):1521–1525.View ArticlePubMedPubMed CentralGoogle Scholar
- Bergmann S, Hammerschmidt S: Versatility of pneumococcal surface proteins. Microbiology 2006, 152(Pt 2):295–303. 10.1099/mic.0.28610-0View ArticlePubMedGoogle Scholar
- Kostiukova NN, Bekhalo VA: [Epidemic process of gonorrhea in modern world]. Zh Mikrobiol Epidemiol Immunobiol 2009, (1):87–93.
- The WHO/IUATLD Global Project on Anti-tuberculosis Drug Resistance Surveillance 2002–2007/[written by Abigail Wright, Matteo Zignol]: Anti-tuberculosis drug resistance in the world: fourth global report Geneva, Switzerland: World Health Organization; 2008. WHO/HTM/TB/2008.394
- Lindahl G, Stalhammar-Carlemalm M, Areschoug T: Surface proteins of Streptococcus agalactiae and related proteins in other bacterial pathogens. Clin Microbiol Rev 2005, 18(1):102–127. 10.1128/CMR.18.1.102-127.2005View ArticlePubMedPubMed CentralGoogle Scholar
- Lin J, Huang S, Zhang Q: Outer membrane proteins: key players for bacterial adaptation in host niches. Microbes Infect 2002, 4(3):325–331. 10.1016/S1286-4579(02)01545-9View ArticlePubMedGoogle Scholar
- Niemann HH, Schubert WD, Heinz DW: Adhesins and invasins of pathogenic bacteria: a structural view. Microbes Infect 2004, 6(1):101–112. 10.1016/j.micinf.2003.11.001View ArticlePubMedGoogle Scholar
- Galperin MY, Koonin EV: Searching for drug targets in microbial genomes. Curr Opin Biotechnol 1999, 10(6):571–578. 10.1016/S0958-1669(99)00035-XView ArticlePubMedGoogle Scholar
- Pizza M, et al.: Identification of vaccine candidates against serogroup B meningococcus by whole-genome sequencing. Science 2000, 287(5459):1816–1820. 10.1126/science.287.5459.1816View ArticlePubMedGoogle Scholar
- Maione D, Margarit I, Rinaudo CD, Masignani V, Mora M, Scarselli M, Tettelin H, Brettoni C, Iacobini ET, Rosini R, D'Agostino N, Miorin L, Buccato S, Mariani M, Galli G, Nogarotto R, Nardi Dei V, Vegni F, Fraser C, Mancuso G, Teti G, Madoff LC, Paoletti LC, Rappuoli R, Kasper DL, Telford JL, Grandi G: Identification of a universal Group B streptococcus vaccine by multiple genome screen. Science 2005, 309(5731):148–150. 10.1126/science.1109869View ArticlePubMedPubMed CentralGoogle Scholar
- Rodriguez-Ortega MJ, Norais N, Bensi G, Liberatori S, Capo S, Mora M, Scarselli M, Doro F, Ferrari G, Garaguso I, Maggi T, Neumann A, Covre A, Telford JL, Grandi G: Characterization and identification of vaccine candidate proteins through analysis of the group A Streptococcus surface proteome. Nat Biotechnol 2006, 24(2):191–197. 10.1038/nbt1179View ArticlePubMedGoogle Scholar
- Ton-That H, Marraffini LA, Schneewind O: Protein sorting to the cell wall envelope of Gram-positive bacteria. Biochim Biophys Acta 2004, 1694(1–3):269–278. 10.1016/j.bbamcr.2004.04.014View ArticlePubMedGoogle Scholar
- Janulczyk R, Rasmussen M: Improved pattern for genome-based screening identifies novel cell wall-attached proteins in gram-positive bacteria. Infect Immun 2001, 69(6):4019–4026. 10.1128/IAI.69.6.4019-4026.2001View ArticlePubMedPubMed CentralGoogle Scholar
- Schulz GE: The structure of bacterial outer membrane proteins. Biochim Biophys Acta 2002, 1565(2):308–317. 10.1016/S0005-2736(02)00577-1View ArticlePubMedGoogle Scholar
- Molloy MP, Herbert BR, Slade MB, Rabilloud T, Nouwens AS, Williams KL, Gooley AA: Proteomic analysis of the Escherichia coli outer membrane. Eur J Biochem 2000, 267(10):2871–2881. 10.1046/j.1432-1327.2000.01296.xView ArticlePubMedGoogle Scholar
- Phadke ND, Molloy MP, Steinhoff SA, Ulintz PJ, Andrews PC, Maddock JR: Analysis of the outer membrane proteome of Caulobacter crescentus by two-dimensional electrophoresis and mass spectrometry. Proteomics 2001, 1(5):705–720. Publisher Full Text 10.1002/1615-9861(200104)1:5<705::AID-PROT705>3.0.CO;2-NView ArticlePubMedGoogle Scholar
- Molloy MP, Phadke ND, Chen H, Tyldesley R, Garfin DE, Maddock JR, Andrews PC: Profiling the alkaline membrane proteome of Caulobacter crescentus with two-dimensional electrophoresis and mass spectrometry. Proteomics 2002, 2(7):899–910. 10.1002/1615-9861(200207)2:7<899::AID-PROT899>3.0.CO;2-YView ArticlePubMedGoogle Scholar
- Nouwens AS, Cordwell SJ, Larsen MR, Molloy MP, Gillings M, Willcox MD, Walsh BJ: Complementing genomics with proteomics: the membrane subproteome of Pseudomonas aeruginosa PAO1. Electrophoresis 2000, 21(17):3797–3809. 10.1002/1522-2683(200011)21:17<3797::AID-ELPS3797>3.0.CO;2-PView ArticlePubMedGoogle Scholar
- Rhomberg TA, Karlberg O, Mini T, Zimny-Arndt U, Wickenberg U, Rottgen M, Jungblut PR, Jeno P, Andersson SG, Dehio C: Proteomic analysis of the sarcosine-insoluble outer membrane fraction of the bacterial pathogen Bartonella henselae. Proteomics 2004, 4(10):3021–3033. 10.1002/pmic.200400933View ArticlePubMedGoogle Scholar
- Sabarth N, Lamer S, Zimny-Arndt U, Jungblut PR, Meyer TF, Bumann D: Identification of surface proteins of Helicobacter pylori by selective biotinylation, affinity purification, and two-dimensional gel electrophoresis. J Biol Chem 2002, 277(31):27896–27902. 10.1074/jbc.M204473200View ArticlePubMedGoogle Scholar
- Delcher AL, Bratke KA, Powers EC, Salzberg SL: Identifying bacterial genes and endosymbiont DNA with Glimmer. Bioinformatics 2007, 23(6):673–679. 10.1093/bioinformatics/btm009View ArticlePubMedPubMed CentralGoogle Scholar
- Delcher AL, Harmon D, Kasif S, White O, Salzberg SL: Improved microbial gene identification with GLIMMER. Nucleic Acids Res 1999, 27(23):4636–4641. 10.1093/nar/27.23.4636View ArticlePubMedPubMed CentralGoogle Scholar
- Salzberg SL, Delcher AL, Kasif S, White O: Microbial gene identification using interpolated Markov models. Nucleic Acids Res 1998, 26(2):544–548. 10.1093/nar/26.2.544View ArticlePubMedPubMed CentralGoogle 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
- Sonnhammer EL, von Heijne G, Krogh A: A hidden Markov model for predicting transmembrane helices in protein sequences. Proc Int Conf Intell Syst Mol Biol 1998, 6: 175–182.PubMedGoogle 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.04625404View ArticlePubMedPubMed CentralGoogle Scholar
- Juncker AS, Willenbrock H, Von Heijne G, Brunak S, Nielsen H, Krogh A: Prediction of lipoprotein signal peptides in Gram-negative bacteria. Protein Sci 2003, 12(8):1652–1662. 10.1110/ps.0303703View ArticlePubMedPubMed CentralGoogle Scholar
- Gardy JL, Laird MR, Chen F, Rey S, Walsh CJ, Ester M, Brinkman FS: PSORTb v.2.0: expanded prediction of bacterial protein subcellular localization and insights gained from comparative proteome analysis. Bioinformatics 2005, 21(5):617–623. 10.1093/bioinformatics/bti057View ArticlePubMedGoogle Scholar
- Gardy JL, Spencer C, Wang K, Ester M, Tusnady GE, Simon I, Hua S, deFays K, Lambert C, Nakai K, Brinkman FS: PSORT-B: Improving protein subcellular localization prediction for Gram-negative bacteria. Nucleic Acids Res 2003, 31(13):3613–3617. 10.1093/nar/gkg602View ArticlePubMedPubMed CentralGoogle Scholar
- Kall L, Krogh A, Sonnhammer EL: An HMM posterior decoder for sequence feature prediction that includes homology information. Bioinformatics 2005, 21(Suppl 1):i251–257. 10.1093/bioinformatics/bti1014View ArticlePubMedGoogle Scholar
- Kall L, Krogh A, Sonnhammer EL: Advantages of combined transmembrane topology and signal peptide prediction--the Phobius web server. Nucleic Acids Res 2007, (35 Web Server):W429–432. 10.1093/nar/gkm256
- Kall 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
- Bagos PG, Tsirigos KD, Liakopoulos TD, Hamodrakas SJ: Prediction of lipoprotein signal peptides in Gram-positive bacteria with a Hidden Markov Model. J Proteome Res 2008, 7(12):5082–5093. 10.1021/pr800162cView ArticlePubMedGoogle Scholar
- Bigelow HR, Petrey DS, Liu J, Przybylski D, Rost B: Predicting transmembrane beta-barrels in proteomes. Nucleic Acids Res 2004, 32(8):2566–2577. 10.1093/nar/gkh580View ArticlePubMedPubMed CentralGoogle Scholar
- Bigelow H, Rost B: PROFtmb: a web server for predicting bacterial transmembrane beta barrel proteins. Nucleic Acids Res 2006, (34 Web Server):W186–188. 10.1093/nar/gkl262Google Scholar
- Billion A, Ghai R, Chakraborty T, Hain T: Augur-a computational pipeline for whole genome microbial surface protein prediction and classification. Bioinformatics 2006, 22(22):2819–2820. 10.1093/bioinformatics/btl466View 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.