- Report
- Open access
- Published:

# Conditional random field approach to prediction of protein-protein interactions using domain information

*BMC Systems Biology*
**volume 5**, Article number: S8 (2011)

## Abstract

### Background

For understanding cellular systems and biological networks, it is important to analyze functions and interactions of proteins and domains. Many methods for predicting protein-protein interactions have been developed. It is known that mutual information between residues at interacting sites can be higher than that at non-interacting sites. It is based on the thought that amino acid residues at interacting sites have coevolved with those at the corresponding residues in the partner proteins. Several studies have shown that such mutual information is useful for identifying contact residues in interacting proteins.

### Results

We propose novel methods using conditional random fields for predicting protein-protein interactions. We focus on the mutual information between residues, and combine it with conditional random fields. In the methods, protein-protein interactions are modeled using domain-domain interactions. We perform computational experiments using protein-protein interaction datasets for several organisms, and calculate AUC (Area Under ROC Curve) score. The results suggest that our proposed methods with and without mutual information outperform EM (Expectation Maximization) method proposed by Deng et al., which is one of the best predictors based on domain-domain interactions.

### Conclusions

We propose novel methods using conditional random fields with and without mutual information between domains. Our methods based on domain-domain interactions are useful for predicting protein-protein interactions.

## Background

Understanding of protein functions and protein-protein interactions is one of important topics in the field of molecular biology and bioinformatics. Recently, many researchers have focused on the investigation of amino acid residues of proteins to reveal interactions and contacts between residues [1–4]. If residues at important sites for interactions between proteins are substituted in one protein, the corresponding residues in interacting partner proteins are expected to be also substituted by selection pressure. Otherwise, such mutated proteins may lose the interactions. Fraser et al. confirmed that interacting proteins evolve at similar evolutionary rates by comparing putatively orthologous protein sequences between *S. cerevisiae* and *C. elegans*[5]. It means that substitutions for contact residues occur in both interacting proteins as long as the proteins keep interacting with each other. Therefore, mutual information (MI) between residues is useful for predicting protein-protein interactions for proteins of unknown function. MI is calculated from multiple sequence alignments for homologous protein sequences. Weigt et al. identified direct residue contacts between sensor kinase and response regulator proteins by message passing, which is an improvement of MI [4]. Burger and van Nimwegen used a dependence tree where a node corresponds to a position of amino acid sequences, and predicted interactions using a Bayesian network method [2]. On the other hand, Markov random field and conditional random field models have been well studied in fields of natural language processing [6, 7]. Also in bioinformatics, protein function prediction methods from protein-protein interaction network and other biological networks were developed using Markov random fields [8, 9]. On the other hand, several prediction methods have been developed based on domain-domain interactions. Deng et al. proposed a domain-based probabilistic model of protein-protein interactions, and developed EM (Expectation Maximization) method [10]. Based on this probabilistic model, LP (Linear Programming)-based methods were developed [11], and Chen et al. improved the accuracy of interaction strength prediction by APM (Association Probabilistic Method) [12]. In this paper, we propose prediction methods based on domain-domain interactions using conditional random fields with and without mutual information. Furthermore, we perform computational experiments for several protein-protein interaction datasets, compare the methods with the EM method proposed by Deng et al. [10], which is one of the best predictors based on domain-domain interactions, and the association method proposed by Sprinzak and Margalit [13] (the APM method for binary interaction data is equivalent to the association method), and show that our methods outperform the EM method and the association method.

## Mutual information between domains

In order to investigate the relationship between two positions of proteins, MI for distributions of amino acids at the positions is used. Such distributions can be obtained from multiple alignments of protein sequences and domain sequences. In this section, we briefly review MI for distributions of amino acids, and explain MI between domains.

We assume that multiple sequence alignments for domains *D*_{
m
} and *D*_{
n
} are obtained, respectively (see Figure 1). In order to calculate MI, we need joint appearance frequencies. However, we cannot see which sequence in the multiple alignment of domain *D*_{
m
} corresponds to a specified sequence in that of *D*_{
n
}. Therefore, we assume that sequences contained in the same organism can be paired. In the example of Figure 1, the second sequence of *D*_{
m
} is paired with the first one of *D*_{
n
}, the third one of *D*_{
m
} is paired with the second one of *D*_{
n
}, and so on. The first sequence of *D*_{
m
} is not counted into the appearance frequencies because it is not paired with any sequence of *D*_{
n
} although it may be paired with sequences of other domains than *D*_{
n
}.

Let ** A** be a set of amino acids,

*f*

_{ i }(

*A*) be the appearance frequency of amino acid

*A*at position

*i*in domains

*D*

_{ m }and

*D*

_{ n }, and

*f*

_{ ij }(

*A*,

*B*) be the joint appearance frequency of a pair of amino acids

*A*at position

*i*in

*D*

_{ m }and

*B*at position

*j*in

*D*

_{ n }, where each frequency is divided by the number of paired sequences

*M*in the multiple alignments such that ∑

_{A∈A}

*f*

_{ i }(

*A*) = ∑

_{A,B∈A}

*f*

_{ ij }(

*A*,

*B*) = 1.

Multiple alignments often include some gaps. Weigt et al. counted the frequencies of gaps as well as amino acids [4]. Therefore, we also consider gaps to be a kind of amino acids, that is, the number of distinct amino acids is |** A**| = 21. Then, mutual information for positions

*i*in

*D*

_{ m }and

*j*in

*D*

_{ n }is defined as the Kullback-Leibler divergence between the multiplication of appearance frequencies,

*f*

_{ i }(

*A*)

*f*

_{ j }(

*B*), and the joint appearance frequencies,

*f*

_{ ij }(

*A*,

*B*), as follows.

If frequency distributions of amino acids at positions *i* and *j* are independent from each other, *f*_{
ij
}(*A*,*B*) ≈ *f*_{
i
}(*A*)*f*_{
j
}(*B*), and *MI*_{
ij
} approaches to zero. This means that the two positions are not related with each other in the evolutionary process. If domains *D*_{
m
} and *D*_{
n
} interact at the positions, it is considered that *MI*_{
ij
} becomes high because the positions have coevolved through the evolutionary process in order to keep the interaction. It should be noted that two positions *i* and *j* do not always directly interact even if *MI*_{
ij
} is high [4]. However, such proteins with high values of MI have a possibility to directly interact with each other at other positions in the proteins.

However, we need to reduce *MI*_{
ij
} because it can be unnecessarily high depending on distributions of *f*_{
i
}(*A*) and *f*_{
j
}(*B*). For that purpose, we make use of , which is the mutual information *MI*_{
ij
} from the joint frequency, *f*_{
ij
}(*A*, *B*), obtained by shuffling at random the combinations of sequences in multiple alignments. In this paper, we repeat the procedure 400 times according to [4], and take the average. For practical uses of MI, *f*_{
i
}(*A*), *f*_{
j
}(*B*) and *f*_{
ij
}(*A*,*B*) should be positive values. Otherwise, we cannot calculate *MI*_{
ij
} by using computers. Therefore, we use the following pseudocount as in [4],

where *η* is a constant value, in this paper we use *η* = 1. It should be noted that the sum over all amino acids ** A**, and because ∑

_{A∈A}

*f*

_{ i }(

*A*) = ∑

_{A,B∈A}

*f*

_{ ij }(

*A*,

*B*) = 1.

In order to investigate interactions between proteins, we need MI between domains included in the proteins. Thus, we define MI between domains *D*_{
m
} and *D*_{
n
}, *M*_{
mn
}, to be the maximum of MI over all positions *i* and *j* as follows.

where 〈*v*〉 means the average of *v*, *i* and *j* are positions of *D*_{
m
} and *D*_{
n
}, respectively. Since *MI*_{
ij
} is calculated to be high for the positions *i* and *j* that include many gaps, we exclude positions that include more than 20% gaps as in [14].

## Conditional random field model for PPI

In this section, we propose a probabilistic model for protein-protein and domain-domain interactions using conditional random fields [6, 7] because it can be considered that two domains *D*_{
m
} and *D*_{
n
} do not always interact even if the mutual information *M*_{
mn
} is large. For example, Weigt et al. improved MI and proposed direct information (DI) because residues do not always contact with each other even if the MI is large [4]. Most proteins contain domains as is well known. If two proteins do not interact with each other, any two domains contained in the proteins must not interact with each other. In the left example of Figure 2, protein *P*_{
i
} consists of domains *D*_{1} and *D*_{2} and protein *P*_{
j
} consists of domain *D*_{3} respectively. If *P*_{
i
} and *P*_{
j
} do not interact, any pair of (*D*_{1}, *D*_{3}) and (*D*_{3}, *D*_{3}) does not interact. Deng et al. proposed a probabilistic model for a pair of proteins as follows [10]. By assuming that proteins *P*_{
i
} and *P*_{
j
} interact if and only if at least a pair of domains included in the proteins interacts, and events that domains interact are independent from each other, they defined

where *P*_{
ij
} = 1 means that proteins *P*_{
i
} and *P*_{
j
} interact, *D*_{
mn
} = 1 means that domains *D*_{
m
} and *D*_{
n
} interact, *D*_{
mn
} ∈ *P*_{
ij
} means that domain *D*_{
m
} is included in protein *P*_{
i
} and *D*_{
n
} is included in *P*_{
j
} and the product in the right hand side is calculated for all domain pairs (*D*_{
m
}, *D*_{
n
}) included in the protein pair (*P*_{
i
}, *P*_{
j
}). By transforming equation (5), we have

where *λ*^{(mn)} = log(1 – *Pr*(*D*_{
mn
} = 1)).

From this equation, we can consider the following Markov random field model for protein pair (*P*_{
i
}, *P*_{
j
}) (see Figure 2).

where *p*_{
ij
} ∈ {0, 1}, ** d** means a set of events on domain-domain interactions,

*D*

_{ mn }=

*d*

_{ mn }(

*d*

_{ mn }∈ {0, 1}), denotes a local feature, is the corresponding weight parameter and related to the joint probability

*Pr*(

*P*

_{ ij }=

*s*,

*D*

_{ mn }=

*t*), and Z

_{ ij }denotes the normalization constant. For instance, equation (8) for

*p*

_{ ij }= 0 is equivalent to equation (7) in the case that for all protein pairs (

*P*

_{ i },

*P*

_{ j }) and if

*s*=

*t*= 0, otherwise 0.

In Markov random fields, random variables have Markov properties represented as an undirected graph [15]. The factor graph for our model is represented to be a bipartite graph *G*(*U*, *V*, *E*) with a set of vertices *U* corresponding to protein-protein interactions *P*_{
ij
}, a set of vertices *V* corresponding to domain-domain interactions *D*_{
mn
}, and a set of edges *E* between *U* and *V* as the right figure of Figure 2. There exists an edge between *P*_{
ij
} ∈ *U* and *D*_{
mn
} ∈ *V* if and only if *D*_{
mn
} ∈ *P*_{
ij
}. For the left example of Figure 2, protein pair (*P*_{
i
}, *P*_{
j
}) includes domain pairs (*D*_{1}, *D*_{3}) and (*D*_{2}, *D*_{3}). Then, in the factor graph, the vertex of *P*_{
ij
} is connected with vertices of *D*_{13} and *D*_{23}, respectively. Although the vertex of *P*_{
ij
} does not have other adjacent vertices than the vertices of *D*_{13} and *D*_{23}, those of *D*_{13} and *D*_{23} can be connected with other vertices than that of *P*_{
ij
}

Since *Pr*(*P*_{
ij
} = 0|*D*_{
mn
} =*t*) = 1 – *Pr*(*P*_{
ij
} = 1|*D*_{
mn
} = *t*), it is redundant to consider both *s* = 0, 1, and it is sufficient to consider only *s* = 1. Therefore, in order to simplify the model, we substitute , , and for all protein pairs (*P*_{
i
}, *P*_{
j
}). Then, we have the following joint probability,

where ** p** means a set of events on protein-protein interactions,

*P*

_{ ij }=

*p*

_{ ij }.

We here introduce mutual information between domains ** M** = {

*M*_{ mn }} as given conditional data in order to combine it with the probabilistic model. Then, equation (9) can be written as

where

*σ*(*x*) = 1/(1 + *e*^{–x}) is an increasing function, and *c* is a positive constant. It should be noted that a negative value, –1, is given to because it is undesired that a pair of domains interact although proteins having the pair do not interact. In this way, the local feature correlates protein-protein interactions *P*_{
ij
} with domain-domain interactions *D*_{
mn
} (see Figure 2).

For a conditional random field model without MI, we use the following local feature instead of .

### Parameter estimation

In this section, we discuss how to estimate the parameters . We assume that protein-protein interaction data ** p** = {

*p*

_{ ij }} are given. Then, the likelihood function is represented by

where *Z*(** M**) = ∏

_{ p }

_{ij∈p}

*Z*

_{ ij }(

**). By taking the logarithm, we have**

*M*We estimate the parameters by maximizing the log-likelihood function, *l*(** λ**). Since log(

*e*

^{x}+

*e*

^{y}) is a convex function for variables

*x*and

*y*, that is,

*l*(

**) is a concave function, we are able to obtain a global maximum. For maximizing such functions, various methods such as the steepest descent method, Newton’s method, and the Broyden-Fletcher-Goldfarb-Shanno (BFGS) [16] method have been developed. Newton’s method calculates the inverse of the Hessian matrix for the objective function. However, the computational cost is high. Therefore, the quasi-Newton method approximates the matrix by some efficient method using the first derivatives, the gradient. In this paper, we use the BFGS method, which is one of the quasi-Newton methods. By differentiating equation (15) partially with respect to each parameter , we have**

*λ*In the BFGS method, this equation is repeatedly applied for updating a solution.

## Computational experiments

### Data and implementation

We used protein-protein interaction data of *H. sapiens*, *D. melanogaster*, and *C. elegans* from the DIP database [17], the file name is ’dip20091230.txt’. We used the UniProt Knowledgebase database (version 15.4) [18] as protein domain inclusion data. We deleted proteins that did not have any domain, and obtained 294 interacting protein pairs as positive data that included 300 distinct proteins and 320 domains for *H. sapiens*, 449 interacting pairs that included 562 proteins and 449 domains for *D. melanogaster*, and 250 interacting pairs that included 602 proteins and 476 domains for *C. elegans*.

We used the Pfam database (version 24.0) [19] to obtain multiple sequence alignments for domains, and calculated MI, *M*_{
mn
}, for each pair of domains. Figure 3 shows the distributions of domain MI *M*_{
mn
} for *H. sapiens*, *D. melanogaster*, and *C. elegans*. We can see from the figure that most domain MIs are distributed in the part of less than about 0.8 for all organisms. It is considered that domains *D*_{
m
} and *D*_{
n
} with *M*_{
mn
} less than 0.8 may not interact, and domains with *M*_{
mn
} more than 0.8 have more possibilities to interact with each other. Therefore, we set the constant *c* in equation (12) to be 0.8. Although we tried several values from 0.6 to 1.0 for *c*, the results were similar to the case of *c* = 0.8.

We selected non-interacting protein pairs as negative data uniformly at random such that negative data did not overlap with the positive data. The number of negative data was the same as that of positive data for each organism.

We used libLBFGS (version 1.9) with default parameters to estimate the parameters , which is a C implementation of the limited memory BFGS method [20], and is available on the web page, http://www.chokkan.org/software/liblbfgs/.

### Result

In order to evaluate our method, we compared the proposed CRF method with MI and that without MI with the EM method by Deng et al. [10] and the association method proposed by Sprinzak and Margalit [13]. The association method and the APM method [12] estimate probabilities *λ*_{
mn
} that domains *D*_{
m
} and *D*_{
n
} interact as and , respectively, where *N*_{
mn
} (*I*_{
mn
}) denotes the number of (interacting) protein pairs that include domain pair (*D*_{
m
}, *D*_{
n
}), and *ρ*_{
ij
} denotes the interaction strength of protein pair (*P*_{
i
}, *P*_{
j
}), 0 ≤ *ρ*_{
ij
} ≤ 1. However, our input interaction data are binary, that is, *ρ*_{
ij
} takes only 0 or 1. Then, the numerator of the APM method becomes *I*_{
mn
}. It means that the APM method for binary interaction data is equivalent to the association method. In the EM method, probabilities *λ*_{
mn
} that domains *D*_{
m
} and *D*_{
n
} interact are estimated by the recursive formula, , where *o*_{
ij
} = 1 denotes that it was observed that proteins *P*_{
i
} and *P*_{
j
} interact with each other, and *fn* = 0.8. In this paper, the solution of the association method was given as the initial value of the EM method.

We performed five-fold cross-validation, that is, split the data into 5 datasets (4 for training and 1 for test), estimated from the training datasets, and calculated *Pr*(*P*_{
ij
} = 1|** M**) of equation (10) for each protein pair in the test dataset and AUC (Area Under ROC Curve) score, where among the test dataset only protein pairs that included at least a parameter estimated from the corresponding training dataset were always used. We repeated 5 times, and took the average. Tables 1, 2, and 3 show the results on AUC for training and test datasets by the CRF method with MI, that without MI, the EM method, and the association method for

*H. sapiens*,

*D. melanogaster*, and

*C. elegans*, respectively. An AUC score is the area under an ROC (Receiver Operating Characteristic) curve, and takes a value between 0 and 1. The ROC curve of a random classifier lies on the diagonal line, and the AUC score is 0.5. The ROC curve of a perfect classifier goes through the point (0 (false positive rate), 1 (true positive rate)), and the AUC score is 1. A classifier with the AUC score closer to 1 has better performance. We can see from these tables that the results by the CRF method with MI are better than those by the CRF method without MI, and that the results by the CRF method without MI are better than those by the EM method and the association method. It is also seen that the results by the EM method are almost the same as those by the association method. It might be because the parameters of the EM method were estimated from the solution of the association method and the solution of the EM method already reached a local optimum. Figures 4, 5, and 6 show the average ROC curves for training and test datasets by the CRF method with MI, that without MI, the EM method, and the association method. For training datasets, the results by all of the methods were almost perfect. For test datasets, the CRF method with MI outperformed that without MI, the EM method, and the association method. It should be noted that the ROC curves of the EM method are almost the same as those of the association method for the same reason discussed above.

## Conclusions

We proposed novel methods which combine conditional random fields with the domain-based model of protein-protein interactions. In order to give better performance, we introduced mutual information to the probabilistic model. In the improved model, mutual information between domains is given as conditions, where MI between domains is defined as the maximum of MIs between residues in the domains. This method was developed based on the fact that amino acid residues at important sites for interactions have coevolved with each other, and MI has been used for identifying contact residues in interactions. We performed five-fold cross-validation experiments, and calculated AUC for probabilities that two proteins interact. The results suggested that our proposed methods, especially the CRF method with mutual information, are useful. However, the results of AUC for training datasets implied that estimated parameters were overfitting to training datasets. For avoiding that problem, we can improve the methods, for instance, by adding regularization terms, *l*_{1}-norm of parameters to the log-likelihood function. Since CRF has an advantage to be able to incorporate large number of features, it remains as a future work to improve the model itself to obtain better accuracy by, for instance, modifying the local feature and adding new features.

## Authors contributions

JS proposed the use of mutual information for predicting protein-protein interactions. Methods were developed and implemented by MH. MK and TA participated in the discussion during development of the methods. The manuscript was prepared by MH, JS, and TA.

## References

White RA, Szurmant H, Hoch JA, Hwa T: Features of protein-protein interactions in two-component signaling deduced from genomic libraries. Methods Enzymol. 2007, 422: 75-101. full_text. full_text

Burger L, van Nimwegen E: Accurate prediction of protein-protein interactions from sequence alignments using a Bayesian method. Molecular Systems Biology. 2008, 4: 165-10.1038/msb4100203.

Halabi N, Rivoire O, Leibler S, Ranganathan R: Protein sectors: Evolutionary units of three-dimensional structure. Cell. 2009, 138: 774-786. 10.1016/j.cell.2009.07.038.

Weigt M, White RA, Szurmant H, Hoch JA, Hwa T: Identification of direct residue contacts in protein-protein interaction by message passing. Proc. Natl. Acad. Sci. USA. 2009, 106: 67-72. 10.1073/pnas.0805923106.

Fraser HB, Hirsh AE, Steinmetz LM, Scharfe C, Feldman MW: Evolutionary rate in the protein interaction network. Science. 2002, 296: 750-752. 10.1126/science.1068696.

Sha F, Pereira F: Shallow parsing with conditional random fields. Proc. HLT-NAACL 2003. 2003, 134-141.

Sutton C, McCallum A: An introduction to conditional random fields for relational learning. Introduction to statistical relational learning. 2006, MIT Press, 93-128.

Deng M, Zhang K, Mehta S, Chen T, Sun F: Prediction of protein function using protein-protein interaction data. Journal of Computational Biology. 2003, 10 (6): 947-960. 10.1089/106652703322756168.

Deng M, Chen T, Sun F: An integrated probabilistic model for functional prediction of proteins. Journal of Computational Biology. 2004, 11: 463-475. 10.1089/1066527041410346.

Deng M, Mehta S, Sun F, Chen T: Inferring domain-domain interactions from protein-protein interactions. Genome Research. 2002, 12: 1540-1548. 10.1101/gr.153002.

Hayashida M, Ueda N, Akutsu T: Inferring strengths of protein-protein interactions from experimental data using linear programming. Bioinformatics. 2003, 19 (suppl 2): ii58-ii65. 10.1093/bioinformatics/btg1061.

Chen L, Wu LY, Wang Y, Zhang XS: Inferring protein interactions from experimental data by association probabilistic method. Proteins. 2006, 62 (4): 833-837. 10.1002/prot.20783.

Sprinzak E, Margalit H: Correlated sequence-signatures as markers of protein-protein interaction. Journal of Molecular Biology. 2001, 311: 681-692. 10.1006/jmbi.2001.4920.

Little DY, Chen L: Identification of coevolving residues and coevolution potentials emphasizing structure, bond formation and catalytic coordination in protein evolution. PLoS One. 2009, 4: e4762-10.1371/journal.pone.0004762.

Moussouri J: Gibbs and Markov random systems with constraints. Journal of Statistical Physics. 1974, 10: 11-33. 10.1007/BF01011714.

Bertsekas DP: Nonlinear Programming. 1999, Athena Scientific

Salwinski L, Miller CS, Smith AJ, Pettit FK, Bowie JU, Eisenberg D: The Database of Interacting Proteins: 2004 update. Nucleic Acids Research. 2004, 32: D449-D451. 10.1093/nar/gkh086.

The UniProt Consortium: The Universal Protein Resource (UniProt) in 2010. Nucleic Acids Research. 2010, 38: D142-D148.

Finn RD, Mistry J, Tate J, Coggill P, Heger A, Pollington JE, Gavin OL, Gunasekaran P, Ceric G, Forslund K, Holm L, Sonnhammer ELL, Eddy SR, Bateman A: The Pfam protein families database. Nucleic Acids Research. 2010, 38: D211-D222. 10.1093/nar/gkp985.

Nocedal J: Updating quasi-Newton matrices with limited storage. Mathematics of Computation. 1980, 35 (151): 773-782. 10.1090/S0025-5718-1980-0572855-7.

## Acknowledgements

This work was partially supported by Grants-in-Aid #22240009 and #21700323 from MEXT, Japan. JS would like to thank the National Health and Medical Research Council of Australia (NHMRC) and the Chinese Academy of Sciences (CAS) for financially supporting this research via the NHMRC Peter Doherty Fellowship and the Hundred Talents Program of CAS.

This article has been published as part of *BMC Systems Biology* Volume 5 Supplement 1, 2011: Selected articles from the 4th International Conference on Computational Systems Biology (ISB 2010). The full contents of the supplement are available online at http://www.biomedcentral.com/1752-0509/5?issue=S1.

## Author information

### Authors and Affiliations

### Corresponding author

## Additional information

### Competing interests

The authors declare that they have no competing interests.

## 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

Hayashida, M., Kamada, M., Song, J. *et al.* Conditional random field approach to prediction of protein-protein interactions using domain information.
*BMC Syst Biol* **5**
(Suppl 1), S8 (2011). https://doi.org/10.1186/1752-0509-5-S1-S8

Published:

DOI: https://doi.org/10.1186/1752-0509-5-S1-S8