- Research article
- Open Access
- Published:

# Topology and weights in a protein domain interaction network – a novel way to predict protein interactions

*BMC Genomics*
**volume 7**, Article number: 122 (2006)

## Abstract

### Background

While the analysis of unweighted biological webs as diverse as genetic, protein and metabolic networks allowed spectacular insights in the inner workings of a cell, biological networks are not only determined by their static grid of links. In fact, we expect that the heterogeneity in the utilization of connections has a major impact on the organization of cellular activities as well.

### Results

We consider a web of interactions between protein domains of the Protein Family database (PFAM), which are weighted by a probability score. We apply metrics that combine the static layout and the weights of the underlying interactions. We observe that unweighted measures as well as their weighted counterparts largely share the same trends in the underlying domain interaction network. However, we only find weak signals that weights and the static grid of interactions are connected entities. Therefore assuming that a protein interaction is governed by a single domain interaction, we observe strong and significant correlations of the highest scoring domain interaction and the confidence of protein interactions in the underlying interactions of yeast and fly.

Modeling an interaction between proteins if we find a high scoring protein domain interaction we obtain 1, 428 protein interactions among 361 proteins in the human malaria parasite *Plasmodium falciparum*. Assessing their quality by a logistic regression method we observe that increasing confidence of predicted interactions is accompanied by high scoring domain interactions and elevated levels of functional similarity and evolutionary conservation.

### Conclusion

Our results indicate that probability scores are randomly distributed, allowing to treat static grid and weights of domain interactions as separate entities. In particular, these finding confirms earlier observations that a protein interaction is a matter of a single interaction event on domain level. As an immediate application, we show a simple way to predict potential protein interactions by utilizing expectation scores of single domain interactions.

## Background

The depiction of interactions between genes, proteins and metabolites as networks has uncovered unexpected similarities in the organization of various biological networks, indicating that generic principles and mechanics give rise to their structure. Although such networks vary extensively in their complexity, corroborative evidence points to a series of simple organizing principles that characterize all complex networks. The most dramatic is the scale-free nature of these networks, a remarkable inhomogeneity that highlights a small number of highly connected nodes which secure the networks integrity [1]. The special role such proteins play for the stability of protein interaction networks is further indicated by their significant propensity to be simultaneously essential as well as evolutionary conserved [2]. Reflecting their inherent cohesive nature, complex networks are characterized by the accumulation of discernible modules. Such clusters of densely interconnected nodes combine in an overlapping manner, share well defined functions and hubs as the modules connectors [1, 3, 4]. Similarly to hubs, cohesively bound motifs of protein networks are frequently conserved as a whole, suggesting their role as evolutionary relevant units [5]. While these findings allowed spectacular insights into the inner workings of a cell, biological networks are generally not only determined by their layout of links. In fact, we expect that the heterogeneity in the utilization of connections has a major impact on the organization of cellular activities as well. Recently, attention turned to weighted scientific collaborations and airways networks [6], allowing a first insight into the intricate interplay between links and their weights. Concluding, analysis of real world networks indicate that the static grid of links and their weights can not be regarded as separate entities. Here, we present a first statistical analysis of a weighted biological network by considering a web of PFAM domain interactions. Each link between domains is weighted by an expectation score, reflecting the probability that a particular domain interaction indeed gives rise to observed protein interactions. Applying metrics that combine the static layout of interactions and their weights, we observe that the patterns of correlations are similar for weighted and unweighted network parameters. In contrast to other real world networks, we find weak signals that do not support an entanglement of static grid and weights of domain interactions, allowing us to confirm that a protein interactions are largely governed by single domain interactions.

Assuming that pairs of interacting proteins in *S. cerevisiae* and *D. melanogaster* are indeed dominated by the highest scoring domain interaction their domain architectures suggest, we find that the confidence score of a protein interaction correlates well with its highest scoring domain interaction. As an application, this observation indicates a simple method to model interactions between proteins of the human malaria parasite *P. falciparum*. Assuming an interaction between proteins if we find at least one high scoring domain interaction we predict 1, 428 novel protein interactions among 321 proteins. The quality of each predicted interaction is assessed by a logistic regression model, allowing us to uncover reliable interactions between proteins that share similar functions and are preferably conserved in evolution.

## Results

As a source of high quality interaction data of protein domains we utilized the results of a recent study by Riley et al. [7]. In this statistical approach, called domain pair exclusion analysis (DPEA), a likelihood ratio test is applied to assess the contribution of each potential PFAM-A and PFAM-B domain [8] interaction to the likelihood of a set of observed protein interactions as of DIP [9]. Applying a statistical framework which evaluates the confidence that domains *i* and *j* indeed interact, the authors obtain a network of 1, 566 domains that are embedded in a web of 2, 767 interactions. Weighting each interaction by its probability score – the expectation value [7] – we are primarily interested in the interplay between topology and the reliability of the underlying interactions.

Allowing a first insight in the weights role, we observe a heavy tail in the cumulative distribution of the expectation value of domain links *E*, which can be roughly approximated by a power-law (*P*(*E*) *~ E*^{-2.7}) (Figure 1a). In real world networks the correlation of the degree product *k*_{
i
}*k*_{
j
}and the weight *w*_{
ij
}follows a power-law shaped curve, potentially indicating an intricate relationship between the static layout and weights of links. In our case, we hardly find such a dependence (Figure 1a, inset). In fact, the mean expectation value is almost constant for more than two decades, indicating a general lack of correlation between weights and the domains number of interaction partners [6]

Investigating further if the topology of the underlying domain interaction network and their weights are indeed independent from each other, we combine both topology and weights by a series of measures that enable a more significant assessment of the impact of weights [6]. In an unweighted domain interaction network, the domains degree is defined as *k*_{
i
}= ∑_{
j
}*a*_{
ij
}where *a*_{
ij
}= 1 if there exists a link between domains *i* and *j*. Extending this definition, the strength of a domain *i* is defined as

{s}_{i}={\displaystyle \sum _{j}{a}_{ij}{E}_{ij},}\left(1\right)

accounting for individual expectation values *E*_{
ij
}as weights of interactions of domain *i*. Comparing the statistical properties of a domains degree *k* and its strength *s* we observe that their frequency distributions follow a generalized Zipf's law *P*(*x*) = *α* × (*β* + *x*)^{-γ}(Figure 1b) [10]. The power-law tail of the degree distribution indicates the presence of scale-free topology [11], suggesting that the integrity of the underlying network basically depends on a small subset of highly connected nodes. Analogously, there exists a majority of nodes having low strength while a minority of nodes reach high levels of strength. A list of highest interacting domains shows prominent protagonists that are responsible for important cellular functions such as signaling and cell-cell contacts (Table 1). In particular, we observe that highly connected domains such as pkinase, rrm1 or Zinc finger C2H2 also pool a lot of strength, indicating a proportionality between high level of interactions and their strength.

Investigating the local cohesiveness of network areas, the unweighted representation of the clustering coeffcient *C*_{
i
}measures the degree of cohesiveness around a particular domain *i* [12]. The dependence of the average clustering coeffcient *C* from the domains degree *k* recovers further information about the structure of the underlying network. In most real world networks *C*(*k*) exhibits a highly nontrivial behavior as exemplified by a power-law decay with increasing degree *k*. Averaging over the clustering coeffcients of domains with a certain degree *k*, we find this particular signature, indicating the presence of a nested hierarchy of modules [1] (Figure 1c). Accounting for weights, Barrat *et al*. [6] extended the initial definition of the clustering coeffcient to combine topological information with weights of network links. Considering the expectation value of each domain interactions *E* as the weight of links, we define the weighted clustering coeffcient as

{C}_{i}^{w}=\frac{1}{{s}_{i}({k}_{i}-1)}{\displaystyle \sum _{j,h}\frac{{E}_{ij}+{E}_{ih}}{2}}{a}_{ij}{a}_{ih}{a}_{jh}.\left(2\right)

Since the structure essentially follows the concept of the original clustering coeffcient, we expect that {C}_{i}^{w} retains its dependence from the degree *k*. Indeed, we find a power-law dependence in both networks (Figure 1c). Considering the mean weighted clustering coeffcient of the whole network as the arithmetic mean over all domains *N*, , we obtain 0.097. Comparing this result to the value of the mean unweighted clustering coeffcient of 0.093, we find that ⟨*C*^{w}⟩/⟨*C*⟩ ≈ 1.0. Since the weighted clustering coeffcient reflects a domain's neighborhood to be connected to domains of similar strength the latter result indicates that local clustering predominately occurs on the level of comparable strength.

Another measure that allows insights in the relationship of network layout and weights are degree-degree correlations. Similarly to *C*^{w}, we define the average weighted nearest-neighbors degree as [6]

In an unweighted network the definition of *k*_{nn,i}recovers the average nearest neighbor degree of a node, where {k}_{nn,i}=\frac{1}{{k}_{i}}{\displaystyle {\sum}_{j}{a}_{ij}{k}_{j}}. In the presence of correlations with connectivity *k*, the behavior of the latter measure *k*_{nn,i}(*k*) identifies two classes of networks. If *k*_{
nn
}(*k*) is an increasing function with *k*, vertices with higher degree have an increased probability to be connected with large-degree vertices, a feature that is known as assortative mixing. If *k*_{
nn
}(*k*) decreases with *k*, the underlying network is disassortative, indicating that high degree vertices predominantly are connected to sparsely linked ones. Similarly to other biological networks [13], we find a weak albeit significant trend toward disassortativity in both the unweighted and weighted domain interaction networks (Figure 1d). Considering the nearest neighbor degree of the whole network as the arithmetic mean over all nodes *N*, \u3008{k}_{nn}^{w}\u3009=\frac{1}{N}{\displaystyle {\sum}_{i=1}^{N}{k}_{nn,i}^{w}}, we obtain 12.81. Comparing this result to the value of the mean unweighted nearest neighbor degree of 12.84, we find that 10216;{k}_{nn}^{w}10217;/κ_nnˆw≈ 1.0, indeed confirming that in both the weighted as well as unweighted representation the disassortative behavior prevails.

The previously introduced topological measures of both unweighted and weighted representations of the same domain interaction network share the same qualitative features, suggesting that weights and topology are entangled entities. However, recalling the observation that the degree product does not correlate with the links underlying weights casts doubt on this assumption. Further insights into a potential interplay of topology and utilization of domain interactions arise from correlations between a domains degree and strength (Figure 2a). Despite the existence of inevitable fluctuations, the dependence of the strength from the degree of a domain in the underlying domain interaction network shows a clear and significant power-law *s*(*k*) *~ k*^{β}with ^{β}= 1.04, allowing us to conclude that topology and utilization of links in domain interaction networks are separate entities since independent weights and connectivities would lead to an exponent ^{β}= 1 [6]. We receive further support of this hypothesis by the disparity value *Y*_{2}, a measure that quantifies biased distributions, defined as

{Y}_{2}(i)={\displaystyle \sum _{j\in \Gamma (i)}\frac{{E}_{ij}^{2}}{{s}_{i}^{2}}}\left(4\right)

where Γ(*i*) is the set of neighbors of domain *i*. In Figure 2b we observe a clear power-law in the dependence of the disparity value *Y*_{2} from the degree *k*, *Y*_{2}(*k*) ~ *k*^{-0.9}. Similarly to the dependence of the strength from the degree (Figure 2a), an exponent close to 1 suggests that the expectation values of domain interactions are distributed in an uncorrelated manner [6, 14].

The absence of any correlations between the structure of the web of domain interactions and their confidence suggests that domain interactions hardly interfere with each other. As a consequence, protein interactions are primarily governed by a single domain interaction. Indeed, a recent survey of protein interactions uncovered a rate of 94% that protein interactions are determined by a single pairwise domain interaction [15] while protein interactions that involve interactions between two or more domains are hardly found. A high *E* reflects the probability that the domains in question indeed interact while a low *E*_{
ij
}suggests that other potential domain interactions are roughly as good at explaining the observed protein interactions [7]. Therefore, we assume that a protein interaction is governed by the domain interaction with the highest expectation value. In order to uncover a potential correlation between the quality of a particular protein interaction and the highest scoring domain interaction, we utilize two well curated sets of protein interactions in *S. cerevisiae* [16] and *D. melanogaster* [17] where each interaction is evaluated by a confidence score. Utilizing information about the domain composition of proteins as of the Integr8 database, we screen each domain pair that is suggested by the domain architectures of the underlying proteins. Provided these pairs indeed map to high scoring domain interactions, each protein interaction is assumed to be governed by the domain interaction with highest expectation score. Applied to the evaluated protein interaction sets of *S. cerevisae* and *D. melanogaster*, we observe a strong and significant correlation between an interactions confidence and the expectation value of the underlying highest scoring domain interaction (Figure 4a). In turn, we can potentially use the previous conclusion that the absence of correlations between interactions and their probability indicates the dominance of single domain interactions as a means to infer protein interactions. As an organism, we chose the human malaria parasite *P. falciparum*. Utilizing domain information from the Integr8 database we annotate Plasmodium proteins with their corresponding PFAM domains. In order to avoid interactions between proteins that appear in different compartments we additionally assign each protein with its cellular component terms as of the GO Slim database [18]. Considering all protein pairs of Plasmodium we select those that share at least one GO Slim term. The domain architectures of candidate protein pairs are screened for domain pairs that have at least one high scoring domain interaction. In case we find more than one high scoring domain interaction, we choose the highest scoring one, according to the statistical argument that domain interactions with higher expectation score have a better chance to explain the underlying protein interaction. In Figure 3a, we give a schematic survey of the procedure. Applying this method to the proteome of *P. falciparum* we find 1, 428 interactions between 361 proteins [see Additional file 1]. In order to evaluate each of these potential protein interactions, we characterize each link by measures that reflect biological significance. In particular, we are interested in parameters that are independent of the initial assumption that the highest scoring domain interaction indeed can be used to predict protein interactions. As such, we choose co-expression correlation values of interacting proteins, since similar expression profiles tend to indicate interacting proteins. For *P. falciparum*, we utilized gene expression data over 48 time points. Compiling gene expression data derived from micro-array analysis [19–21], we determine Pearson's correlation coeffcients *r*_{
P
}of each interaction (see Materials and Methods). In addition, we calculated hypergeometric clustering coeffcients *C*_{
vw
}for each interaction, a topological measure that reflects local cohesiveness around a certain link and strongly correlates with the quality of the underlying protein interaction [22] (see Material and Methods). Combining these measures, we utilized a logistic regression method (see Material and Methods) trained by carefully selected sets of 213 true positive and 173 negative interactions, allowing us to assess the quality of each interaction by a confidence score between 0 and 1 (Figure 3b). As a quality measure of the utilized training sets, we performed a leave-one-out strategy, allowing us to obtain 95% accuracy.

Binning interactions according to their confidence value, we observe that about half of the interactions have an elevated degree of confidence (Figure 4b). In each bin, we averaged the expectation score of the domain interactions and observe that high quality of protein interactions – as exemplified by high confidence – are strongly linked to high expectation scores of the underlying domain interaction (Figure 5a). Supported by significant correlation values, this observation is a confirmation of our original assumption that protein interactions are dominated by the highest scoring domain interactions, while high scoring domain interactions indicate the presence of a potential protein interaction. As additional measures of quality, we make use of the well known fact that protein interactions occur between proteins of similar function [23]. As a measure of functional homogeneity of interacting proteins, we apply a hypergeometric test (see Materials and Methods) of the distributions of the proteins GO terms [18]. In particular, this statistical measure reflects the probability that GO terms of interacting proteins have been distributed randomly. Averaging over all interaction specific values in each bin, we find a strong and significant correlation, confirming that protein interactions of increasing confidence tend to occur between functionally related proteins (Figure 5b). As a final test, we wondered if the predicted protein interactions in *P. falciparum* have an evolutionary signature. In particular, we utilized three protein interaction sets of the organisms *S. cerevisiae* [16], *D. melanogaster* [17] and *H. sapiens* [23, 24]. Utilizing orthologous protein information from the InParanoid database [25], we sampled all protein interactions in each organisms that have a fully conserved counterpart – an interolog [26] – in the predicted set of interactions of *P. falciparum*. In Figure 5c, we observe that especially predictions with high confidence pool most of the found interologs in each organism, strongly indicating the reliability of our predictions.

We compared the predicted sets of interactions to a recently published set of experimentally determined protein interactions of *P. falciparum* [27]. Although many interactions of this set have been assigned potential protein domain interactions, the utilized domain information does not overlap strongly with PFAM, restricting the overlap with our predicted set to only 2 interactions. In particular, we find self interactions of the hypothetical Plasmodium proteins PFL0275w and PF10_0232. In the first case, a self interaction of the FHA domain gives rise to the observed protein interaction, while a self interaction of chromo domain determines the latter one. In both cases, the interacting proteins are hypothetical, meaning that their function is unclear. However, the fact that we found domain interactions suggests a role for these proteins. In particular, the forkhead-associated FHA domain is a phosphopeptide recognition domain found in many regulatory proteins, while the chromo (CHRromatin Organization MOdifier) domain is a conserved region of around 60 amino acids involved in the alteration of the structure of chromatin. Putatively, PFL0275w is involved in regulatory activities while PF10_0232 might play a role in chromatin remodeling. In general, our predictions show a prevalence of functions revolving around the proteasome, spliceosome and ribosome. In particular, Table 2 ranks the domain interactions that gave rise to the highest number of predictions in *P. falciparum*. In particular, we observe that domain interactions between the RNA recognition motif rrm1, proteasome and LSM domains appear among the most prevalent domain interactions. As the previous examples illustrates, many interactions are related to self interactions of the underlying domains. As such, we observe a total of 154 self interactions. Indeed, it is well known that multi-protein complexes contain homo-dimers including proteasome [28], ribosome [29] and spliceosome [30]. In particular, rrm's are found in a variety of RNA binding proteins, including various hnRNP proteins, proteins implicated in regulation of alternative splicing, and protein components of snRNPs. The LSM domain contains Sm proteins as well as other related LSM (Like Sm) proteins. The U1, U2, U4/U6, and U5 small nuclear ribonucleoprotein particles (snRNPs) involved in pre-mRNA splicing contain seven Sm proteins in common, which assemble around the Sm site present in four of the major spliceosomal small nuclear RNAs. The U6 snRNP binds to the LSM (Like Sm) proteins. The proteasome is a multicatalytic proteinase complex that is involved in an ATP/ubiquitin-dependent proteolytic pathway. In eukaryotes, the proteasome is composed of about 28 distinct subunits, which form a highly ordered ring-shaped structure (20S ring). Concluding, in the proteasome, ribosome and spliceosome proteins which carry those domains tend to shape stable structures which are mostly governed by self domain interactions, validating the presence of self interactions in our predictions.

## Discussion & conclusion

Assessing the statistical characteristics of a weighted domain interaction network we show that the confidence as exemplified by the expectation value of domain interactions is far from being evenly distributed. Characterizing the underlying weighted domain interactions network, we observe that weighted and unweighted measures of topology follow the same trends. Despite these observations we do not find any significant proof that topology and weights in the domain interaction network are necessarily dependent from each other. In fact, correlations between strength and connectivity as well as disparity suggest that weights as exemplified by the expectation score of each domain interaction are randomly distributed, allowing us to (i) treat the static layout of links and their weights as separate entities and (ii) conclude that protein interactions are indeed governed by a single protein domain interaction [15].

The presence of highly reliable domain interactions offers potential new ways for the prediction and evaluation of protein interactions. In particular, we observe a correlation between an elevated confidence level of a protein interaction in yeast and fly and an increase in the reliability of the underlying domain interactions. As an application, we propose a novel method for the inference of potential protein interactions. While this method can be applied to the prediction of protein interactions in any organism for which PFAM annotation of the organisms proteome is available, we chose the human malaria parasite *P. falciparum*. Screening through all pairs of proteins that provide at least one high scoring domain interaction, we sample potential candidates. Here, we stress that the determination of a high scoring domain interaction has been used as a preselection step of potential protein interaction candidates. In order to evaluate each interaction we resort to interaction specific parameters that are independent from the underlying domain interactions. We find interactions between proteins, that not only show an elevated degree of functional similarity and evolutionary conservation, but also validate our assumption that high scoring domain interactions indeed give rise to reliable interactions. Predominately, we find an enrichment of protein interactions caused by domain interactions that represent functions in the ribosome, proteasome and spliceosome. As reported in protein complexes in other eukaryotes, these functions emphasize a considerable amount of self interactions, we also find in our predictions.

Comparing with existing experimental data sets, we only find a minimal overlap, caused by the fact that many proteins of *P. falciparum* currently are not annotated with PFAM domains. On the other hand, experimental determination of protein interactions in *P. falciparum* is in its starting phase covering about a quarter of known proteins. As such, our predictions can help focus experimental studies on specific interactions unique to this pathogen.

## Methods

### Domain-domain interactions

As a source of high quality interaction data of protein domains we utilized the results of a recent study by Riley et al. [7]. In this statistical approach called domain pair exclusion analysis (DPEA), a likelihood ratio test is applied to assess the contribution of each potential PFAM-A and PFAM-B domain [8] interaction to the likelihood of a set of observed protein interactions. DPEA consists of three steps: (i) Utilizing protein interaction data from DIP [9], the frequency *S*_{
ij
}of an interaction between *i* and *j* in relation to their abundance in the data is computed. (ii) Using *S*_{
ij
}as an initial guess, an expectation maximization algorithm is applied to obtain a maximum likelihood estimate of Θ_{
ij
}which stands for the probability of domain interaction *ij* among all the possible domain interactions which are suggested by the domain architectures of the interacting protein pairs where domain *i* and *j* co-occur. In a third step, all possible interactions of domains *i* and *j* are excluded from the mixture of competing hypotheses for the presence of corresponding protein interactions, EM is rerun, and the change in likelihood is expressed as a log odds score, *E*_{
ij
}, reflecting the confidence that domains *i* and *j* indeed interact. As such, a high value of *E*_{
ij
}indicates that there is extensive evidence in protein interaction data that domains *i* and *j* interact while a low *E*_{
ij
}suggests that other potential domain interactions are roughly as good at explaining the observed protein interactions [7]. As a proof of concept, domain pairs inferred to interact with high *E* are significantly enriched among domain pairs known to interact in the Protein Data Bank (PDB). The domain interaction network thus obtained comprises 1, 566 domains which are embedded in 2, 767 interactions that score *E*_{
ij
}≥ 3.

### Protein interactions

We utilized a large scale compilation of yeast protein interactions. In particular, this data set combines 47, 783 experimentally obtained protein interactions among 4, 175 proteins in *S. cerevisiae* [16] obtained from sources as diverse as mRNA expression studies and yeast2hybrid screens. Each interaction was characterized by a confidence score obtained by the application of a logistic regression model. Analogously, the quality of experimentally protein interactions in *D. melanogaster* was assessed, allowing for 6, 222 proteins and 16, 914 links [17]. As for direct experimental observations of protein interactions in *P. falciparum*, we utilized a set of 2, 475 interactions among 1, 304 proteins that have been obtained by the modification of a yeast2hybrid method [27]. Additionally, we utilized a large-scale compilation of human interactions totaling 89, 572 interactions among 9, 018 proteins [23, 24].

### Protein domain data

The advent of fully sequenced genomes of various organisms has facilitated the investigation of proteomes. The Integr8 database has been set up to provide comprehensive statistical and comparative analyzes of complete proteomes of fully sequenced organisms. The initial version of the application contained data for genomes and proteomes of 182 sequenced organisms (including 19 archae, 150 bacteria and 13 eukaryotes) and proteome analyzes derived through the integration of UniProt [31], InterPro [32], CluSTr [33], GO/GOA [34], EMSD, Genome Reviews and IPI [35]. In particular, we utilized IPI (International Protein Index) files to elucidate the domain architecture of the corresponding proteins. For our analysis, we focused on domain data retrieved from the PFAM database, a reliable collection of multiple sequence alignments of protein families and profile hidden Markov models [36].

### Orthologous protein data

The InParanoid database [25] provides putative orthologous sequence information for the complete proteomes of organism pairs *S. cerevisiae*, *D. melanogaster*, *H. sapiens* and *P. falciparum*. The algorithm for detecting orthologous relationships is based on pairwise similarity scores which are by default calculated with the BLASTP program. InParanoid detects mutual best hits between sequences from two different species. These are two main orthologs that form an orthologous group. Other sequences are added to this group if they are closely related to one of the main orthologs. These members of the orthologous group are called in-paralogs. A confidence value provided by a standard bootstrap procedure for each in-paralog shows how closely related it is to the main ortholog. In our study, we only selected the main sequence pairs of each orthologous group allowing us to obtain 2, 319 yeast proteins, 1, 351 in *D. melanogaster* and 1, 525 in *H. sapiens* with putative orthologs in *P. falciparum*.

### Co-expression data

Genes with similar expression profiles are likely encoding interacting proteins. For *P. falciparum*, we utilized gene expression data, compiling 5, 156 genes over 48 time points as of Winzeler et al., [19, 21] and of Bozdech et al. collecting 4, 318 genes over 48 time points [37]. As a gene similarity metric we calculated Pearson's correlation coeffcient for every protein interaction over *m* time points defined as

{r}_{p}=\frac{\frac{1}{m}{\displaystyle {\sum}_{i=1}^{m}{x}_{i}{y}_{i}-\u3008x\u3009\u3008y\u3009}}{{\sigma}_{i}{\sigma}_{j}}\left(5\right)

where ⟨*x*⟩ and ⟨*y*⟩ are the sample means of expression values *x*_{
i
}and *x*_{
j
}, and *σ*_{
i
}and *σ*_{
j
}are their standard deviations.

### Logistic regression

In order to get an estimate of an interactions reliability, we employed a logistic regression model. According to the logistic regression, the probability of a true interaction *T*_{
vw
}given the two input variables, hypergeometric clustering coeffcient *x*_{1} = *C*_{
vw
}and co-expression correlation coeffcient *x*_{2} = *r*_{
P
}, *X* = (*x*_{1}*, x*_{2})

Pr({T}_{vw}|X)=\frac{exp({\beta}_{0}+{\beta}_{1}{x}_{1}+{\beta}_{2}{x}_{2})}{1+exp({\beta}_{0}+{\beta}_{1}{x}_{1}+{\beta}_{2}{x}_{2})}\left(6\right)

where *β*_{
n
}are the parameters of the distribution. Given training data we optimized the distribution parameters by maximizing the likelihood of the data. Here, we applied the corresponding routines as of the Biopython package [38]. As a training set for true positives we choose 213 high scoring protein-interactions in yeast [16] that are fully conserved in Plasmodium. In the same way, we selected 173 low scoring interactions as true negative training set. Applying a leave-one-out analysis to determine the prediction accuracy, our model is recalculated from the training data after removing the interaction to be predicted (leave-one-out), allowing us to obtain the correct result in 95% of cases.

### Hypergeometric clustering coeffcient

Recently, a network topology based approach uncovered a remarkable correlation between enhanced quality of protein interactions and the degree of clustering of their immediate network neighborhood [22]. Considering a network with *N* nodes, we define the hypergeometric clustering coeffcient as

{C}_{vw}=-\mathrm{log}{\displaystyle \sum _{i=|N(v)\cap N(w)|}^{\mathrm{min}(|N(v)|,|N(w)|)}\frac{\left(\begin{array}{c}|N(v)|\\ i\end{array}\right)\left(\begin{array}{c}N-|N(v)|\\ |N(w)|-i\end{array}\right)}{\left(\begin{array}{c}N\\ |N(w)|\end{array}\right)}\left(7\right)}

where *N*(*x*) represents the neighborhood of a vertex *x*. Given fixed neighborhood sizes *N*(*v*) and *N*(*w*) of nodes *v* and *w*, the hypergeometric clustering coeffcient increases with elevated overlap between the nodes neighborhoods. Provided that the neighborhoods are independent, the summation can be interpreted as a *p* value, reflecting the probability of obtaining a number of mutual neighbors between nodes *v* and *w* at or above the observed number by chance.

### GO annotation data and functional homogeneity

Similarly to the hypergeometric clustering coeffcient, we define the functional homogeneity of a domain pair *ij*

f{h}_{ij}=-\mathrm{log}{\displaystyle \sum _{i=|GO(v)\cap GO(w)|}^{\mathrm{min}(|GO(v)|,|GO(w)|)}\frac{\left(\begin{array}{c}|GO(v)|\\ i\end{array}\right)\left(\begin{array}{c}T-|GO(v)|\\ |GO(w)|-i\end{array}\right)}{\left(\begin{array}{c}T\\ |GO(w)|\end{array}\right)}\left(8\right)}

where *GO*(*i*) is the set of GO Terms of protein *i*, and *T* is the total number of different GO terms [18]. In analogy, the summation can be interpreted as a *p* value, reflecting the probability that a protein pair shares a certain number of GO terms at or above the observed number by chance.

## References

Barabaśi A, Oltvai Z: Network Biology: Understanding the Cell's Functional Organization. Nat Rev Genet. 2004, 101-113. 10.1038/nrg1272. 5

Wuchty S: Topology and Evolution in Yeast Interaction Networks. Genome Res. 2004, 14: 1310-1314. 10.1101/gr.2300204.

Han J, Bertin N, Hao T, Goldberg DS, Berriz G, Zhang L, Dupuy D, Walhout A, Cusick M, Roth F, Vidal M: Evidence for dynamically organized modularity in the yeast protein-protein interaction network. Nature. 2004, 430: 88-93. 10.1038/nature02555.

Guimera R, Amaral L: Functional cartography of complex metabolic networks. Nature. 2005, 433: 895-900. 10.1038/nature03288.

Wuchty S, Oltvai Z, Barabási AL: Evolutionary conservation of motif constituents within the yeast protein interaction network. Nature Genetics. 2003, 35: 176-179. 10.1038/ng1242.

Barrat A, Barthélemy M, Pastor-Satorras R, Vespignani A: The architecture of complex weighted networks. Proc Natl Acad Sci USA. 2004, 101 (11): 3747-3752. 10.1073/pnas.0400087101.

Riley R, C Lee CS, Eisenberg D: Inferring protein domain interactions from databases of interacting proteins. Genome Biol. 2005, 6 (10): R89-10.1186/gb-2005-6-10-r89.

Bateman A, Coin L, Durbin R, Finn R, Hollich V, Griffths-Jones S, Khanna A, Marshall M, Moxon S, Sonnhammer E, Studholme D, Yeats C, Eddy S: The Pfam protein families database. Nucl Acids Res. 2004, 32: D138-D141. 10.1093/nar/gkh121.

Xenarios I, Salwinski L, Duan X, Higney P, Kim SM, Eisenberg D: DIP, the Database of Interacting Proteins: a research tool for studying cellular networks of protein interactions. Nucl Acids Res. 2002, 30: 303-305. 10.1093/nar/30.1.303.

Park J, Lappe M, Teichmann S: Mapping Protein Family Interactions: Intramolecular and Intermolecular Protein Family Interaction Repertoires in the PDB and Yeast. J Mol Biol. 2001, 307: 929-938. 10.1006/jmbi.2001.4526.

Albert R, Barabási AL: Statistical mechanics of complex networks. Rev Mod Phys. 2002, 74: 47-10.1103/RevModPhys.74.47.

Watts D, Strogatz S: Collective dynamics of 'small-world' networks. Nature. 1998, 393: 440-442. 10.1038/30918.

Newman M: Assortative mixing in networks. Phys Rev Lett. 2002, 89: 208701-10.1103/PhysRevLett.89.208701.

Barthelemy M, Gondran B, Guichard E: Spatial structure of the internet traffc. Physica A. 2003, 319: 633-642. 10.1016/S0378-4371(02)01382-1.

Aloy P, Böttcher B, Ceulemans H, Leutwein C, Mellwig C, Fischer S, Gavin AC, Bork P, Superti-Furga G, Serrano L, Russell R: Structure-Based Assembly of Protein Complexes in Yeast. Science. 2004, 303: 2026-2029. 10.1126/science.1092645.

Bader JS, Chaudhuri JRD, Chant J: Gaining confidence in high-throughput protein interaction networks. Nat Biotechnol. 2004, 22: 78-85. 10.1038/nbt924.

Giot L, Bader J, Brouwer C, Chaudhuri A, Kuang B, Li Y, Hao Y, Ooi C, Godwin B, Vitols E, Vijayadamodar G, Pochart P, Machineni H, Welsh M, Kong Y, Zerhusen B, Malcolm R, Varrone Z, Collis A, Minto M, Burgess S, McDaniel L, Stimpson E, Spriggs F, Williams J, Neurath K, Ioime N, Agee M, Voss E, Furtak K, Renzulli R, Aanensen N, Carrolla S, Bickelhaupt E, Lazovatsky Y, DaSilva A, Zhong J, Stanyon C, Finley R, White K, Braverman M, Jarvie T, Gold S, Leach M, Knight J, Shimkets R, McKenna M, Chant J, Rothberg J: A Protein Interaction Map of Drosophila melanogaster. Science. 2003, 302: 1727-1736. 10.1126/science.1090289.

Consortium GO: The Gene Ontology (GO) database and informatics resource. Nucl Acids Res. 2004, 32: D258-D261. 10.1093/nar/gkh036.

Le Roch K, Zhou Y, Blair P, Grainger M, Moch J, Haynes J, De la Vega P, Holder A, Batalov S, Carucci D, Winzeler E: Discovery of Gene Function by Expression Profiling of the Malaria Parasite Life Cycle. Science. 2003, 301: 1503-1508. 10.1126/science.1087025.

Johnson KLJ, Florens L, Zhou Y, Santrosyan A, Grainger M, Yan S, Williamson K, Holder A, Carucci D, Yates III J, Winzeler E: Global analysis of transcript and protein levels across the Plasmodium falciparum life cycle. Genome Res. 2004, 14: 2308-2318. 10.1101/gr.2523904.

Winzeler E: Applied systems biology and malaria. Nat Rev Microbiol. 2006, 4: 145-151. 10.1038/nrmicro1327.

Goldberg D, Roth F: Assessing experimentally derived interactions in a small world. Proc Natl Acad Sci USA. 2003, 100: 4372-4376. 10.1073/pnas.0735871100.

Lehner B, Fraser A: A first-draft human protein-interaction map. Genome Biol. 2004, 5 (9): R63-10.1186/gb-2004-5-9-r63.

Ramani A, Bunescu R, Mooney R, Marcotte E: Consolidating the set of known human protein-protein interactions in preparation for large-scale mapping of the human interactome. Genome Biol. 2005, 6 (5): R40-10.1186/gb-2005-6-5-r40.

Remm M, Storm C, Sonnhammer E: Automatic clustering of orthologs and in-paralogs from pairwise species comparisons. J Mol Biol. 2001, 314: 1041-1052. 10.1006/jmbi.2000.5197.

Ge H, Liu Z, Church G, Vidal M: Correlation between transcriptome and interactome mapping data from Saccharomyces cerevisiae. Nature Genetics. 2001, 29: 482-486. 10.1038/ng776.

LaCount D, Vignali M, Chettier R, Phansalkar A, Bell R, Hesselberth J, Schoenfeld L, I Ota SS, Kurschner C, Fields S, Hughes R: A protein interaction network of the malaria parasite Plasmodium falciparum. Nature. 2005, 438: 103-107. 10.1038/nature04104.

Bochtler M, Ditzel L, Groll M, Hartmann C, Huber R: The Proteasome. Annu Rev Biophys Biomol Struct. 1999, 28: 295-317. 10.1146/annurev.biophys.28.1.295.

Matadeen R, Patwardhan A, Gowen B, Orlova E, Pape T, Cuf M, Mueller F, Brimacombe R, van Heel M: The Escherichia coli large ribosomal subunit at 7.5A resolution. Structure Fold Des. 1999, 7: 1575-1583. 10.1016/S0969-2126(00)88348-3.

Mura C, Cascio D, Sawaya M, Eisenberg D: The crystal structure of a heptameric archaeal SM protein: implications for the eukaryotic snRNP core. Proc Natl Acad Sci USA. 2001, 98: 5532-5537. 10.1073/pnas.091102298.

Apweiler R, Bairoch A, Wu C, Barker W, Boeckmann B, Ferro S, Gasteiger E, Huang H, Lopez R, Magrane M, Martin M, Natale D, O'Donovan C, Redaschi N, Yeh L: Uniprot: the universal protein knowledgebase. Nucl Acids Res. 2004, 32: D115-D119. 10.1093/nar/gkh131.

Mulder N, Apweiler R, Attwood T, Bairoch A, Barrell D, Bateman A, Binns D, Biswas M, Bradley P, Bork P, Bucher P, Copley R, Courcelle E, Das U, Durbin R, Falquet L, Fleischmann W, Griffths-Jones S, Haft D, Harte N, Hulo N, Kahn D, Kanapin A, Krestyaninova M, Lopez R, Letunic I, Lonsdale D, Silventoinen V, Orchard S, Pagni M, Peyruc D, Ponting C, Selengut J, Servant F, Sigrist C, Vaughan R, Zdobnov E: The InterPro Database, 2003 brings increased coverage and new features. Nucl Acids Res. 2003, 31: 315-318. 10.1093/nar/gkg046.

Kriventseva E, Fleischmann W, Zdobnov E, Apweiler R: CluSTr: a database of clusters of SWISS-PROT+TrEMBL proteins. Nucl Acids Res. 2001, 29: 33-36. 10.1093/nar/29.1.33.

Consortium GO: The gene ontology (go) database and information resource. Nucl Acids Res. 2004, 32: D258-D261. 10.1093/nar/gkh036.

Kersey P, Duarte J, Williams A, Apweiler R, Karavidopoulou Y, Birney E: The international protein index: An integrated database for proteomics experiments. Proteomics. 2004, 4 (7): 1985-1988. 10.1002/pmic.200300721.

Bateman A, Coin L, Durbin R, Finn R, Hollich V, Griffths-Jones S, Khanna A, Marshall M, Moxon S, Sonnhammer E, Studholme D, Yeats C, Eddy S: The Pfam protein families database. Nucl Acids Res. 2004, 32: D138-D141. 10.1093/nar/gkh121.

Bozdech Z, Llinas M, Pulliam B, Wong E, Zhu J, DeRisi J: The Transcriptome of the Intraerythrocytic Developmental Cycle of Plasmodium falciparum. PLoS Biol. 2003, 1: 1-16. 10.1371/journal.pbio.0000005.

The Biopython package. [http://www.biopython.org]

## Acknowledgements

This project was entirely funded by the Northwestern Institute on Complexity (NICO).

## Author information

### Authors and Affiliations

### Corresponding author

## Electronic supplementary material

### 12864_2006_505_MOESM1_ESM.TXT

Additional File 1: **Predicted protein interactions in** ** P. falciparum**. This file contains the full set of predicted protein interactions in

*P. falciparum*. Each column is tab-delineated, presenting the interacting proteins (column 1&2), the underlying domain interactions (columns 3 & 4), their expectation score (column 5), the interactions confidence score (column 6) and the proteins functional annotation (columns 7 & 8) (TXT 173 KB)

## Authors’ original submitted files for images

Below are the links to the authors’ original submitted files for images.

## Rights and permissions

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.

## About this article

### Cite this article

Wuchty, S. Topology and weights in a protein domain interaction network – a novel way to predict protein interactions.
*BMC Genomics* **7, **122 (2006). https://doi.org/10.1186/1471-2164-7-122

Received:

Accepted:

Published:

DOI: https://doi.org/10.1186/1471-2164-7-122

### Keywords

- Protein Interaction
- Domain Interaction
- Real World Network
- Domain Pair
- Human Malaria Parasite