Subcellular Localization of Extracytoplasmic Proteins in Monoderm Bacteria: Rational Secretomics-Based Strategy for Genomic and Proteomic Analyses (original) (raw)

Open Access

Peer-reviewed

Research Article

Subcellular Localization of Extracytoplasmic Proteins in Monoderm Bacteria: Rational Secretomics-Based Strategy for Genomic and Proteomic Analyses

PLOS

x

Figures

Abstract

Genome-scale prediction of subcellular localization (SCL) is not only useful for inferring protein function but also for supporting proteomic data. In line with the secretome concept, a rational and original analytical strategy mimicking the secretion steps that determine ultimate SCL was developed for Gram-positive (monoderm) bacteria. Based on the biology of protein secretion, a flowchart and decision trees were designed considering (i) membrane targeting, (ii) protein secretion systems, (iii) membrane retention, and (iv) cell-wall retention by domains or post-translocational modifications, as well as (v) incorporation to cell-surface supramolecular structures. Using Listeria monocytogenes as a case study, results were compared with known data set from SCL predictors and experimental proteomics. While in good agreement with experimental extracytoplasmic fractions, the secretomics-based method outperforms other genomic analyses, which were simply not intended to be as inclusive. Compared to all other localization predictors, this method does not only supply a static snapshot of protein SCL but also offers the full picture of the secretion process dynamics: (i) the protein routing is detailed, (ii) the number of distinct SCL and protein categories is comprehensive, (iii) the description of protein type and topology is provided, (iv) the SCL is unambiguously differentiated from the protein category, and (v) the multiple SCL and protein category are fully considered. In that sense, the secretomics-based method is much more than a SCL predictor. Besides a major step forward in genomics and proteomics of protein secretion, the secretomics-based method appears as a strategy of choice to generate in silico hypotheses for experimental testing.

Citation: Renier S, Micheau P, Talon R, Hébraud M, Desvaux M (2012) Subcellular Localization of Extracytoplasmic Proteins in Monoderm Bacteria: Rational Secretomics-Based Strategy for Genomic and Proteomic Analyses. PLoS ONE 7(8): e42982. https://doi.org/10.1371/journal.pone.0042982

Editor: Ivo G. Boneca, Institut Pasteur Paris, France

Received: February 10, 2012; Accepted: July 13, 2012; Published: August 9, 2012

Copyright: © Renier et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

Funding: This work was supported in part by the “Institut National de la Recherche Agronomique (INRA)” and the European Framework Programme 6 (FP6) with the ProSafeBeef (Advancing Beef Safety and Quality through Research and Innovation) research consortium. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Competing interests: The authors have declared that no competing interests exist.

Introduction

All living cells interface with their surrounding through proteins that are located in the cell envelope, displayed on the cell surface or released into the extracellular milieu, and even beyond when injected into a host cell. Such proteins are translocated in the first instance through biological membranes by protein-conducting channels, i.e. translocons. In bacteria, several secretion systems enabling protein transport from the inside to the outside of the cell have been characterised. Distinction must be made between bacteria possessing (i) one biological membrane (the cytoplasmic membrane), the so-called monodermata, and (ii) two biological membranes (the inner membrane and outer membrane), the so-called didermata [1], [2], [3]. Until recently and in contrast to diderm bacteria [4], [5], [6], [7], [8], comprehensive knowledge on protein secretion systems in monoderm species was restricted to the non-pathogenic Bacillus subtilis species [9], [10], [11] or scattered among different species for specific systems [12], [13], [14], [15].

Because of the absence of an outer membrane, the numerical classification of protein secretion systems does not apply to monodermata and export across the cytoplasmic membrane actually corresponds to a secretion event [16]. As in didermata, the Sec (secretion) and Tat (twin-arginine translocation) machineries are found in the cytoplasmic membrane but additional secretion systems can be present in monoderms, i.e. the FPE (fimbrilin-protein exporter), ABC (ATP-binding cassette) transporters, FEA (flagellum export apparatus), holins (hole-forming) and Wss (WXG100 secretion system) [17], [18]. As thoroughly explain by several specialists in the field of bacterial protein secretion [2], [19], [20], [21], [22], [23], [24], we will abstain to use the “T7SS” terminology to describe the Wss in monoderms, which is actually ascribed to the chaperone-usher pathway in diderm-LPS and at best only apply to diderm-mycolate, which is restricted to bacteria of the genus Mycobacterium.

As in any living cell, extracytoplasmic proteins cover a vast variety of functions, including nutrient uptake, chemosensing, motility, adhesion or cell envelope biogenesis. Moreover, their subcellular localization (SCL) and biological functions provide distinguishing clues regarding the physiology, lifestyle, position and interactions of the bacterial cell in an ecological niche and more generally in an ecosystem, as revealed for example by extracellular degradative enzymes in saprophytes or cell-surface virulence factors in pathogens. The final SCL of a protein results from a series of molecular mechanisms, involving post-translational and/or post-translocational modifications. The secretome concept is very useful for considering these different steps as it includes both the proteins secreted across the cytoplasmic membrane to the membrane–cell wall interface, the cell wall or the extracellular environment, and importantly, the secretory machineries themselves [10], [25], [26]. As such and contrary to what it sometimes misconceptualized by some authors, the secretome is not a proteome per se, let alone the subset of extracellular proteins, that is actually the exoproteome [1]. In bacteria with a Gram-positive cell-envelope architecture, the proteins actively transported via these secretion systems, the so-called secreted proteins, can have radically different final destinations and be either (i) anchored to the cytoplasmic membrane, (ii) associated with the cell wall, (iii) released into the extracellular milieu, or even (iv) injected into a host cell [1]. Description of SCL now follows the Gene Ontology (GO) recommendations for describing “Cellular component”, one of the three structured controlled vocabularies [27].

Because experimental investigation of the membrane and cell wall proteomes is hindered by technical limitation of protein extraction from their subcellular fractions, genomic prediction of SCL has been the subject of intense research effort. Numerous localization predictors have been developed for predicting the final destination of proteins. These bioinformatic tools can be divided into (i) specialized prediction tools, essentially based on the identification of signal peptides or retention sequences to the membrane or cell wall, e.g. SignalP [28], LipoP [29], TMHMM [30] or CW-PRED [31], and (ii) global prediction tools indicating the protein final SCL, e.g. PSORTb [32], LocTree [33], CELLO [34] or Gpos-mPLoc [35]. Such ensemble classifiers based on support vector machine (SVM) or neural network (NN) have been constructed on algorithms with a rationale somehow disconnected from the biology of the system investigated. Each of these tools having its own prediction limits, though, an alternative and powerful strategy consists in combining predictions [36]. For Gram-positive bacteria, different pipelines have been developed to predict final location of protein, e.g. Augur [37], LocateP [38] or SurfG+ [39], but none of them is comprehensive. A momentous limitation is that, by essence, their workflows are not evolutive but established once and for all and cannot be willingly adjusted in light of new findings in the field. Consequently, new specialized prediction tools cannot be swiftly implemented since any modifications remain at the discretion of their designers. In addition, results of some of these tools are frozen in databases and interrogations cannot be readily performed on demand or on newly available bacterial genomes.

In order to provide a much more comprehensive, reliable, flexible and adjustable prediction of protein SCL in monoderm bacteria, we aimed at developing a strategy for analysing genomic and proteomic data for secreted proteins. The originality of the method resides in the use of an extensive number of readily available bioinformatic tools organized in a workflow and decision trees mimicking very closely the molecular steps encountered by a protein in the course of secretion. Compared to all available predictions tools to date, this methodology embraces the secretome concept as it does not only consider the presence of various domains that can target or retain the secreted proteins within the cell-envelope of monoderm bacterium but most importantly considers the presence/absence of protein secretion pathways. Taking into account (i) it potentially expresses a prodigious amount of extracellular proteins and cell-surface proteins, which are membrane or cell-wall attached [4], [40], [41], (ii) it is a pathogenic Gram-positive bacterium contrary to the paradigm B. subtilis, and (iii) numerous experimental proteomic approaches have been dedicated to its extracytoproteome, including the surface and extracellular proteomes [42], [43], [44], [45], [46], [47], it prompts us to focus on the extracytoproteome of Listeria monocytogenes as a case study.

Results

Design of the secretomics-based method for genomic analysis of secreted proteins in monoderm bacteria

The rationale for analyzing the final SCL of secreted proteins is based on the biology of protein secretion, i.e. the secretome. This concept provides an integrated and global view by considering protein routing, transport systems, post-translocational (translational) modifications and subcellular location [10], [25], [26]. Based on the secretome in monoderm bacteria (Figure 1), the steps considered are (i) the targeting of protein to the membrane by an N-terminal signal peptide, (ii) the protein secretion systems present, i.e. Sec, Tat, ABC, FPE, FEA, holin and/or Wss, (iii) the membrane retention of secreted proteins by domains or post-translational modifications, (iv) the cell-wall retention of secreted proteins by domains or post-translational modifications, and (v) the incorporation to cell-surface supramolecular structures. A summary of the different abbreviations in use in relation to protein secretion, especially protein categories and subcellular localizations, is provided in Table 1.

thumbnail

Figure 1. Synoptic view of the secretomics-based method for monoderm bacterium.

Based on the biology of protein secretion, the coding sequences (CDS) are sequentially analysed in a workflow for (i) signal peptide (SP), (ii) secretory pathway, (iii) membrane retention, (iv) cell-wall retention, and (v) surface appendage. For each step, a combination of different tools allows defining different databases, as indicated in the detailed flowchart (Figure 2). From there, the resulting databases are analysed as depicted in the detailed decision trees (Figure 3). In the end, proteins are discriminated into different categories and different SCL are predicted. Sec-SP: Sec-dependent SP; Unc-SP: uncleaved SP; TMD: transmembrane domain; PGBD: peptidoglycan-binding domain; CWBD: cell-wall binding domain; SLHD: S-layer homology domain; SL-Prot: S-layer protein; Doc/Coh: dockerin/cohesin domain; IMP: integral membrane protein; GO: gene ontology.

https://doi.org/10.1371/journal.pone.0042982.g001

As one of the cornerstones of the secretomics-based method, the type of signal peptide (SP) is defined (Figure 2). Besides, the search in parallel of all known protein secretion pathways in monoderm bacteria constitutes the keystone of this analytical strategy. Secretion pathways not only focus on translocases of the secretion systems but include the associated post-translocation maturation pathways, namely the respective signal peptidases, the lipoprotein maturation pathway, and the covalent cell-wall anchoring sortase pathway (Table 2). It is important to stress that the absence of a SP cannot rule out protein secretion by alternative systems as detailed below. Presence of the different types of SP is associated with the respective secretion systems, whereas substrates to alternative secretion systems lacking SP are subjected to similarity search (Figure 3). Together with pseudopilus and flagellum secreted via FPE and FEA respectively, other surface supramolecular structures are considered, namely S-layer, pilus and/or cellulosome. For integral membrane protein (IMP) and besides the presence of transmembrane domain (TMD), the presence of uncleaved N-terminal SP, which could serve as a signal anchor, is utterly considered. Secreted proteins, which are neither predicted as retained to the membrane, cell-wall nor subunits of surface appendages, are predicted as exoproteins. Once the different protein categories are defined, SCL can be predicted (Figure 3). Cytoproteins potentially secreted by non-classical pathway (NC) are also predicted as localized in the extracellular milieu. Position +2 of the cleavage site (C+2) is checked as an indicator that the lipoprotein can be potentially released into the extracellular milieu [26], [48]. The number of GW modules is carefully considered for protein SCL in the extracellular milieu rather than in the cell wall. Besides discriminating the protein category from the SCL to avoid any confusing statement or misleading interpretation, the secretomics-based method provides a detailed description of protein routing and permits considering multiple SCL for a given protein.

thumbnail

Figure 2. Comprehensive flowchart of the secretomics-based method in a monoderm bacterium.

The analysis considered the (i) signal peptide (SP), (ii) type of SP (Signal anchor, FPE-SP, Tat-SP, ABC-S and, lipobox), (iv) protein secretion systems, (v) exported proteins lacking a SP (Export, FEA-, Holin, and Wss-substrates), (iv) transmembrane domain (TMD), and (iv) relevant conserved domains (LPXTG, WXL, LysM, CWBD1,...etc). Details of the prediction tools used for the analysis and definition of the databases are provided inthe Materials & Methods section.

https://doi.org/10.1371/journal.pone.0042982.g002

thumbnail

Figure 3. Detailed decision trees for prediction of protein category and SCL of secreted proteins.

(1) Proteins exhibiting N-terminal SP are extracted from TMD database (Y). (2) The different types of SP (Tat-SP, ABC-SP and FPE-SP) are extracted (Y). (3) Absence of a signal anchor (N) and export (Y) define Sec-dependent SP (Sec-SP). (1) Proteins with TMD but no predicted SP (N), (4) are checked for uncleaved SP (Unc-SP), i.e. TMD of at least 7 amino acid within the first 100 N-terminal residues and with Nin-Cout topology (Type II signal) (Y). (3) Unc-SP also comprises proteins with signal anchor (Y) and SP categorised as non-exported (N). (5) From the types of SP are clearly defined. Together with (4) proteins with TMD but no SP (N) and (6) protein substrates of holins, Wss and FEA, (7) the presence of the respective protein secretion systems is checked (Y). (8) When the respective protein secretion system is absent (N) or (9) proteins are not predicted as secreted (N), proteins are considered as cytoproteins and located in the CP. (10) Cytoproteins predicted as exported by NC (Y) are further considered as located extracellularly. (7) Secreted proteins and their respective secretion system are defined from there. (11) Translocated proteins with Unc-SP (Y) are IMPs. (11) Translocated proteins without Unc-SP (N), and (12) with TMD (Y) but no SP (N) are IMPs. (13) Remaining translocated proteins with a cleavable SP (Y) and a single predicted TMD (TMD = 1) (Y) cannot be IMP, and are checked for (16) the presence of a lipobox. (13) Remaining translocated proteins with more than one TMD (N) are checked for (14) the absence of overlap (N) with SP region and LPXTG domain respectively (TMD = 2 AND TMD = LPXTG) to be IMP, otherwise (Y) are checked for (16) the presence of a lipobox. (15) From TMD topology prediction (Figure 2), IMPs are further subcategorised and considered as integral to CM. (16) Presence of a lipobox (Y) define lipoproteins anchored to the CM. (17) The presence of glycine residue at position C+2 (C+2 = G) (Y) indicates potential release into the EM [26], [48]. (18) Presence of cell-wall retention domains define (19) parietal proteins (CW-protein) that are further subcategorised (Figure 2) and considered as located at the CW. (20) Proteins with less tha 2 GW modules is not defined as CW-protein located at the CW [57], [58]. (21) Proteins part of S-layer, pilus and cellulosome, as well as (22) pseudopilus and flagellum are defined. (23) Secreted proteins with none of the cell-envelope retention are as exoproteins located in the EM. N: No, Y: Yes.

https://doi.org/10.1371/journal.pone.0042982.g003

SCL prediction of secreted proteins in L. monocytogenes following the secretomics-based method

Refining the analysis of genome-encoded proteins exhibiting a SP in L. monocytogenes EGD-e and consolidating the results with the prediction of uncleaved SP, it appears that 723 proteins exhibit a SP, i.e. 1 protein with a Tat-SP, 3 an ABC-SP and 5 a FPE-SP as well as 224 with a SP of Type I (SP I), 74 with a SP II and 416 with an uncleaved SP (Unc-SP) (Table S1). The protein secretion pathways (Table 2) are associated with the identification of their respective protein substrates. Altogether and following majority vote approach, 741 proteins appear secreted, including 714 proteins targeted to Sec, 1 to Tat, 4 to ABC (including a leaderless bacteriocin), 3 to holin, 3 to Wss, 5 to FPE and 11 to FEA (Table S1). In addition, 162 IMPs lack a SP but could nonetheless be targeted and translocated within the membrane via YidC. Moreover, 108 proteins primarily predicted as cytoproteins, were further predicted as putatively secreted by NC. Altogether, proteins secreted via these different pathways are predicted as either located to the cytoplasmic membrane (CM; GO:0031226), cell wall (CW; GO:0009275), cell surface (CS; GO:0009986) and/or extracellular milieu (EM; GO:0005576).

Location intrinsic to the cytoplasmic membrane (GO:0031226).

This location splits into two further classes, i.e. either integral to the cytoplasmic membrane (GO:0005887) or anchored to the cytoplasmic membrane (GO:0046658).

Among the 686 predicted IMPs (GO:0005887) as revealed by majority-vote scheme, 524 bear an N-terminal SP including 100 with a SP I, 8 with a SP II and 416 with an Unc-SP (Table S2). From the most recent review [40], the number of protein with a C-terminal “hydrophobic tail” was estimated at only 10, no estimate have ever been provided for N-terminal “hydrophobic tail” proteins [49]. It is also worth noting that the view that only single-spanning IMP (ssIMP) with N-terminal or C-terminal TMD can be considered as surface exposed is simplistic. On the one hand, it is not only the TMD position but rather its orientation that should be considered. On the other hand, even a short strand of amino acid interacting with the external side, such as a loop in multi-spanning IMP (msIMP) can have significant biological function. Based on topogenic elements defined by TMD orientation, where Type I and Type II modules have Nout–Cin and Nin–Cout orientations respectively [50], ssIMPs are further classified into (i) Type I where the ssIMP remains membrane-integrated by a Type I module following cleavage of an N-terminal SP, (ii) Type II where the ssIMP remains membrane-integrated by a Type I module, such an Unc-SP for instance, and (iii) Type III where the ssIMP remains membrane-integrated by a Type I module but is deprivated of an N-terminal SP [51], [52], [53]. Following the standard nomenclature and among the 157 ssIMP, 20 are of Type I (ssIMP I), i.e. exhibiting a cleavable SP, 120 are ssIMP II (including 87 exhibiting an Unc-SP), and 17 ssIMP III (Table S2). Among the 529 ms IMP (msIMP), 88 exhibit a cleavable Sec-SP and 323 an Unc-SP (Table S2). From literature survey, the number of IMP in L. monocytogenes was estimated at 267 [4], proteins with a C-terminal “hydrophobic tail” was estimated at only 10 [40] and no estimate had ever been provided for N-terminal “hydrophobic tail” proteins [49]. Besides the fact that the terminology in usage for IMPs in L. monocytogenes is quite inadequate, the estimate is here considerably changed and qualitatively improved.

Among the 74 lipoproteins (GO:0046658) predicted through majority voting (Table S3), 41 exhibited a glycine residue at position +2 of the cleavage site, conferring them a double possible location in the CM and the EM. Besides, 8 lipoproteins are also predicted as IMP (2 ssIMP I and 6 msIMP); 5 of them (Lmo0269, Lmo0641, Lmo0821, Lmo2687 and Lmo2793) are predicted as lipoproteins for the first time. While one-third of them have unknown function, even after PSI-BLAST search, the majority of them would be substrate-binding protein from ABC transport system. The combination of comprehensive prediction tools allowed here correcting the previous prediction of 68 lipoproteins [40].

Location at the cell wall (GO:0009275).

Parietal proteins (CW-protein) covalently anchored to the cell-wall bear a C-terminal LPXTG domain. Of note, the LPXTG-domain refers to the original regular expression of the motif now surpassed by HMM profiles encompassing the diversity of the cell-wall covalent-anchoring domain, i.e. the LPXTG-protein family [54]. In L. monocytogenes, the 43 LPXGT-proteins previously mentioned have been here predicted [40], [49], including two LPXTG-proteins substrate to SrtB, a lipoprotein (Lmo1136) and a LXPTG protein bearing LysM domains (Lmo880) (Table S4). No YSIRK motif (IPR005877) could be identified within Sec-dependent signal peptides of these proteins; when present, this motif is systematically associated with a LPXTG domain for efficient protein secretion [55] and/or specific final localization within the bacterial cell wall [56].

Besides proteins covalently linked to the cell wall, several CW-proteins are predicted attached by alternative means (Table S4). Contrary to the most recent review [40], only 5 instead of 9 GW-proteins are estimated as CW-proteins. Indeed, the 4 other proteins only exhibit one GW module, which is insufficient for binding to lipoteichoic acids in a noncovalent manner, and are consequently considered here as located extracellularly (Table S6) [57], [58], [59]. In addition, 4 WXL-proteins and one PGBD1-protein (Lmo1851) have been also identified. Previously, the PGBD1-protein had been announced as a IMP because of the presence of a TMD [40]; however, the only TMD predicted is located in the cleavable N-terminal signal peptide of Type I.

Location at cell surface (GO:0009986).

According to GO, this term is intended to annotate gene products that are attached to the CM (intrinsic or loosely bound, GO:0031226 and GO:0031232, respectively) or to the CW (integral or loosely bound, GO:0009275 and GO:0010339, respectively), and thus not only covers IMP, lipoprotein and CW-protein categories but also includes protein complexes (GO:0043234).

Concerning secreted proteins belonging to supramolecular cell-surface appendages (GO:0043234 AND 0009986), neither S-layer proteins, cellulosome components nor prepilins with LPXTG domains could be identified in L. monocytogenes. Nonetheless, 11 flagellar components are predicted secreted and assembled by the FEA (Table S5). In addition, 5 genes encoding for pseudo-prepilin exported by the FPE are identified. Surprisingly, these supramolecular structures had rarely been reported as part of the surface proteins in L. monocytogenes [40], [60].

Location into the extracellular milieu (GO:0005576).

80 exoproteins were predicted as localized extracellularly, including 69 proteins transported via Sec, 1 protein secreted via Tat, 4 bacteriocins translocated via ABC transporters, 3 proteins exported via holins and 3 proteins secreted via Wss (Table S6). While three proteins were predicted with a Tat-SP in the first place (Table S1), only the Dyp (dye decolorising peroxydase)-type peroxydase Lmo0367 was finally predicted as a Tat substrate (Table S6). Among the five predicted bacteriocins, four of them would be secreted via ABC transporters, i.e. Lmo0335, Lmo0615, Lmo2574, and Lmo2753, which is leaderless. As a lactococcin 972 homolog Lmo2776 exhibits an N-terminal SP is most certainly exported in Sec-dependent manner.

Besides, some lipoproteins (Table S3) and primarily cytoplasmic proteins were also predicted with extracellular localization following NC (Table S1). For example, it can be noted that some proteins lacking a SP can be secreted in a route involving the alternative cytosolic ATPase SecA2, a paralogue of SecA, that most certainly converge to the Sec translocon [4], [61].

Comparison with SCL predictors applied to the extracytoproteome of L. monocytogenes

With 8 different predictable SCL in agreement with GO terms in Gram-positive bacteria, i.e. EM (GO:0005876), CS (GO:0009986), cell-surface protein complex (GO:0043234 AND 0009986), CW (GO:0009275), intrinsic to the CM (GO:0031226), integral to the CM (GO:0005887), anchored to the CM (GO:0046658), and cytoplasm (GO:0005737), other predictors never reach such a level of discrimination (Tables 3 and S7). Usually, 4 SCLs are considered as in CELLO [34], PSORTb [62] and Gpos-mPLoc [35] with the cytoplasm, membrane, cell wall and extracellular milieu. SubLoc [63] or LocTree [33] only differentiate cytoplasmic from extracellular location, and Augur the surface from the extracellular milieu [37]. Of note, results of some of these tools are not strictly in accordance with the GO and suffer of several misconceptions regarding the field of protein secretion [1], e.g. (i) SurfG+ [39] and LocateP [38] misleadingly refer to “secreted” for EM location, (ii) they mixes up SCLs with protein categories using descriptors such as LPXTG, multi-transmembrane or N-terminally anchored, which are further not in agreement with standart nomenclature for IMP for instance, (iii) SurfG+ considers “potentially surface exposed (PSE)” needs a certain amino acid length to cross the cell wall, whereas first, this is protein folding that would rather matter, and second, there is no need for a protein to poke out the confine of the cell wall to interact with their environment [60]. In addition and contrary to SubLoc, LocTree, Augur, PSORTb, SurfG+ or LocateP, the secretomics-based method allows attributing multiple SCL for a single protein.

A major and original advance from the secretomics-based analysis resides in the clear and systematic differentiation of the protein category from the protein SCL. To objectively evaluate the performance to predict the protein SCL, the secretomics-based method and other SCL predictors were tested against a dataset of 337 distinct L. monocytogenes proteins, which actual location is defined and referenced (Table S7). SubLoc and LocTree exhibit the lowest overall performance as indicated by MCC values lower than 0.5 (Table 3). While they both show quite similar percentages of overall and single accuracy for SCL prediction, they exhibit an imbalance in their specificity with regard to the sensivity for protein prediction into the cytoplasm. With quite similar percentages of overall and single accuracy for SCL prediction, CELLO and Gpos-mPloc have MCC values just above 0.5. They both present an imbalance in their sensivity with regard to the specificity for SCL prediction to the EM and CW, especially CELLO, which appears inappropriate for EM and CW prediction as its sensitivity is null in the latter case.

With correct overall performance, Augur and PSORTb also present an imbalance in their sensivity/specificity for SCL prediction to the EM (Table 3). While the Augur database identifies the so-called “cell anchor” responsible for cell surface location [37], it suffers from numerous inaccuracies. The presence of LRR and NLPC/P60 domains are considered as cell anchors, which are not. Other well-known cell-wall binding domains, namely PGBD, CWBD, SLHD or WXL domain, are not taken into account. The absence of SP when a “cell anchor” is present does not underdetermine the surface prediction. From there, considering a “single secretome analysis” is performed, when just identifying all secreted proteins with SP and without “cell anchor”, is misleading. In the end, the proportion of genome encoded L. monocytogenes proteins, which SCL prediction is in agreement with the secretomics-based method is quite low, i.e. 67% (Table S7); with a number of Sec-secreted exoproteins estimated at 112, only 49% are in agreement with our present investigation as the remaining proteins exhibit cell-envelope retention domains. Conscious of the ambiguities in the Gram nomenclature especially in the field of protein secretion [1], [2], PSORTb proposes a confusing description to differentiate the analysis for bacteria with Gram positive or negative staining and those “positive with outer membrane” and “negative without outer membrane”. In line with a phylum level perspective on bacterial cell envelope architecture, however, the terminology of monoderm and diderm bacteria, which can even be discriminated between diderm-LPS and diderm-mycolate bacteria, is much more appropriate especially in the field of secretomics [1], [2], [64].

Among tested SCL predictors, SurfG+ and LocateP have the highest overall performances (Table 3). At a genome scale level, the agreement in SCL predictions of L. monocytogenes proteins between these tools and the secretomics-based approach is quite high, i.e. 95% and 92% respectively (Table S7). Back to the performance evaluation metrics, however, LocateP present an imbalance in its sensivity with regard to the specificity for SCL prediction to the EM, as well as to the CW for SurfG+ (Table 3). While the attempt by LocateP to differentiate the IMPs between multi-transmembrane, C-terminally and N-terminally anchored proteins with or without cleavable SP can be acknowledged, it also results in misprediction to the cytoplasm for several ssIMP II or III where the TMD is not located in the N-terminal or C-terminal region and then do not stand within the classification above (Table S7). Similarly in SurfG+ the differentiation of IMP topology into Loop-out, Nterm-out or Cterm-out tails does not compile with the standard nomenclature for IMP. Moreover, LocateP only considers LPXTG motif for SCL to the cell wall; as a result several proteins predicted as IMPs or exoproteins in L. monocytogenes EGD-e are in fact CW-proteins with LysM, GW, PGBD or WXL domains. Finally, the possibility for a single protein to belong to different protein categories is not fully considered by either of these two tools.

Compared to all other SCL predictors, the secretomics-based methods clearly outperform them both considering single and overall performances with accuracy reaching or close to 100% and MCC reaching or close to 1 for all of the 8 SCLs considered by the approach (Table 3). The secretomics-based method is the sole to take into consideration protein cell-surface protein complexes (GO:0043234 AND 0009986), i.e. S-layer (GO:0030115), cellulosome (GO:0043263), flagellum (GO:0019861), pseudo-pilus and pilus (GO:0009289). Actually, PSORTb considers flagellar components but it has not been especially designed to identify substrates to the FEA [62]; as a result, it could only find 3 out of the 11 components secreted by FEA in L. monocytogenes EGD-e and misleadingly predict them as located into the extracellular milieu (Table S7). Protein substrates of FPE or FEA are systematically mispredicted by other SCL predictors, most often as located in the cytoplasm or release in the extracellular milieu.

As the keystone of the secretomics-based method, none of the SCL predictors consider the secretion pathways, including protein routing, transport mechanisms and post-translocational modifications. In fact, LocateP was the best attempt to mimic secretion process but it is restricted to only 3 protein secretion systems, Sec, Tat and ABC protein exporter, plus NC secretion. Another major difference is that the pipeline misses some treads and does not follow the sequential steps of the biology of protein secretion, e.g. TMD prediction is considered before SP, processing of Type 4 prepilin is ruled out, or non-covalent anchoring of proteins to the cell wall is overlooked.

Comparison with experimental proteomic analyses of the extracytoproteome of L. monocytogenes

Globally, the SCL from the secretomics-based approach is in good agreement with experimental data as revealed from the analyses of different extracytoplasmic fractions in L. monocytogenes EGD-e (Table S8). Information on protein trafficking has been experimentally confirmed for only a few proteins and concerns the (i) secretion of proteins dependent on the SPases I (SipX, SipY or SipZ) or Spase II LspA, (ii) SecA2-dependent protein secretion, (iii) protein anchoring to the CM by the prolipoprotein diacylglyceryltransferase Lgt, (iv) non-covalent anchoring to cell wall via GW repeats, or (v) covalent anchoring of proteins to the cell wall via sortases, i.e. SrtA and SrtB. Protein secretion system in used was demonstrated for none of them.

Concerning exoproteins, 26 were experimentally identified, mostly in culture supernatant but also in multi-fractions for some of them, namely membrane, cell wall and/or cell surface exposed fractions (Table S8). SP cleavage site was confirmed for only 3 of them, the involvement of SPase I SipZ for 2 of them, and the SecA2-dependent secretion for one of them. Most exoproteins experimentally identified are predicted secreted via Sec and two of them would be substrates of holin or ABC exporter, respectively.

Among CW-proteins, 35 were experimentally identified, essentially in the cell wall but some of them were also found in the membrane and/or supernatant fractions (Table S8). While the presence in the supernatant of CW-proteins covalently attached to the cell wall can be more surprising than those linked by weak interactions, it must considered that LPXTG-protein SCL is strictly associated with cell-wall biogenesis. Cell wall anchoring via sortases, SrtA or SrtB, was confirmed for 15 LPXTG-proteins but SP cleavage site as well as the involvement of SPase I SipZ and/or SipX was confirmed for only 2 CW-proteins, and the SecA2-dependent secretion for one of them.

With 117 IMP experimentally identified and categorized into ssIMP I, ssIMP II, ssIMP III and msIMP (Table S8), 36 lipoproteins essentially originate from membrane fraction (GO:005624). Again, some of them were also found in the cell wall, supernatant and/or cell surface exposed fractions. For IMP, the involvement of SPase I was confirmed for only one of them as well as SecA2-dependent translocation. Interestingly, most lipoproteins identified in the supernatant exhibited a G at C+2. SecA2-dependent secretion was demonstrated for 4 of them and the involvement of SPase II LspA for only 2 of them. SP cleavage site was confirmed for none. Anchoring to the CM by Lgt was actually confirmed for 26 lipoproteins.

With 192 representatives, the cytoproteins represent the largest category of experimentally identified proteins in extracytoplamic fractions of L. monocytogenes (Table S8). While only 6 of them could be predicted as secreted via NC in the first place, SecA2-dependent translocation was demonstrated for 12 of them. This category stresses the importance to discriminate the GO terms for SCL from the GO terms for the different cell fractions. For example, the presence of ribosomal proteins in membrane fractions is certainly related to SRP (Signal Recognition Particule)-dependent pathway, where ribosomes interact closely with the Sec translocon in the course of co-translational translocation. The SCL of such proteins would then be extrinsic to the CM (GO:0019897) and more precisely on the internal side (GO:0031234), still they would be identified in the membrane fraction (GO:0005624) and even more precisely peripheral to membrane of membrane fraction (GO:0000300) depending on the experimental protocol applied. Besides being secreted in a SecA2-dependent manner, the primary glycolytic enzyme enolase was demonstrated to moonlight when located extracytoplasmically since, together with DnaK, EF-TU and GAPDH, it binds to human plasminogen [44].

Discussion

Compared with other SCL predictors for Gram-positive bacteria available to date, the secretomics-based method clearly outperforms them in term of single and overall performances. Accordingly, (i) protein routing is described and detailed in terms of secretion system, post-translocational modification, cell-envelope retention and/or incorporation to cell-surface supramolecular structures, (ii) the number of distinct SCL, as well as the different protein categories is the most comprehensive for monoderm bacterial cell, (iii) detailed description of protein type and topology is provided, (iv) it unambiguously differentiates the SCL from the protein category, and (v) the possibility for a single protein to belong to different protein categories as well as to have multiple SCL is fully considered. In that sense, the secretomics-based method is much more than a SCL predictor. In L. monocytogenes, this approach conveyed a highly significant quantity of new information, e.g. i) by considering for the first time ABC transporters in protein secretion, ii) presenting the first ever comprehensive estimate of IMPs with careful concern about their topology (msIMP, ssIMP I, II, or III) and double affiliation (to lipoprotein for instance), iii) reporting cell-surface supramolecular structures (flagellum and pseudo-pilus), iv) correcting the prediction of the number of lipoproteins or v) of non-covalently cell-wall anchored proteins by including PGBD1 and excluding proteins exhibiting only one GW domain. While protein attachment to the lipoteichoic acids of the cell wall is clearly modulated by the number of GW modules, i.e. complete release with one GW module and complete retention with 8 GW modules [57], [58], the proportion of retention/release from the bacterial cell for proteins with an in-between number of GW modules remain to be established. As a boundary, this combinatorial approach is based on the biology of protein secretion, which means it relies on the veracity of the knowledge in the field at a given time. This is also the strength of the secretomics-based method as revealed by its performance and its ability to be readily adjusted in light of new findings in the field.

As the keystone of the method, the protein secretion pathways in Gram-positive bacteria are fully taken into consideration and provided for each secreted proteins. Such an in-depth integration of the secretome concept for genomic analysis had never been achieved before. LocateP is the only SCL predictor attempting to reconsider protein secretion systems in monoderm bacteria but it misses some points, e.g. non-covalent cell-wall anchoring of secreted proteins, cell-surface supramolecular structures, or alternative secretion systems to Sec, Tat and ABC pathways [38]. The recently developed CoBaltDB, which combines at once prediction outputs from bioinformatic tools related to protein SCL is a valuable tool to get data prior to the secretomics-based analysis following the decision tree here described for monoderm bacteria [65]. While the basic limit of these tools resides in their lack of flexibility as new predictions tools, new basic findings in the field of protein secretion or genomes are made available, the secretomics-based approach can on the contrary be readily implemented within an Excel data table (Table S1). Nonetheless, we are currently working to make also available in the near future an online tool to ease and fasten the use of the secretomics-based method. Though, as for all pipelines combining results of individual prediction tools and developed to date, the secretomics-based method is based on majority-voting scheme for predictions; for example, LocateP applies unique, consensus and/or unanimously scheme at different stages of the decision trees [38]. Taking lipoprotein prediction tools for instance, they display different overall performances but also distinct precision and recall values [66], which should even be further considered in their single performances for the cleavage site prediction for instance. An awaited development in the field is the attempt to benchmark each individual tool for the subsequent weightings of each prediction.

Besides genomic analysis, this rational strategy is also useful for categorizing, analyzing, and supporting proteomic data. As an example, the prediction of a lipoprotein even by only one tool can be of high relevance if it backups some experimental proteomic data. An interesting aspect is the possibility to generate a database containing mature protein sequences, where any cleavable parts of pre(pro)proteins are removed [43]. The root of proteomics is indeed the protein identification based on expected peptide sizes deduced from in silico trypsinization of amino acid sequence. Confrontation of experimental data with the mature database can then allow experimental determination of the cleavage site in SP or sorting signal [43].

The secretomics-based method allows rigorous interpretation of experimental data by clearly differentiating the predicted SCL GO terms from the experimental cell fraction GO terms. While GO terms are available for proteins identified in membrane fractions (GO:0005624, 0000299, 0000300), they are still awaited for cell wall, surface exposed and supernatant fractions. This further pinpointed some limitations for SCL prediction. Despite very recent improvement in SignalP for example [67], straight prediction in bacteria of uncleaved SP serving as a signal anchor are still unavailable and thus requires intricate analyses (Figure 3). Due to limited experimental data on the IWZ (inner wall zone) [68], for which no GO term is as yet available, prediction of IWZ protein is not reliable. No bioinformatic tools have been developed to predict secreted peripheral proteins (GO:0019897, 0010339). Similarly, no SCL prediction into the cytosol of a host cell for substrates to cytolysin-meditated translocation could be performed.

By considering all molecular mechanisms known to date, the present secretomics-based analysis constitutes the ultimate step in mimicking the protein secretion steps for prediction of the final SCL. Contrary to all other SCL predictors, it does not only supply a static snapshot of protein SCL but also offers the full picture of the secretion process dynamics. Such insight provides strong basis to generate a wealth of in silico hypotheses that further fuel experimental work, e.g. demonstrating the secretion of a protein via a specific secretion system, its post-translational modification by specific sorting pathway, or the biological significance of peripheral proteins. Hopefully, this analytical strategy should be inspirational for the development of rational secretomics-based approaches dedicated to diderm-LPS (Gram-negative) bacteria. As a leitmotif in the field of bacterial protein secretion, the reason for the presence of primarily cytoproteins in the supernatant is an important issue. Besides autolysis, several alternative hypotheses can be formulated (i) allolysis, (ii) uncovered secretion systems, (iii) piggybacking, or (iv) leakage [69]. As exemplified with L. monocytogenes, the secretomics-based method provides the most updated and comprehensive genomic analysis of the extracytoproteome in a monoderm bacterium and should promote further experimental testing based on generated in silico hypotheses.

Materials and Methods

Bioinformatic analyses were performed from web-based servers or under Linux environment with Sun Grid Engine (SGE) server hosted at INRA MIGALE Bioinformatics Platform from (INRA, Jouy-en-Josas, France) with a Intel Quad Core W3520 2.67 GHz. The cluster is organized around 4 computing groups, (i) AMD Quad Core 8354 2.8 GHz, (ii) Intel Quad Core E5520 2.27 GHz, (iii) Intel Quad Core 5340 2.33 GHz, and (iv) Intel Dual Core 5140 2.33 GHz. The complete genome, coding sequences (CDS) and annotation files for L. monocytogenes EGD-e were downloaded from GenBank (ftp://ftp.ncbi.nih.gov/genbank/genomes/Bacteria/Listeria_monocytogenes/).

Searches against various databases were performed using different tools, namely RPS-BLAST v2.2.19 (Reverse Position-Specific BLAST) [70], HMMER v2.3.2 for hidden Markov models (HMM) [71], InterProScan v4.3 [72], or ScanProsite v1.0 [73]. Interrogated databases included InterPro (IPR) v32.0 [74], Pfam (PF) v24.0 [75], SMART (SM) v6.1 [76], TIGRfam (TIGR) v10.1 [77], SuperFamily (SSF) SCOP v1.73 [78], [79], PIRSF v2.74 [80], PRK v3.0 [81], COG v1.0 [82] and Prosite (PS) v20.7 [83]. Position-specific iterated BLAST (PSI-BLAST v2.2.25) [70] searches were executed against UniProtKB v2011_07 [84] until convergence with unlimited number of database sequences for matrix building.

Identification of N-terminal signal peptide

N-terminal signal peptide (SP) were predicted combining results from (i) SignalP v2.0 and v3.0 using both neural network (NN) and hidden Markov model (HMM) [85] with truncation set either at 35, 70, 140 or disabled and predictions with best scores and propensity were considered, (ii) PrediSi v1.0 [86], (iii) Phobius v1.0 [87], (iv) SOSUIsignal v1.0 [88], (v) Signal-3L v1.0 [89], (vi) SIG-PRED v1.0, (vii) Signal-CF v1.0 [90], (vii) SPScan v1.0, an implementation of von Heijne's weight matrix approach with McGeoch criteria where prompted parameters and optional parameter -ADJustscores were used and predictions with best scores and propensity were considered [91], [92], (viii) SPOctopus v1.0 [93], (ix) RPSP v1.0 [94], and (x) Signal-BLAST v1.0 [95]. Tat SP prediction was performed from (i) TatP v1.0 [96], (ii) TATFIND v1.4 [97], (iii) PRED-TAT v1.0 [98], and (iv) InterPro (IPR006311, IPR019546, PS51318, PF10518, TIGR01409). Whenever possible, all these previously cited tools were trained on prokaryotes, bacteria or Gram-positive bacteria. Pseudopilin-like SP were searched from PilFind v1.0 and with ScanProsite syntax [73] for consensus motif [AG]-F-[TSI]-[LTY]-x-[EAF] located between the N- and H-domains of the SP [4]. SP for protein substrates of ABC transporters were identified using BAGEL v1.0 [99]. Uncleaved signal peptides, i.e. N-terminal signal anchors, were predicted from SignalP using HMM trained on eukaryotes [85]. Identification of proteins with no signal peptide and secreted via alternative systems was based on genomic proximity with the genes encoding the respective secretion pathway as described below, coupled to similarity searches, namely substrates to holin and Wss (IPR018921, IPR010310, TIGR03930, PF10663, PF06013). Prediction of non-classical secreted proteins, i.e. lacking a signal peptide and translocated via unknown secretion system, was performed from (i) SecretomeP v2.0 [100], (2) SecretP v2.0 [101] and (3) NClassG+ v1.0 [102] trained on Gram-positive bacteria.

Identification of protein secretion systems

Protein secretion systems were identification performing BLAST search against TC-DB v2011_07 [103], namely Sec (TC#3.A.5), Tat (TC#2.A.64), ABC (TC#3.A.1), FPE (TC#3.A.14), FEA (TC#3.A.6.2), holins (TC#1.E), and Wss (TC#9.A.44) pathways. To discriminate ABC transporters involved in protein secretion, results were compiled with BAGEL [99]. Results were refined by similarity searches, also to identify signal peptidases of Type I (SPase I, IPR000223, COG0681, PRK10861, TIGR02227, SSF51306), SPaseII (COG0597, PRK00376, PRK01574, IPR001872, TIGR00077, PF01252) and prepilin peptidase (COG1989, IPR000045, PF01478), the lipoprotein maturation pathway for some Sec- or Tat substrates, namely in addition to SPase II, the prolipoprotein diacylglyceryltransferase Lgt (COG0682, PRK00052, PRK12437, PRK13108, IPR001640, TIGR00544, PF01790) and the apolipoprotein N-acyltransferase (Lnt) (COG0815, IPR004563, TIGR00546) and the covalent cell-wall anchoring sortase (COG3764, IPR005754, IPR009835, IPR015986, IPR23365, IPR022445, TIGR01076, TIGR03064, TIGR03784, SSF63817, PF04203, PF07170, PIRSF029877, PIRSF030150) pathway for some Sec-substrates.

Identification of transmembrane hydrophobic domains and lipoproteins

Transmembrane hydrophobic α-helice domain (TMD) were predicted combining (i) TMHMM v2.0 [30], (ii) SVMtm v1.0 [104] (iii) THUMBUP v1.0 [105], (iv) SOSUI v1.10 [106], (v) HMMTOP v2.0 [107], (vi) PHDhtm v8.94_69 [108], (vii) UMDHMMTMHP v1.0 [105], (vii) MEMSAT v3.0 [109], (ix) MemBrain v1.0 [110], (x) OCTOPUS v1.0 [111], (xi) MINNOU v1.0 [112], (xii) Philius v1.0 [113], (xiii) SCAMPI v1.0 [114], (xiv) TMPred v1.0 [115], (xv) SPLIT v4.0 [116], and (xvi) TOPCONS v1.0 [117].

For the identification of lipobox, sequences were submitted to (i) DOLOP v2.0 [118], (ii) LipoP v1.0 [29], (iii) SPEPLip v1.0 [119], (iv) LipPred v1.0 [120], (v) ScanProsite for a scan with both PS51257 profile and G+LPP v2.0 pattern [121], (vi) LIPO [122], and (vii) PRED-LIPO [123].

Identification of cell-wall attachment domains

Besides similarity search from InterPro (IPR001899, IPR017502, IPR019931, IPR017503, IPR019948, TIGR01167, TIGR03063, TIGR03065, PF00746, PS50847), LPXTG domain were specifically identified by LPXTG-HMM profile [54] and CW-PRED v2.0 [31]. The concomitant presence of YSIRK domain (IPR005877, PF04650, TIGR01168) within Sec-SP was also checked [55], [56]. WXL domain were found scanning for the motif [LI]-[TE]-W-[TS]-L with ScanProsite where C-terminal location of the motif was also taken into account [124]. Similarity searches were performed for other cell-wall binding motifs, namely (i) LysM (IPR018392, IPR002482, PF01476, SM00257, SSF54106), (ii) GW (SSF82057), (iii) choline-binding domain, also called cell-wall binding domain of type 1 (CWBD1; IPR018337, PF01473, PS51170, SSF69360), (iv) CWBD2 (IPR007253, PF04122), (v) peptidoglycan-binding domain of type 1 (PGBD1; IPR002477, PF01471, SSF47090), (vi) PGBD2 (IPR014927, PF08823), (vii) PGBD3 (IPR018537, PF09374, SSF53955), PGBD4 (IPR022029, PF12229, SSF143985), and S-layer homology domain (SLHD; IPR001119, PF00395, PS51272).

Identification of cell-surface supramolecular structure

Besides the identification of pseudo-pilus components based on the identification of proteins with pseudopilin-like signal peptides as described above, pilus in Gram-positive bacteria was identified from LOCP [125]. Identification of proteins substrates to FEA followed similarity search (IPR001029, IPR00149, IPR010809, IPR003481, IPR013384, IPR002371, IPR020013, IPR005648, IPR021136, IPR001444, IPR006300, IPR006299, IPR001624, COG1345, COG1334, COG1344, COG1256, COG1749, COG1843, COG4786, COG4787, COG1815, COG1558, COG1677, SSF64518, SSF117143), as well as proteins substrates to Wss (IPR018921, IPR010310, TIGR03930, PF10663, PF06013). Identification of cellulosome components was based on dockerin/cohesin domains (IPR002102, IPR018452, IPR016134, IPR002105, IPR018242, IPR009034, IPR002883) and for S-layer on SLHD (IPR001119, COG1361) coupled to PSI-BLAST search.

The secretomics-based method

The method is truly based on the biology of protein secretion in monoderm bacterium (Figure 1). The different types of N-terminal SP targeting a protein to the membrane are discriminated into SP specific for Sec, Tat, ABC and FPE, as well as uncleaved SP (Unc-SP). The protein secretion systems considered are Sec, Tat, ABC, holins, Wss, FPE, FEA, and non-classical secretion (NC). Both transmembrane domain (TMD), including uncleaved SP, or post-translational modifications, i.e. lipobox, are considered for membrane retention. Domains (PGBD, CWBD, WXL, LysM, SLHD, GW) or post-translational modification (LPXTG) are considered for cell-wall retention. Besides flagellum and pseudopilus secreted via the FPE and FEA, S-layer, pilus and cellulosome are also considered as cell-surface supramolecular structures.

The coding sequences (CDS) are screened simultaneously following a flowchart (Figure 2), which provides results from the prediction tools described above and organizes to define the different databases as exemplified in Table S1, namely (i) signal peptide (SP), (ii) type of SP, (iv) protein secretion systems, (v) export (proteins lacking a SP), (iv) transmembrane domain (TMD), and (iv) relevant conserved motifs. Results of different tools for similar prediction as combined into a majority vote approach, i.e. SP, FPE-SP, Tat-SP, ABC-SP, Lipobox, Export and TMD. Using specific tools, proteins exhibiting SP specific for Tat, FPE and ABC transporters are defined as well as proteins with signal anchor, i.e. uncleaved SP, or lipobox in SP of Type II (SP II); combined with the 14 distinct predictions resulting from 11 tools, proteins exhibiting a N-terminal SP are delineated after checking the compatibility with functional annotation verified by similarity search with InterProScan and/or PSI-BLAST. Besides secretion systems, substrates of FEA, holin and Wss are identified following BLAST search against TC-DB. Exported proteins are identified using specific tools for NC. Prediction of proteins with TMD combined 16 tools, 12 of them provide further information on protein topology. Topological information allows to categorised proteins into msIMP, ssIMP I, ssIMP II or ssIMP III. Proteins with a lipobox in SP II are extracted from 8 distinct predictions. Finally, identification of conserved motifs for cell wall anchoring or surface appendages results from interrogation of 10 distinct databases. Parietal proteins are discriminated between LPXTG-, WXL-, SLHD-, LysM-, GW-, CWBD1/2- and/or PGBD1/2/3/4-protein.

From there and as a key step of the the secretomics-based method, the resulting databases are analysed as depicted in the detailed decision trees (Figure 3). As a result of the decision trees, proteins are discriminated into 10 distinct primary categories; IMPs can be further discriminated into multi-spanning IMP (msIMP), single-spanning IMP of Type I (ssIMP I), II (ssIMP II) or III (ssIMP III), whereas parietal protein (CW-protein) can be dicrimnated into 12 subtypes (LPXTG-protein, LysM protein, CWBD1-protein,...etc). In agreement with GO terms, the secretomics-based method provides 4 primary SCL, CP (GO:0005737), CM (GO:0005886), CW (GO:0009275) and EM (GO:0005576); protein complex (GO:0043234), cell surface (GO:0009986), integral (GO:0005887) and anchored to CM (GO:0046658) SCL are further discriminated.

Subcellular localization predictors

Predictions for SCL were performed from pipelines, i.e. Augur v2010_07 [37], LocateP v1.0 [38] and SurfG+ [39], as well as support vector machines (SVMs), i.e. SubLoc v1.0 [63] and LocTree v1.0 [33] trained on prokaryotes where only extracellular and cytoplasmic SCL are considered, whereas CELLO v2.5 [34], PSORTb v3.0.2 [62] and Gpos-mPLoc v1.0 [35] trained on Gram-positive bacteria considered extracellular, cell wall, membrane and cytoplasmic subcellular compartments as well as multiple localization sites.

Performance evaluation metrics

Sensitivity, specificity, accuracy and MCC (Matthews Correlation Coefficient) of the SCL methods were calculated from the four basic statistics, i.e. true positives (TP), false negatives (FN), false positives (FP) and true negatives (TN) [126]. Sensitivity was calculated as TP/(TP+FN), specificity as TN/(TN+FP), accuracy as (TP+TN)/(TP+FN+FP+TN), and MCC as [(TP×TN)-(FN×FP)]/√[(TP+FN)×(TP+FP)×(TN+FP)×(TN+FN)]. Sensitivity, specificity and accuracy were expressed as percentage, whereas MMC values can vary between −1 and 1. Taking the extracellular milieu (EM) as an example of SCL prediction, (i) TP corresponds to instances both predicted and actually located in the EM, (ii) FN to instances actually located in the EM but not predicted in the EM, (iii) FP to instances predicted in the EM but actually not located in the EM, and (iv) TN to instances both predicted and actually not located in the EM. Besides the overall performance of a given tool, single performances were evaluated for each of the different SCL considered. A dataset was constructed from L. monocytogenes proteins, which location was experimentally confirmed and/or acknowledged from literature survey (Table S7). It consisted of 337 distinct proteins localized in the (i) extracellular milieu (EM; 49 proteins), (ii) cell surface (CS; 199 proteins), (iii) protein complex (PC; 16 proteins), (iv) cell wall (CW; 30 proteins), (v) intrinsic to the cytoplasmic membrane (CM; 148 proteins), (vi) anchored to the CM (aCM; 37 proteins), (vii) integral to the CM (iCM; 114 proteins), and/or (viii) cytoplasm (CP; 117 proteins) (Table S7).

Supporting Information

Table S2.

Summarised information about protein categories, secretion pathways and GO terms for IMPs, lipoproteins, cell-wall proteins, subunits of supramolecular cell-surface appendages and exoproteins, respectively, as predicted by the secretomics-based method in L. monocytogenes EGD-e.

https://doi.org/10.1371/journal.pone.0042982.s002

(PDF)

Table S3.

Summarised information about protein categories, secretion pathways and GO terms for IMPs, lipoproteins, cell-wall proteins, subunits of supramolecular cell-surface appendages and exoproteins, respectively, as predicted by the secretomics-based method in L. monocytogenes EGD-e.

https://doi.org/10.1371/journal.pone.0042982.s003

(PDF)

Table S4.

Summarised information about protein categories, secretion pathways and GO terms for IMPs, lipoproteins, cell-wall proteins, subunits of supramolecular cell-surface appendages and exoproteins, respectively, as predicted by the secretomics-based method in L. monocytogenes EGD-e.

https://doi.org/10.1371/journal.pone.0042982.s004

(PDF)

Table S5.

Summarised information about protein categories, secretion pathways and GO terms for IMPs, lipoproteins, cell-wall proteins, subunits of supramolecular cell-surface appendages and exoproteins, respectively, as predicted by the secretomics-based method in L. monocytogenes EGD-e.

https://doi.org/10.1371/journal.pone.0042982.s005

(PDF)

Table S6.

Summarised information about protein categories, secretion pathways and GO terms for IMPs, lipoproteins, cell-wall proteins, subunits of supramolecular cell-surface appendages and exoproteins, respectively, as predicted by the secretomics-based method in L. monocytogenes EGD-e.

https://doi.org/10.1371/journal.pone.0042982.s006

(PDF)

Acknowledgments

This work was supported in part by INRA (French National Institute for Agronmical Research), and the European Framework Programme 6 (FP6) with the ProSafeBeef (Advancing Beef Safety and Quality through Research and Innovation) research consortium (http://www.prosafebeef.eu/). Sandra Renier is a PhD research fellow granted by the “Ministère de l′Education Nationale, de l′Education Supérieure et de la Recherche”. The authors are obliged to Valentin Loux (INRA, Jouy-en-Josas, France), Menglong Li (Sichuan University, Chengdu, China) and Manuel A. Patarroyo (Universidad del Rosario, Bogota, Colombia) for batch query against SurfG+, SecretP and NClassG+, respectively. We are grateful to the INRA MIGALE bioinformatics platform (http://migale.jouy.inra.fr) for providing computational resources. The authors are thankful to anonymous reviewers for their very constructive comments.

Author Contributions

Conceived and designed the experiments: SR MD. Performed the experiments: SR PM MD. Analyzed the data: SR RT MH MD. Contributed reagents/materials/analysis tools: SR PM MD. Wrote the paper: SR RT MH MD.

References

  1. 1.Desvaux M, Hébraud M, Talon R, Henderson IR (2009) Secretion and subcellular localizations of bacterial proteins: a semantic awareness issue. Trends Microbiol 17: 139–145.
  2. 2.Sutcliffe IC (2010) A phylum level perspective on bacterial cell envelope architecture. Trends Microbiol 18: 464–470.
  3. 3.Sutcliffe IC (2011) Cell envelope architecture in the Chloroflexi: a shifting frontline in a phylogenetic turf war. Environ Microbiol.
  4. 4.Desvaux M, Hébraud M (2006) The protein secretion systems in Listeria: inside out bacterial virulence. FEMS Microbiol Rev 30: 774–805.
  5. 5.Sibbald MJ, Ziebandt AK, Engelmann S, Hecker M, de Jong A, et al. (2006) Mapping the pathways to staphylococcal pathogenesis by comparative secretomics. Microbiol Mol Biol Rev 70: 755–788.
  6. 6.Kleerebezem M, Hols P, Bernard E, Rolain T, Zhou M, et al. (2010) The extracellular biology of the lactobacilli. FEMS Microbiol Rev 34: 199–230.
  7. 7.DiGiuseppe Champion PA, Cox JS (2007) Protein secretion systems in Mycobacteria. Cell Microbiol 9: 1376–1384.
  8. 8.Chater KF, Biro S, Lee KJ, Palmer T, Schrempf H (2010) The complex extracellular biology of Streptomyces. FEMS Microbiol Rev 34: 171–198.
  9. 9.Tjalsma H, Antelmann H, Jongbloed JD, Braun PG, Darmon E, et al. (2004) Proteomics of protein secretion by Bacillus subtilis: separating the “secrets” of the secretome. Microbiol Mol Biol Rev 68: 207–233.
  10. 10.Tjalsma H, Bolhuis A, Jongbloed JD, Bron S, van Dijl JM (2000) Signal peptide-dependent protein transport in Bacillus subtilis: a genome-based survey of the secretome. Microbiol Mol Biol Rev 64: 515–547.
  1. 11.Simonen M, Palva I (1993) Protein secretion in Bacillus species. Microbiol Rev 57: 109–137.
  1. 12.Tan KS, Wee BY, Song KP (2001) Evidence for holin function of tcdE gene in the pathogenicity of Clostridium difficile. J Med Microbiol 50: 613–619.
  1. 13.Harboe M, Oettinger T, Wiker HG, Rosenkrands I, Andersen P (1996) Evidence for occurrence of the ESAT-6 protein in Mycobacterium tuberculosis and virulent Mycobacterium bovis and for its absence in Mycobacterium bovis BCG. Infect Immun 64: 16–22.
  1. 14.Schaerlaekens K, Schierova M, Lammertyn E, Geukens N, Anne J, et al. (2001) Twin-arginine translocation pathway in Streptomyces lividans. J Bacteriol 183: 6727–6732.
  1. 15.Chen I, Provvedi R, Dubnau D (2006) A macromolecular complex formed by a pilin-like protein in competent Bacillus subtilis. J Biol Chem 281: 21720–21727.
  1. 16.Desvaux M, Parham NJ, Scott-Tucker A, Henderson IR (2004) The general secretory pathway: a general misnomer? Trends Microbiol 12: 306–309.
  1. 17.Desvaux M, Hebraud M, Henderson IR, Pallen MJ (2006) Type III secretion: what's in a name? Trends Microbiol 14: 157–160.
  1. 18.Desvaux M, Khan A, Scott-Tucker A, Chaudhuri RR, Pallen MJ, et al. (2005) Genomic analysis of the protein secretion systems in Clostridium acetobutylicum ATCC824. Biochim Biophys Acta-Mol Cell Res 1745: 223–253.
  1. 19.Sutcliffe IC (2011) New insights into the distribution of WXG100 protein secretion systems. Antonie Van Leeuwenhoek.
  1. 20.Desvaux M, Hébraud M, Talon R, Henderson IR (2009) Outer membrane translocation: numerical protein secretion nomenclature in question in mycobacteria. Trends Microbiol 17: 338–340.
  1. 21.Bleves S, Viarre V, Salacha R, Michel GP, Filloux A, et al. (2010) Protein secretion systems in Pseudomonas aeruginosa: A wealth of pathogenic weapons. Int J Med Microbiol 300: 534–543.
  1. 22.Fagerlund A, Lindback T, Granum PE (2010) Bacillus cereus cytotoxins Hbl, Nhe and CytK are secreted via the Sec translocation pathway. BMC Microbiol 10: 304.
  1. 23.Romine MF (2011) Genome-wide protein localization prediction strategies for Gram negative bacteria. BMC Genomics 12 Suppl 1: S1.
  1. 24.Desvaux M, Khan A, Beatson S, Scott-Tucker A, Henderson IR (2005) Protein secretion systems in Fusobacterium nucleatum: genomic identification of Type 4 piliation and complete Type V pathways brings new insight in mechanisms of pathogenesis. Biochim Biophys Acta-Biomembr 1713: 92–112.
  1. 25.Sarvas M, Harwood CR, Bron S, van Dijl JM (2004) Post-translocational folding of secretory proteins in Gram-positive bacteria. Biochim Biophys Acta-Mol Cell Res 1694: 311–327.
  1. 26.Antelmann H, Tjalsma H, Voigt B, Ohlmeier S, Bron S, et al. (2001) A proteomic view on genome-based signal peptide predictions. Genome Res 11: 1484–1502.
  1. 27.Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, et al. (2000) Gene ontology: tool for the unification of biology. The Gene Ontology Consortium. Nat Genet 25: 25–29.
  1. 28.Nielsen H, Brunak S, von Heijne G (1999) Machine learning approaches for the prediction of signal peptides and other protein sorting signals. Protein Eng 12: 3–9.
  1. 29.Juncker AS, Willenbrock H, Von Heijne G, Brunak S, Nielsen H, et al. (2003) Prediction of lipoprotein signal peptides in Gram-negative bacteria. Protein Sci 12: 1652–1662.
  1. 30.Sonnhammer EL, von Heijne G, Krogh A (1998) A hidden Markov model for predicting transmembrane helices in protein sequences. Proc Int Conf Intell Syst Mol Biol 6: 175–182.
  1. 31.Litou ZI, Bagos PG, Tsirigos KD, Liakopoulos TD, Hamodrakas SJ (2008) Prediction of cell wall sorting signals in Gram-positive bacteria with a hidden Markov model: application to complete genomes. J Bioinform Comput Biol 6: 387–401.
  1. 32.Gardy JL, Laird MR, Chen F, Rey S, Walsh CJ, et al. (2005) PSORTb v.2.0: expanded prediction of bacterial protein subcellular localization and insights gained from comparative proteome analysis. Bioinformatics 21: 617–623.
  1. 33.Nair R, Rost B (2005) Mimicking cellular sorting improves prediction of subcellular localization. J Mol Biol 348: 85–100.
  1. 34.Yu CS, Chen YC, Lu CH, Hwang JK (2006) Prediction of protein subcellular localization. Proteins 64: 643–651.
  1. 35.Shen HB, Chou KC (2009) Gpos-mPLoc: a top-down approach to improve the quality of predicting subcellular localization of Gram-positive bacterial proteins. Protein Pept Lett 16: 1478–1484.
  1. 36.Dzeroski S, Zenko B (2004) Is combining classifiers with stacking better than selecting the best one? Machine Learning 54: 255–273.
  1. 37.Billion A, Ghai R, Chakraborty T, Hain T (2006) Augur – a computational pipeline for whole genome microbial surface protein prediction and classification. Bioinformatics 22: 2819–2820.
  1. 38.Zhou M, Boekhorst J, Francke C, Siezen RJ (2008) LocateP: genome-scale subcellular-location predictor for bacterial proteins. BMC Bioinformatics 9: 173.
  1. 39.Barinov A, Loux V, Hammani A, Nicolas P, Langella P, et al. (2009) Prediction of surface exposed proteins in Streptococcus pyogenes, with a potential application to other Gram-positive bacteria. Proteomics 9: 61–73.
  1. 40.Bierne H, Cossart P (2007) Listeria monocytogenes surface proteins: from genome predictions to function. Microbiol Mol Biol Rev 71: 377–397.
  1. 41.Renier S, Hebraud M, Desvaux M (2011) Molecular biology of surface colonization by Listeria monocytogenes: an additional facet of an opportunistic Gram-positive foodborne pathogen. Environ Microbiol 13: 835–850.
  1. 42.Pucciarelli MG, Calvo E, Sabet C, Bierne H, Cossart P, et al. (2005) Identification of substrates of the Listeria monocytogenes sortases A and B by a non-gel proteomic analysis. Proteomics. 10.1002/pmic.200402075.
  1. 43.Desvaux M, Dumas E, Chafsey I, Chambon C, Hebraud M (2010) Comprehensive appraisal of the extracellular proteins from a monoderm bacterium: theoretical and empirical exoproteomes of Listeria monocytogenes EGD-e by secretomics. J Proteome Res 9: 5076–5092.
  1. 44.Schaumburg J, Diekmann O, Hagendorff P, Bergmann S, Rohde M, et al. (2004) The cell wall subproteome of Listeria monocytogenes. Proteomics 4: 2991–3006.
  1. 45.Calvo E, Pucciarelli MG, Bierne H, Cossart P, Albar JP, et al. (2005) Analysis of the Listeria cell wall proteome by two-dimensional nanoliquid chromatography coupled to mass spectrometry. Proteomics 5: 433–443.
  1. 46.Garcia-Del Portillo F, Calvo E, D'Orazio V, Pucciarelli MG (2011) Association of ActA to the peptidoglycan revealed by cell wall proteomics of intracellular Listeria monocytogenes. J Biol Chem.
  1. 47.Trost M, Wehmhoner D, Karst U, Dieterich G, Wehland J, et al. (2005) Comparative proteome analysis of secretory proteins from pathogenic and nonpathogenic Listeria species. Proteomics 5: 1544–1557.
  1. 48.Tjalsma H, van Dijl JM (2005) Proteomics-based consensus prediction of protein retention in a bacterial membrane. Proteomics 5: 4472–4482.
  1. 49.Glaser P, Frangeul L, Buchrieser C, Rusniok C, Amend A, et al. (2001) Comparative genomics of Listeria species. Science 294: 849–852.
  1. 50.White SH, von Heijne G (2004) The machinery of membrane protein assembly. Curr Opin Struct Biol 14: 397–404.
  1. 51.Goder V, Spiess M (2001) Topogenesis of membrane proteins: determinants and dynamics. FEBS Lett 504: 87–93.
  1. 52.Desvaux M, Hébraud M (2008) Chapter 12: Analysis of cell envelope proteins. In: Handbook of Listeria monocytogenes. Ed: Liu D., CRC Press, Taylor and Francis Group, Boca Raton, Florida, USA. 359–393.
  2. 53.Higy M, Junne T, Spiess M (2004) Topogenesis of membrane proteins at the endoplasmic reticulum. Biochemistry 43: 12716–12722.
  1. 54.Boekhorst J, de Been MW, Kleerebezem M, Siezen RJ (2005) Genome-wide detection and analysis of cell wall-bound proteins with LPXTG-like sorting motifs. J Bacteriol 187: 4928–4934.
  1. 55.Bae T, Schneewind O (2003) The YSIRK-G/S motif of staphylococcal protein A and its role in efficiency of signal peptide processing. J Bacteriol 185: 2910–2919.
  1. 56.DeDent A, Bae T, Missiakas DM, Schneewind O (2008) Signal peptides direct surface proteins to two distinct envelope locations of Staphylococcus aureus. EMBO J 27: 2656–2668.
  1. 57.Braun L, Dramsi S, Dehoux P, Bierne H, Lindahl G, et al. (1997) InlB: an invasion protein of Listeria monocytogenes with a novel type of surface association. Mol Microbiol 25: 285–294.
  1. 58.Jonquieres R, Bierne H, Fiedler F, Gounon P, Cossart P (1999) 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 34: 902–914.
  1. 59.Marino M, Banerjee M, Jonquieres R, Cossart P, Ghosh P (2002) GW domains of the Listeria monocytogenes invasion protein InlB are SH3-like and mediate binding to host ligands. EMBO J 21: 5623–5634.
  1. 60.Desvaux M, Dumas E, Chafsey I, Hébraud M (2006) Protein cell surface display in Gram-positive bacteria: from single protein to macromolecular protein structure. FEMS Microbiol Lett 256: 1–15.
  1. 61.Rigel NW, Braunstein M (2008) A new twist on an old pathway – accessory Sec systems. Mol Microbiol 69: 291–302.
  1. 62.Yu NY, Wagner JR, Laird MR, Melli G, Rey S, et al. (2010) PSORTb 3.0: improved protein subcellular localization prediction with refined localization subcategories and predictive capabilities for all prokaryotes. Bioinformatics 26: 1608–1615.
  1. 63.Hua S, Sun Z (2001) Support vector machine approach for protein subcellular localization prediction. Bioinformatics 17: 721–728.
  1. 64.Gupta RS (2011) Origin of diderm (Gram-negative) bacteria: antibiotic selection pressure rather than endosymbiosis likely led to the evolution of bacterial cells with two membranes. Antonie Van Leeuwenhoek 100: 171–182.
  1. 65.Goudenège D, Avner S, Lucchetti-Miganeh C, Barloy-Hubler F (2010) CoBaltDB: Complete bacterial and archaeal orfeomes subcellular localization database and associated resources. BMC Microbiol 10: 88.
  1. 66.Rahman O, Cummings SP, Harrington DJ, Sutcliffe IC (2008) Methods for the bioinformatic identification of bacterial lipoproteins encoded in the genomes of Gram-positive bacteria. World Journal of Microbiology & Biotechnology 24: 2377–2382.
  1. 67.Petersen TN, Brunak S, von Heijne G, Nielsen H (2011) SignalP 4.0: discriminating signal peptides from transmembrane regions. Nat Methods 8: 785–786.
  1. 68.Matias VR, Beveridge TJ (2005) Cryo-electron microscopy reveals native polymeric cell wall structure in Bacillus subtilis 168 and the existence of a periplasmic space. Mol Microbiol 56: 240–251.
  1. 69.Desvaux M (2012) Contribution of holins to protein trafficking: secretion, leakage or lysis? Trends Microbiol 20: 259–261.
  1. 70.Altschul SF, Madden TL, Schaffer AA, Zhang J, Zhang Z, et al. (1997) Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res 25: 3389–3402.
  1. 71.Eddy SR (1996) Hidden Markov models. Curr Opin Struct Biol 6: 361–365.
  1. 72.Zdobnov EM, Apweiler R (2001) InterProScan – an integration platform for the signature-recognition methods in InterPro. Bioinformatics 17: 847–848.
  1. 73.de Castro E, Sigrist CJ, Gattiker A, Bulliard V, Langendijk-Genevaux PS, et al. (2006) ScanProsite: detection of PROSITE signature matches and ProRule-associated functional and structural residues in proteins. Nucleic Acids Res 34: W362–365.
  1. 74.Mulder NJ, Apweiler R, Attwood TK, Bairoch A, Bateman A, et al. (2007) New developments in the InterPro database. Nucleic Acids Res 35: D224–D228.
  1. 75.Bateman A, Coin L, Durbin R, Finn RD, Hollich V, et al. (2004) The Pfam protein families database. Nucleic Acids Res 32: D138–141.
  1. 76.Schultz J, Milpetz F, Bork P, Ponting CP (1998) SMART, a simple modular architecture research tool: identification of signaling domains. Proc Natl Acad Sci U S A 95: 5857–5864.
  1. 77.Selengut JD, Haft DH, Davidsen T, Ganapathy A, Gwinn-Giglio M, et al. (2007) TIGRFAMs and Genome Properties: tools for the assignment of molecular function and biological process in prokaryotic genomes. Nucleic Acids Res 35: D260–D264.
  1. 78.Murzin AG, Brenner SE, Hubbard T, Chothia C (1995) SCOP: a structural classification of proteins database for the investigation of sequences and structures. J Mol Biol 247: 536–540.
  1. 79.Wilson D, Madera M, Vogel C, Chothia C, Gough J (2007) The SUPERFAMILY database in 2007: families and functions. Nucleic Acids Res 35: D308–D313.
  1. 80.Wu CH, Nikolskaya A, Huang H, Yeh LS, Natale DA, et al. (2004) PIRSF: family classification system at the Protein Information Resource. Nucleic Acids Res 32: D112–D114.
  1. 81.Oneill K, Klimke W, Tatusova T (2007) Protein clusters: a collection of proteins grouped by sequence similarity and function.
  1. 82.Tatusov RL, Natale DA, Garkavtsev IV, Tatusova TA, Shankavaram UT, et al. (2001) The COG database: new developments in phylogenetic classification of proteins from complete genomes. Nucleic Acids Res 29: 22–28.
  1. 83.Hulo N, Bairoch A, Bulliard V, Cerutti L, De Castro E, et al. (2006) The PROSITE database. Nucleic Acids Res 34: D227–230.
  1. 84.Wu CH, Apweiler R, Bairoch A, Natale DA, Barker WC, et al. (2006) The Universal Protein Resource (UniProt): an expanding universe of protein information. Nucleic Acids Res 34: D187–D191.
  1. 85.Bendtsen JD, Nielsen H, von Heijne G, Brunak S (2004) Improved prediction of signal peptides: SignalP 3.0. J Mol Biol 340: 783–795.
  1. 86.Hiller K, Grote A, Scheer M, Munch R, Jahn D (2004) PrediSi: prediction of signal peptides and their cleavage positions. Nucleic Acids Res 32: W375–W379.
  1. 87.Käll L, Krogh A, Sonnhammer EL (2004) A combined transmembrane topology and signal peptide prediction method. J Mol Biol 338: 1027–1036.
  1. 88.Gomi M, Sonoyama M, Mitaku S (2004) High performance system for signal peptide prediction: SOSUIsignal. Chem-Bio Info J 4: 142–147.
  1. 89.Shen HB, Chou KC (2007) Signal-3L: a 3-layer approach for predicting signal peptides. Biochem Biophys Res Commun 363: 297–303.
  1. 90.Chou KC, Shen HB (2007) Signal-CF: A subsite-coupled and window-fusing approach for predicting signal peptides. Biochemical and Biophysical Research Communications 357: 633–640.
  1. 91.von Heijne G (1986) A new method for predicting signal sequence cleavage sites. Nucleic Acids Res 14: 4683–4690.
  1. 92.McGeoch DJ (1985) On the predictive recognition of signal peptide sequences. Virus Res 3: 271–286.
  1. 93.Viklund H, Bernsel A, Skwark M, Elofsson A (2008) SPOCTOPUS: a combined predictor of signal peptides and membrane protein topology. Bioinformatics 24: 2928–2929.
  1. 94.Plewczynski D, Slabinski L, Tkacz A, Kajan L, Holm L, et al. (2007) The RPSP: Web server for prediction of signal peptides. Polymer 48: 5493–5496.
  1. 95.Frank K, Sippl MJ (2008) High-performance signal peptide prediction based on sequence alignment techniques. Bioinformatics 24: 2172–2176.
  1. 96.Bendtsen JD, Nielsen H, Widdick D, Palmer T, Brunak S (2005) Prediction of twin-arginine signal peptides. BMC Bioinformatics 6: 167.
  1. 97.Rose RW, Bruser T, Kissinger JC, Pohlschroder M (2002) Adaptation of protein secretion to extremely high-salt conditions by extensive use of the twin-arginine translocation pathway. Mol Microbiol 45: 943–950.
  1. 98.Bagos PG, Nikolaou EP, Liakopoulos TD, Tsirigos KD (2010) Combined prediction of Tat and Sec signal peptides with hidden Markov models. Bioinformatics 26: 2811–2817.
  1. 99.de Jong A, van Hijum SA, Bijlsma JJ, Kok J, Kuipers OP (2006) BAGEL: a web-based bacteriocin genome mining tool. Nucleic Acids Res 34: W273–279.
  1. 100.Bendtsen JD, Kiemer L, Fausboll A, Brunak S (2005) Non-classical protein secretion in bacteria. BMC Microbiol 5: 58.
  1. 101.Yu L, Guo Y, Li Y, Li G, Li M, et al. (2010) SecretP: identifying bacterial secreted proteins by fusing new features into Chou's pseudo-amino acid composition. J Theor Biol 267: 1–6.
  1. 102.Restrepo-Montoya D, Pino C, Nino LF, Patarroyo ME, Patarroyo MA (2011) NClassG+: A classifier for non-classically secreted Gram-positive bacterial proteins. BMC Bioinformatics 12: 21.
  1. 103.Saier MH Jr, Yen MR, Noto K, Tamang DG, Elkan C (2009) The Transporter Classification Database: recent advances. Nucleic Acids Res 37: D274–278.
  1. 104.Yuan Z, Mattick JS, Teasdale RD (2004) SVMtm: support vector machines to predict transmembrane segments. J Comput Chem 25: 632–636.
  1. 105.Zhou H, Zhou Y (2003) Predicting the topology of transmembrane helical proteins using mean burial propensity and a hidden-Markov-model-based method. Protein Sci 12: 1547–1555.
  1. 106.Hirokawa T, Boon-Chieng S, Mitaku S (1998) SOSUI: classification and secondary structure prediction system for membrane proteins. Bioinformatics 14: 378–379.
  1. 107.Tusnády GE, Simon I (2001) The HMMTOP transmembrane topology prediction server. Bioinformatics 17: 849–850.
  1. 108.Rost B, Casadio R, Fariselli P, Sander C (1995) Transmembrane helices predicted at 95% accuracy. Protein Sci 4: 521–533.
  1. 109.Jones DT (2007) Improving the accuracy of transmembrane protein topology prediction using evolutionary information. Bioinformatics 23: 538–544.
  1. 110.Shen H, Chou JJ (2008) MemBrain: improving the accuracy of predicting transmembrane helices. PLoS One 3: e2399.
  1. 111.Viklund H, Elofsson A (2008) OCTOPUS: improving topology prediction by two-track ANN-based preference scores and an extended topological grammar. Bioinformatics 24: 1662–1668.
  1. 112.Cao B, Porollo A, Adamczak R, Jarrell M, Meller J (2006) Enhanced recognition of protein transmembrane domains with prediction-based structural profiles. Bioinformatics 22: 303–309.
  1. 113.Reynolds SM, Kall L, Riffle ME, Bilmes JA, Noble WS (2008) Transmembrane topology and signal peptide prediction using dynamic bayesian networks. PLoS Comput Biol 4: e1000213.
  1. 114.Bernsel A, Viklund H, Falk J, Lindahl E, von Heijne G, et al. (2008) Prediction of membrane-protein topology from first principles. Proc Natl Acad Sci U S A 105: 7177–7181.
  1. 115.Hofmann K, Stoffel W (1993) TMbase – A database of membrane spanning proteins segments. Biol Chem Hoppe-Seyler 374: 166.
  1. 116.Juretic D, Zoranic L, Zucic D (2002) Basic charge clusters and predictions of membrane protein topology. Journal of Chemical Information and Computer Sciences 42: 620–632.
  1. 117.Bernsel A, Viklund H, Hennerdal A, Elofsson A (2009) TOPCONS: consensus prediction of membrane protein topology. Nucleic Acids Res 37: W465–468.
  1. 118.Babu MM, Priya ML, Selvan AT, Madera M, Gough J, et al. (2006) A database of bacterial lipoproteins (DOLOP) with functional assignments to predicted lipoproteins. J Bacteriol 188: 2761–2773.
  1. 119.Fariselli P, Finocchiaro G, Casadio R (2003) SPEPlip: the detection of signal peptide and lipoprotein cleavage sites. Bioinformatics 19: 2498–2499.
  1. 120.Taylor PD, Toseland CP, Attwood TK, Flower DR (2006) LipPred: A web server for accurate prediction of lipoprotein signal sequences and cleavage sites. Bioinformation 1: 335–338.
  1. 121.Hutchings MI, Palmer T, Harrington DJ, Sutcliffe IC (2009) Lipoprotein biogenesis in Gram-positive bacteria: knowing when to hold 'em, knowing when to fold 'em. Trends Microbiol 17: 13–21.
  1. 122.Berven FD, Karlsen OS, Straume AH, Flikka K, Murrell JC, et al. (2006) Analysing the outer membrane subproteome of Methylococcus capsulatus (Bath) using proteomics and novel biocomputing tools. Archives of Microbiology 184: 362–377.
  1. 123.Bagos PG, Tslrigos KD, Liakopoulos TD, Hamodrakas SJ (2008) Prediction of lipoprotein signal peptides in Gram-positive bacteria with a hidden Markov model. Journal of Proteome Research 7: 5082–5093.
  1. 124.van Pijkeren JP, Canchaya C, Ryan KA, Li Y, Claesson MJ, et al. (2006) Comparative and functional analysis of sortase-dependent proteins in the predicted secretome of Lactobacillus salivarius UCC118. Appl Environ Microbiol 72: 4143–4153.
  1. 125.Plyusnin I, Holm L, Kankainen M (2009) LOCP – locating pilus operons in Gram-positive bacteria. Bioinformatics 25: 1187–1188.
  1. 126.Gardy JL, Brinkman FS (2006) Methods for predicting bacterial protein subcellular localization. Nat Rev Microbiol 4: 741–751.