- Open Access
Modularity analysis based on predicted protein-protein interactions provides new insights into pathogenicity and cellular process of Escherichia coli O157:H7
© Wang et al; licensee BioMed Central Ltd. 2011
Received: 27 October 2011
Accepted: 22 December 2011
Published: 22 December 2011
With the development of experimental techniques and bioinformatics, the quantity of data available from protein-protein interactions (PPIs) is increasing exponentially. Functional modules can be identified from protein interaction networks. It follows that the investigation of functional modules will generate a better understanding of cellular organization, processes, and functions. However, experimental PPI data are still limited, and no modularity analysis of PPIs in pathogens has been published to date.
In this study, we predict and analyze the functional modules of E. coli O157:H7 systemically by integrating several bioinformatics methods. After evaluation, most of the predicted modules are found to be biologically significant and functionally homogeneous. Six pathogenicity-related modules were discovered and analyzed, including novel modules. These modules provided new information on the pathogenicity of O157:H7. The modularity of cellular function and cooperativity between modules are also discussed. Moreover, modularity analysis of O157:H7 can provide possible candidates for biological pathway extension and clues for discovering new pathways of cross-talk.
This article provides the first modularity analysis of a pathogen and sheds new light on the study of pathogens and cellular processes. Our study also provides a strategy for applying modularity analysis to any sequenced organism.
Most cellular processes are carried out by groups of physically interacting proteins. Protein-protein interactions (PPIs) are at the heart of biological activities. A complete and reliable interaction map representing the specific binary interactions within a cell would provide a significant platform for understanding many biochemically relevant processes.
Several high-throughput experimental methods - such as pull down , immunoprecipitation , two-hybrid system, protein chips [3, 4] - have been developed to detect PPIs among all the proteins encoded by a genome. While the data from these experimental approaches have been useful to biologists, there are several shortcomings. In particular, the results from high-throughput interaction mappings have low accuracy, and even reliable techniques can generate many false positives when applied genome-wide. Estimated error rates of high-throughput interaction results range from 41 to 90% [5, 6]. Detecting experimental interactions is also labor-intensive and costly, in part because the number of possible PPIs is very large.
Computational methods provide a complementary approach to detecting PPIs and extending protein interactomes. A variety of computational methods have been applied to observe or predict the PPI networks in biological systems. These methods enable one to discover novel putative interactions and often provide information for designing new experiments for specific protein sets. The computational approaches to in-silico prediction can be classified into several general categories: genomic scale approaches , sequence-based approaches [8, 9], structure-based approaches [10, 11], protein domain-based approaches , learning-based approaches , and network-topology-based approaches . Among these, certain domain-based prediction methods have shown great sensitivity and specificity to experimental PPIs, including MLE (Maximum Likelihood Estimation)  and MSSC (Maximum Specificity Set Cover) . These approaches infer potential DDIs by relying on a training set of PPIs, and then use these potential DDIs to predict PPIs in testing sets.
However, obtaining networks of PPIs is not the final target. A major challenge is how to manage and analyze the huge number of data on PPIs. It has been reported that the metabolic networks of 43 distinct organisms are organized into many small, highly connected topological modules that combine in a hierarchical manner into larger, less cohesive units . A module of a PPI network may represent a protein complex, or a group of proteins participating in the same cellular process. The prediction and analysis of PPI modules will aid us in elucidating the basic mechanisms of biological activities, while modularity analysis of the PPIs of pathogens could give us a better understanding of their pathogenicity.
Cluster analysis is an obvious choice of method for extracting functional modules from networks of PPIs. Clustering can be defined as the grouping of objects based on their shared discrete, measurable properties. A variety of clustering algorithms have been developed and successfully used in diverse fields. Recently, a systematic quantitative evaluation of the four most important clustering algorithms has been presented by Brohee and Van Helden . The four algorithms were RNSC, MCODE, SPC, and Markov Cluster algorithm (MCL). Their results showed that the MCL algorithm was both remarkably robust to graph alterations and superior in the extraction of complexes from interaction networks.
The bacterium Escherichia coli O157:H7, which causes diarrhea and hemolytic uremic syndrome (HUS), is a worldwide threat to public health and has been implicated in many outbreaks of hemorrhagic colitis. The death rate for infected populations is between 5 and 10 percent. Currently, there is no effective method for curing or preventing infection. In 2001, the U.S.A. and Japan published the genome sequences of the EHEC 0157:H7 EDL933 and Sakai strains, respectively, which made genome-scale research on O157:H7 possible [19, 20].
Several papers concerning modularity analysis of E. coli have been published in recent years. Von Mering et al. integrated techniques of conserved gene neighborhood, gene fusion events, and common phylogenic distribution to find functional modules of E. coli K12. By comparing to known metabolic pathways, they discussed pathway extension and functional links among pathways . Gerdes et al. gave a system level analysis of essential genes in E. coli MG1655 based on experimental results, and discussed these essential genes in topological modules . Li et al. developed a four-step approach for genome-wide discovery of parallel modules from protein functional linkages. This approach recovers known parallel complexes and pathways and discovers new ones . To date, however, there are no published papers referring to modularity analysis of pathogens, whose PPIs are under-studied.
In this paper, the EHEC O157:H7 Sakai strain is selected for further research. Our aim is to analyze the modularity of the pathogen E. coli O157:H7 PPI network, without any known experimental PPI data. We also want to see what can be interpreted about pathogenicity and cellular processes by the modularity analysis. First, a domain-based method was used to predict the PPIs of O157:H7. Then we used the Markov Cluster algorithm (MCL) and separated 172 modules out of the predicted O157 PPIs. After evaluation, we found that most of these modules were functionally homogeneous and biologically significant. One hundred and twenty-one modules were considered highly reliable and may provide directions for experimental research. Six pathogenicity-related modules were analyzed, some of which are new and deserve further experimental validation. After investigation of the relationships among modules, the modularity of cellular function and cooperative effects are discussed. In view of these modules, our analysis can provide a better understanding of cell function. Moreover, the predicted modules can provide possible candidates for biological pathway extension and clues for discovering new modes of cross-talk between pathways. Overall, these results provide the first modularity analysis of a pathogen and shed new light on the study of pathogenicity and cellular processes.
Prediction of E. coli O157:H7 PPIs
Predicted data of O157 protein interactions
Total protein of O157
Average PPIs per protein
Evaluation of predicted PPIs
The reliability of our predicted PPIs is the basis of this research, as the module prediction and further analysis are all founded on these data. Although MLE and MSSC were proved in their original papers to show good sensitivity and specificity, we still adopted several other methods (below) to validate the reliability of prediction.
Comparison of predicted PPIs with STRING database
Currently, there are too few experimental PPIs from E. coli O157 to evaluate our predictions. So we chose the STRING database, which collects experimental and predicted PPIs. The PPIs in the STRING database include direct and indirect (functional) associations, which are derived from four sources: genomic context, high-throughput experiments, co-expression, and previous knowledge.
We downloaded the PPI data from the STRING database and 548,828 PPIs of EHEC O157:H7 Sakai strain were selected, relating to 5201 proteins. After comparing our predicted PPIs with the STRING datasets, we found that 2478 (20.4%) of them overlapped. PPIs are dynamic, so the results of different assays or even the same assays at different times or under different conditions can vary. The overlap percentage here is not far from the ~25% overlap of the high-throughput experimental yeast interactome by Gavin  and Krogan  respectively. So we can conclude that our predicted PPIs overlapped well with the existing data in the STRING database.
Topological analysis of predicted PPIs network
To evaluate the predicted PPIs, we also analyzed the topological structure of the predicted PPI network to see if it had the same characteristics as the PPI networks generated by experiments.
Yook et al. compared four available databases that represent the protein interaction network of S. cerevisiae and found that the yeast protein interaction network in each database shows scale-free topology and hierarchical modularity . Li et al. analyzed three of the largest protein interaction networks of S. cerevisiae, C. elegans and D. melanogaster and also confirm the scale-free, small-world property .
Clustering coefficients can reflect the network's modularity property. It is defined as Ci = 2ni/ki(ki-1), where ni represents the number of direct links between the ki neighbors of node i. The clustering coefficient of a network is defined as the average value of all individual nodes. For the random network model, all the nodes share the same clustering coefficient: , where < k > is the average degree of the network and N is the total number of nodes. The relative clustering coefficient of the predicted PPI network is calculated and normalized by a random network of similar size . The result is 54.3, which indicates that the nature of the clustering is far from random.
Topological characteristics of experimental PPI networks  and predicted dataset of O157
No. of proteins
No. of interactions
Diameter of network
Relative clustering coefficient
Predicted dataset of O157
GO distance analysis of predicted PPIs
The GO annotations  of E. coli O157:H7 Sakai strain are included in the Uniprot GOA annotation and thus are downloaded from the GO database. For each pair of proteins in our predicted dataset and random dataset, we computed the GO distance. The smaller the value, the more specific a category the two proteins belong to, so they are more likely to interact. The average GO distance of the predicted PPI dataset is 0.50, while that of the random dataset is 0.92. Student's t test shows that the predicted PPI dataset has a significantly low GO distance towards the random dataset (t = 82.5, P < < 10-15).
Semantic similarity analysis of predicted PPIs
Besides the GO distance, the semantic similarity of GO terms is calculated to evaluate the predicted PPI dataset. For pairs of interacting proteins, GO terms in the "Biological Process" category are used to compute the semantic similarity score. The biological process aspect can be related to protein interactions, both physical and indirect (involved in same process). Of the 1651 proteins in the predicted PPI dataset, 1356 (82.1%) had the GO Biological Process annotation. The average semantic similarity of our predicted PPI dataset is 3.49, while the score of random set is 1.66. Student's t test shows that the predicted PPI dataset has a significantly high semantic similarity score (t = 52.6, P < < 10-15). The average scores of both datasets may be somehow higher, as we ignored interacting pairs lacking GO annotations in the "Biological Process" category.
Combining the comparison of STRING database, topological analysis, GO distance and semantic similarity analysis, we can tell that the predicted dataset is reliable enough for further analysis.
Prediction of modules from PPI network
The Markov cluster algorithm (MCL)  is used to predict modules from the PPIs obtained above, with an inflation coefficient of I = 1.8. However, the modules predicted by the MCL algorithm have no overlapping components, while in real organisms some proteins exist in multiple complexes or participate in several cellular processes at the same time. We therefore identified the proteins shared between modules by a post-processing step. The shared proteins are shown in Additional File 4 Table S1; 172 modules were separated and the size distribution of the final predicted modules dataset (Additional File 5) is displayed in Additional File 6, Figure S1.
Assessing the quality of derived modules
In order to determine whether the predicted protein modules are biologically significant and deserving of further research, we assessed the derived modules in three ways, described below.
GO annotation analysis of the modules
We used the GO Biological Process annotation to evaluate the functional coherence of the modules predicted above. As we mentioned earlier, 1356 (82.1%) of the 1651 proteins in predicted PPI dataset had GO Biological Process annotation. We used these proteins as the background set. For each module and each GO term, we computed the enrichment of the term in the module versus the background set using a hypergeometric test and derived a P-value for the module. These P-values were further FDR-corrected for multiple testing. For each module, we chose the term that yielded the highest level of significance. Finally, the results showed that 121 modules (70.3%) had an enriched GO annotation (P < 0.01), which means that our predicted modules have good functional coherence. Details of the enriched GO terms and P-values of the modules can be seen in Additional File 7, Table S1.
Comparison with known conserved protein complexes
Although there is no known database of the protein complexes of the E. coli O157:H7 Sakai strain, those complexes related to survival or reproduction will probably be conserved in other related organisms. We can find these modules from our predicted data, and then compare them with conserved protein complexes validated experimentally in other bacteria to evaluate the reliability of our predicted modules. To achieve this, we searched our modules in BOND (Biomolecular Object Network Databank) and in published references in PubMed. The results show that 55 of our predicted modules (32.0%) have complexes conserved in other bacteria, details of which can be seen in Additional File 8, Table S1. Three examples are given below.
Module 62 consisted of eight proteins (Additional File 9, Figure S1A), seven of which are components of the F1F0-type ATP synthase. This enzyme is of crucial importance almost ubiquitously because ATP is the common "energy currency" of cells [33, 34], so it is conserved among almost all organisms. Module 114 (Additional File 9, Figure S1B) consisted of three proteins, RecO, RecF, and RecR. The RecFOR complex is also conserved in other organisms; it modulates the assembly and disassembly of RecA filaments on DNA and is essential in DNA recombination, repair and replication [35, 36]. Module 126 (Additional File 9, Figure S1C) consisted of HypC, HypD, and HybG, which are also conserved in other bacteria. The complex functions in the assembly of the active site of the [NiFe]-hydrogenase enzymes .
Comparison with KEGG pathway database
The KEGG pathway database  collects pathway data from metabolic processes, genetic information processes, environmental information processes and cellular processes. As we mentioned earlier, a protein module may not only represent a protein complex, but could also refer to a group of proteins participating in the same biological pathway or cellular process. So comparing our predicted module with KEGG pathways is useful for the evaluation. We adopted the following method to achieve this goal. First, we used Uniprot GOA annotation to annotate the GO Term (Biological Process) for each KEGG pathway and an enriched GO Term was obtained for each KEGG pathway. If more than two proteins in our predicted modules belonged to KEGG pathways, we selected the most overlapping KEGG pathway for further comparison. If the enriched GO term of both the module and the most overlapping KEGG pathway were the same, we deemed that the module was reliable; otherwise, we decided that it was not. Eighty modules have at least two proteins that belong to 92 KEGG pathways. Among these, 51 (63.8%) modules have an enriched GO term identical with the KEGG pathway. Because of the incompleteness and bias of the KEGG pathway data, we can infer that other modules that have no overlap with a KEGG pathway may have the same rate of reliability.
Combining GO annotation, comparison with conserved protein complexes and KEGG pathways, we can conclude that the predicted modules have good functional homogeneity and biological significance, which gives us confidence for further discussion. In all of the 172 predicted modules, 121 functionally significantly enriched modules identified by GO annotation analysis were considered highly reliable, which can give clues for further research.
Predicted modules provide new information on pathogenicity
Two main factors contribute to the pathogenicity of EHEC O157:H7. First, the strain's ability to adhere and colonize, involving invasion, proliferation and resistance to the host immune system, makes it an effective pathogen. Second, the strain produces toxins, including endotoxin and exotoxin . Among our predicted modules, we found six functional modules that contain known virulence factors and may relate to E. coli O157's pathogenicity.
Novel module in cell adhesion
NCBI annotation, SignalP and TMpred prediction result of the seven central proteins
D-alanyl-D-alanine carboxypeptidase fraction A
glycyl-tRNA synthetase beta subunit
penicillin-binding protein 7
Landscape of O157:H7 iron acquisition system
Importing both siderophores and heme proteins into pathogens requires outer membrane receptors and TonB-ExbB-ExbD import systems , which is the case in module 35 in Figure 4. Besides iron, the uptake of vitamin B12 also depends on the TonB system, and BtuB in Figure 4 is the outer membrane receptor for B12 . The TolA, TolR, and TolQ elements in the yellow module share similarities with the TonB system and may contribute to the integrity of the cellular envelope .
Two proteins are annotated as "hypothetical protein" in the green module; the Blast result shows that they have high similarity to a TonB-dependent outer membrane receptor. We also used Signal P and Tmpred  to see whether they have signal peptides or transmembrane domains. Both of the proteins were predicted to have signal peptides, meaning that they are probably secreted proteins. One of them was predicted to have four possible transmembrane domains, while no transmembrane domains were predicted for the other. Further experiments are suggested to validate their roles as outer membrane receptors.
We can see from Figure 4 that there are multiple outer membrane receptors for iron, probably because iron is essential for O157; in this way the pathogen can increase its tolerance to mutations and ensure the continued import of iron. It has been reported that mutations in one of the outer membrane receptors do not affect iron acquisition, while a TonB mutant failed to use heme as an iron source or to utilize the siderophores, and showed reduced virulence [50, 51]. Our predicted interaction map can explain this well, as TonB plays a hub-like role, such that mutation or deletion of it will destroy the network and thereby cause inefficient iron acquisition.
Shiga toxin: the most important virulence factor of O157
Other pathogenicity-related modules
Module 88, displayed in Additional File 10, Figure S1A, contains Tir (Translocated Intimin Receptor), which is secreted by a type III secretion system (TTSS) and inserted into the host cell membrane as a receptor for intimin . CesT is Tir's chaperone, and it directs Tir to TTSS ATPase EscN . EscN functions in substrate recognition and in chaperone release from, and unfolding of, type III secreted proteins . Module 115 has three proteins, shown in Additional File 10, Figure S1B: the α, β, and γ subunits of urease. The complex catalyzes the hydrolysis of urea into ammonia and carbon dioxide. Expression of urease could modify internal and/or surrounding anion concentrations, enabling O157 to survive in acidic conditions and perhaps contributing to its low infectious dose . Module 130 consists of two periplasmic proteins, RseA and RseB, together with rpoE. RpoE is a sigma factor responsible for the transcription of genes for the cell envelope stress response. Under normal conditions, it is kept inactive by its interaction with the periplasmic proteins RseA and RseB. Under stress conditions, a protease in the periplasm cleaves the interaction between rpoE and the RseAB complex. Then rpoE is free to regulate the expression of genes associated with the stress response. Inactivation of rpoE diminishes bacterial survival and growth inside host macrophages, as it also regulates genes required for oxidative stress resistance .
From our study of the pathogenicity-related modules above, we can see that our modularity analysis of O157:H7 provided new information about the pathogen and clues for further experimental validation. This new information not only includes novel modules and speculations about the functions of new proteins, but also pathogenicity from the point of view of the module level. However, owing to false negatives in PPI prediction, other pathogenicity-related modules have not been shown here. With more complete and reliable experimental PPI data and more accurate in silico prediction methods in the future, the modularity analysis of pathogens will help us generate a better understanding of pathogenicity.
Investigation of relationships among predicted modules
Predicted modules demonstrate modularity of cell function
In the 21st century, research in cell biology is changing from molecular to modular, and from studying the function of only one gene or protein to investigating how a group of biological molecules functions as a module. The conception of modularity of cell function is that cell functions are carried out by different modules, which comprise many types of molecules. Through studying the relationships among modules, we found that our predicted modules demonstrate the modularity of cell function.
Numbers of modules related to other cellular functions in our predicted data
Number of related modules(module ID)
9 (66, 67, 81, 82, 86, 98, 108, 114, 119)
6 (1, 13, 25, 44, 56, 83)
Signal transducer activity
4 (8, 17, 28, 29)
Structural constituent of ribosome
7 (91, 92, 97, 122, 127, 128, 151)
8 (2, 4, 14, 15, 75, 35, 94, 96)
Transcription factor (regulator)
7 (5, 7, 24, 27, 48, 74, 112)
There are still many problems to solve regarding the modularity of cell function, such as how other cell functions are modular, or what the relationship is among modules fulfilling the same cellular function or how they coordinate and assemble. Once these problems are solved, the era of synthetic biology will really arrive. We will be able to design and synthesize a functional module in the same way that we design complex electronic modules or chips for a computer. Moreover, we will be able to arrange different functional modules properly, so as to simulate the cellular environment or synthesize new life.
Predicted modules show cooperative effects
Cooperativity among biological molecules, for example enzymes, is established beyond doubt. We investigate the relationships among our predicted modules and find that cooperativity also obtains at the module level. The effect of different modules with similar functions cooperating to fulfill a cellular process is described as a positive cooperative effect; if these modules perform opposite functions, then we call it a negative cooperative effect.
Positive cooperative effect
We can see clearly that these four modules of two systems carry out similar functions and cooperate to fulfill the cell's protein export functions.
Negative cooperative effect
In Figure 7, we have mentioned that the green module plays a role in forming a dynamic Z-ring structure at the midcell boundary and stabilizing this structure. Since the components in the yellow module play an opposing role, they are destabilizing factors of the Z-ring. These two modules together shift the equilibrium of FtsZ expression between an unassembled cytoplasmic pool and the assembled ring. The negative cooperative effect between these two modules ensures that the Z-ring only forms in the midcell instead of other places in the cell, so cell division is accomplished properly.
Of course, besides the cooperative effect, there might be other relationships among modules such as regulation, activation, deactivation and so on. But investigating these relationships requires other information and remains a great challenge.
Predicted modules provide possible candidates for biological pathway extension
Biological pathways can be regarded as a collection of known relationships or reactions among biological objects (i.e., genes or gene products). However, knowledge about biological pathways is not complete, and known pathways are insufficient to cover all genes or gene products. In the case of humans, for example, about 3,000 genes are covered by the major biological pathway databases. The rest, and their relationships with other genes or gene products, remain to be explored .
In the previous assessment, we found modules that overlapped with the KEGG pathway and had the same enriched GO terms. In these modules, proteins not overlapping may be possible candidates for KEGG and may provide information for pathway extension. Two examples are given below.
Module 4 consists of 53 proteins, 51 of which overlap with KEGG pathway ecs02010: ATP-binding cassette (ABC) transporters (Additional File 11, Figure S1). ABC transporters play an important role in bacteria, importing various nutrients required for survival in different niches and exporting substances toxic to the cell. Characteristically, ABC transporters have three components: a substrate-binding protein, a permease protein, and an ATP-binding protein. The NCBI annotations of the remaining two proteins in module 4 were "putative binding-protein dependent transport protein" and "putative transport system permease protein", respectively. This suggests that these two proteins are possible candidates for KEGG pathway extension and deserve further experimental validation.
From the two examples above, we can see that modularity analysis of O157:H7 has the capacity to provide possible candidates for, and facilitate research into, biological pathway extension.
Predicted modules gave clues for discovering new cross-talk
Few pathways are isolated. Cross-talk among pathways links them into a complex network . Such a network, with cross-talk linking distinct pathways, would confer on the cell a sophisticated ability to sense multiple environmental signals impinging upon it, thus providing a means of adapting or regulating its response to a particular range of effectors . Major new challenges have arisen from attempts to identify cross-talk among pathways.
Modularity analysis of PPIs has become an important and challenging topic in life science. It will be instructive for scientists exploring the basic rules of cell processes or studying the mechanisms of pathogens. In this research, we integrated several bioinformatics methods for modularity analysis of E. coli O157:H7. We predicted PPIs inside the bacterium, then we derived predicted modules. Evaluation showed that these predicted modules were functional, enriched, and biologically reliable. One hundred and forty-four highly reliable modules were provided as directions for experimental research. Through our modularity analysis, pathogenicity-related modules can be found, which provide new information regarding the mechanisms of pathogenicity and enable us to study pathogens in a new light. Our predicted modules also suggest that some cellular functions are modular and give a more comprehensive understanding of cell functions. Cooperativity among modules was discussed. Moreover, our modularity analysis of O157:H7 provides possible candidates for the extension of biological pathways and clues for discovering new cross-talk between pathways. The method for modularity analysis provided in this study can be applied to other pathogens or any other organism of interest that has been sequenced.
The bottleneck of this method is its reliance on the prediction of PPIs. Because there are no available, experimentally validated data on PPIs for E. coli O157:H7, we adopted a domain-based bioinformatics method to predict such PPIs. Although this method has been successful and has some advantages compared to other predictive methods, we can see from the results that less than 1/3 of the proteins of O157 were predicted to have PPIs, which caused high false negative rates. False positives also exist in this incomplete set of PPIs. These problems limited later modularity analysis. We believe that with the increase of experimental data and the development of analytical methods, modularity analysis of all biological interactions will greatly facilitate research in life science and the development of synthetic biology.
Materials and methods
All datasets used in this study were downloaded in November 2007. DIP  provides 3722 creditable protein interactions, which are validated by two or more experimental methods. These 3722 PPIs were selected from 54,511 PPIs deposited in the DIP database referring to more than 200 organisms, including both non-pathogens and pathogens. Both direct (physical) and indirect (functional) associations are included among the 548,828 PPIs of E. coli O157:H7 used for overlap analysis and downloaded from the String database. The PPIs in String are derived from four sources: genomic context, high-throughput experiments, coexpression, and previous knowledge. Sequences for the proteins of E. coli O157:H7 are downloaded from the NCBI Refseq database in Fasta format. One hundred and twelve pathways of O157 are downloaded from the KEGG pathway database .
Prediction of E. coli O157 PPIs
Maximum Likelihood Estimation (MLE) and Maximum Specificity Set Cover (MSSC) are both based on the Association Method (AM). These methods use currently available PPI data, and estimate the probabilities of domain-domain interaction observed in PPIs. Using the inferred domain-domain interaction, these methods can then predict previously unknown protein interactions. As MLE and MSSC modify AM in different and independent ways to improve accuracy, we combined these two methods to achieve a better result. In this work, the program InterProScan version 12.0 is adopted to scan domains of the creditable protein interactions from DIP. Then MLE and MSSC are used to build a domain interaction matrix for the prediction of PPIs. InterProScan was used again to scan domains of proteins of E. coli O157:H7 Sakai. The total matched pairs of the O157:H7 proteins were compared with the domain interaction matrix, and raw predicted PPI data were obtained. Two post-processing steps were applied to the raw PPI data. First, we eliminated directional repeats from the PPIs. Because the prediction program cannot predict weighted directional PPIs, directional PPIs are actually the same. Second, we eliminated self-interactions. The existence of self-interactions will generate single protein modules when the MCL algorithm is used for prediction. Although simple protein modules may represent homogeneous multimers, these would be difficult to analyze without additional information.
Computation of GO distance
In this formula, T GO are the sets of the proteins' associated GO terms, while |T GO | stands for their number of elements and Δ is the symmetrical difference between the two sets. This distance formula emphasizes the importance of the shared GO terms by giving more weight to similarities than to differences. Therefore, for two genes that share no GO terms, the distance is 1, while for two proteins sharing exactly the same set of GO terms, the distance is 0. Based on this algorithm, some Perl programs were written to implement the calculation and analysis of GO distance in our study.
Calculation of GO semantic similarity
We used an R based package, csbl.go http://csbi.ltdk.helsinki.fi/csbl.go, to calculate GO semantic similarity. The package can compute similarities for arbitrary numbers of genes and supports custom GO annotation tables. For pairs of interacting proteins, the semantic similarity of GO terms in the "Biological process" taxonomy is calculated on the basis of the Resnik measure.
Prediction of PPI modules
The Markov Cluster algorithm (MCL)  was downloaded at http://micans.org/mcl. The algorithm simulates flow on the PPI graph by constructing its adjacency matrix and computing its successive powers to increase the contrast between regions with high flow and regions with low flow. The graph was then partitioned into high-flow regions corresponding to protein modules, separated by regions of no flow.
In the MCL algorithm, the inflation coefficient (I) is the main value that affects cluster granularity. A big I value will tend to result in fine-grained clusters. We selected five I values to predict protein modules: 1.4, 1.8, 2.2, 2.6, and 3.0. Since there is no known protein complex database available for E. coli O157 that would allow us to identify the optimal value, we used the MCL evaluation program combined with references and finally chose I = 1.8 for further research.
The modules predicted by the MCL algorithm have no overlapping components, but in real organisms some proteins exist in multiple complexes or participate in several cellular processes at the same time. Shuye Pu et al.  solved this problem by a step of simple post-processing. In that step, proteins assigned to each cluster (donor cluster) were scanned for interaction partners in other clusters, and proteins interacting with a sufficiently large fraction of partners in another cluster (acceptor cluster) were also assigned to that cluster.
A map of predicted PPIs was drawn using Pajek software downloaded at http://vlado.fmf.uni-lj.si/pub/networks/pajek. Other PPI maps of modules were drawn by Cytoscape 2.4.0, which can be downloaded at http://www.cytoscape.org/download.php
The TIGR tools we used to categorize proteins can be obtained at http://www.jcvi.org/cms/research/projects/cmr/overview/. The TMpred program makes a prediction of membrane-spanning regions and their orientation. The program can be used at http://www.ch.embnet.org/software/TMPRED_form.html. SignalP predicts the presence and location of signal peptide cleavage sites in amino acid sequences from different organisms: Gram-positive prokaryotes, Gram-negative prokaryotes, and eukaryotes. It is also an online tool and can be obtained at http://www.cbs.dtu.dk/services/SignalP/
We are grateful to Chao Niu, Yan Ling, and Yuxia Li for helpful comments and discussion of the work.
- Martzen MR, McCraith SM, Spinelli SL, Torres FM, Fields S, Grayhack EJ, Phizicky EM: A biochemical genomics approach for identifying genes by the activity of their products. Science. 1999, 286: 1153-1155. 10.1126/science.286.5442.1153.View ArticlePubMedGoogle Scholar
- Phizicky EM, Fields S: Protein-Protein Interaction - Method for Detection and Analysis. Microbiological Reviews. 1995, 59: 94-123.PubMed CentralPubMedGoogle Scholar
- Chen Y, Xu D: Computational analyses of high-throughput protein-protein interaction data. Current Protein & Peptide Science. 2003, 4: 159-180. 10.2174/1389203033487225.View ArticleGoogle Scholar
- Pelletier J, Sidhu S: Mapping protein-protein interactions with combinatorial biology methods. Current Opinion in Biotechnology. 2001, 12: 340-347. 10.1016/S0958-1669(00)00225-1.View ArticlePubMedGoogle Scholar
- von Mering C, Krause R, Snel B, Cornell M, Oliver SG, Fields S, Bork P: Comparative assessment of large-scale data sets of protein-protein interactions. Nature. 2002, 417: 399-403.View ArticlePubMedGoogle Scholar
- Yu J, Fotouhi F: Computational approaches for predicting protein-protein interactions: a survey. J Med Syst. 2006, 30: 39-44. 10.1007/s10916-006-7402-3.View ArticlePubMedGoogle Scholar
- Sun JC, Sun Y, Ding GH, Liu Q, Wang C, He YY, Shi TL, Li YX, Zhao ZM: InPrePPI: an integrated evaluation method based on genomic context for predicting protein-protein interactions in prokaryotic genomes. Bmc Bioinformatics. 2007, 8: 10-10.1186/1471-2105-8-10.View ArticleGoogle Scholar
- Shen J, Zhang J, Luo X, Zhu W, Yu K, Chen K, Li Y, Jiang H: Predicting protein-protein interactions based only on sequences information. Proc Natl Acad Sci USA. 2007, 104: 4337-4341. 10.1073/pnas.0607879104.PubMed CentralView ArticlePubMedGoogle Scholar
- Bock JR, Gough DA: Predicting protein--protein interactions from primary structure. Bioinformatics. 2001, 17: 455-460. 10.1093/bioinformatics/17.5.455.View ArticlePubMedGoogle Scholar
- Aytuna AS, Gursoy A, Keskin O: Prediction of protein-protein interactions by combining structure and sequence conservation in protein interfaces. Bioinformatics. 2005, 21: 2850-2855. 10.1093/bioinformatics/bti443.View ArticlePubMedGoogle Scholar
- Sprinzak E, Margalit H: Correlated sequence-signatures as markers of protein-protein interaction. J Mol Biol. 2001, 311: 681-692. 10.1006/jmbi.2001.4920.View ArticlePubMedGoogle Scholar
- Singhal M, Resat H: A domain-based approach to predict protein-protein interactions. Bmc Bioinformatics. 2007, 8: 19-10.1186/1471-2105-8-19.View ArticleGoogle Scholar
- Neugebauer A, Hartmann RW, Klein CD: Prediction of protein-protein interaction inhibitors by chemoinformatics and machine learning methods. Journal of Medicinal Chemistry. 2007, 50: 4665-4668. 10.1021/jm070533j.View ArticlePubMedGoogle Scholar
- Yu H, Paccanaro A, Trifonov V, Gerstein M: Predicting interactions in protein networks by completing defective cliques. Bioinformatics. 2006, 22: 823-829. 10.1093/bioinformatics/btl014.View ArticlePubMedGoogle Scholar
- Deng M, Mehta S, Sun F, Chen T: Inferring domain-domain interactions from protein-protein interactions. Genome Res. 2002, 12: 1540-1548. 10.1101/gr.153002.PubMed CentralView ArticlePubMedGoogle Scholar
- Huang CB, Morcos F, Kanaan SP, Wuchty S, Chen DZ, Izaguirre JA: Predicting protein-protein interactions from protein domains using a set cover approach. Ieee-Acm Transactions on Computational Biology and Bioinformatics. 2007, 4: 78-87.View ArticlePubMedGoogle Scholar
- Ravasz E, Somera AL, Mongru DA, Oltvai ZN, Barabasi AL: Hierarchical organization of modularity in metabolic networks. Science. 2002, 297: 1551-1555. 10.1126/science.1073374.View ArticlePubMedGoogle Scholar
- Brohee S, van Helden J: Evaluation of clustering algorithms for protein-protein interaction networks. BMC Bioinformatics. 2006, 7: 488-10.1186/1471-2105-7-488.PubMed CentralView ArticlePubMedGoogle Scholar
- Hayashi T, Makino K, Ohnishi M, Kurokawa K, Ishii K, Yokoyama K, Han CG, Ohtsubo E, Nakayama K, Murata T: Complete genome sequence of enterohemorrhagic Escherichia coli O157:H7 and genomic comparison with a laboratory strain K-12. DNA Res. 2001, 8: 11-22. 10.1093/dnares/8.1.11.View ArticlePubMedGoogle Scholar
- Perna NT, Plunkett G, Burland V, Mau B, Glasner JD, Rose DJ, Mayhew GF, Evans PS, Gregor J, Kirkpatrick HA: Genome sequence of enterohaemorrhagic Escherichia coli O157:H7. Nature. 2001, 409: 529-533. 10.1038/35054089.View ArticlePubMedGoogle Scholar
- von Mering C, Zdobnov EM, Tsoka S, Ciccarelli FD, Pereira-Leal JB, Ouzounis CA, Bork P: Genome evolution reveals biochemical networks and functional modules. Proceedings of the National Academy of Sciences of the United States of America. 2003, 100: 15428-15433. 10.1073/pnas.2136809100.PubMed CentralView ArticlePubMedGoogle Scholar
- Gerdes SY, Scholle MD, Campbell JW, Balazsi G, Ravasz E, Daugherty MD, Somera AL, Kyrpides NC, Anderson I, Gelfand MS: Experimental determination and system level analysis of essential genes in Escherichia coli MG1655. Journal of Bacteriology. 2003, 185: 5673-5684. 10.1128/JB.185.19.5673-5684.2003.PubMed CentralView ArticlePubMedGoogle Scholar
- Li HY, Pellegrini M, Eisenberg D: Detection of parallel functional modules by comparative analysis of genome sequences. Nature Biotechnology. 2005, 23: 253-260. 10.1038/nbt1065.View ArticlePubMedGoogle Scholar
- Xenarios I, Fernandez E, Salwinski L, Duan XJ, Thompson MJ, Marcotte EM, Eisenberg D: DIP: The Database of Interacting Proteins: 2001 update. Nucleic Acids Res. 2001, 29: 239-241. 10.1093/nar/29.1.239.PubMed CentralView ArticlePubMedGoogle Scholar
- Quevillon E, Silventoinen V, Pillai S, Harte N, Mulder N, Apweiler R, Lopez R: InterProScan: protein domains identifier. Nucleic Acids Res. 2005, 33: W116-120. 10.1093/nar/gki442.PubMed CentralView ArticlePubMedGoogle Scholar
- Gavin AC, Aloy P, Grandi P, Krause R, Boesche M, Marzioch M, Rau C, Jensen LJ, Bastuck S, Dumpelfeld B: Proteome survey reveals modularity of the yeast cell machinery. Nature. 2006, 440: 631-636. 10.1038/nature04532.View ArticlePubMedGoogle Scholar
- Krogan NJ, Cagney G, Yu H, Zhong G, Guo X, Ignatchenko A, Li J, Pu S, Datta N, Tikuisis AP: Global landscape of protein complexes in the yeast Saccharomyces cerevisiae. Nature. 2006, 440: 637-643. 10.1038/nature04670.View ArticlePubMedGoogle Scholar
- Yook SH, Oltvai ZN, Barabasi AL: Functional and topological characterization of protein interaction networks. Proteomics. 2004, 4: 928-942. 10.1002/pmic.200300636.View ArticlePubMedGoogle Scholar
- Li D, Li J, Ouyang S, Wang J, Wu S, Wan P, Zhu Y, Xu X, He F: Protein interaction networks of Saccharomyces cerevisiae, Caenorhabditis elegans and Drosophila melanogaster: large-scale organization and robustness. Proteomics. 2006, 6: 456-461. 10.1002/pmic.200500228.View ArticlePubMedGoogle Scholar
- Barabasi AL, Albert R: Emergence of scaling in random networks. Science. 1999, 286: 509-512. 10.1126/science.286.5439.509.View ArticlePubMedGoogle Scholar
- Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, Cherry JM, Davis AP, Dolinski K, Dwight SS, Eppig JT: Gene ontology: tool for the unification of biology. The Gene Ontology Consortium. Nat Genet. 2000, 25: 25-29. 10.1038/75556.PubMed CentralView ArticlePubMedGoogle Scholar
- Enright AJ, Van Dongen S, Ouzounis CA: An efficient algorithm for large-scale detection of protein families. Nucleic Acids Res. 2002, 30: 1575-1584. 10.1093/nar/30.7.1575.PubMed CentralView ArticlePubMedGoogle Scholar
- Senior AE, Nadanaciva S, Weber J: The molecular mechanism of ATP synthesis by F1F0-ATP synthase. Biochim Biophys Acta. 2002, 1553: 188-211. 10.1016/S0005-2728(02)00185-8.View ArticlePubMedGoogle Scholar
- Weber J, Senior AE: ATP synthesis driven by proton transport in F1F0-ATP synthase. FEBS Lett. 2003, 545: 61-70. 10.1016/S0014-5793(03)00394-6.View ArticlePubMedGoogle Scholar
- Kuzminov A: Recombinational repair of DNA damage in Escherichia coli and bacteriophage lambda. Microbiol Mol Biol Rev. 1999, 63: 751-813. table of contentsPubMed CentralPubMedGoogle Scholar
- Madiraju MV, Templin A, Clark AJ: Properties of a mutant recA-encoded protein reveal a possible role for Escherichia coli recF-encoded protein in genetic recombination. Proc Natl Acad Sci USA. 1988, 85: 6592-6596. 10.1073/pnas.85.18.6592.PubMed CentralView ArticlePubMedGoogle Scholar
- Butland G, Zhang JW, Yang WH, Sheung A, Wong P, Greenblatt JF, Emili A, Zamble DB: Interactions of the Escherichia coli hydrogenase biosynthetic proteins: HybG complex formation. Febs Letters. 2006, 580: 677-681. 10.1016/j.febslet.2005.12.063.View ArticlePubMedGoogle Scholar
- Kanehisa M, Goto S, Kawashima S, Nakaya A: The KEGG databases at GenomeNet. Nucleic Acids Res. 2002, 30: 42-46. 10.1093/nar/30.1.42.PubMed CentralView ArticlePubMedGoogle Scholar
- Kaper JB, Nataro JP, Mobley HLT: Pathogenic Escherichia coli. Nature Reviews Microbiology. 2004, 2: 123-140. 10.1038/nrmicro818.View ArticlePubMedGoogle Scholar
- Torres AG, Giron JA, Perna NT, Burland V, Blattner FR, Avelino-Flores F, Kaper JB: Identification and characterization of lpfABCC'DE, a fimbrial operon of enterohemorrhagic Escherichia coli O157:H7. Infect Immun. 2002, 70: 5416-5427. 10.1128/IAI.70.10.5416-5427.2002.PubMed CentralView ArticlePubMedGoogle Scholar
- Langermann S, Palaszynski S, Barnhart M, Auguste G, Pinkner JS, Burlein J, Barren P, Koenig S, Leath S, Jones CH, Hultgren SJ: Prevention of mucosal Escherichia coli infection by FimH-adhesin-based systemic vaccination. Science. 1997, 276: 607-611. 10.1126/science.276.5312.607.View ArticlePubMedGoogle Scholar
- Schaeffer AJ, Schwan WR, Hultgren SJ, Duncan JL: Relationship of type 1 pilus expression in Escherichia coli to ascending urinary tract infections in mice. Infect Immun. 1987, 55: 373-380.PubMed CentralPubMedGoogle Scholar
- Bendtsen JD, Nielsen H, von Heijne G, Brunak S: Improved prediction of signal peptides: SignalP 3.0. J Mol Biol. 2004, 340: 783-795. 10.1016/j.jmb.2004.05.028.View ArticlePubMedGoogle Scholar
- Crosa JH: Genetics and Molecular-Biology of Siderophore-Mediated Iron Transport in Bacteria. Microbiological Reviews. 1989, 53: 517-530.PubMed CentralPubMedGoogle Scholar
- Law D, Kelly J: Use of Heme and Hemoglobin by Escherichia-Coli O157 and Other Shiga-Like-Toxin-Producing Escherichia-Coli Serogroups. Infection and Immunity. 1995, 63: 700-702.PubMed CentralPubMedGoogle Scholar
- Braun V: Energy-Coupled Transport and Signal-Transduction through the Gram-Negative Outer-Membrane Via Tonb-Exbb-Exbd-Dependent Receptor Proteins. Fems Microbiology Reviews. 1995, 16: 295-307. 10.1111/j.1574-6976.1995.tb00177.x.View ArticlePubMedGoogle Scholar
- Taylor R, Burgner JW, Clifton J, Cramer WA: Purification and characterization of monomeric Escherichia coli vitamin B12 receptor with high affinity for colicin E3. J Biol Chem. 1998, 273: 31113-31118. 10.1074/jbc.273.47.31113.View ArticlePubMedGoogle Scholar
- Lloubes R, Cascales E, Walburger A, Bouveret E, Lazdunski C, Bernadac A, Journet L: The Tol-Pal proteins of the Escherichia coli cell envelope: an energized system required for outer membrane integrity?. Res Microbiol. 2001, 152: 523-529. 10.1016/S0923-2508(01)01226-8.View ArticlePubMedGoogle Scholar
- Hofmann K, Stoffel W: TMbase - A database of membrane spanning proteins segments. Biol Chem Hoppe-Seyler. 1993, 374: 166-Google Scholar
- Perry RD, Shah J, Bearden SW, Thompson JM, Fetherston JD: Yersinia pestis TonB: Role in iron, heme, and hemoprotein utilization. Infection and Immunity. 2003, 71: 4159-4162. 10.1128/IAI.71.7.4159-4162.2003.PubMed CentralView ArticlePubMedGoogle Scholar
- Russo TA, McFadden CD, Carlino-MacDonald UB, Beanan JM, Barnard TJ, Johnson JR: IroN functions as a siderophore receptor and is a urovirulence factor in an extraintestinal pathogenic isolate of Escherichia coli. Infection and Immunity. 2002, 70: 7156-7160. 10.1128/IAI.70.12.7156-7160.2002.PubMed CentralView ArticlePubMedGoogle Scholar
- Louise CB, Obrig TG: Specific Interaction of Escherichia-Coli O157/H7-Derived Shiga-Like Toxin-Ii with Human Renal Endothelial-Cells. Journal of Infectious Diseases. 1995, 172: 1397-1401. 10.1093/infdis/172.5.1397.View ArticlePubMedGoogle Scholar
- Siegler RL, Obrig TG, Pysher TJ, Tesh VL, Denkers ND, Taylor FB: Response to Shiga toxin 1 and 2 in a baboon model of hemolytic uremic syndrome. Pediatric Nephrology. 2003, 18: 92-96.PubMedGoogle Scholar
- Ostroff SM, Kobayashi JM, Lewis JH: Infections with Escherichia-Coli O157-H7 in Washington State - the 1st Year of Statewide Disease Surveillance. Jama-Journal of the American Medical Association. 1989, 262: 355-359. 10.1001/jama.1989.03430030043031.View ArticleGoogle Scholar
- Kenny B, DeVinney R, Stein M, Reinscheid DJ, Frey EA, Finlay BB: Enteropathogenic E-coli (EPEC) transfers its receptor for intimate adherence into mammalian cells. Cell. 1997, 91: 511-520. 10.1016/S0092-8674(00)80437-7.View ArticlePubMedGoogle Scholar
- Abe A, de Grado M, Pfuetzner RA, Sanchez-SanMartin C, DeVinney R, Puente JL, Strynadka NCJ, Finlay BB: Enteropathogenic Escherichia coli translocated intimin receptor, Tir, requires a specific chaperone for stable secretion. Molecular Microbiology. 1999, 33: 1162-1175.View ArticlePubMedGoogle Scholar
- Gauthier A, Finlay BB: Translocated intimin receptor and its chaperone interact with ATPase of the type III secretion apparatus of enteropathogenic Escherichia coli. Journal of Bacteriology. 2003, 185: 6747-6755. 10.1128/JB.185.23.6747-6755.2003.PubMed CentralView ArticlePubMedGoogle Scholar
- Heimer SR, Welch RA, Perna NT, Posfai G, Evans PS, Kaper JB, Blattner FR, Mobley HL: Urease of enterohemorrhagic Escherichia coli: evidence for regulation by fur and a trans-acting factor. Infect Immun. 2002, 70: 1027-1031. 10.1128/IAI.70.2.1027-1031.2002.PubMed CentralView ArticlePubMedGoogle Scholar
- Missiakas D, Mayer MP, Lemaire M, Georgopoulos C, Raina S: Modulation of the Escherichia coli sigma(E) (RpoE) heat-shock transcription-factor activity by the RseA, RseB and RseC proteins. Molecular Microbiology. 1997, 24: 355-371. 10.1046/j.1365-2958.1997.3601713.x.View ArticlePubMedGoogle Scholar
- Hartwell LH, Hopfield JJ, Leibler S, Murray AW: From molecular to modular cell biology. Nature. 1999, 402: C47-52. 10.1038/35011540.View ArticlePubMedGoogle Scholar
- Bi E, Lutkenhaus J: Ftsz Ring Structure Associated with Division in Escherichia-Coli. Nature. 1991, 354: 161-164. 10.1038/354161a0.View ArticlePubMedGoogle Scholar
- Lu CL, Stricker J, Erickson HP: FtsZ from Escherichia coli, Azotobacter vinelandii, and Thermotoga maritima - Quantitation, GTP hydrolysis, and assembly. Cell Motility and the Cytoskeleton. 1998, 40: 71-86. 10.1002/(SICI)1097-0169(1998)40:1<71::AID-CM7>3.0.CO;2-I.View ArticlePubMedGoogle Scholar
- Hale CA, deBoer PAJ: Direct binding of FtsZ to ZipA, an essential component of the septal ring structure that mediates cell division in E-coli. Cell. 1997, 88: 175-185. 10.1016/S0092-8674(00)81838-3.View ArticlePubMedGoogle Scholar
- Margolin W: Spatial regulation of cytokinesis in bacteria. Current Opinion in Microbiology. 2001, 4: 647-652. 10.1016/S1369-5274(01)00264-8.View ArticlePubMedGoogle Scholar
- Ebersbach G, Gerdes K: Plasmid segregation mechanisms. Annu Rev Genet. 2005, 39: 453-479. 10.1146/annurev.genet.38.072902.091252.View ArticlePubMedGoogle Scholar
- Ikeda M, Wachi M, Jung H, Ishino F, Matsuhashi M: Homology among MurC, MurD, MurE and MurF proteins in Escherichia coli and that between E. coli MurG and a possible MurG protein in Bacillus subtilis. J Genet Appl Microbiol. 1990, 36: 179-187. 10.2323/jgam.36.179.View ArticleGoogle Scholar
- Goffin C, Ghuysen JM: Multimodular penicillin binding proteins: An enigmatic family of orthologs and paralogs. Microbiology and Molecular Biology Reviews. 1998, 62: 1079-+.PubMed CentralPubMedGoogle Scholar
- Hanada M, Nishiyama KI, Mizushima S, Tokuda H: Reconstitution of an efficient protein translocation machinery comprising SecA and the three membrane proteins, SecY, SecE, and SecG (p12). J Biol Chem. 1994, 269: 23625-23631.PubMedGoogle Scholar
- Fekkes P, van der Does C, Driessen AJ: The molecular chaperone SecB is released from the carboxy-terminus of SecA during initiation of precursor protein translocation. Embo J. 1997, 16: 6105-6113. 10.1093/emboj/16.20.6105.PubMed CentralView ArticlePubMedGoogle Scholar
- Powers T, Walter P: Co-translational protein targeting catalyzed by the Escherichia coli signal recognition particle and its receptor. Embo J. 1997, 16: 4880-4886. 10.1093/emboj/16.16.4880.PubMed CentralView ArticlePubMedGoogle Scholar
- Robinson C, Bolhuis A: Protein targeting by the twin-arginine translocation pathway. Nat Rev Mol Cell Biol. 2001, 2: 350-356. 10.1038/35073038.View ArticlePubMedGoogle Scholar
- Chung TS, Kim J, Kim K, Kim JH: Biological Pathway Extension Using Microarray Gene Expression Data. Genomics & Informatics. 2008, 6: 202-209. 10.5808/GI.2008.6.4.202.View ArticleGoogle Scholar
- Schwartz MA, Ginsberg MH: Networks and crosstalk: integrin signalling spreads. Nature Cell Biology. 2002, 4: E65-E68. 10.1038/ncb0402-e65.View ArticlePubMedGoogle Scholar
- Houslay MD: 'Crosstalk': a pivotal role for protein kinase C in modulating relationships between signal transduction pathways. Eur J Biochem. 1991, 195: 9-27. 10.1111/j.1432-1033.1991.tb15671.x.View ArticlePubMedGoogle Scholar
- Brun C, Chevenet F, Martin D, Wojcik J, Guenoche A, Jacq B: Functional classification of proteins for the prediction of cellular function from a protein-protein interaction network. Genome Biology. 2004, 5:Google Scholar
- Ovaska K, Laakso M, Hautaniemi S: Fast gene ontology based clustering for microarray experiments. BioData Min. 2008, 1: 11-10.1186/1756-0381-1-11.PubMed CentralView ArticlePubMedGoogle Scholar
- Pu S, Vlasblom J, Emili A, Greenblatt J, Wodak SJ: Identifying functional modules in the physical interactome of Saccharomyces cerevisiae. Proteomics. 2007, 7: 944-960. 10.1002/pmic.200600636.View 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.