- Methodology article
- Open Access
LocateP: Genome-scale subcellular-location predictor for bacterial proteins
BMC Bioinformatics volume 9, Article number: 173 (2008)
In the past decades, various protein subcellular-location (SCL) predictors have been developed. Most of these predictors, like TMHMM 2.0, SignalP 3.0, PrediSi and Phobius, aim at the identification of one or a few SCLs, whereas others such as CELLO and Psortb.v.2.0 aim at a broader classification. Although these tools and pipelines can achieve a high precision in the accurate prediction of signal peptides and transmembrane helices, they have a much lower accuracy when other sequence characteristics are concerned. For instance, it proved notoriously difficult to identify the fate of proteins carrying a putative type I signal peptidase (SPIase) cleavage site, as many of those proteins are retained in the cell membrane as N-terminally anchored membrane proteins. Moreover, most of the SCL classifiers are based on the classification of the Swiss-Prot database and consequently inherited the inconsistency of that SCL classification. As accurate and detailed SCL prediction on a genome scale is highly desired by experimental researchers, we decided to construct a new SCL prediction pipeline: LocateP.
LocateP combines many of the existing high-precision SCL identifiers with our own newly developed identifiers for specific SCLs. The LocateP pipeline was designed such that it mimics protein targeting and secretion processes. It distinguishes 7 different SCLs within Gram-positive bacteria: intracellular, multi-transmembrane, N-terminally membrane anchored, C-terminally membrane anchored, lipid-anchored, LPxTG-type cell-wall anchored, and secreted/released proteins. Moreover, it distinguishes pathways for Sec- or Tat-dependent secretion and alternative secretion of bacteriocin-like proteins. The pipeline was tested on data sets extracted from literature, including experimental proteomics studies. The tests showed that LocateP performs as well as, or even slightly better than other SCL predictors for some locations and outperforms current tools especially where the N-terminally anchored and the SPIase-cleaved secreted proteins are concerned. Overall, the accuracy of LocateP was always higher than 90%. LocateP was then used to predict the SCLs of all proteins encoded by completed Gram-positive bacterial genomes. The results are stored in the database LocateP-DB http://www.cmbi.ru.nl/locatep-db.
LocateP is by far the most accurate and detailed protein SCL predictor for Gram-positive bacteria currently available.
In bacteria, secreted proteins are involved in stress sensing, substrate binding, cell communication, microbe-host interaction, adhesion, and other essential processes relevant to the environment and life style of the organisms. The secreted proteins are exported via various mechanisms and are retained by the bacterial cell via various interactions or released to the medium (Figure 1A). To identify the "secretome"  on a genome scale, subcellular proteomic studies have been carried out [3–7]. Although these experimental methods have contributed greatly to our knowledge of the subcellular location (SCL) of a variety of proteins, until now their scope has remained limited. In contrast, high-throughput computational methods for prediction of SCL sequence characteristics can be easily applied to every species whose genome has been sequenced.
Computational methods have gained considerable precision in the past decades. Initial tools focussed on detecting the presence, type and location of protein transmembrane segments, including signal peptides for targeting and translocation of proteins. One of the very first SCL prediction methods was introduced by Kyte and Doolittle  in 1982 with their amino acid hydropathy index. Since the late 90's machine-learning methods became more prominent, including neural networks [9–11], hidden Markov models (HMM) [12–14], support vector machines [2, 15–23], Bayesian networks [24, 25], and combined algorithms [26–33]. Moreover, present studies tend to combine different resources and methods [34–37]. For example, Chou et al.  combined gene ontology and functional domain databases, Shatkay et al.  combined text search and sequence data, and Marcotte et al.  combined protein homology and phylogenetic profiles in their studies.
Unfortunately, as a result of the trade-off between specificity and accuracy, computational methods will always be prone to error. Moreover, the number of false predictions increases even further when the SCL-related sequence characteristics have not been properly identified. For instance, among the Sec-dependent exported proteins, current predictors have severe difficulties to distinguish the proteins that are cleaved from the cell membrane by the type I signal peptidase (SPIase) – in this paper we will refer to these proteins as "secreted" – from a relatively large group of membrane-anchored proteins that also contain a putative SPIase-cleavage site but are not cleaved – in this paper we will lump these proteins in the category "N-anchored" [41–43].
As knowledge on the precise SCL of a protein is especially important to judge the biological nature and role of its activity, we constructed a new SCL prediction pipeline called LocateP. Our pipeline is geared to identify the detailed SCL of bacterial proteins by combining existing and novel prediction tools. Special effort was made to increase the accuracy of the prediction of N-anchored proteins. The version of LocateP presented here focuses on SCL prediction of proteins from Gram-positive bacteria.
The construction of the SCL-prediction pipeline LocateP
A major drawback of most current sub-cellular location (SCL) predictors is that they are not aimed at the prediction of very specific SCLs but merely at the rather broad locations intracellular, membrane bound/associated and extracellular, in line with the Swiss Prot classification system. We therefore constructed a SCL predictor pipeline LocateP, that distinguishes 7 SCLs and 3 targeting pathways that can be identified in Gram-positive bacteria, with a focus on extracellular SCLs (see Figure 1A).
The LocateP pipeline was designed such that it mimics the protein secretion process in Gram-positive bacteria. The pipeline structure can be categorized as follows: (1) secretion pathway prediction, (2) transmembrane-segment detection, (3) signal peptide identification, and (4) cleavage and retention signal recognition. The LocateP pipeline employs existing SCL prediction tools (Table 1) as well as our own new and more accurate methods for the prediction of lipoproteins, Tat-secreted, N-terminally anchored, C-terminally anchored and secreted proteins (see Methods). LocateP uses at least 2 prediction methods for each SCL, in order to increase prediction accuracy. The selection criteria imposed on these methods were derived from literature. The LocateP pipeline is depicted in Figure 2; its construction is described in more detail in the "Methods" section and in the legend of Figure 2. A detailed flow chart is presented in Additional file 1.
Making the distinction between N-anchored and secreted proteins containing a SPI-cleavage site
In the past, the sequence corresponding to the signal peptide has been subdivided into three distinct regions: the N, H and C regions [28, 44–46] (Figure 1B). Most of the membrane proteins with a single N-terminal TM anchor are easily identified as they do not have a predicted cleavage site for signal peptidases. However, as mentioned above, the prediction of SCL of proteins containing a putative signal peptidase type I (SPIase) cleavage site appears particularly difficult for current SCL predictors. Although many Sec-exported proteins are cleaved by the SPIase, a considerable number of proteins is not cleaved and remains membrane-anchored via the N-terminus .
To identify the features that determine cleavage, the multiple sequence alignments of the signal peptides from experimentally validated N-anchored and secreted proteins  containing a putative SPI cleavage site in Bacillus subtilis were analyzed. To enhance the signal, orthologous sequences from other Bacilli were added in the analysis (see Materials and Methods). The Weblogos  of the two collections of sequences are given in Figure 3A. No distinguishing pattern could be detected by eye. Therefore, a series of HMMs were constructed based on the sequence alignments of the N-anchored and secreted proteins. Nine pairs of HMMs were built for sequences surrounding the putative SPI cleavage site. Different numbers of residues on either side of the putative cleavage site were included in the models in order to investigate the roles of the H-region and the C-region in cleavage-site recognition. When the HMM pairs were applied to the two respective sets of sequences, it appeared that the HMM pair containing an equal number of residues on either side of the putative cleavage site performed best in predicting correctly whether the cleavage site was genuine or not (Figure 3B). The individual HMMs were not mutually exclusive for either of the two sets of sequences. However, when the scoring of the two HMMs of the pair was combined into a scoring matrix, the experimentally determined non-cleaved and cleaved sequences could be distinguished almost perfectly. The scoring matrix was included in LocateP; details of the matrix are described in the legend of Figure 3 and in the "Methods" section.
Initial validation of LocateP
Ideally, the performance of the LocateP pipeline should be checked with large experimentally validated data sets. Unfortunately, the availability of such large data sets is rather limited. LocateP was tested first with the experimental data set of Tjalsma et al.  which was used to create the HMM pair that distinguishes the N-anchored and secreted proteins containing a putative SPI-cleavage site. LocateP was able to distinguish these proteins with an accuracy of >90%. Then the performance of LocateP was tested on ten other data sets. These sets were extracted from literature describing other SCL prediction tools. LocateP performed extremely well on these sets, as the prediction accuracy was always higher than 90% (Table 2).
A second check was done with data collected from TransportDB . Based on expert knowledge on the composition and location of various transport systems and their functional components, the SCL prediction of 1336 transport-related proteins from Bacillus subtilis 168, Bacillus cereus ATCC14579 and Lactobacillus plantarum WCFS1 was verified (Table 3). For a difficult group like the substrate-binding proteins of ABC transport systems, LocateP identified 113 of 124 proteins in a correct SCL for substrate binding (96 lipoproteins, 6 secreted and 11 N-anchored proteins) [49–59]. For the other groups of transport-related proteins, the predicted SCL fitted the biological role of the transport proteins in ~98% of the cases.
A third less quantitative check included a comparison of the LocateP predictions for all N-anchored and secreted proteins of the Bacillus subtilis genome with their NCBI functional annotations (Additional file 2, 3 and 4). Nearly all of the predictions appeared to make biological sense according to literature: most of the predicted N-anchored proteins were annotated to be involved in processes that are related to the cell-envelope, such as cell division, transport, cell-envelope biogenesis, mobility, competence, signal transduction, protein turnover, etc; most predicted secreted proteins were indeed known to be secreted enzymes such as extracellular carbohydrases , alkaline phosphatases [61, 62], metalloproteases , neutral proteases, and subtilisin-family proteases .
Further validation and comparison of LocateP with other tools and pipelines
Recently, Gardy et al.  have compared most of the current SCL classifiers, and some tools showed excellent performance. We compared the performance of LocateP to a selection of these tools, including the individual SCL predictors and other general integrative SCL classifiers that were considered best (Table 2).
On N-anchored, secreted, lipid-anchored and multi-transmembrane protein prediction
LocateP and several individual SCL prediction tools were applied to the same collection of reference data sets. LocateP showed similar or higher recall to PrediSi, Phobius and SignalP 3.0 at signal peptide detection, respectively (Table 2, test sets 1, 2, 3, 4). LocateP performed clearly better than all other tools at predicting lipoproteins and multi-transmembrane proteins (Table 2, test sets 6, 7, 8). For the group of N-anchored proteins, LocateP clearly outperformed all other tools with a much higher prediction specificity and accuracy (Table 2, test set 5).
As has been noted by others [12, 14, 42], the N-anchored membrane proteins form an ambiguous group with respect to the location of their biological activity, i.e. outside or inside the cell. Various N-anchored proteins are actually active at the cytoplasmic side of the bacterial cell membrane [65–75]. Due to the lack of reliable distinguishing algorithms and experimental data, no reliable prediction methods for these "outside-in" proteins are available yet [12, 29, 43, 76–78]. As a result, in the current version of LocateP, proteins are only annotated as "N-terminally anchored"; most are presumed to function outside the cell, while some might have intracellular activity. A few of the known intracellular cases are indicated in Additional file 2.
On Tat-secreted protein prediction
Recent research pointed out that the Tat-export pathway plays an important role as a parallel protein secretion pathway to the Sec-pathway in some Gram-positive organisms [79–85]. Unfortunately, Sec-signal peptide detectors have a high false-negative prediction rate on Tat-substrates . Therefore, we considered it necessary to include a Tat-secreted protein prediction tool in the LocateP pipeline, and we combined two newly created Tat-secreted protein-specific HMMs (see Methods) with Tat-find v.1.2 for the SCL prediction of these proteins. Tat-signal peptides are known to have an almost invariable double Arg or Lys+Arg motif (RR-motif) [87–90] upstream of the transmembrane segment. It appeared important that the Tat-secreted protein predictors can discriminate the Tat-signal peptides from sequences (especially transmembrane helices) that contain consecutive positively charged residues.
We compared the performance of LocateP, TatP and Tat-find v1.2 on the proteins containing a RR/RK pattern in their N-terminus (test sets 3, 4 and 9). LocateP clearly performed better than the other two specific tools when tested with intracellular and membrane proteins sets, and thus showed an excellent capability of Tat-signal peptide detection (Table 2). Moreover, it appeared that TatP and Tat-find v1.2 predicted several proteins to be secreted via the Tat-pathway in 22 species that apparently lack the relevant pathway genes , whereas LocateP did not find any Tat-pathway substrates in those species. Thus, LocateP showed the best overall performance among the Tat-pathway prediction tools for gram positive bacteria.
Comparing LocateP and other integrative SCL classifiers
According to the comparative study of Gardy et al. , CELLO  is one of the best SCL classification pipelines. We therefore evaluated the performance of LocateP as an integrative SCL classifier by comparing it to CELLO and the widely used pipeline Psortb.v.2.0 . Other pipelines like SubLoc , LOCtree , Proteome Analyst , P-CLASSIFIER  and PSLpred  were not selected because they either do not provide prediction of membrane proteins, or are tailored for Gram-negative bacteria, or in the best case showed similar performance to Psortb.v.2.0 or CELLO. Recently, a SCL prediction tool called Gpos-Ploc  was published that classifies Gram-positive proteins. LocateP was not compared to Gpos-Ploc because its web server accepts only one sequence per search. Moreover, the overall accuracy of the tool is reported to be only ~85% .
LocateP had an accuracy lower than CELLO (Table 2) when tested with data extracted from the Swiss-Prot database (test set 10 ). However, when compared using experimental data (test sets 3, 4, 5, 6, 7), CELLO and Psortb v2.0 showed dramatically poor prediction rates (Table 2). This poor performance relates to the fact that the training data of CELLO and Psortb v.2.0 were from the Swiss-Prot database (i.e., part of test set 10 and 11). This database does not distinguish between N-anchored, secreted and lipoproteins, and at the same time the members of these groups are distributed over two general classes: "membrane" and "extracellular". Thus, in essence the poor performance of CELLO and Psortv.b.2.0 is a consequence of the less-specific classification in Swiss-Prot (Table 2). Vice versa, the lower accuracy of LocateP on the Swiss-Prot data is related to the inconsistency in the classification.
TBPred  is a SCL classifier that was especially designed for mycobacteria, based on the idea that organism-specific methods might have higher accuracy [96, 97]. We compared LocateP with TBPred using the training data of TBPred (test set 11). Surprisingly, LocateP showed considerably higher accuracy than TBPred, especially on lipoprotein and secreted protein prediction, even though no mycobacterial proteins were involved in the lipoprotein prediction training process of LocateP.
Finally, the performance of LocateP was compared to Augur , a computational pipeline that also combines many existing tools. Augur detects signal peptides and transmembrane helices using only SignalP and TMHMM, and consequently the accuracy of N-anchored protein prediction of Augur is much lower than with LocateP. Augur also falsely predicted 8 lipoproteins out of a test-set of 114 non-lipoproteins (test sets 4 and 5), which implied a higher false-positive rate than LocateP on lipoprotein prediction.
Comparative analysis of protein subcellular location in Gram-positive bacteria
LocateP was applied to the encoded proteins of all complete Gram-positive bacterial genomes available in the NCBI database. The average distribution of proteins grouped by predicted SCL was calculated for each genome. Despite the different genome sizes, Gram-positive bacteria tend to have a similar distribution of proteins over certain SCLs independent of class or family, and this independency also holds for individual Gram-positive bacterial genomes (Table 4). We note that the fractions of intracellular and membrane proteins predicted by LocateP in Gram-positive genomes were consistent with what was previously estimated by other tools [12, 18, 24, 25]. The complete genome predictions can be viewed in our database LocateP-DB .
Although the early SCL-prediction tools performed rather poorly, current tools perform rather well on specific categories of signal-peptide containing proteins and membrane proteins , reaching an accuracy of 96%. Nevertheless, for other groups like secreted, N-anchored and lipoproteins these tools still perform rather poorly. As the latter groups represent a considerable part of the secretome, we decided to design a new SCL-identification pipeline called LocateP.
The performance of LocateP was checked against the best current tools and it outperformed all of them, particularly when difficult groups of proteins and SCLs were concerned. The outstanding performance was achieved though the generation of specific HMMs based on protein sequences whose cellular fate had been experimentally tested. For instance, it has long been a problem to identify secreted and N-anchored proteins from the group of proteins carrying a putative SPI-cleavage site motif. Formerly, the H-region together with the cleavage site were considered to be the key elements of SPIase-substrate recognition. Therefore, previous signal-peptide predictors were constructed focusing on the H-region and/or on the cleavage site [14, 18, 26, 98–100]. However, Carlos et al.  found that the H-region of the SPIase substrate was not critical for peptidase-cleavage capability but that, in contrast, mutations in the C-region of originally non-cleaved proteins caused alternative cleavage. They therefore claimed that specific substrate-enzyme interactions around the C-region should be decisive for SPIase-cleavage site recognition. Indeed, our analysis of the signal sequences of a group of secreted and N-anchored proteins indicated that the C-region is important, but that at the same time also the H-region carries properties that determine the protein's fate (i.e. to be or not to be cleaved). The fact that the performance of the dedicated HMMs became worse when the sequence was extended beyond 30 residues implies that the decisive information is present in this stretch of sequence. LocateP improved the separation of N-anchored and secreted proteins from ~40% (by Phobius ) to > 90% without disturbing the SCL prediction of the other types of proteins.
LocateP was designed as a pipeline, and hence could have performed less well on specific categories than specialized tools. In particular, the performance would have been considerably lower if the flow scheme had been chosen wrongly. However, a comparison of the performance on lipoproteins, membrane proteins, Sec-secreted and Tat-secreted proteins with the specialized tools LipoP 1.0, TMHMM 2.0, Phobius, SignalP 3.0, Predisi, Tat-find v.1.2 and TatP shows that LocateP does not suffer from being a pipeline tool. Apparently, our choice to mimic the order in the bacterial secretion process was a correct one. In fact, it has been shown by others that the SCL prediction can be improved considerably by simulating the protein sorting processes [93, 102]. Overall, LocateP performed very well, with an accuracy higher than 95% for nearly all categories, and only slightly lower in one case (91% for N-anchored proteins), but still considerably better than all other tools. LocateP could be used to distinguish 7 SCLs and 3 sorting pathways and avoided the inconsistent SCL classification which most SCL classifiers inevitably inherited from Swiss-Prot.
Because of the high prediction accuracy of LocateP on proteins of known biological function (see e.g. Additional file 2), we expect that the SCL prediction of proteins of unknown function should also be equally reliable. In principle, the genome-scale SCL predictions made by LocateP provide an excellent starting point for functional annotation and experimental analysis of encoded proteins of unknown function, as they provide numerous clues about where to look for a certain biological activity.
Although LocateP already performs quite well, there is inevitably room for improvement. For instance, in the Swiss-Prot database, many of the annotated cell-wall proteins are secreted proteins bound to the cell surface via non-covalent interactions. Known elements of non-covalent binding include choline-binding domains, LysM domains, type 2 cell-wall binding domains, GW-modules, Lysin-binding motifs, ChW-binding motifs, WxL domains, LPP-region binding, S-layer proteins, and others [103–117]. The current version of LocateP was designed to predict only the covalent cell-wall (peptidoglycan) binding mechanism of proteins by dedicated sortases. For instance, among the 14 cell-wall proteins in test set 10d , 13 are non-covalently cell-wall bound secreted proteins. LocateP correctly predicted 10 of these as "secreted proteins", but does not allow for the fact that these proteins could be cell-wall bound via non-covalent mechanisms after secretion (Table 2). Future versions of LocateP will include non-covalent binding to the cell wall.
Not all mechanisms of protein secretion or modification are known to date and not all have been included in the LocateP pipeline yet [118–123]. This is the case for proteins that have been shown to occur at various locations or those that are secreted via minor pathways. Examples are proteins that are either cleaved multi-domain proteins , auto-transporters found in both cytoplasmic and extracellular locations , or proteins with various SCL depending on growth phase and/or specific environment [124, 126]. The multi-compartment proteins and minor-pathway secreted proteins appear to be rare in most bacteria, and their sorting mechanisms are not completely understood yet. Therefore, the current version of LocateP predicts only one SCL for such proteins, which may be only partially correct. In contrast, Psortb v2.0 and CELLO were claimed to be capable of multi-location prediction [126, 127]. Both tools employ machine-learning methods and the predicted multiple locations should represent a certain statistical significance even without large-scale experimental evidence. However, both tools inevitably generate a considerable number of false positives. Similarly, SecretomeP 2.0 , which was made to predict non-classically secreted proteins, was not included in LocateP because of its high false-prediction rate. The recent predictors Euk-mPloc  and Hum-mPloc  incorporated up-to-date Eukaryotic proteins that were found to have multiple compartments and the tools achieved rather satisfying accuracies. Similar tools will be included or constructed for LocateP when more experimental data on multiple locations of bacterial proteins are available.
Another group of proteins that is not treated separately by LocateP is the group that is exported by unknown mechanisms and is known as the Gram-positive periplasmic proteins [130–133]. Carlsson et al.  recently reported that in Gram-positive bacteria the secreted proteins could be directed to different extracellular regions including a periplasmic space. In fact, the prediction of a subcellular location "periplasmic" in Gram-positive was not included in any published SCL prediction tools for Gram-positive bacteria, except in Gpos-PLoc . However, the Gpos-PLoc prediction algorithm was based on only 5 proteins which were extracted from the Swiss-Prot database. Indeed, among these 5 proteins, four were expressed in the E. coli periplasmic space, but no evidence exists that they are also expressed in the periplasm of a Gram-positive organism [135–138]. Moreover, one protein (P29166) was proven to be located inside the cell . LocateP predicted correctly that 4 of the 5 proteins should be secreted and are located outside of the plasma membrane, while P29166 was predicted to be cytoplasmic, in line with the experimental evidence.
LocateP was first tailored for the SCL prediction of Gram-positive bacterial proteins; therefore prediction of the Gram-negative specific proteins, such as β-barrel membrane proteins, was not yet included in the pipeline. LocateP was compared to SigTree , a signal-peptide detector based on sets of experimentally verified E. coli proteins, using the same data set from E. coli (data not shown). LocateP showed slightly lower accuracy than SigTree did, which suggests that the sequence composition of signal peptides from Gram-positive and Gram-negative bacterial proteins could be different. Future versions of LocateP will be improved and extended to Gram-negative bacterial protein SCL prediction by incorporating Gram-negative specific subcellular-location prediction tools.
Finally, we must emphasize that in several cases an automatic SCL prediction of a protein will inevitably give an incorrect prediction using LocateP or any other tool: (i) when the start codon of a gene encoding a protein with a signal peptide has been wrongly identified (either too far upstream or too far downstream), (ii) when a frame shift in the open-reading frame leads to different fragments of encoded proteins, and (iii) when an intracellular protein contains a signal peptide-like hydrophobic helix near the N-terminus; in this case such helices generally fold into the interior of the globular protein [141–144].
As detailed and accurate genome-scale SCL prediction of encoded proteins is highly desired by scientists in various biological research areas, numerous existing and newly developed tools were combined into one pipeline: LocateP. To date, LocateP is the most detailed protein SCL predictor for Gram-positive bacterial proteins among all tools that have been reported, in that it presently distinguishes 7 different SCLs and 3 sorting pathways, with focus on extracellular SCLs. Moreover, it is also the most accurate SCL predictor, especially on distinguishing N-anchored and secreted proteins. LocateP was applied on all completed Gram-positive bacterial genomes from the NCBI sequence database. The results are updated synchronously with Genbank updates and are publicly available via the database LocateP-DB . The present version contains SCL predictions for 436,771 proteins in 148 genomes of Gram-positive bacteria. These genome-scale SCL predictions provide an excellent starting point for experimentalists to improve the functional annotation of proteins.
Sources of sequence information and location data
The genome sequences of Gram-positive bacteria were extracted from GenBank on May 30th, 2007 ("ORGANISM" annotation fields: 'Firmicutes' or 'Actinobacteria'), and were continuously updated since then. Protein sequences of Bacillus species were collected both from GenBank on April 1st, 2007 and from the ERGO database  on November 15th, 2006.
Eight different protein data sets of known subcellular location were selected from literature describing other tools and describing proteome studies (Table 2 legend). In order to check the performance of the LocateP pipeline, the SCL predictions were checked against an expert evaluation of the functional location of transport-related proteins from several Gram-positive bacterial genomes in TransportDB  on May 30th, 2007, and an expert evaluation against the protein function annotation as retrieved from GenBank on August, 10th, 2007.
Sequence analysis and evaluation of performance
Multiple sequence alignments were built with MUSCLE . HMMs were built with HMMER . Wherever appropriate, HMMs of varying length and different regions of the aligned N-terminal sequences of proteins were made, and the HMM that performed best was selected. Performance was evaluated using the statistical measure recall (or sensitivity) which is the number of true positives divided by the sum of the true positives and the false negatives.
Bioinformatics tools included in the LocateP pipeline
Many studies have compared and evaluated currently available transmembrane segment and signal peptide predictors [41, 44, 76, 98–100, 148]. Based on those studies and our own preliminary trials the following tools were selected to be included in our SCL prediction pipeline LocateP: TMHMM 2.0 , Phobius , SignalP 3.0 , PrediSi , and Bagel  (Table 1). Of these, TMHMM 2.0 and SignalP 3.0 are the most popular ones in the field; Phobius was selected for its high specificity on transmembrane segment identification; PrediSi was selected because it was trained with comparatively recent experimental data, and because it slightly outperformed SignalP 3.0 when applied to Gram-positive bacterial proteins . We also included the predictor Bagel for non-classically secreted bacteriocin-like proteins . The membrane protein predictor MemType-2L  includes topology prediction of N-anchored proteins but showed rather low accuracy with our experimental datasets; therefore this tool was not included in LocateP. Some other tools were not incorporated either because of a high false-prediction rate (e.g. HMMTOP  and SecretomeP 2.0 ), a low specificity for Gram-positive bacteria (e.g. LipoP 1.0 ), or simply the lack of stand-alone installable software packages (e.g. TatP , Signal-3L, Signal-CF and Tat-pred ).
Signal peptide detection
LocateP detects signal peptides by scanning the protein N-terminus, which was defined as the initial 60 amino acids of the protein, using SignalP 3.0, Phobius and PrediSi. Some proteins have a signal peptide shortly after these 60 amino acids. These proteins were predicted as "intracellular", but we added the extra remark of "TMH start AFTER 60" to the annotation indicating that these proteins could be secreted. No attempt was made to choose alternative start codons of incorrectly predicted start sites of ORFs.
Specific HMMs to determine the SCL of proteins with a putative SPI-cleavage site
Recently, Tjalsma et al. have experimentally determined the SCL of a large number of Bacillus subtilis proteins . The experimental set contained 66 proteins with a putative SPI-cleavage site. Of these 36 appeared to be cleaved and thus secreted, whereas 30 were shown to remain N-anchored. We named these sets "EXP-secreted" and "EXP-anchored", respectively, and used them to construct set-specific HMMs. To enhance the inherent signal, both sets were expanded by adding orthologous sequences from other Bacilli. First, homologs were searched with BLASTP  in the ERGO genome database  using full-length sequences. Only the three best BLAST hits were considered orthologs, when they also showed conserved gene context and functional annotation, high similarity and similar protein length. In this way, after removing orthologs containing identical N-terminal sequences, 27 secreted and 23 N-anchored orthologs could be added to the "EXP-" sets.
Pairs of HMMs were built to separate the group of proteins with a putative SPI cleavage-site into those that are cleaved (i.e. secreted/released) and those that are not (i.e. N-anchored). The sequences were aligned around the putative cleavage site and the length of the HMMs was varied (length >8). All HMM pairs were applied to both "EXP-" sets; the E-value was set at 10,000 to assure each protein gained an HMM score. For each pair the separation between truly cleaved and truly N-anchored proteins was analyzed and it appeared that the HMM pair containing equivalent amounts of H and C region residues achieved the highest specificity in distinguishing the two sets (see Figure 3B). The most specific HMM pair had a length of 25 amino acids and ran from residue -14 to +10 relative to the cleavage site Alanine (see Figure 3A).
The individual HMMs of the selected pair (HMMnon-cleaved, HMMcleaved) each displayed a relatively high specificity, but this was increased significantly by combining the two HMMs. A generic scoring scheme was derived via the following procedure: i) The HMM scores were rounded to discrete integers and the score distribution for the EXP-anchored and EXP-secreted protein sets was used to determine a first cut-off. The discrete HMM scores related to the HMMnon-cleaved ranged from -19 to +20 with all non-cleaved (i.e. N-anchored) proteins scoring higher then 3, those related to the HMMcleaved ranged from -29 to +20 with all cleaved (i.e. secreted) proteins scoring higher then 0. In fact, for both HMMs only in a small scoring range the two protein groups overlapped. Therefore, all sequences with a score ≤2 using the HMMnon-cleaved were attributed the SCL: SEC-secreted, and those with a score ≤-1 using the HMMcleaved were attributed the SCL: N-anchored. ii) For those sequences that scored >2 with the HMMnon-cleaved and >-1 with the HMMcleaved, the score with both models was compared. In case HMMcleaved > HMMnon-leaved score, the sequence was considered SEC-secreted, whereas, in case HMMnon-cleaved ≥ HMMcleaved score, the sequence was considered N-anchored.
The creation and selection of a specific HMM for lipoprotein prediction
The experimental data of Tjalsma et al.  indicated that at least 42 distinct proteins of Bacillus subtilis are lipoproteins. This set of proteins was taken and expanded with orthologs from 18 closely related Bacillus species using an Inparanoid  search for best bi-directional hits. After removing the sequences which contain identical initial 50 residues, 219 putative orthologous lipoproteins could be added. As all lipoproteins are anchored to the cell membrane by thioether linkage of the conserved lipobox cysteine to a diglyceride [41, 56, 151, 157, 158], the sequences were aligned around the lipobox. Eight HMMs were built based on different N-terminal regions from these proteins varying in length between 5 and 30 residues. Each HMM was applied to the original dataset of Tjalsma et al. and the performance was evaluated. The HMM with a length of 21 residues (-20 residues to the lipobox Cysteine) showed the highest specificity when the T-score was set to 3. Gaps were allowed in this model except in the region of the lipobox (residue -5 to the lipobox Cysteine)
The creation and selection of a specific HMM for Tat-secreted protein prediction
The 105 putative Tat-secreted proteins (according to Swiss-Prot) from the TatP-positive training set  were taken as the initial set for generating Tat-specific HMMs. The sequences were aligned either around the double Arg or Lys+Arg motif (RR-motif) [87–90] upstream of the transmembrane helix or the putative AxA triplet cleavage site  predicted by TatP downstream of the transmembrane helix. Eleven HMMs with different lengths were generated. A combination of two HMMs was found to be most specific with the training data, together with the restriction of an E-value smaller than 10: one HMM contained 2 residues in front and 16 after the twin-arginine motif, and the other HMM contained 17 amino acids in front and 1 residue after the triplet cleavage site. Interestingly, these two HMMs partly overlapped each other by the transmembrane (H) region. According to Taylor et al. , the -3 to +7 residues surrounding the twin-arginine should be the most characteristic for Tat-secreted protein identification. This conclusion was reaffirmed by our HMM model. The current tools Tat-find v1.2 and TatP both focus solely on the twin-arginine motif and consecutive transmembrane helix detection. In the prediction of the Tat-secretion signal our HMMs were combined with the Tat-find v1.2 program (in a scoring matrix) and therefore more weight was given to the prediction of the twin-arginine motif and its following hydrophobic region as Tat identifiers. This combined method was tested with 22 independent experimentally verified Tat-secreted proteins (20 of them are from E. coli [159–162], while PhoD and YwbN were from Bacillus subtilis [163, 164]). The SCL of 20 of these proteins was correctly identified by LocateP, including PhoD and YwbN (these 2 proteins were not in the HMM training set). Using this procedure the false prediction rate was significantly decreased compared to Tat-find v1.2 and TatP (Table 2). It was suggested that the Tat-pathway in Gram-positive bacteria is structurally different from Gram-negative bacteria [80, 86, 90, 120, 163, 164]. Although the Tat-secreted prediction of LocateP outperformed current tools, this part of the tool was trained with Gram-negative bacterial proteins due to the lack of experimental data from Gram-positive bacteria (see above). In order to avoid potential errors, LocateP also scans all proteins assuming them to be Sec-secreted, except for the bacteriocin-like secreted proteins. If the Tat-secreted possibility score of a protein was significant, the final subcellular location of this protein was marked "Possibly Tat-" as an extra reference.
Specific criteria for LPxTG-anchored and C-anchored protein prediction
The following topological criteria were used to identify LPxTG-type cell-wall anchored and C-anchored membrane proteins. For the selection of LPxTG-anchored proteins, the criteria were : (i) the protein has only one N-terminal signal peptide/TM segment and only one C-terminal TM segment, (ii) the C-terminus of the protein contains an LPxTG-type motif; (iii) the LPxTG-type motif is followed by the C-terminal transmembrane helix and a positively charged C-terminal tail. These criteria were validated with 85 experimentally verified LPxTG-anchored proteins [166–171] and 83 of them were correctly identified.
The criteria used for predicting C-anchored proteins were: (i) the protein has only 2 predicted TM helices, one situated at the N-terminus and one at the C-terminus, (ii) the protein has a cleaved N-terminal signal peptide, (iii) the protein has a C-terminal transmembrane helix and a positively charged C-terminal tail but no LPxTG motif, (iv) the distance between the N-terminal and C-terminal helices is larger than 45 residues.
Tjalsma H, Bolhuis A, Jongbloed JD, Bron S, van Dijl JM: Signal peptide-dependent protein transport in Bacillus subtilis: a genome-based survey of the secretome. Microbiol Mol Biol Rev 2000, 64(3):515–547.
Huang F, Parmryd I, Nilsson F, Persson AL, Pakrasi HB, Andersson B, Norling B: Proteomics of Synechocystis sp. strain PCC 6803: identification of plasma membrane proteins. Mol Cell Proteomics 2002, 1(12):956–966.
Molloy MP, Phadke ND, Maddock JR, Andrews PC: Two-dimensional electrophoresis and peptide mass fingerprinting of bacterial outer membrane proteins. Electrophoresis 2001, 22(9):1686–1696.
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.
Murakami Y, Imai M, Nakamura H, Yoshimura F: Separation of the outer membrane and identification of major outer membrane proteins from Porphyromonas gingivalis. Eur J Oral Sci 2002, 110(2):157–162.
Bumann D, Aksu S, Wendland M, Janek K, Zimny-Arndt U, Sabarth N, Meyer TF, Jungblut PR: Proteome analysis of secreted proteins of the gastric pathogen Helicobacter pylori. Infect Immun 2002, 70(7):3396–3403.
Kyte J, Doolittle RF: A simple method for displaying the hydropathic character of a protein. J Mol Biol 1982, 157(1):105–132.
Rost B, Fariselli P, Casadio R: Topology prediction for helical transmembrane proteins at 86% accuracy. Protein Sci 1996, 5(8):1704–1718.
Sidhu A, Yang ZR: Prediction of signal peptides using bio-basis function neural networks and decision trees. Appl Bioinformatics 2006, 5(1):13–19.
Boden M, Hawkins J: Prediction of subcellular localization using sequence-biased recurrent networks. Bioinformatics 2005, 21(10):2279–2286.
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.
Tusnady GE, Simon I: The HMMTOP transmembrane topology prediction server. Bioinformatics 2001, 17(9):849–850.
Kall L, Krogh A, Sonnhammer EL: An HMM posterior decoder for sequence feature prediction that includes homology information. Bioinformatics 2005, 21 Suppl 1: i251–7.
Chou KC, Elrod DW: Prediction of membrane protein types and subcellular locations. Proteins 1999, 34(1):137–153.
Hua S, Sun Z: Support vector machine approach for protein subcellular localization prediction. Bioinformatics 2001, 17(8):721–728.
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.
Bendtsen JD, Nielsen H, von Heijne G, Brunak S: Improved prediction of signal peptides: SignalP 3.0. J Mol Biol 2004, 340(4):783–795.
Emanuelsson O: Predicting protein subcellular localisation from amino acid sequence information. Brief Bioinform 2002, 3(4):361–376.
Yu CS, Lin CJ, Hwang JK: Predicting subcellular localization of proteins for Gram-negative bacteria by support vector machines based on n-peptide compositions. Protein Sci 2004, 13(5):1402–1406.
Hoglund A, Donnes P, Blum T, Adolph HW, Kohlbacher O: MultiLoc: prediction of protein subcellular localization using N-terminal targeting sequences, sequence motifs and amino acid composition. Bioinformatics 2006, 22(10):1158–1165.
Su CY, Lo A, Chiu HS, Sung TY, Hsu WL: Protein subcellular localization prediction based on compartment-specific biological features. Comput Syst Bioinformatics Conf 2006, 325–330.
Lei Z, Dai Y: An SVM-based system for predicting protein subnuclear localizations. BMC Bioinformatics 2005, 6: 291.
Taylor PD, Attwood TK, Flower DR: Toward bacterial protein sub-cellular location prediction: single-class discrimminant models for all gram- and gram+ compartments. Bioinformation 2006, 1(8):276–280.
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.
Hawkins J, Boden M: Detecting and sorting targeting peptides with neural networks and support vector machines. J Bioinform Comput Biol 2006, 4(1):1–18.
Billion A, Ghai R, Chakraborty T, Hain T: Augur - a computational pipeline for whole genome microbial surface protein prediction and classification. Bioinformatics 2006.
Shen HB, Chou KC: Gpos-PLoc: an ensemble classifier for predicting subcellular localization of Gram-positive bacterial proteins. Protein Eng Des Sel 2007.
Zhou H, Zhou Y: Predicting the topology of transmembrane helical proteins using mean burial propensity and a hidden-Markov-model-based method. Protein Sci 2003, 12(7):1547–1555.
Lorena AC, de Carvalho AC: Protein cellular localization prediction with Support Vector Machines and Decision Trees. Comput Biol Med 2007, 37(2):115–125.
Lu Z, Szafron D, Greiner R, Lu P, Wishart DS, Poulin B, Anvik J, Macdonell C, Eisner R: Predicting subcellular localization of proteins using machine-learned classifiers. Bioinformatics 2004, 20(4):547–556.
Bulashevska A, Eils R: Predicting protein subcellular locations using hierarchical ensemble of Bayesian classifiers based on Markov chains. BMC Bioinformatics 2006, 7: 298.
Wang J, Sung WK, Krishnan A, Li KB: Protein subcellular localization prediction for Gram-negative bacteria using amino acid subalphabets and a combination of multiple support vector machines. BMC Bioinformatics 2005, 6: 174.
Doderer M, Yoon K, Salinas J, Kwek S: Protein subcellular localization prediction using a hybrid of similarity search and error-correcting output code techniques that produces interpretable results. In Silico Biol 2006, 6(5):419–433.
Guo J, Pu X, Lin Y, Leung H: Protein subcellular localization based on PSI-BLAST and machine learning. J Bioinform Comput Biol 2006, 4(6):1181–1195.
Bhasin M, Garg A, Raghava GP: PSLpred: prediction of subcellular localization of bacterial proteins. Bioinformatics 2005, 21(10):2522–2524.
Wang M, Li A, Xie D, Jiang Z, Feng H, Fan Z: Improving prediction of protein subcellular localization using evolutionary information and sequence-order information. Conf Proc IEEE Eng Med Biol Soc 2005, 4: 4434–4436.
Chou KC, Cai YD: A new hybrid approach to predict subcellular localization of proteins by incorporating gene ontology. Biochem Biophys Res Commun 2003, 311(3):743–747.
Shatkay H, Hoglund A, Brady S, Blum T, Donnes P, Kohlbacher O: SherLoc: high-accuracy prediction of protein subcellular localization by integrating text and protein sequence data. Bioinformatics 2007, 23(11):1410–1417.
Marcotte EM, Xenarios I, van Der Bliek AM, Eisenberg D: Localizing proteins in the cell from their phylogenetic profiles. Proc Natl Acad Sci U S A 2000, 97(22):12115–12120.
Tjalsma H, van Dijl JM: Proteomics-based consensus prediction of protein retention in a bacterial membrane. Proteomics 2005, 5(17):4472–4482.
Tjalsma H: Feature-based reappraisal of the Bacillus subtilis exoproteome. Proteomics 2007, 7(1):73–81.
Lao DM, Okuno T, Shimizu T: Evaluating transmembrane topology prediction methods for the effect of signal peptide in topology prediction. In Silico Biol 2002, 2(4):485–494.
Gardy JL, Brinkman FS: Methods for predicting bacterial protein subcellular localization. Nat Rev Microbiol 2006, 4(10):741–751.
Nilsson I, Witt S, Kiefer H, Mingarro I, von Heijne G: Distant downstream sequence determinants can control N-tail translocation during protein insertion into the endoplasmic reticulum membrane. J Biol Chem 2000, 275(9):6207–6213.
Emanuelsson O, Nielsen H, Brunak S, von Heijne G: Predicting subcellular localization of proteins based on their N-terminal amino acid sequence. J Mol Biol 2000, 300(4):1005–1016.
Crooks GE, Hon G, Chandonia JM, Brenner SE: WebLogo: a sequence logo generator. Genome Res 2004, 14(6):1188–1190.
Ren Q, Kang KH, Paulsen IT: TransportDB: a relational database of cellular membrane transport systems. Nucleic Acids Res 2004, 32(Database issue):D284–8.
Koide A, Hoch JA: Identification of a second oligopeptide transport system in Bacillus subtilis and determination of its role in sporulation. Mol Microbiol 1994, 13(3):417–426.
Sa-Nogueira I, Nogueira TV, Soares S, de Lencastre H: The Bacillus subtilis L-arabinose (ara) operon: nucleotide sequence, genetic organization and expression. Microbiology 1997, 143 ( Pt 3): 957–969.
Coppee JY, Auger S, Turlin E, Sekowska A, Le Caer JP, Labas V, Vagner V, Danchin A, Martin-Verstraete I: Sulfur-limitation-regulated proteins in Bacillus subtilis: a two-dimensional gel electrophoresis study. Microbiology 2001, 147(Pt 6):1631–1640.
Leskela S, Kontinen VP, Sarvas M: Molecular analysis of an operon in Bacillus subtilis encoding a novel ABC transporter with a role in exoprotein production, sporulation and competence. Microbiology 1996, 142 ( Pt 1): 71–77.
Miethke M, Klotz O, Linne U, May JJ, Beckering CL, Marahiel MA: Ferri-bacillibactin uptake and hydrolysis in Bacillus subtilis. Mol Microbiol 2006, 61(6):1413–1427.
Ollinger J, Song KB, Antelmann H, Hecker M, Helmann JD: Role of the Fur regulon in iron transport in Bacillus subtilis. J Bacteriol 2006, 188(10):3664–3673.
Sebulsky MT, Heinrichs DE: Identification and characterization of fhuD1 and fhuD2, two genes involved in iron-hydroxamate uptake in Staphylococcus aureus. J Bacteriol 2001, 183(17):4994–5000.
Schneider R, Hantke K: Iron-hydroxamate uptake systems in Bacillus subtilis: identification of a lipoprotein as part of a binding protein-dependent transport system. Mol Microbiol 1993, 8(1):111–121.
Gomez A, Ramon D, Sanz P: The Bacillus subtilis lipoprotein LplA causes cell lysis when expressed in Escherichia coli. Microbiology 1994, 140 ( Pt 8): 1839–1845.
Kim MS, Shin J, Lee W, Lee HS, Oh BH: Crystal structures of RbsD leading to the identification of cytoplasmic sugar-binding proteins with a novel folding architecture. J Biol Chem 2003, 278(30):28173–28180.
Stentz R, Zagorec M: Ribose utilization in Lactobacillus sakei: analysis of the regulation of the rbs operon and putative involvement of a new transporter. J Mol Microbiol Biotechnol 1999, 1(1):165–173.
Steinmetz M, Le Coq D, Aymerich S, Gonzy-Treboul G, Gay P: The DNA sequence of the gene for the secreted Bacillus subtilis enzyme levansucrase and its genetic control sites. Mol Gen Genet 1985, 200(2):220–228.
Kobayashi T, Hakamada Y, Adachi S, Hitomi J, Yoshimatsu T, Koike K, Kawai S, Ito S: Purification and properties of an alkaline protease from alkalophilic Bacillus sp. KSM-K16. Appl Microbiol Biotechnol 1995, 43(3):473–481.
Shirai T, Suzuki A, Yamane T, Ashida T, Kobayashi T, Hitomi J, Ito S: High-resolution crystal structure of M-protease: phylogeny aided analysis of the high-alkaline adaptation mechanism. Protein Eng 1997, 10(6):627–634.
Kunst F, Ogasawara N, Moszer I, Albertini AM, Alloni G, Azevedo V, Bertero MG, Bessieres P, Bolotin A, Borchert S, Borriss R, Boursier L, Brans A, Braun M, Brignell SC, Bron S, Brouillet S, Bruschi CV, Caldwell B, Capuano V, Carter NM, Choi SK, Codani JJ, Connerton IF, Danchin A, et al.: The complete genome sequence of the gram-positive bacterium Bacillus subtilis. Nature 1997, 390(6657):249–256.
Wu XC, Lee W, Tran L, Wong SL: Engineering a Bacillus subtilis expression-secretion system with a strain deficient in six extracellular proteases. J Bacteriol 1991, 173(16):4952–4958.
Helmann JD: Deciphering a complex genetic regulatory network: the Bacillus subtilis sigmaW protein and intrinsic resistance to antimicrobial compounds. Sci Prog 2006, 89(Pt 3–4):243–266.
Katis VL, Wake RG: Membrane-bound division proteins DivIB and DivIC of Bacillus subtilis function solely through their external domains in both vegetative and sporulation division. J Bacteriol 1999, 181(9):2710–2718.
Katis VL, Harry EJ, Wake RG: The Bacillus subtilis division protein DivIC is a highly abundant membrane-bound protein that localizes to the division site. Mol Microbiol 1997, 26(5):1047–1055.
Noirclerc-Savoye M, Le Gouellec A, Morlot C, Dideberg O, Vernet T, Zapun A: In vitro reconstitution of a trimeric complex of DivIB, DivIC and FtsL, and their transient co-localization at the division site in Streptococcus pneumoniae. Mol Microbiol 2005, 55(2):413–424.
Kramer N, Hahn J, Dubnau D: Multiple interactions among the competence proteins of Bacillus subtilis. Mol Microbiol 2007, 65(2):454–464.
Yoshimura M, Asai K, Sadaie Y, Yoshikawa H: Interaction of Bacillus subtilis extracytoplasmic function (ECF) sigma factors with the N-terminal regions of their potential anti-sigma factors. Microbiology 2004, 150(Pt 3):591–599.
Helmann JD: The extracytoplasmic function (ECF) sigma factors. Adv Microb Physiol 2002, 46: 47–110.
Coque JJ, Liras P, Martin JF: Genes for a beta-lactamase, a penicillin-binding protein and a transmembrane protein are clustered with the cephamycin biosynthetic genes in Nocardia lactamdurans. Embo J 1993, 12(2):631–639.
Fleischmann RD, Alland D, Eisen JA, Carpenter L, White O, Peterson J, DeBoy R, Dodson R, Gwinn M, Haft D, Hickey E, Kolonay JF, Nelson WC, Umayam LA, Ermolaeva M, Salzberg SL, Delcher A, Utterback T, Weidman J, Khouri H, Gill J, Mikula A, Bishai W, Jacobs Jr WR Jr., Venter JC, Fraser CM: Whole-genome comparison of Mycobacterium tuberculosis clinical and laboratory strains. J Bacteriol 2002, 184(19):5479–5490.
Chirakkal H, O'Rourke M, Atrih A, Foster SJ, Moir A: Analysis of spore cortex lytic enzymes and related proteins in Bacillus subtilis endospore germination. Microbiology 2002, 148(Pt 8):2383–2392.
Rotanova TV, Botos I, Melnikov EE, Rasulova F, Gustchina A, Maurizi MR, Wlodawer A: Slicing a protease: structural features of the ATP-dependent Lon proteases gleaned from investigations of isolated domains. Protein Sci 2006, 15(8):1815–1828.
Moller S, Croning MD, Apweiler R: Evaluation of methods for the prediction of membrane spanning regions. Bioinformatics 2001, 17(7):646–653.
Shen HB, Yang J, Chou KC: Fuzzy KNN for predicting membrane protein types from pseudo-amino acid composition. J Theor Biol 2006, 240(1):9–13.
Bernsel A, Von Heijne G: Improved membrane protein topology prediction by domain assignments. Protein Sci 2005, 14(7):1723–1728.
McDonough JA, Hacker KE, Flores AR, Pavelka MS Jr., Braunstein M: The twin-arginine translocation pathway of Mycobacterium smegmatis is functional and required for the export of mycobacterial beta-lactamases. J Bacteriol 2005, 187(22):7667–7679.
Widdick DA, Dilks K, Chandra G, Bottrill A, Naldrett M, Pohlschroder M, Palmer T: The twin-arginine translocation pathway is a major route of protein export in Streptomyces coelicolor. Proc Natl Acad Sci U S A 2006, 103(47):17927–17932.
Posey JE, Shinnick TM, Quinn FD: Characterization of the twin-arginine translocase secretion system of Mycobacterium smegmatis. J Bacteriol 2006, 188(4):1332–1340.
De Keersmaeker S, Vrancken K, Van Mellaert L, Lammertyn E, Anne J, Geukens N: Evaluation of TatABC overproduction on Tat- and Sec-dependent protein secretion in Streptomyces lividans. Arch Microbiol 2006, 186(6):507–512.
De Keersmaeker S, Van Mellaert L, Lammertyn E, Vrancken K, Anne J, Geukens N: Functional analysis of TatA and TatB in Streptomyces lividans. Biochem Biophys Res Commun 2005, 335(3):973–982.
Caldelari I, Mann S, Crooks C, Palmer T: The Tat pathway of the plant pathogen Pseudomonas syringae is required for optimal virulence. Mol Plant Microbe Interact 2006, 19(2):200–212.
Schaerlaekens K, Van Mellaert L, Lammertyn E, Geukens N, Anne J: The importance of the Tat-dependent protein secretion pathway in Streptomyces as revealed by phenotypic changes in tat deletion mutants and genome analysis. Microbiology 2004, 150(Pt 1):21–31.
Bendtsen JD, Nielsen H, Widdick D, Palmer T, Brunak S: Prediction of twin-arginine signal peptides. BMC Bioinformatics 2005, 6: 167.
Berks BC: A common export pathway for proteins binding complex redox cofactors? Mol Microbiol 1996, 22(3):393–404.
Stanley NR, Palmer T, Berks BC: The twin arginine consensus motif of Tat signal peptides is involved in Sec-independent protein targeting in Escherichia coli. J Biol Chem 2000, 275(16):11591–11596.
Lee PA, Tullman-Ercek D, Georgiou G: The bacterial twin-arginine translocation pathway. Annu Rev Microbiol 2006, 60: 373–395.
Jongbloed JD, Antelmann H, Hecker M, Nijland R, Bron S, Airaksinen U, Pries F, Quax WJ, van Dijl JM, Braun PG: Selective contribution of the twin-arginine translocation pathway to protein secretion in Bacillus subtilis. J Biol Chem 2002, 277(46):44068–44078.
Dilks K, Rose RW, Hartmann E, Pohlschroder M: Prokaryotic utilization of the twin-arginine translocation pathway: a genomic survey. J Bacteriol 2003, 185(4):1478–1483.
Guo T, Hua S, Ji X, Sun Z: DBSubLoc: database of protein subcellular localization. Nucl Acids Res 2004, 32(suppl_1):D122–124.
Nair R, Rost B: Mimicking cellular sorting improves prediction of subcellular localization. J Mol Biol 2005, 348(1):85–100.
Szafron D, Lu P, Greiner R, Wishart DS, Poulin B, Eisner R, Lu Z, Anvik J, Macdonell C, Fyshe A, Meeuwis D: Proteome Analyst: custom predictions with explanations in a web-based tool for high-throughput proteome annotations. Nucleic Acids Res 2004, 32(Web Server issue):W365–71.
Rashid M, Saha S, Raghava GP: Support Vector Machine-based method for predicting subcellular localization of mycobacterial proteins using evolutionary information and motifs. BMC Bioinformatics 2007, 8(1):337.
Chou KC, Shen HB: Hum-PLoc: a novel ensemble classifier for predicting human protein subcellular localization. Biochem Biophys Res Commun 2006, 347(1):150–157.
Shen HB, Chou KC: Hum-mPLoc: an ensemble classifier for large-scale human protein subcellular location prediction by incorporating samples with multiple sites. Biochem Biophys Res Commun 2007, 355(4):1006–1011.
Hiller K, Grote A, Scheer M, Munch R, Jahn D: PrediSi: prediction of signal peptides and their cleavage positions. Nucleic Acids Res 2004, 32(Web Server issue):W375–9.
Zhang Z, Henzel WJ: Signal peptide prediction based on analysis of experimentally verified cleavage sites. Protein Sci 2004, 13(10):2819–2824.
Menne KM, Hermjakob H, Apweiler R: A comparison of signal sequence prediction methods using a test set of signal peptides. Bioinformatics 2000, 16(8):741–742.
Carlos JL, Paetzel M, Brubaker G, Karla A, Ashwell CM, Lively MO, Cao G, Bullinger P, Dalbey RE: The role of the membrane-spanning domain of type I signal peptidases in substrate cleavage site selection. J Biol Chem 2000, 275(49):38813–38822.
Taylor PD, Attwood TK, Flower DR: Combining algorithms to predict bacterial protein sub-cellular location: Parallel versus concurrent implementations. Bioinformation 2006, 1(8):285–289.
Joris B, Englebert S, Chu CP, Kariyama R, Daneo-Moore L, Shockman GD, Ghuysen JM: Modular design of the Enterococcus hirae muramidase-2 and Streptococcus faecalis autolysin. FEMS Microbiol Lett 1992, 70(3):257–264.
Waligora AJ, Hennequin C, Mullany P, Bourlioux P, Collignon A, Karjalainen T: Characterization of a cell surface protein of Clostridium difficile with adhesive properties. Infect Immun 2001, 69(4):2144–2153.
Furst P, Mosch HU, Solioz M: A protein of unusual composition from Enterococcus faecium. Nucleic Acids Res 1989, 17(16):6724.
Aubert-Pivert E, Davies J: Biosynthesis of butirosin in Bacillus circulans NRRL B3312: identification by sequence analysis and insertional mutagenesis of the butB gene involved in antibiotic production. Gene 1994, 147(1):1–11.
Margot P, Pagni M, Karamata D: Bacillus subtilis 168 gene lytF encodes a gamma-D-glutamate-meso-diaminopimelate muropeptidase expressed by the alternative vegetative sigma factor, sigmaD. Microbiology 1999, 145 ( Pt 1): 57–65.
Noback MA, Holsappel S, Kiewiet R, Terpstra P, Wambutt R, Wedler H, Venema G, Bron S: The 172 kb prkA-addAB region from 83 degrees to 97 degrees of the Bacillus subtilis chromosome contains several dysfunctional genes, the glyB marker, many genes encoding transporter proteins, and the ubiquitous hit gene. Microbiology 1998, 144 ( Pt 4): 859–875.
Steiner B, Romero-Steiner S, Cruce D, George R: Cloning and sequencing of the hyaluronate lyase gene from Propionibacterium acnes. Can J Microbiol 1997, 43(4):315–321.
Mesnage S, Fontaine T, Mignot T, Delepierre M, Mock M, Fouet A: Bacterial SLH domain proteins are non-covalently anchored to the cell surface via a conserved mechanism involving wall polysaccharide pyruvylation. Embo J 2000, 19(17):4473–4484.
Fernandez-Tornero C, Lopez R, Garcia E, Gimenez-Gallego G, Romero A: A novel solenoid fold in the cell wall anchoring domain of the pneumococcal virulence factor LytA. Nat Struct Biol 2001, 8(12):1020–1024.
Jonquieres R, Bierne H, Fiedler F, Gounon P, Cossart P: Interaction between the protein InlB of Listeria monocytogenes and lipoteichoic acid: a novel mechanism of protein association at the surface of gram-positive bacteria. Mol Microbiol 1999, 34(5):902–914.
Cabanes D, Dehoux P, Dussurget O, Frangeul L, Cossart P: Surface proteins and the pathogenic potential of Listeria monocytogenes. Trends Microbiol 2002, 10(5):238–245.
Sara M: Conserved anchoring mechanisms between crystalline cell surface S-layer proteins and secondary cell wall polymers in Gram-positive bacteria? Trends Microbiol 2001, 9(2):47–9; discussion 49–50.
Chauvaux S, Matuschek M, Beguin P: Distinct affinity of binding sites for S-layer homologous domains in Clostridium thermocellum and Bacillus anthracis cell envelopes. J Bacteriol 1999, 181(8):2455–2458.
Foster SJ: Cloning, expression, sequence analysis and biochemical characterization of an autolytic amidase of Bacillus subtilis 168 trpC2. J Gen Microbiol 1991, 137(8):1987–1998.
Nolling J, Breton G, Omelchenko MV, Makarova KS, Zeng Q, Gibson R, Lee HM, Dubois J, Qiu D, Hitti J, Wolf YI, Tatusov RL, Sabathe F, Doucette-Stamm L, Soucaille P, Daly MJ, Bennett GN, Koonin EV, Smith DR: Genome sequence and comparative analysis of the solvent-producing bacterium Clostridium acetobutylicum. J Bacteriol 2001, 183(16):4823–4838.
Desvaux M, Khan A, Beatson SA, Scott-Tucker A, Henderson IR: Protein secretion systems in Fusobacterium nucleatum: genomic identification of Type 4 piliation and complete Type V pathways brings new insight into mechanisms of pathogenesis. Biochim Biophys Acta 2005, 1713(2):92–112.
Desvaux M, Hebraud M, Henderson IR, Pallen MJ: Type III secretion: what's in a name? Trends Microbiol 2006, 14(4):157–160.
Desvaux M, Hebraud M: The protein secretion systems in Listeria: inside out bacterial virulence. FEMS Microbiol Rev 2006, 30(5):774–805.
Desvaux M, Dumas E, Chafsey I, Hebraud M: Protein cell surface display in Gram-positive bacteria: from single protein to macromolecular protein structure. FEMS Microbiol Lett 2006, 256(1):1–15.
Desvaux M, Scott-Tucker A, Turner SM, Cooper LM, Huber D, Nataro JP, Henderson IR: A conserved extended signal peptide region directs posttranslational protein translocation via a novel mechanism. Microbiology 2007, 153(Pt 1):59–70.
Ghosh J, Caparon MG: Specificity of Streptococcus pyogenes NAD(+) glycohydrolase in cytolysin-mediated translocation. Mol Microbiol 2006, 62(4):1203–1214.
Henderson IR, Cappello R, Nataro JP: Autotransporter proteins, evolution and redefining protein secretion. Trends Microbiol 2000, 8(12):529–532.
Henderson IR, Navarro-Garcia F, Desvaux M, Fernandez RC, Ala'Aldeen D: Type V protein secretion pathway: the autotransporter story. Microbiol Mol Biol Rev 2004, 68(4):692–744.
Rey S, Gardy JL, Brinkman FS: Assessing the precision of high-throughput computational and laboratory approaches for the genome-wide identification of protein subcellular localization in bacteria. BMC Genomics 2005, 6: 162.
Yu CS, Chen YC, Lu CH, Hwang JK: Prediction of protein subcellular localization. Proteins 2006, 64(3):643–651.
Bendtsen JD, Kiemer L, Fausboll A, Brunak S: Non-classical protein secretion in bacteria. BMC Microbiol 2005, 5: 58.
Chou KC, Shen HB: Euk-mPLoc: a fusion classifier for large-scale eukaryotic protein subcellular location prediction by incorporating multiple sites. J Proteome Res 2007, 6(5):1728–1734.
Zuber B, Haenni M, Ribeiro T, Minnig K, Lopes F, Moreillon P, Dubochet J: Granular layer in the periplasmic space of gram-positive bacteria and fine structures of Enterococcus gallinarum and Streptococcus gordonii septa revealed by cryo-electron microscopy of vitreous sections. J Bacteriol 2006, 188(18):6652–6660.
Matias VR, Beveridge TJ: Cryo-electron microscopy reveals native polymeric cell wall structure in Bacillus subtilis 168 and the existence of a periplasmic space. Mol Microbiol 2005, 56(1):240–251.
Matias VR, Beveridge TJ: Native cell wall organization shown by cryo-electron microscopy confirms the existence of a periplasmic space in Staphylococcus aureus. J Bacteriol 2006, 188(3):1011–1021.
Pooley HM, Merchante R, Karamata D: Overall protein content and induced enzyme components of the periplasm of Bacillus subtilis. Microb Drug Resist 1996, 2(1):9–15.
Carlsson F, Stalhammar-Carlemalm M, Flardh K, Sandin C, Carlemalm E, Lindahl G: Signal sequence directs localized secretion of bacterial surface proteins. Nature 2006, 442(7105):943–946.
Rothe B, Rothe B, Roggentin P, Schauer R: The sialidase gene from Clostridium septicum: cloning, sequencing, expression in Escherichia coli and identification of conserved sequences in sialidases and other proteins. Mol Gen Genet 1991, 226(1–2):190–197.
Bahl H, Burchhardt G, Spreinat A, Haeckel K, Wienecke A, Schmidt B, Antranikian G: alpha-Amylase of Clostridium thermosulfurogenes EM1: nucleotide sequence of the gene, processing of the enzyme, and comparison of other alpha-amylases. Appl Environ Microbiol 1991, 57(5):1554–1559.
Shen SH, Chretien P, Bastien L, Slilaty SN: Primary sequence of the glucanase gene from Oerskovia xanthineolytica. Expression and purification of the enzyme from Escherichia coli. J Biol Chem 1991, 266(2):1058–1063.
Rothe B, Roggentin P, Frank R, Blocker H, Schauer R: Cloning, sequencing and expression of a sialidase gene from Clostridium sordellii G12. J Gen Microbiol 1989, 135(11):3087–3096.
Adams MW: The structure and mechanism of iron-hydrogenases. Biochim Biophys Acta 1990, 1020(2):115–145.
Jean-Luc Falcone RK Dominique Belin and Bastien Chopard: Understanding the Signal Sequences with Machine Learning. In EvoBIO 2007 proceedings. Volume LNCS 4447. Springer; 2007:pp. 57–67.
Legler PM, Cai M, Peterkofsky A, Clore GM: Three-dimensional solution structure of the cytoplasmic B domain of the mannitol transporter IImannitol of the Escherichia coli phosphotransferase system. J Biol Chem 2004, 279(37):39115–39121.
Pollastri G, Martin AJ, Mooney C, Vullo A: Accurate prediction of protein secondary structure and solvent accessibility by consensus combiners of sequence and structure information. BMC Bioinformatics 2007, 8(1):201.
Gatlin CL, Pieper R, Huang ST, Mongodin E, Gebregeorgis E, Parmar PP, Clark DJ, Alami H, Papazisi L, Fleischmann RD, Gill SR, Peterson SN: Proteomic profiling of cell envelope-associated proteins from Staphylococcus aureus. Proteomics 2006, 6(5):1530–1549.
Sanyal SC, Pal S, Chowdhury S, DasGupta C: 23S rRNA assisted folding of cytoplasmic malate dehydrogenase is distinctly different from its self-folding. Nucleic Acids Res 2002, 30(11):2390–2397.
Overbeek R, Larsen N, Walunas T, D'Souza M, Pusch G, Selkov E Jr., Liolios K, Joukov V, Kaznadzey D, Anderson I, Bhattacharyya A, Burd H, Gardner W, Hanke P, Kapatral V, Mikhailova N, Vasieva O, Osterman A, Vonstein V, Fonstein M, Ivanova N, Kyrpides N: The ERGO genome analysis and discovery system. Nucleic Acids Res 2003, 31(1):164–171.
Edgar RC: MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res 2004, 32(5):1792–1797.
Eddy SR: Profile hidden Markov models. Bioinformatics 1998, 14(9):755–763.
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.
de Jong A, van Hijum SA, Bijlsma JJ, Kok J, Kuipers OP: BAGEL: a web-based bacteriocin genome mining tool. Nucleic Acids Res 2006, 34(Web Server issue):W273–9.
Chou KC, Shen HB: MemType-2L: a web server for predicting membrane proteins and their types by incorporating evolution information through Pse-PSSM. Biochem Biophys Res Commun 2007, 360(2):339–345.
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.
Shen HB, Chou KC: Signal-3L: A 3-layer approach for predicting signal peptides. Biochem Biophys Res Commun 2007, 363(2):297–303.
Chou KC, Shen HB: Signal-CF: a subsite-coupled and window-fusing approach for predicting signal peptides. Biochem Biophys Res Commun 2007, 357(3):633–640.
Taylor PD, Toseland CP, Attwood TK, Flower DR: TATPred: a Bayesian method for the identification of twin arginine translocation pathway signal sequences. Bioinformation 2006, 1(5):184–187.
Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ: Basic local alignment search tool. J Mol Biol 1990, 215(3):403–410.
O'Brien KP, Remm M, Sonnhammer EL: Inparanoid: a comprehensive database of eukaryotic orthologs. Nucleic Acids Res 2005, 33(Database issue):D476–80.
Sutcliffe IC, Harrington DJ: Pattern searches for the identification of putative lipoprotein genes in Gram-positive bacterial genomes. Microbiology 2002, 148(Pt 7):2065–2077.
Bengtsson J, Tjalsma H, Rivolta C, Hederstedt L: Subunit II of Bacillus subtilis cytochrome c oxidase is a lipoprotein. J Bacteriol 1999, 181(2):685–688.
DeLisa MP, Tullman D, Georgiou G: Folding quality control in the export of proteins by the bacterial twin-arginine translocation pathway. Proc Natl Acad Sci U S A 2003, 100(10):6115–6120.
Palmer T, Sargent F, Berks BC: Export of complex cofactor-containing proteins by the bacterial Tat pathway. Trends Microbiol 2005, 13(4):175–180.
Robinson C, Bolhuis A: Tat-dependent protein targeting in prokaryotes and chloroplasts. Biochim Biophys Acta 2004, 1694(1–3):135–147.
Berks BC, Palmer T, Sargent F: Protein targeting by the bacterial twin-arginine translocation (Tat) pathway. Curr Opin Microbiol 2005, 8(2):174–181.
Tjalsma H, Antelmann H, Jongbloed JD, Braun PG, Darmon E, Dorenbos R, Dubois JY, Westers H, Zanen G, Quax WJ, Kuipers OP, Bron S, Hecker M, van Dijl JM: Proteomics of protein secretion by Bacillus subtilis: separating the "secrets" of the secretome. Microbiol Mol Biol Rev 2004, 68(2):207–233.
Jongbloed JD, Grieger U, Antelmann H, Hecker M, Nijland R, Bron S, van Dijl JM: Two minimal Tat translocases in Bacillus. Mol Microbiol 2004, 54(5):1319–1325.
Boekhorst J, de Been MW, Kleerebezem M, Siezen RJ: Genome-wide detection and analysis of cell wall-bound proteins with LPxTG-like sorting motifs. J Bacteriol 2005, 187(14):4928–4934.
Marraffini LA, Dedent AC, Schneewind O: Sortases and the art of anchoring proteins to the envelopes of gram-positive bacteria. Microbiol Mol Biol Rev 2006, 70(1):192–221.
Marraffini LA, Schneewind O: Targeting proteins to the cell wall of sporulating Bacillus anthracis. Mol Microbiol 2006.
Gaspar AH, Marraffini LA, Glass EM, Debord KL, Ton-That H, Schneewind O: Bacillus anthracis sortase A (SrtA) anchors LPXTG motif-containing surface proteins to the cell wall envelope. J Bacteriol 2005, 187(13):4646–4655.
Mazmanian SK, Ton-That H, Schneewind O: Sortase-catalysed anchoring of surface proteins to the cell wall of Staphylococcus aureus. Mol Microbiol 2001, 40(5):1049–1057.
Navarre WW, Schneewind O: Surface proteins of gram-positive bacteria and mechanisms of their targeting to the cell wall envelope. Microbiol Mol Biol Rev 1999, 63(1):174–229.
Calvo E, Pucciarelli MG, Bierne H, Cossart P, Albar JP, Garcia-Del Portillo F: Analysis of the Listeria cell wall proteome by two-dimensional nanoliquid chromatography coupled to mass spectrometry. Proteomics 2005, 5(2):433–443.
van Dijl JM, de Jong A, Vehmaanpera J, Venema G, Bron S: Signal peptidase I of Bacillus subtilis: patterns of conserved amino acids in prokaryotic and eukaryotic type I signal peptidases. Embo J 1992, 11(8):2819–2828.
Berks BC, Sargent F, Palmer T: The Tat protein export pathway. Mol Microbiol 2000, 35(2):260–274.
Rose RW, Bruser T, Kissinger JC, Pohlschroder M: Adaptation of protein secretion to extremely high-salt conditions by extensive use of the twin-arginine translocation pathway. Mol Microbiol 2002, 45(4):943–950.
Rey S, Acab M, Gardy JL, Laird MR, deFays K, Lambert C, Brinkman FSL: PSORTdb: a protein subcellular localization database for bacteria. Nucl Acids Res 2005, 33(suppl_1):D164–168.
We thank the authors of the tools used in the LocateP pipeline, especially Anne de Jong, Anders Krogh, Karsten Hiller, and Andre Billion for Bagel, LipoP 1.0, PrediSi and Augur, repectively. We thank Antoinette Killian for stimulating discussion and Qinghu Ren for his help with TransportDB. Special thanks to the BioRange programme members Jack Leunissen for PatScan software support, Tim Hulsen for the help in database construction, Elena Marchiori and Jaap Heringa for suggestions on machine-learning methods. This work is part of the BioRange programme of the Netherlands Bioinformatics Centre (NBIC), which is supported by a BSIK grant through the Netherlands Genomics Initiative (NGI). Christof Francke is supported by a grant from the Kluyver Centre for Genomics of Industrial Fermentation.
RS conceived the project and acquired the funding, MM performed the research and constructed the LocateP pipeline and database, JB helped in tool selection, CF and RS supervised the project, while MM, CF and RS wrote the manuscript. All authors read and approved the manuscript.
Electronic supplementary material
Additional file 1: . The different SCL tools used at each decision step are indicated. The different SCLs distinguished are boxed in the middle; "Bacteriocin" signifies bacteriocin-like proteins secreted by non-classical pathways, identified by Bagel. a) all tools agreed, b) all possible hits, c) majority vote, d) 2–3 TM segments and C-terminus detected by LPxTG HMM. (PDF 77 KB)
About this article
Cite this article
Zhou, M., Boekhorst, J., Francke, C. et al. LocateP: Genome-scale subcellular-location predictor for bacterial proteins. BMC Bioinformatics 9, 173 (2008). https://doi.org/10.1186/1471-2105-9-173
- Hide Markov Model
- Transmembrane Helix
- Putative Cleavage Site
- Hide Markov Model Pair
- Hide Markov Model Score