Email updates

Keep up to date with the latest news and content from BMC Bioinformatics and BioMed Central.

This article is part of the supplement: Proceedings of the Fourth Annual MCBIOS Conference. Computational Frontiers in Biomedicine

Open Access Proceedings

Structure and function predictions of the Msa protein in Staphylococcus aureus

Vijayaraj Nagarajan and Mohamed O Elasri*

Author Affiliations

Department of Biological Sciences, The University of Southern Mississippi, Hattiesburg, MS, 39406, USA

For all author emails, please log on.

BMC Bioinformatics 2007, 8(Suppl 7):S5  doi:10.1186/1471-2105-8-S7-S5

The electronic version of this article is the complete one and can be found online at: http://www.biomedcentral.com/1471-2105/8/S7/S5


Published:1 November 2007

© 2007 Nagarajan and Elasri; licensee 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.

Abstract

Background

Staphylococcus aureus is a human pathogen that causes a wide variety of life-threatening infections using a large number of virulence factors. One of the major global regulators used by S. aureus is the staphylococcal accessory regulator (sarA). We have identified and characterized a new gene (modulator of sarA: msa) that modulates the expression of sarA. Genetic and functional analysis shows that msa has a global effect on gene expression in S. aureus. However, the mechanism of Msa function is still unknown. Function predictions of Msa are complicated by the fact that it does not have a homologous partner in any other organism. This work aims at predicting the structure and function of the Msa protein.

Results

Preliminary sequence analysis showed that Msa is a putative membrane protein. It would therefore be very difficult to purify and crystallize Msa in order to acquire structure information about this protein. We have used several computational tools to predict the physico-chemical properties, secondary structural features, topology, 3D tertiary structure, binding sites, motifs/patterns/domains and cellular location. We have built a consensus that is derived from analysis using different algorithms to predict several structural features. We confirm that Msa is a putative membrane protein with three transmembrane regions. We also predict that Msa has phosphorylation sites and binding sites suggesting functions in signal transduction.

Conclusion

Based on our predictions we hypothesise that Msa is a novel signal transducer that might be involved in the interaction of the S. aureus with its environment.

Background

Introduction

Staphylococcus aureus is an important human pathogen that causes several diseases ranging from superficial skin infections to life-threatening diseases such as osteomyelitis and endocarditis. S. aureus is capable of infecting a wide range of tissues in humans because of the large number of virulence factors and the complex regulatory networks that control them [1]. In addition, S. aureus is increasingly resistant to multiple antibiotics thus becoming a growing threat to public health. There is an urgent need to understand the complex regulatory networks used by S. aureus to cause disease. Regulatory networks are attractive therapeutic targets for future treatment of antibiotic resistant infections.

Modulator of sarA (msa)

One of the important global regulators of virulence in S. aureus is the Staphylococcal accessory regulator (sarA) [2]. sarA regulates over 100 genes in S. aureus several of which are associated with virulence [3]. sarA plays an important role in disease [4]. sarA itself is regulated by several loci that modulate its function. We recently identified a novel gene, msa, that modulate the function of sarA [5]. We showed that msa is essential for full expression of sarA and that mutation of msa affected the expression of several virulence factors in both sarA-dependent and sarA-independent manners [5]. Microarray analyses of the msa mutant show that Msa has a global effect on genes in S. aureus (unpublished data). These studies indicate that msa is an important locus in S. aureus and that the characterization of the Msa protein would be very useful in understanding staphylococcal regulatory networks.

Computational tools

Several bioinformatics tools have been developed to predict the structure and functional properties of bio molecules. These tools use a wide variety of algorithms to predict the properties of proteins at different levels [6,7]. The accuracy of these bioinformatics tools has been improving; however, each tool has its own advantages and disadvantages. A particular algorithm has its own characteristic specificity, sensitivity, robustness, computational cost, etc. These characteristics can be tested against benchmarks of known datasets (e. g., Critical Assessment of Techniques for Protein Structure Prediction – CASP). In order to make the most accurate predictions, several methods should be used to build a consensus.

The aim of this work is to predict the structure and functional properties of the Msa protein of S. aureus to the highest possible accuracy. Our prediction results show that the Msa is a putative integral membrane protein with three probable transmembrane regions. We also predict that the Msa contains phosphorylation sites in the loop regions (both inside and outside the membrane). The 3-D structure analysis of the Msa also predicts the presence of putative binding sites. Thus, based on this computational analysis, and previous experimental data [5] we hypothesise that Msa might play a role in signal transduction. The fact that Msa has no known homolog means that it would be a novel signal transducer.

Results and discussion

Primary sequence analysis

The conceptually translated Msa protein is made of 133 amino acids with a predicted molecular weight of 15.6571 kDa and an isoelectric point (pI) of 6.71. The GRAVY index value 1.021 shows that Msa is probably an insoluble protein. The Codon adaptation index (CAI) value predicts the Msa as a highly expressed protein. This is consistent with experimental results described previously by our group [5].

Homology and similarity

The Msa is highly conserved among the different strains of S. aureus (RF122, MRSA252, MSSA476, MW2, COL, Mu50, N315, and NCTC 8325). Even though there were several variations in the nucleotide sequences, we observed good conservation at the amino acid level. Multiple sequence alignment and phylogenetic analysis of both nucleotide sequences (SAUSA300_1294, SACOL1436, SAOUHSC_01402, SAV1401, msa, SAS1342, MW1289, SAR1413, and SAB1257c) and protein sequences (YP_493991, YP_186288, YP_499929, NP_371925, NP_374514, YP_043463, NP_646106) from different strains show that they are identical. The only two exceptions were strains RF122 and MRSA252 which showed slight variations in the Msa sequences. In RF122, the protein sequence (YP_416734) was 97% similar to the Msa sequence from N315 while in MRSA252, the protein sequence (YP_040815) was 98% similar to the Msa from N315. The phylogeny of the Msa protein closely resembled that of the phylogeny of these organisms as determined by Multi Locus Sequence Typing (MLST) [8]. The position and effect of mutations in the Msa protein sequence of the strains MRSA252 and RF122 are discussed in the "3-D structure prediction and analysis" section.

Our similarity search results against several sequence and structure databases, using different BLAST programs, showed that there were no significant closely related homologs for the Msa protein, except for one in S. epidermidis. Even though there were no significant (based on E-value and score) homolog for Msa, BLAST also listed several membrane proteins with remote similarities (alignment Score of 32–35 and E values scores from 0.91–10) only to the first few amino acids of the Msa protein (that corresponds to the predicted signal peptide region).

Localization predictions

All the tools used to predict the cellular location of the protein indicated that Msa is a putative membrane protein. This prompted us to examine the sequence for presence of signal peptide and potential cleavage sites in the Msa protein sequence. Seven out of eight signal peptide prediction tools indicated the presence of a potential signal peptide in the Msa protein (Table 1). The majority of the programs also predicted an N-terminal cleavage site between the amino acid 19 and 20.

Table 1. Signal peptide and cleavage position prediction for the Msa protein

Topology predictions

We performed topology analysis on the Msa sequence using several prediction programs that yielded widely discrepant results (Table 2). Even though most programs failed to recognize the signal peptide, a consensus topology emerged (Table 3). The predicted topology of the Msa is IN-OUT with three putative transmembrane segments (from amino acid positions 27–47, 54–75, 108–125). The N-terminal is predicted as present in the cytoplasmic side of the membrane while the C-terminal is predicted as outside the membrane. Our consensus topology also passed the positive-inside rule and charge bias test [9], with a charge bias of +1 towards the inside of the membrane.

Table 2. Topology predictions for the Msa protein

Table 3. Consensus topology for the Msa protein including the N-terminal signal peptide prediction

Secondary structure prediction results indicated the presence of four distinct helical regions (Figure 1). One helical region corresponds to the cytoplasmic helix while the other three correspond to the integral membrane helices. These results are consistent with the predicted topology.

thumbnailFigure 1. Consensus secondary structure predictions for the Msa protein. Three transmembrane segments (TMS) and a cytoplasmic helix are predicted.

Domains/patterns/motifs

We searched for the presence of domains, patterns and motifs in the Msa protein sequence, to gain insight into its functions and structure. The SMART results showed the presence of all the structural domains that we earlier identified using topology prediction programs and signal peptide prediction programs, viz. an N-terminal signal peptide and three transmembrane regions. In addition, SMART also predicted the presence of a PreATP-grasp domain (d1gsa_1) from the SCOP database. Even though this result had an E-value of 1.5, it was interesting because the predicted domain is a putative binding domain and falls in the predicted cytoplasmic region of Msa (residues 85–116). Our pattern search in the Msa protein sequence, using different programs against the PROSITE database, gave similar results (except for PPSearch, which did not predict the Tyrosine kinase site at position 48), showing the presence of three putative phosphorylation sites (Table 4). All of the predicted sites were found in the exposed regions of the Msa. Analysis of the location of these putative phosphorylation sites showed that two of the putative phosphorylation sites are outside the membrane while one of them is predicted in the cytoplasmic region. We also observed that these putative phosphorylation sites are highly conserved among different strains of S. aureus. This suggests that Msa might be phosphorylated by kinases in the cytoplasm as well as kinases on the outside of the membrane (e.g. from the host cells). These predictions further suggest that Msa might function as a signal transducer and provides important targets for mutagenesis experiments to test this hypothesis.

Table 4. Prosite patterns predicted in the Msa protein

Membrane bound receptors are important components of signal transduction in all living systems. The major class of receptors in eukaryotes contain seven transmembrane segments (7 TM). Prokaryotes use 7 TM class receptors also, however, a recent study showed that prokaryotes carry novel receptor classes that have transmembrane segments ranging from one to eight [10,11]. The Msa protein sequence did not have significant homology with any of the known receptors and experimental studies are underway to evaluate its function as a signal transducer.

3-D structure prediction and analysis

Homology based tertiary structure prediction for the Msa protein failed, because of the lack of homologous structures. We used fold recognition based structure prediction server Phyre to model the tertiary structure of the Msa protein. Visualization and analysis of the predicted structure using Swiss-PDB Viewer (SPDBV) showed that the predicted structure correlated with the other predicted structural features of Msa in terms of the number and positions of the transmembrane helices (Figure 2). We refined the predicted structure by fixing side chains, fixing problematic loops, removal of amino acid clashes (bumps) and energy minimization. The refinements did not yield any drastic change in the initial predicted structure. This was confirmed by visually inspecting the structure and verifying the backbone structure using Ramachandran plot (Figure 3) and computing the total energy difference between the initial model and the refined model.

thumbnailFigure 2. Predicted tertiary structure of the Msa protein showing the three transmembrane helices. Arrow indicates the predicted cleavage site for the putative signal peptide. N, N-terminus; C, C-terminus

thumbnailFigure 3. Ramachandran plot for the predicted tertiary structure of the Msa protein pre (A) and post (B) refinement.

We analysed the predicted tertiary structure for clefts and binding sites using ProFunc server and found putative binding sites in the cytoplasmic region between the second and the third transmembrane helices (Figure 4A). We also used PINUP to predict putative interface residues in the similar region (Figure 4B). Another binding site prediction server Q-SiteFinder also predicted similar binding site and binding site residues (Figure 4C).

thumbnailFigure 4. Binding site predictions for the Msa protein. (A) ProFunc predicted binding site (red); (B) PINUP predicted binding site (interface in green); (C) Q-SiteFinder predicted binding site and binding residues (pink)

ProFunc also predicted a "nest" near the putative phosphorylation site (residues 47–50) which was predicted outside the membrane [12]. The Msa has all the conserved residues that make up the predicted "nest". The predicted "nest" in Msa shows features of an anion-binding site. Such "nests" are characteristic functional motifs, which are found in ATP- or GTP binding proteins.

Multiple sequence alignment of the Msa protein sequence from 11 different strains of S. aureus revealed 12 mutations in strain RF122 and seven mutations in strain MRSA252 relative to consensus. Mutations at amino acid positions 111, 131 and 133 were found in both MRSA252 and RF122 strains. None of these mutations were found in the predicted phosphorylation sites, predicted signal peptide sites or in the predicted anion-binding "nest". But many of the mutations were found both in the integral membrane segments as well as in the other parts of the loop regions. Only one out of the 12 mutations had the replacement (functionally different amino acid), while others were substitutions (functionally similar amino acids), in the strain RF122. In the strain MRSA152, two out of seven mutations were replacements, while others were substitutions. MRSA strain had three mutations in the predicted pre-ATP grasp domain, out of which one had an amino acid replacement. RF122 strain had only one amino acid substitution in the pre-ATP grasp domain. This indicates that the predicted functional sites are constrained from mutation.

Conclusion

We predict that Msa is a membrane protein with a cleavable N-terminal signal peptide sequence, followed by three integral transmembrane regions. The Msa is also predicted to have an IN-OUT topology with at least two putative phosphorylation sites, one outside the membrane and one in the cytoplasmic region. A putative binding site is also predicted in the cytoplasmic region of Msa. Based on these predictions we put forward a model for the Msa protein (Figure 5). This model also prompted us to hypothesise that Msa might function as a novel signal transducer between the environment and the cytoplasm. This model will be used to design and execute experiments to confirm the functions and topology of Msa and further our understanding of its role in the pathogenesis of S. aureus.

thumbnailFigure 5. Predicted model for the Msa protein showing structural and functional features.

Methods

For a complete list of online tools used, see additional file 1.

Additional file 1. URL's of the tools and applications used in this study. URL's of the tools and applications used in this study. This file contains the list of URL's of all the online tools used to predict the structure and functional aspects of Msa protein.

Format: TXT Size: 3KB Download fileOpen Data

Primary sequence analysis

We used the protein sequence (Accession ID: NP_374514) obtained by conceptual translation of the msa open reading frame from the S. aureus N315 genome (NCBI database). The primary sequence analysis was performed using ProtParam, ProtScale [13] and SAPS [14]. ProtScale was used to predict the Msa profile based on several amino acid scales. ProtParam computes properties like molecular weight, theoretical pI, instability index and grand average of hydropathicity (GRAVY). SAPS predicts significant features of protein sequences like charge-clusters, hydrophobic regions, compositional domains etc.

Similarity searching

Similarity searching was done using different programs at NCBI like BLASTP, PSI-BLAST [15] and CDART [16] against several different databases like NR, SWISSPROT and PDB. Multiple sequence alignment and phylogenetic analysis were done using Accelrys Gene v2.5 (Accelrys Inc., San Diego, CA).

Sub-cellular localization

The Sub-cellular localization and the functional categorization of Msa were predicted using ProtFun 2.2 [17], PSORT [18], ProtCompB V-3 [19], PRED-CLASS [20] and SVMProt [21].

ProtFun uses ab initio methods to predict the cellular role category. PSORT uses a rule-based method to predict protein localization sites. ProtCompB combines several methods such as linear discriminant function-based predictions, direct comparison with homologous proteins of known localization, prediction of functional peptide sequences etc., to identify the sub-cellular localization of proteins. PRED-CLASS uses cascading neural networks to classify proteins in to different classes like membrane, globular, fibrous and mixed. SVMProt uses a support vector machine based approach to functionally classify protein sequences.

Signal peptide prediction

Signal peptide prediction was done using SignalP [22], PrediSi [23], sigcleave [24], PSORT [18], Phobius [25], SIG-Pred [26], SOSUIsignal [27] and iPSORT [28].

SignalP 3.0 uses artificial neural networks and hidden Markov models to predict signal peptides and their cleavage sites. PrediSi predicts signal peptide sequences and their cleavage sites based on a position weight matrix that also takes into consideration the amino acid bias present in the proteins. Sigcleave is one of the early tools to predict the signal cleavage sites based on weight matrices. Sigcleave is distributed as part of the EMBOSS package. Phobius is a combined transmembrane protein topology and signal peptide predictor that uses a well trained hidden Markov model. SIG-Pred predicts signal peptides and their cleavage position based on weight matrices. SOSUIsignal uses a high performance system to predict signal peptides, using a three module software system that recognises the three-domain structure of signal peptides. iPSORT predicts the signal peptides based on a rule based system.

Topology prediction

The Topology of Msa protein was predicted using TopPred [29], TMpred [30], PHDhtm [31], TMHMM [32], SPLIT [33], HMMTOP [34], MEMSAT [35], DAS [36] and TSEG [37]. We also computed the charge bias of the generated models, based on the positive-inside rule [9].

TopPred II predicts the topology of a protein based on its hydrophobicity profile and positive-inside rule. TMpred algorithm is based on the statistical analysis of TMbase, a database of naturally occurring transmembrane proteins, using a combination of several weight-matrices for scoring. PHDhtm uses a neural network based approach with the evolutionary information to predict the locations of the transmembrane helices. TMHMM predicts transmembrane regions based on the hidden Markov model. SPLIT 4.0 predicts location of transmembrane helices by performing an automatic selection of optimal amino acid attribute and corresponding preference functions. HMMTOP 2.0 prediction is based on the hypothesis that the difference in the amino acid distributions in various structural parts determines the localization of the transmembrane segments. MEMSAT applies a novel dynamic programming algorithm to recognize membrane topology models by expectation maximization. DAS uses dense alignment surface method to predict transmembrane regions. TSEG uses a discriminant function to predict the transmembrane segments.

Secondary structure prediction

We used the NPS (Network Protein Sequence Analysis) consensus secondary structure server [38]. This server runs the input sequence against several different secondary structure prediction tools and generates a consensus secondary structure out of them.

Domains/patterns/motifs prediction

SMART (Simple Modular Architecture Research Tool) [39] was used to identify the presence of any domains in the Msa protein. We used different pattern searching applications (PPSearch [40], PSITE [41] and ScanProsite [42]), that use PROSITE [43] database, to predict functionally relevant patterns in Msa protein.

3-D Structure prediction and analysis

Initial attempts to predict the tertiary structure of Msa were done using different approaches like homology modelling, threading and ab initio. Automated homology modelling servers Swiss-Model [44] and ModWeb [45] were used for homology modelling. Predictions described in this study were done using fold recognition tools 123D+ [46], GenThreader [47], a new version of 3-D PSSM (Phyre) [48].

The quality of the predicted structure was examined using an online version of the WHATIF [49] program. Structure refinement was done using both WHATIF and Swiss-PDB Viewer [50]. Structure visualization was done using Swiss-PDB Viewer.

The 3-D structure of the Msa protein was analysed for clefts and binding surfaces using ProFunc [51], Q-SiteFinder [52], PINUP [53] and SuMo [54].

Meta servers

We also used meta servers like SCRATCH [55], ProSAL [56] and MetaPP [57] for predicting structure and functional properties of the Msa protein.

Competing interests

The authors declare that they have no competing interests.

Authors' contributions

VN did the analysis and drafted the manuscript. MOE directed the whole research and critically revised the manuscript.

Acknowledgements

This study was supported by a grant (1R15AI062727-01A1) from the National Institute of Allergy and Infectious Diseases (NIAID) to MOE and by The Mississippi Functional Genomics Network (NIH/NCRR P20 RR016476). We extend our sincere thanks to the anonymous reviewers of this manuscript for their valuable suggestions.

This article has been published as part of BMC Bioinformatics Volume 8 Supplement 7, 2007: Proceedings of the Fourth Annual MCBIOS Conference. Computational Frontiers in Biomedicine. The full contents of the supplement are available online at http://www.biomedcentral.com/1471-2105/8?issue=S7.

References

  1. Novick RP: Autoinduction and signal transduction in the regulation of staphylococcal virulence.

    Mol Microbiol 2003, 48(6):1429-1449. PubMed Abstract | Publisher Full Text OpenURL

  2. Cheung AL, Projan SJ: Cloning and sequencing of sarA of Staphylococcus aureus, a gene required for the expression of agr.

    J Bacteriol 1994, 176(13):4168-4172. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  3. Dunman PM, Murphy E, Haney S, Palacios D, Tucker-Kellogg G, Wu S, Brown EL, Zagursky RJ, Shlaes D, Projan SJ: Transcription profiling-based identification of Staphylococcus aureus genes regulated by the agr and/or sarA loci.

    J Bacteriol 2001, 183(24):7341-7353. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  4. Blevins JS, Elasri MO, Allmendinger SD, Beenken KE, Skinner RA, Thomas JR, Smeltzer MS: Role of sarA in the pathogenesis of Staphylococcus aureus musculoskeletal infection.

    Infect Immun 2003, 71(1):516-523. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  5. Sambanthamoorthy K, Smeltzer MS, Elasri MO: Identification and characterization of msa (SA1233), a gene involved in expression of SarA and several virulence factors in Staphylococcus aureus.

    Microbiology 2006, 152(Pt 9):2559-2572. PubMed Abstract | Publisher Full Text OpenURL

  6. Rost B: Review: protein secondary structure prediction continues to rise.

    J Struct Biol 2001, 134(2–3):204-218. PubMed Abstract | Publisher Full Text OpenURL

  7. Pandey G, Kumar V, Steinbach M: Computational Approaches for Protein Function Prediction: A Survey. Twin Cities: Department of Computer Science and Engineering, University of Minnesota; 2006.

  8. Holden MT, Feil EJ, Lindsay JA, Peacock SJ, Day NP, Enright MC, Foster TJ, Moore CE, Hurst L, Atkin R, et al.: Complete genomes of two clinical Staphylococcus aureus strains: evidence for the rapid evolution of virulence and drug resistance.

    Proc Natl Acad Sci USA 2004, 101(26):9786-9791. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  9. von Heijne G: Membrane protein structure prediction. Hydrophobicity analysis and the positive-inside rule.

    J Mol Biol 1992, 225(2):487-494. PubMed Abstract | Publisher Full Text OpenURL

  10. Galperin MY: A census of membrane-bound and intracellular signal transduction proteins in bacteria: bacterial IQ, extroverts and introverts.

    BMC Microbiol 2005, 5:35. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  11. D'Souza M, Glass EM, Syed MH, Zhang Y, Rodriguez A, Maltsev N, Galperin MY: Sentra: a database of signal transduction proteins for comparative genome analysis.

    Nucleic Acids Res 2007, (35 Database):D271-273. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  12. Watson JD, Milner-White EJ: A novel main-chain anion-binding site in proteins: the nest. A particular combination of phi, psi values in successive residues gives rise to anion-binding sites that occur commonly and are found often at functionally important regions.

    J Mol Biol 2002, 315(2):171-182. PubMed Abstract | Publisher Full Text OpenURL

  13. Gasteiger E, Hoogland C, Gattiker A, Duvaud S, Wilkins MR, Appel RD, Bairoch A: Protein Identification and Analysis Tools on the ExPASy Server. In The Proteomics Protocols Handbook. Edited by Walker JM. Humana Press; 2005:571-607. OpenURL

  14. Brendel V, Bucher P, Nourbakhsh IR, Blaisdell BE, Karlin S: Methods and algorithms for statistical analysis of protein sequences.

    Proc Natl Acad Sci USA 1992, 89(6):2002-2006. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  15. Altschul SF, Madden TL, Schaffer AA, Zhang J, Zhang Z, Miller W, Lipman DJ: Gapped BLAST and PSI-BLAST: a new generation of protein database search programs.

    Nucleic Acids Res 1997, 25(17):3389-3402. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  16. Geer LY, Domrachev M, Lipman DJ, Bryant SH: CDART: protein homology by domain architecture.

    Genome Res 2002, 12(10):1619-1623. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  17. Jensen LJ, Gupta R, Staerfeldt HH, Brunak S: Prediction of human protein function according to Gene Ontology categories.

    Bioinformatics 2003, 19(5):635-642. PubMed Abstract | Publisher Full Text OpenURL

  18. Nakai K, Horton P: PSORT: a program for detecting sorting signals in proteins and predicting their subcellular localization.

    Trends Biochem Sci 1999, 24(1):34-36. PubMed Abstract | Publisher Full Text OpenURL

  19. Softberry [http://www.softberry.com] webcite

  20. Pasquier C, Promponas VJ, Hamodrakas SJ: PRED-CLASS: cascading neural networks for generalized protein classification and genome-wide applications.

    Proteins 2001, 44(3):361-369. PubMed Abstract | Publisher Full Text OpenURL

  21. Cai CZ, Han LY, Ji ZL, Chen X, Chen YZ: SVM-Prot: Web-based support vector machine software for functional classification of a protein from its primary sequence.

    Nucleic Acids Res 2003, 31(13):3692-3697. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  22. Bendtsen JD, Nielsen H, von Heijne G, Brunak S: Improved prediction of signal peptides: SignalP 3.0.

    J Mol Biol 2004, 340(4):783-795. PubMed Abstract | Publisher Full Text OpenURL

  23. Hiller K, Grote A, Scheer M, Munch R, Jahn D: PrediSi: prediction of signal peptides and their cleavage positions.

    Nucleic Acids Res 2004, (32 Web Server):W375-379. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  24. von Heijne G: A new method for predicting signal sequence cleavage sites.

    Nucleic Acids Res 1986, 14(11):4683-4690. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  25. Kall L, Krogh A, Sonnhammer EL: A combined transmembrane topology and signal peptide prediction method.

    J Mol Biol 2004, 338(5):1027-1036. PubMed Abstract | Publisher Full Text OpenURL

  26. SIG-Pred: Signal Peptide Prediction [http://www.bioinformatics.leeds.ac.uk/prot_analysis/Signal.html] webcite

  27. Gomi M, Sonoyama M, Mitaku S: High performance system for signal peptide prediction: SOSUIsignal.

    Chem-Bio Informatics Journal 2004, 4(4):142-147. Publisher Full Text OpenURL

  28. Bannai H, Tamada Y, Maruyama O, Nakai K, Miyano S: Extensive feature detection of N-terminal protein sorting signals.

    Bioinformatics 2002, 18(2):298-305. PubMed Abstract | Publisher Full Text OpenURL

  29. Claros MG, von Heijne G: TopPred II: an improved software for membrane protein structure predictions.

    Comput Appl Biosci 1994, 10(6):685-686. PubMed Abstract OpenURL

  30. Hofmann K, Stoffel W: TMBASE – A database of membrane spanning protein segments.

    Biol Chem Hoppe-Seyler 1993, 374:166. OpenURL

  31. Rost B, Casadio R, Fariselli P, Sander C: Transmembrane helices predicted at 95% accuracy.

    Protein Sci 1995, 4(3):521-533. PubMed Abstract | Publisher Full Text OpenURL

  32. Krogh A, Larsson B, von Heijne G, Sonnhammer EL: Predicting transmembrane protein topology with a hidden Markov model: application to complete genomes.

    J Mol Biol 2001, 305(3):567-580. PubMed Abstract | Publisher Full Text OpenURL

  33. Juretic D, Zoranic L, Zucic D: Basic charge clusters and predictions of membrane protein topology.

    J Chem Inf Comput Sci 2002, 42(3):620-632. PubMed Abstract | Publisher Full Text OpenURL

  34. Tusnady GE, Simon I: The HMMTOP transmembrane topology prediction server.

    Bioinformatics 2001, 17(9):849-850. PubMed Abstract | Publisher Full Text OpenURL

  35. Jones DT, Taylor WR, Thornton JM: A model recognition approach to the prediction of all-helical membrane protein structure and topology.

    Biochemistry 1994, 33(10):3038-3049. PubMed Abstract | Publisher Full Text OpenURL

  36. Cserzo M, Wallin E, Simon I, von Heijne G, Elofsson A: Prediction of transmembrane alpha-helices in prokaryotic membrane proteins: the dense alignment surface method.

    Protein Eng 1997, 10(6):673-676. PubMed Abstract | Publisher Full Text OpenURL

  37. Kihara D, Shimizu T, Kanehisa M: Prediction of membrane proteins based on classification of transmembrane segments.

    Protein Eng 1998, 11(11):961-970. PubMed Abstract | Publisher Full Text OpenURL

  38. Deleage G, Blanchet C, Geourjon C: Protein structure prediction. Implications for the biologist.

    Biochimie 1997, 79(11):681-686. PubMed Abstract | Publisher Full Text OpenURL

  39. Letunic I, Copley RR, Schmidt S, Ciccarelli FD, Doerks T, Schultz J, Ponting CP, Bork P: SMART 4.0: towards genomic data integration.

    Nucleic Acids Res 2004, (32 Database):D142-144. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  40. PPSearch [http://www.ebi.ac.uk/ppsearch/] webcite

  41. Solovyev VV, Kolchanov NA: Search for functional sites using consensus. In Computer analysis of Genetic macromolecules. Edited by Kolchanov NA, Lim HA. World Scientific; 1994:16-21. OpenURL

  42. de Castro E, Sigrist CJ, Gattiker A, Bulliard V, Langendijk-Genevaux PS, Gasteiger E, Bairoch A, Hulo N: ScanProsite: detection of PROSITE signature matches and ProRule-associated functional and structural residues in proteins.

    Nucleic Acids Res 2006, (34 Web Server):W362-365. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  43. Hulo N, Bairoch A, Bulliard V, Cerutti L, De Castro E, Langendijk-Genevaux PS, Pagni M, Sigrist CJ: The PROSITE database.

    Nucleic Acids Res 2006, (34 Database):D227-230. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  44. Schwede T, Kopp J, Guex N, Peitsch MC: SWISS-MODEL: An automated protein homology-modeling server.

    Nucleic Acids Res 2003, 31(13):3381-3385. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  45. Sali A, Blundell TL: Comparative protein modelling by satisfaction of spatial restraints.

    J Mol Biol 1993, 234(3):779-815. PubMed Abstract | Publisher Full Text OpenURL

  46. Alexandrov NN, Nussinov R, Zimmer RM: Fast protein fold recognition via sequence to structure alignment and contact capacity potentials.

    Pac Symp Biocomput 1996, 53-72. PubMed Abstract OpenURL

  47. Bryson K, McGuffin LJ, Marsden RL, Ward JJ, Sodhi JS, Jones DT: Protein structure prediction servers at University College London.

    Nucleic Acids Res 2005, (33 Web Server):W36-38. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  48. Kelley LA, MacCallum RM, Sternberg MJ: Enhanced genome annotation using structural profiles in the program 3D-PSSM.

    J Mol Biol 2000, 299(2):499-520. PubMed Abstract | Publisher Full Text OpenURL

  49. Vriend G: WHAT IF: a molecular modeling and drug design program.

    J Mol Graph 1990, 8(1):52-56.

    29.

    PubMed Abstract | Publisher Full Text OpenURL

  50. Guex N, Peitsch MC: SWISS-MODEL and the Swiss-PdbViewer: an environment for comparative protein modeling.

    Electrophoresis 1997, 18(15):2714-2723. PubMed Abstract | Publisher Full Text OpenURL

  51. Laskowski RA, Watson JD, Thornton JM: ProFunc: a server for predicting protein function from 3D structure.

    Nucleic Acids Res 2005, (33 Web Server):W89-93. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  52. Laurie AT, Jackson RM: Q-SiteFinder: an energy-based method for the prediction of protein-ligand binding sites.

    Bioinformatics 2005, 21(9):1908-1916. PubMed Abstract | Publisher Full Text OpenURL

  53. Liang S, Zhang C, Liu S, Zhou Y: Protein binding site prediction using an empirical scoring function.

    Nucleic Acids Res 2006, 34(13):3698-3707. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  54. Jambon M, Imberty A, Deleage G, Geourjon C: A new bioinformatic approach to detect common 3D sites in protein structures.

    Proteins 2003, 52(2):137-145. PubMed Abstract | Publisher Full Text OpenURL

  55. SCRATCH [http://www.ics.uci.edu/%7Ebaldig/scratch/] webcite

  56. ProSAL [http://xray.bmc.uu.se/sbnet/prosal.html] webcite

  57. MetaPP [http://www.predictprotein.org/newwebsite/meta/submit3.php] webcite