P2RP: a web-based framework for the identification and analysis of regulatory proteins in prokaryotic genomes
© Barakat et al.; licensee BioMed Central Ltd. 2013
Received: 25 January 2013
Accepted: 19 March 2013
Published: 20 April 2013
Regulatory proteins (RPs) such as transcription factors (TFs) and two-component system (TCS) proteins control how prokaryotic cells respond to changes in their external and/or internal state. Identification and annotation of TFs and TCSs is non-trivial, and between-genome comparisons are often confounded by different standards in annotation. There is a need for user-friendly, fast and convenient tools to allow researchers to overcome the inherent variability in annotation between genome sequences.
We have developed the web-server P2RP (Predicted Prokaryotic Regulatory Proteins), which enables users to identify and annotate TFs and TCS proteins within their sequences of interest. Users can input amino acid or genomic DNA sequences, and predicted proteins therein are scanned for the possession of DNA-binding domains and/or TCS domains. RPs identified in this manner are categorised into families, unambiguously annotated, and a detailed description of their features generated, using an integrated software pipeline. P2RP results can then be outputted in user-specified formats.
Biologists have an increasing need for fast and intuitively usable tools, which is why P2RP has been developed as an interactive system. As well as assisting experimental biologists to interrogate novel sequence data, it is hoped that P2RP will be built into genome annotation pipelines and re-annotation processes, to increase the consistency of RP annotation in public genomic sequences. P2RP is the first publicly available tool for predicting and analysing RP proteins in users’ sequences. The server is freely available and can be accessed along with documentation at http://www.p2rp.org.
In prokaryotes biological processes tend to be regulated at the level of transcription, with subsets of genes/operons being up/down-regulated by specific DNA-binding proteins known as transcription factors (TFs). TFs can be divided into a few major categories, including sigma factors (SFs), one-component systems (OCSs) and response regulators (RRs), and the DNA-binding activity of these proteins is often regulated. SFs are the specificity-conferring sub-units of RNA polymerase holoenzymes, and they direct the transcription machinery towards particular promoter sequences . The activity of SFs is often regulated by accessory proteins such as anti-SFs, which bind to and inhibit specific SFs. In addition to DNA-binding domains, OCSs possess sensory domains, which modulate DNA-binding activity according to the presence/absence of a particular stimulus . Finally, the DNA-binding activity of RRs is regulated by the phosphorylation-state of their receiver domains, which can be phosphorylated by partner receptor kinase proteins called histidine kinases (HKs). Together HKs and their partner RRs (including non-DNA-binding RRs), form two-component systems (TCSs), which are the dominant phosphorylation-dependent signal transduction pathways of prokaryotes .
A typical prokaryotic genome encodes around 5% TFs  and 1.5% TCS proteins , and for most regulatory proteins (RPs), multiple homologues are usually found in each genome. Therefore for RPs, sequence similarity does not necessarily imply a similar functional role, and annotation of RPs by sequence similarity has resulted in many erroneous annotations.
Over-specific annotation is a common problem. For example, the E. coli PhoB/OmpR family of RRs regulate diverse processes, including potassium homeostasis (KdpE), copper tolerance (CusR) and trimethylamine N-oxide respiration (TorR), in addition to phosphatase expression (PhoB) and osmoregulation (OmpR) [6–10]. However multiple PhoB/OmpR family members in a genome are sometimes ascribed the same role. For example, Clostridium botulinum B str. Eklund 17B encodes 28 OmpR family RRs, of which seven are annotated as regulating phosphatase expression, and 11 are annotated as being VanR, which regulates vancamycin resistance .
Due to intrinsic problems in defining the physiological function of regulatory proteins by sequence homology, functional annotation by sequence similarity has now largely been superseded by categorisation on the basis of domain architecture [12–14]. In this manner, RPs can be divided into families, and family membership then correlates with mechanism of action rather than biological function. Several on-line databases are now available which provide the results of such classification approaches as applied to whole genomes. For example, P2CS , P2TF , MiST2  and DBD .
However there are remaining problems with RP annotation. Many RPs contain multiple domains, and some domains are found in multiple categories of RP. This has led historically to the mis-annotation of many proteins. For instance, SAB1964 is an RR from Staphylococcus aureus RF122, yet it is annotated as a ‘two component sensor protein’, while YPA_3835 is a HK from Yersinia pestis Antiqua, which is annotated as an ‘ATPase-like ATP-binding protein’. Currently, 1.5% of all proteins now classified as RRs in the P2CS database , were originally described in some way as ‘sensor kinase’ proteins in the annotated genome files as retrieved from Refseq/Genbank. This problem is exacerbated by the current lack of a community-defined consensus set of categorisation criteria, or even a consensus naming system, for multi-domain RPs. However, this has been accomplished for a subset of RPs (RRs), by Galperin [13, 14].
Due to their multiplicity within genomes and their multi-domain architectures, RPs are non-trivial to identify and annotate. Currently, the annotation of regulatory genes/proteins in individual genomes and databases is often idiosyncratic, misleading or wrong, confounding between-genome comparisons, and naming conventions are also typically different between genomes/databases. There is consequently a profound need for the adoption of a consistent and harmonised categorisation and annotation system for RPs, which can be applied to any sequence dataset, whether newly derived sequences needing annotation, or previously annotated sequences which might benefit from re-annotation .
We have therefore developed P2RP (Predicted Prokaryotic Regulatory Proteins) – primarily to help increase (re-)annotation consistency of RPs in published genomes, and for experimental biologists who wish to investigate regulatory genes in their novel sequence data. P2RP accepts two types of input – DNA and protein sequences. For nucleotide queries there is an initial gene prediction step (using MED-Start) to generate a proteome, although. gbk (GenBank) files can also be inputted. Predicted and supplied proteomes are then screened for the presence of particular TF/TCS domains, and proteins categorised and annotated according to their domain architecture [15, 18]. Every user query is given an ID, which allows later retrieval of results, and results of the P2RP process can be viewed as a web server interface page, or downloaded in a variety of user-specified formats. P2RP can be accessed at http://www.p2rp.org and is free and open to all users, with no login requirement.
The identification of RP candidates in protein sequences is accomplished by domain analysis of each predicted protein, using RPSBLAST, as previously developed for P2CS and P2TF [4, 15, 18]. The pool of domains used to search for RP proteins was manually selected from the literature and extracted from within Pfam  and SMART  databases. P2RP is a set of PHP scripts (PHP: Hypertext Preprocessor, a server-side scripting language designed for Web development), designed to search the numerous combinations of RP modules and to categorize RP proteins into families based on similarity and/or domain architecture. To circumvent the prediction of false-positive RPs, a post-analysis process is implemented. For instance, enzymes erroneously classified as RP proteins and enzymes containing DNA-binding domains (for instance transposases) are discarded or categorized as ‘Other DNA-binding Proteins’ (ODPs) respectively. The secondary structure of RP proteins is computed using the PSIPRED method . The result of the protein structure prediction is presented as a summary of the number of strands and helices and their location on the protein sequence, on which the identified domains are highlighted. Finally, the cellular localization of each TCS protein is determined by the presence or absence of transmembrane (TM) segments, using the HMMTOP predictor . The server time required for execution of these procedures is less than 100 seconds for genomes of up to 4 megabases in size (Intel Xeon 6-Core 2x2.66 GHz).
P2RP results can be outputted for download in a user-specified format, as an Excel worksheet (see Additional file 1 and Additional file 2 for examples), FASTA file or as tab-delimited text. To keep each user session private, every user query is given an ID, which allows later retrieval of results, using the ‘Search’ menu. The data are stored on the server for one month. For easy navigation, in addition to the homepage, each P2RP page contains a navigation bar that allows users to execute a new RP analysis, to search for a previous job, to access a help page or to contact the authors.
Results and discussion
Annotation and classification of RPs
The classification strategy for TCS proteins was based on that implemented by the P2CS database of genomic and metagenomic TCSs [15, 18]. Possession of one of an established list of TCS domains defines proteins as being TCS components. Identified TCS proteins are then classified based on domain composition , for example the presence of HATPase and HisKA domains leads to classification of a protein as a HK. At this stage proteins are classified as RR, HK or phosphotransfer proteins (which usually shuttle phosphoryl groups between receiver domains).
RRs are then assigned to families according to the nature of any output domain present on the protein, for instance any RRs containing a HTH_8 DNA-binding domain are classified as PrrA family members, while HKs are sub-classified as classic, hybrid, unorthodox, CheA or ‘possible incomplete’ (those lacking an obvious phosphorylatable His residue).
In addition to TCS proteins, P2RP also identifies and categorises DNA-binding proteins (classing them as TFs), if they give hits to a compiled pool of DNA-binding domains from the Pfam and SMART libraries, in a scheme developed initially for P2TF . TFs are then divided into families according to their domain architecture as proposed previously [4, 17, 24, 25], and as implemented in the P2TF database . The P2RP analysis then identifies TF proteins also containing a CheY-like receiver (phosphoacceptor) domains and annotates them as RRs. SFs were divided into 3 sub-groupings; RpoN, RpoD and ECF (extra-cytoplasmic function) SFs. OCSs were defined as proteins that contain sensory ‘input’ and DNA-binding ‘output’ domains but lack histidine kinase and receiver domains characteristic of two-component systems . TFs with only a DNA-binding domain are named transcriptional regulators (TRs). OCSs, TRs, RRs and other non-SF DNA-binding proteins are then divided into 76 families depending upon which domains are present in the proteins. For instance MerR family members contain MerR DNA-binding and B12-binding domains (a protein domain which binds to cobalamin (vitamin B12)).
Validation of the categorisation process was achieved by comparing P2RP output with manually curated datasets, as has been described for P2TF and P2CS [4, 15, 18]. Since P2RP implements the same algorithms developed for P2CS and P2TF, validation statistics (including sensitivity and specificity) for TCS prediction and TF prediction can be obtained from the help pages at http://www.p2cs.org and http://www.p2tf.org respectively. Hundreds of inputs were analysed during testing (amino acid and nucleotide sequences), and beyond the authors, fourteen individuals were involved in the testing and validation (including undergraduate students of Aberystwyth University registered on module BS33120 Molecular Genetics of Microbes).
We have thus developed a high-quality automated analysis system, which builds homogenized genome annotations and increases the consistency of RP prediction in publicly available genomic sequences. As well as assisting experimental biologists, P2RP could be built into genome annotation pipelines, and could thereby generate a significant proportion of all genome annotation. For instance, in a rifamycin SV-producing actinomycete, Amycolatopsis mediterranei U32, over 14% of proteins are predicted and annotated as RPs.
Regulatory processes are fundamental to how microbes alter gene expression in response to environmental changes such as those encountered during infection. Regulation is achieved by RPs such as transcription factors and two-component system proteins. RPs are relatively difficult to annotate because of their multi-domain nature and their paralogy within most prokaryotic genomes. P2RP provides users with the opportunity to investigate and consistently annotate RPs within novel sequence data, or to re-annotate published sequences. It is hoped that this will prove a useful resource to experimental biologists, in addition to increasing consistency in the annotation of RPs in public databases - potentially being used routinely within annotation pipelines.
Availability and requirements
Project name: P2RP.
Project home page:http://www.p2rp.org.
Operating system(s): Platform independent.
License: This website is free and open to all users and there is no login requirement.
Other DNA-binding protein
We are very grateful to the following students at Aberystwyth University studying module BS33120 (Molecular Genetics of Microbes), for testing the server functionalities as part of their studies: Aimee Abbiss, Danielle Aspland, Jonathon Challinor, Jenna Bowen, Simon Davis, Anya-May Hope, Christine Lockey and Thomas Savage. We would also like to thank Wafa Achouak, Marie Bertrand, Gilles De Luca, Sylvain Fochesato, David Lalaouna and Catherine Santaella for useful comments and extensive testing of the web server. We are grateful to DSV/IBITEC-S/GIPSI team and particularly Arnaud Martel and Jean-Marc Le Failler for the hosting server installation.
- Österberg S, del Peso-Santos T, Shingler V: Regulation of alternative sigma factor use. Annu Rev Microbiol. 2011, 65: 37-55. 10.1146/annurev.micro.112408.134219.View ArticlePubMedGoogle Scholar
- Ulrich LE, Koonin EV, Zhulin IB: One-component systems dominate signal transduction in prokaryotes. Trends Microbiol. 2005, 13 (2): 52-56. 10.1016/j.tim.2004.12.006.PubMed CentralView ArticlePubMedGoogle Scholar
- Whitworth DE: Two-component regulatory systems in prokaryotes. Bacterial regulatory networks. Edited by: Filloux A. 2012, Norfolk: Caister Academic Press, 191-222.Google Scholar
- Ortet P, De Luca G, Whitworth DE, Barakat M: P2TF: a comprehensive resource for analysis of prokaryotic transcription factors. BMC Genomics. 2012, 13: 628-10.1186/1471-2164-13-628.PubMed CentralView ArticlePubMedGoogle Scholar
- Whitworth DE: Genomes and knowledge - a questionable relationship?. Trends Microbiol. 2008, 16 (11): 512-519. 10.1016/j.tim.2008.08.001.View ArticlePubMedGoogle Scholar
- Walderhaug MO, Polarek JW, Voelkner P, Daniel JM, Hesse JE, Altendorf K, Epstein W: KdpD and KdpE, proteins that control expression of the kdpABC operon, are members of the two-component sensor-effector class of regulators. J Bacteriol. 1992, 174 (7): 2152-2159.PubMed CentralPubMedGoogle Scholar
- Munson GP, Lam DL, Outten FW, O'Halloran TV: Identification of a copper-responsive two-component system on the chromosome of Escherichia coli K-12. J Bacteriol. 2000, 182 (20): 5864-5871. 10.1128/JB.182.20.5864-5871.2000.PubMed CentralView ArticlePubMedGoogle Scholar
- Simon G, Méjean V, Jourlin C, Chippaux M, Pascal MC: The torR gene of Escherichia coli encodes a response regulator protein involved in the expression of the trimethylamine N-oxide reductase genes. J Bacteriol. 1994, 176 (18): 5601-5606.PubMed CentralPubMedGoogle Scholar
- Yagil E, Bracha M, Lifshitz Y: The regulatory nature of the phoB gene for alkaline phosphatase synthesis in Escherichia coli. Mol Gen Genet. 1975, 137 (1): 11-16.PubMedGoogle Scholar
- Wurtzel ET, Chou MY, Inouye M: Osmoregulation of gene expression. I. DNA sequence of the ompR gene of the ompB operon of Escherichia coli and characterization of its gene product. J Biol Chem. 1982, 257 (22): 13685-13691.PubMedGoogle Scholar
- Arthur M, Molinas C, Courvalin P: The VanS-VanR two-component regulatory system controls synthesis of depsipeptide peptidoglycan precursors in Enterococcus faecium BM4147. J Bacteriol. 1992, 174 (8): 2582-2591.PubMed CentralPubMedGoogle Scholar
- Galperin MY: A census of membrane-bound and intracellular signal transduction proteins in bacteria: bacterial IQ, extroverts and introverts. BMC Microbiol. 2005, 5: 35-10.1186/1471-2180-5-35.PubMed CentralView ArticlePubMedGoogle Scholar
- Galperin MY: Structural classification of bacterial response regulators: diversity of output domains and domain combinations. J Bacteriol. 2006, 188 (12): 4169-4182. 10.1128/JB.01887-05.PubMed CentralView ArticlePubMedGoogle Scholar
- Galperin MY: Diversity of structure and function of response regulator output domains. Curr Opin Microbiol. 2010, 13 (2): 150-159. 10.1016/j.mib.2010.01.005.PubMed CentralView ArticlePubMedGoogle Scholar
- Barakat M, Ortet P, Whitworth DE: P2CS: a database of prokaryotic two-component systems. Nucleic Acids Res. 2011, 39: D771-776. 10.1093/nar/gkq1023.PubMed CentralView ArticlePubMedGoogle Scholar
- Ulrich LE, Zhulin IB: The MiST2 database: a comprehensive genomics resource on microbial signal transduction. Nucleic Acids Res. 2010, 38: D401-407. 10.1093/nar/gkp940.PubMed CentralView ArticlePubMedGoogle Scholar
- Wilson D, Charoensawan V, Kummerfeld SK, Teichmann SA: DBD--taxonomically broad transcription factor predictions: new content and functionality. Nucleic Acids Res. 2008, 36: D88-92. 10.1093/nar/gkn386.PubMed CentralView ArticlePubMedGoogle Scholar
- Barakat M, Ortet P, Jourlin-Castelli C, Ansaldi M, Méjean V, Whitworth DE: P2CS: a two-component system resource for prokaryotic signal transduction research. BMC Genomics. 2009, 10: 315-10.1186/1471-2164-10-315.PubMed CentralView ArticlePubMedGoogle Scholar
- Zhu HQ, Hu GQ, Ouyang ZQ, Wang J, She ZS: Accuracy improvement for identifying translation initiation sites in microbial genomes. Bioinformatics. 2004, 20 (18): 3308-3317. 10.1093/bioinformatics/bth390.View ArticlePubMedGoogle Scholar
- Finn RD, Tate J, Mistry J, Coggill PC, Sammut SJ, Hotz HR, Ceric G, Forslund K, Eddy SR, Sonnhammer EL: The Pfam protein families database. Nucleic Acids Res. 2008, 36: D281-288. 10.1093/nar/gkn226.PubMed CentralView ArticlePubMedGoogle Scholar
- Letunic I, Doerks T, Bork P: SMART 7: recent updates to the protein domain annotation resource. Nucleic Acids Res. 2012, 40: D302-305. 10.1093/nar/gkr931.PubMed CentralView ArticlePubMedGoogle Scholar
- McGuffin LJ, Bryson K, Jones DT: The PSIPRED protein structure prediction server. Bioinformatics. 2000, 16 (4): 404-405. 10.1093/bioinformatics/16.4.404.View ArticlePubMedGoogle Scholar
- Tusnády GE, Simon I: The HMMTOP transmembrane topology prediction server. Bioinformatics. 2001, 17 (9): 849-850. 10.1093/bioinformatics/17.9.849.View ArticlePubMedGoogle Scholar
- Minezaki Y, Homma K, Nishikawa K: Genome-wide survey of transcription factors in prokaryotes reveals many bacteria-specific families not found in archaea. DNA Res. 2005, 12 (5): 269-280.View ArticlePubMedGoogle Scholar
- Sierro N, Makita Y, de Hoon M, Nakai K: DBTBS: a database of transcriptional regulation in Bacillus subtilis containing upstream intergenic conservation information. Nucleic Acids Res. 2008, 36: D93-96. 10.1093/nar/gkn421.PubMed CentralView ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.