Skip to main content

Advertisement

WebNetCoffee: a web-based application to identify functionally conserved proteins from Multiple PPI networks

Article metrics

Abstract

Background

The discovery of functionally conserved proteins is a tough and important task in system biology. Global network alignment provides a systematic framework to search for these proteins from multiple protein-protein interaction (PPI) networks. Although there exist many web servers for network alignment, no one allows to perform global multiple network alignment tasks on users’ test datasets.

Results

Here, we developed a web server WebNetcoffee based on the algorithm of NetCoffee to search for a global network alignment from multiple networks. To build a series of online test datasets, we manually collected 218,339 proteins, 4,009,541 interactions and many other associated protein annotations from several public databases. All these datasets and alignment results are available for download, which can support users to perform algorithm comparison and downstream analyses.

Conclusion

WebNetCoffee provides a versatile, interactive and user-friendly interface for easily running alignment tasks on both online datasets and users’ test datasets, managing submitted jobs and visualizing the alignment results through a web browser. Additionally, our web server also facilitates graphical visualization of induced subnetworks for a given protein and its neighborhood. To the best of our knowledge, it is the first web server that facilitates the performing of global alignment for multiple PPI networks.

Availability

http://www.nwpu-bioinformatics.com/WebNetCoffee

Background

Proteins are involved in almost all life processes. The discovering of protein function is of significance for understanding the underlying molecular mechanism in organisms and diseases. Thanks to the development of high-throughput technologies, computational approaches become a major force to uncover molecular structures and functions [1, 2]. Since then, a large amount to genomics, proteins, interactions, pathways and functional annotation data have been accumulated and freely available in public databases, such as Genbank [3], Uniprot [4], KEGG [5] and RSCB PDB [6]. So far, many computational tools [710] have been developed to understand molecular function by using genomic sequences, pathways and molecular networks. However, there is still a large gap to obtain a comprehensive knowledge of protein function for various species [11].

Global network alignment is an efficient framework to systematically identify functionally conserved proteins from different species. It aims to search for an optimal global node map for all nodes in different PPI networks. These proteins matched in one group are thought to be functionally conserved, which are also called function-oriented ortholog (FO) groups. Therefore, one can predict the function of an uncharacterized protein according to the functional annotation of another protein in its FO group.

The results of a global network alignment consist of a series of matchsets for two or more PPI networks. Each matchset of the node map represents a putative FO group. IsoRank [12] was firstly developed to identify FO groups between two given PPI networks. The alignment was intuitively guided by an idea that two proteins matched if and only if their neighbors can also be well matched. The algorithm was later improved in IsoRankN [13] to find a global alignment for multiple networks (GAMN). IsoRankN takes spectral partitioning method to find alignment clusters on the induced graphs of pairwise alignment scores. Since then, a bunch of multiple global alignment tools have been developed one after another, most notably Graemlin 2.0 [14], SMETANA [15], NETAL [16], NetCoffee [17], MAGNA [18] and MAGNA++ [19].

However, most of these alignment tools depend on some existing libraries, environment configuration or many input datasets with a particularly designed format, which is difficult to prepare. And the command line interface (CLI) of these computational tools cannot provide a graphical user interface to manage the data, visualize the processing status and annotate the alignment results using public databases. Several web servers have been developed to perform local alignments, to query pathways and patterns, including NetworkBlast [20], PINALOG [21], NetAligner [22], PathBlast [23], NetalieQ [24], etc. However, as far as we know, very few of web servers facilitate the aligning of multiple PPI networks. To make the task of GAMN easier to be done for non-expert users, here, we present a web server WebNetCoffee based on the NetCoffee algorithm, which can fast and accurately search for a global node map for multiple PPI networks. The web server is available at http://www.nwpu-bioinformatics.com/WebNetCoffee/, which enables users to upload their own test datasets or select three or more species from four well-known PPI databases IntAct [25], STRING [26], DIP [27] and BioGRID [28].

Implementation

We implemented the web server WebNetCoffee in several programming languages, which include C++, html, php, css and mysql. The graph library LEMON version 1.2.3 [29] was used in the implementation of NetCoffee. The Apache HTTP Server and MySQL provide the fundamental web server environment. The server runs on a CPU of Intel(R) Xeon(R) CPU E5-2682 v4 @ 2.50GHz.

The workflow of NetCoffee

Given a set of networks G1,G2,,Gk, k≥3, each network can be modeled as a graph Gi=(Vi,Ei), where Vi and Ei represent proteins and interactions appearing in networks. Proteins aligned in one group is a matchset, which is a subset of \(\cup _{i=1}^{k}V_{i}\). The global network alignment problem is to search for a set of mutually disjoint matchsets for two or more PPI networks. We assumed that the sequence similarity and topology similarity can imply the functional conservation of proteins in different species.

An integrated model was adopt to measure the similarity of a given pair of nodes by using both topology and sequence information in the NetCoffee algorithm [17]. It employs simulated annealing to optimize a target function, aiming to search for an optimal global one-to-one map based on similarity of both network topology and protein sequences. There are four major steps: 1) building PPI networks and a library of bipartite graphs; 2) calculating integrated weight using triplet extension; 3) collecting candidate edges with maximum matching; 4) optimization with the simulated annealing approach. NetCoffee was distinguished itself with other existing algorithms by its fast speed and biologically meaningful alignment. It can perform a GAMN job on three or more PPI networks. The alignment result consists of a lot of matchsets, each of which represents a putative functional ortholog group.

Using WebNetCoffee

The WebNetCoffee provides a simple web interface for performing GAMN tasks. The home page of WebNetCoffee briefly introduces the foundation of NetCoffee and the resources of online datasets. The help page can quickly guide a new user to perform a GAMN task and query the results through a jobid step by step. Besides, it also gives more detailed description for each panel in the result page.

Users can launch a WebNetCoffee job on both the online datasets and users’ own datasets (see in Fig. 1a,b). Each job can be assigned a user-specified job title. The default parameter of alpha is 0.5, which is used to balance the contribution of topology and sequence score in the alignment result. To launch a job on the online datasets, totally, 15 species from four databases are available for users’ options. One can choose three or more species for each GAMN task. To avoid users’ very large computational tasks, WebNetCoffee allows performing a job on 3-5 networks in BioGRID and STRING, 3-11 networks in IntAct and DIP. Each file of users’ datasets uploaded to the server is restricted to be less than 200M.

Fig. 1
figure1

A workflow of WebNetCoffee. a Networks of up to 15 species are available for users’ options, which were manually collected from four commonly used databases. b Users can perform WebNetCoffee on their own datasets. c Statistics of test datasets. d Convergence curves of simulated annealing. e Table of the functionally conserved protein groups. f Visualization of induced sub-networks

In case many tasks were submitted to the server simultaneously, we designed a job queue to manage all the jobs. Each submitted job would firstly go into the job queue, waiting for a time slot in the server. There is a watchdog managing the submitted jobs with the principle of first come first serve (FCFS). It checks the status of the job queue in the background at regular intervals. The earliest job will start to run when a time slot was assigned to it. Users can query their results through a jobid within one week after the job finished. In one week, it will expire automatically. For the protection of privacy, each user can only see these jobs submitted by themselves (with the same IP address) in the job list. It also allows users to set a password before launching a job, which can avoid privacy leaks when multiple users shares a same IP address.

In the result page, we present statistics of test datasets in the first part, which include nodes and edges of each input network, the final alignment score and input parameters etc (see in Fig. 1c). To visualize the process of simulated annealing, the convergence curve is plotted in the second part. From Fig. 1d, users can see how fast the computation can converge to a stable score. In the third part, there is a large table separated in many pages, each page contains at most ten matchsets (see in Fig. 1e). Each matchset implies a group of functionally conserved proteins, which can be used in the “annotation transfer”. Additionally, our web server can provide information from open accessible databases to annotate the alignment results, such as the Uniprot ID, and gene ontology annotations. Each protein is linked to its GenPept page in the open accessible database NCBI protein. Users can easily check their sequence similarity and download function annotations (GO terms) by a simple click. Besides, users can also search for specific matchsets with a pattern in a search box. For example, if a substring “P535” was queried in the search box, the protein accession identifiers matched to the pattern “%P355%” will be extracted from the result table. In the column of visualization, the graphical view of induced sub-networks would be extracted from its PPI networks (see in Fig. 1f). The test dataset and alignment results of each task are available in the result page, which makes it very easy to run other methods on our online datasets and to compare the alignment quality.

Results

As network topology and protein sequences were required for scoring the alignment, we manually collected a set of data from a series of openly accessible online resources storing protein-protein interactions, protein sequences, and protein annotation data (see in Fig. 2). Totally, 218,339 proteins and 4,009,541 interactions were extracted from six freely available databases, including IntAct, STRING, DIP, BioGRID, UniprotKB [30] and Ensembl [31] (see details in Fig. 3). Since proteins in different sources labeled by different identifiers, we convert all different sources identifiers into Uniprot identifier in the format acession:version, which is commonly used in many famous databases such as IntAct, QuikGO [32], UniProtKB/Swiss-Prot [33] and the NCBI protein databases [34]. This identifiers were further used to query the GO annotation in our web tools.

Fig. 2
figure2

Overview of the WebNetCoffee resources. We manually collected our testing datasets from eight openly accessible databases for up to 15 species. Protein-protein interactions (PPIs) and protein sequences were extracted from IntAct, STRING, DIP, BioGRID, UniProt, and Ensembl. Protein annotation data such as gi number and gene ontology annotations were extracted from the NCBI Protein and Uniprot-GOA project. Using the BLAST package, we performed pairwise sequence alignments to search for similar protein sequences

Fig. 3
figure3

Statistics of proteins and molecular interactions in our online datasets. All these 215,002 proteins and 4,005,485 interactions of 15 species are openly accessible to researchers

Proteins of up to 15 species and their associated annotation data including taxonomy information, gene ontology were collected to build our online datasets. To improve the data quality, interactions generated by co-complex expansion (in IntAct) and these with a small combined score (less than 900) in STRING were filtered away. The package of BLASTP [35] was performed to calculate the pairwise sequence similarity. These pairs of homologous proteins with sufficient common region (evalue <1e−7) are likely to be in a same matchset. Our online datasets including interactions and homologous proteins are openly accessible to researchers for performance comparison in the Download page. All our online datasets are regularly updated once every three months if there exist latest updates in the corresponding databases. Hopefully, it can promote the development of more advanced network alignment tools, facilitate the build of benchmark datasets, and lead us to a better understanding of molecular evolution and functions of these uncharacterized proteins.

Conclusion

In this paper, we present a fast and versatile web server based on the network alignment algorithm of NetCoffee to search for functional orthologous groups from multiple PPI networks. To the best of our knowledge, it is the first web server for globally aligning multiple PPI networks. Compared to NetCoffee, it also provides a friendly graphical user interface for easily performing GAMN tasks, managing submitted jobs and visualizing the alignment results. In the following work, we will support network alignment of multiple networks from more commonly used databases, and keep WebNetCoffee up-to-date with the latest data from these databases. Furthermore, the algorithm will be extended to take into account GO terms and pathways, which is expected to improve the prediction accuracy and facilitate more advanced applications such as prediction of disease-related genes.

Availability and requirements

Project name: WebNetcoffee

Project home page: http://www.nwpu-bioinformatics.com/WebNetCoffee/

Operating system(s): Platform independent

Programming language: C++, html, css, PHP

Other requirements: None

License: GNU GPL v3

Any restrictions to use by non-academics: license needed

Abbreviations

CLI:

Command line interface

FO:

Function-oriented ortholog

FCFS:

First come first serve

GAMN:

Global alignment for multiple networks

GO:

Gene ontology

PPI:

Protein-protein interactions

References

  1. 1

    Gupta PB, Onder TT, Jiang G, et al. Identification of selective inhibitors of cancer stem cells by high-throughput screening. Cell. 2009; 138(4):645–59.

  2. 2

    Hu J, Gao Y, Zheng Y, Shang X. Kf-finder: Identification of key factors from host-microbial networks in cervical cancer. BMC Syst Biol. 2018; 12(S4):54.

  3. 3

    Clark K, Karschmizrachi I, Lipman DJ, et al. Genbank. Nucleic Acids Res. 2016; 44(Database issue):67–72.

  4. 4

    Consortium U. Uniprot: a hub for protein information. Nucleic Acids Res. 2015; 43(Database issue):204–12.

  5. 5

    Ogata H, Goto S, Fujibuchi W, et al. Computation with the kegg pathway database. Biosystems. 1998; 47(1-2):119–28.

  6. 6

    Sussman J, Lin D, Jiang J, et al. Protein data bank (pdb): database of three-dimensional structural information of biological macromolecules. Acta Crystallogr. 2010; 54(6-1):1078–84.

  7. 7

    Sharan R, Ulitsky I, Shamir R. Network-based prediction of protein function. Mol Syst Biol. 2007; 3(88). https://doi.org/10.1038/msb4100129.

  8. 8

    Hu J, Shang X. Detection of network motif based on a novel graph canonization algorithm from transcriptional regulation networks. Molecules. 2017; 22(12):2194.

  9. 9

    Hu J, Reinert K. LocalAli: an evolutionary-based local alignment approach to identify functionally conserved modules in multiple networks. Bioinformatics. 2015; 31(3):363–72.

  10. 10

    Song L, Li D, Zeng X, Wu Y, Guo L, Zou Q. ndna-prot: identification of dna-binding proteins based on unbalanced classification. BMC Bioinforma. 2014; 15(1):298. https://doi.org/10.1186/1471-2105-15-298.

  11. 11

    Radivojac P, Clark WT, Oron TR, et al. A large-scale evaluation of computational protein function prediction. Nat Med. 2013; 10(3):221–7.

  12. 12

    Singh R, Xu J, Berger B. Global alignment of multiple protein interaction networks with application to functional orthology detection. Natl Acad Sci. 2008; 105:12763–8.

  13. 13

    Liao C-S, Lu K, Baym M, et al. Isorankn: spectral methods for global alignment of multiple protein networks. Bioinformatics. 2009; 25(12):253–8.

  14. 14

    Flannick J, Novak A, Srinivasan BS, et al. Graemlin: general and robust alignment of multiple large interaction networks. Genome Res. 2006; 16(9):1169.

  15. 15

    Sahraeian SME, Yoon B-J. Smetana: accurate and scalable algorithm for probabilistic alignment of large-scale biological networks. PLoS ONE. 2013; 8(7):67995.

  16. 16

    Neyshabur B, Khadem A, Hashemifar S, et al. Netal: a new graph-based method for global alignment of protein-protein interaction networks. Bioinformatics. 2013; 29(13):1654–62.

  17. 17

    Hu J, Kehr B, Reinert K. Netcoffee: a fast and accurate global alignment approach to identify functionally conserved proteins in multiple networks. Bioinformatics. 2014; 30(4):540–8.

  18. 18

    Saraph V, Milenković T. Magna: Maximizing accuracy in global network alignment. Bioinformatics. 2014; 30(20):2931–40.

  19. 19

    Vijayan V, Saraph V, Milenković T. Magna++: Maximizing accuracy in global network alignment via both node and edge conservation. Bioinformatics. 2015; 31(14):2409–11.

  20. 20

    Kalaev M, Smoot M, Ideker T, et al. Networkblast: comparative analysis of protein networks. Bioinformatics. 2008; 24(4):594–596. https://doi.org/10.1093/bioinformatics/btm630.

  21. 21

    Phan HTT, Sternberg M. Pinalog: a novel approach to align protein interaction networks–implications for complex detection and function prediction. Bioinformatics. 2012; 28(9):1239–45.

  22. 22

    Pache RA, Céol A, Aloy P. Netaligner–a network alignment server to compare complexes, pathways and whole interactomes. Nucleic Acids Res. 2012; 40(Web Server issue):157–61.

  23. 23

    Kelley BP, Yuan B, Lewitter F, et al. Pathblast: a tool for alignment of protein interaction networks. Nucleic Acids Res. 2004; 32(Web Server issue):83.

  24. 24

    Elkebir M, Brandt BW, Heringa J, et al. Natalieq: A web server for protein-protein interaction network querying. BMC Syst Biol. 2014; 8(1):1–7.

  25. 25

    Kerrien S, Aranda B, Breuza L, et al. The intact molecular interaction database in 2012. Nucleic Acids Res. 2012; 40(Database issue):841.

  26. 26

    Szklarczyk D, Franceschini A, Kuhn M, et al. The string database in 2011: functional interaction networks of proteins, globally integrated and scored. Nucleic Acids Res. 2011; 39(Database issue):561–8.

  27. 27

    Xenarios I, Salwínski L, Duan XJ, et al.Dip, the database of interacting proteins: a research tool for studying cellular networks of protein interactions. Nucleic Acids Res. 2002; 30(1):303.

  28. 28

    Chatr-Aryamontri A, Breitkreutz B-J, Heinicke S, et al.The biogrid interaction database: 2013 update. Nucleic Acids Res. 2013; 41(Database issue):816.

  29. 29

    Dezso B, Juttner A, Kovacs P. Lemon – an open source c++ graph template library. Electron Notes Theor Comput Sci. 2011; 264(5):23–45. https://doi.org/10.1016/j.entcs.2011.06.003. Proceedings of the Second Workshop on Generative Technologies (WGT) 2010.

  30. 30

    Boutet E, Lieberherr D, Tognolli MT, et al.Uniprotkb/swiss-prot. Methods Mol Biol. 2007; 406:89.

  31. 31

    Flicek P, Amode MR, Barrell D, et al. Ensembl 2014. Nucleic Acids Res. 2014; 42(Database issue):749–55.

  32. 32

    Binns D, Dimmer E, Huntley R, et al.Quickgo: a web-based tool for gene ontology searching. Bioinformatics. 2009; 25(22):3045.

  33. 33

    Gane PJ, Bateman A, Mj M, et al. Uniprot: A hub for protein information. Nucleic Acids Res. 2015; 43(Database issue):204–12.

  34. 34

    Coordinators NR. Database resources of the national center for biotechnology information. Nucleic Acids Res. 2018; 46(D1):8–13. https://doi.org/10.1093/nar/gkx1095.

  35. 35

    Camacho C, Coulouris G, Avagyan V, et al.Blast+: architecture and applications. BMC Bioinforma. 2009; 10(1):421.

Download references

Acknowledgements

Not applicable.

Funding

Publication costs were funded by the National Natural Science Foundation of China (Grant No. 61702420); This project has been funded by the National Natural Science Foundation of China (Grant No. 61332014, 61702420 and 61772426); the China Postdoctoral Science Foundation (Grant No. 2017M613203); the Natural Science Foundation of Shaanxi Province (Grant No. 2017JQ6037); the Fundamental Research Funds for the Central Universities (Grant No. 3102018zy032).

Availability of data and materials

All our datasets and material can be found at http://www.nwpu-bioinformatics.com/WebNetCoffee/download.php.

Author information

JLH designed the computational framework and workflow, developed the web-tool and wrote the manuscript; YQG developed the web-based tool and performed all data analyses; JHH and YZ collected all necessary data and preprocessed the raw datasets; XQS is the major coordinator, who contributed a lot of time and efforts in the discussion of this project. All authors read and approved the final manuscript.

Correspondence to Xuequn Shang.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

Not applicable.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License(http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver(http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Hu, J., Gao, Y., He, J. et al. WebNetCoffee: a web-based application to identify functionally conserved proteins from Multiple PPI networks. BMC Bioinformatics 19, 422 (2018) doi:10.1186/s12859-018-2443-4

Download citation

Keywords

  • Multiple network alignment
  • Webserver
  • PPI networks
  • Protein databases
  • Gene ontology