Skip to main content

Extract interaction detection methods from the biological literature



Considerable efforts have been made to extract protein-protein interactions from the biological literature, but little work has been done on the extraction of interaction detection methods. It is crucial to annotate the detection methods in the literature, since different detection methods shed different degrees of reliability on the reported interactions. However, the diversity of method mentions in the literature makes the automatic extraction quite challenging.


In this article, we develop a generative topic model, the Correlated Method-Word model (CMW model) to extract the detection methods from the literature. In the CMW model, we formulate the correlation between the different methods and related words in a probabilistic framework in order to infer the potential methods from the given document. By applying the model on a corpus of 5319 full text documents annotated by the MINT and IntAct databases, we observe promising results, which outperform the best result reported in the BioCreative II challenge evaluation.


From the promising experiment results, we can see that the CMW model overcomes the issues caused by the diversity in the method mentions and properly captures the in-depth correlations between the detection methods and related words. The performance outperforming the baseline methods confirms that the dependence assumptions of the model are reasonable and the model is competent for the practical processing.


Interaction detection method extraction

The study of protein interactions is one of the most pressing biological problems. In the literature mining community, considerable efforts have been made to automatically extract the protein-protein interactions (PPI) from the literature [13] and some practical systems have been put into use [4, 5].

Nevertheless, little work has been done to automatically extract the interaction detection methods from the literature. The detection methods available to identify protein interactions vary in their level of resolution and the confidence of reliability. Therefore, it is important to identify such detection methods in order to validate the reported interactions. Some interaction databases, such as MINT [6] and IntAct [7], require the interaction entries to be experimentally confirmed. However, manually annotating the detection methods in the literature is time-consuming: on average, the curation of a manuscript takes up 2–3 hours of an expert curator [8]. Therefore, there is great practical demand of automatically extracting the detection methods from the literature.

The first critical assessment of detection method extraction was carried out by the BioCreative II challenge evaluation [9]. But only two groups (out of sixteen) submitted their results.

The diversity of method mentions in the literature is the major obstacle precluding the automatic extraction. In the real situation, different authors prefer different words and phrases to describe the same methods. For example, the detection method "two hybrid" (MI:0018) has 7 related synonyms, e.g. "2-hybrid", "2 H ", "2 h", "classical two hybrid", "Gal4 transcription regeneration", "two-hybrid", "yeast two hybrid", and one exact synonym, e.g. "2 hybrid", in the MI ontology [10] definition (it includes the terms describing the interaction detection methods). Although the ontology has already included so many different descriptions, biologists would just mention "yeast 2-h", which is not included in the ontology, in their manuscripts.

To validate the diversity, we apply a string matching algorithm with all the names/synonyms from the MI ontology on a set of 740 documents, annotated with 96 methods and provided by the BioCreative II challenge evaluation. The matching performance is demonstrated in Table 1.

Table 1 String matching performance.

As Table 1 illustrates, the poor recall performance confirms the serious diversity, and the inferior precision stems from the simple matching algorithm, which does not take the context into consideration, since most of the matched names are not the exact methods applied in the document but the background knowledge. In this sense, the rigid dictionary-based matching strategy fails to address the practical problem.

Another straightforward solution is to treat the extraction issue as a classification problem – for each detection method in the ontology definition, a set of binary classifiers are built to make yes/no decisions [11, 12]. But the traditional discriminative classifiers make little attempt to uncover the probabilistic structure and the correlation within both input and output spaces. In the biological domain, ignoring the correlation within both methods and words would hinder the performance since there are intrinsic relations.

In another point of view, from the perspective of involvement of domain experts, some approaches achieved acceptable results on the small data set. In Rinaldi's work [13], they invited the biologists to summarize the keywords and patterns for the extraction task and manually refined the patterns according to the performance. Obviously, this manner is not suitable for the large-scale data processing and its flexibility is not desirable.

Generative topic model

Nowadays, in the machine learning community, the generative topic model is receiving more and more attentions. Latent Dirichlet Allocation (LDA) [14] is one of the most typical models. LDA reduces the complex process of producing a document into a small number of simple probabilistic steps and thus specifies a probability distribution over all possible documents. Using standard statistical techniques, one can invert the process and infer the set of latent topics responsible for generating a given set of documents [15].

LDA-like topic models are rapidly developed into quite different domains. Xing Wei [16] introduced the LDA model into information retrieval system and improved the retrieval performance; David Mimno [17] proposed the Author-Persona-Topic model to formulate the expertise of authors based on their publications; Fei-Fei Li [18] advanced a hierarchical generative model to classify natural scene in an unsupervised manner.

The advantages of the generative topic models are: 1) it would be easy to postulate complex latent structures responsible for a set of observations; 2) the correlation between different factors could be easily exploited by introducing the latent topic variables.

In this article, in order to extract the detection methods from the biological literature, we propose to formulate the correlation between the detection methods and related word occurrences in a probabilistic framework. In particular, we assume the applied methods are governed by a set of latent topics and the corresponding word descriptions are also influenced by the same topic factors, which characterize the correlation between the methods and related words. Under this setting, we appeal to the generative topic model to capture such latent correlations and infer the potential methods from the observed words by the statistic inference technique.

The intuitive notion behind the proposed model is that: different documents contain informative commonality in the descriptions of the same methods, therefore we propose to discover the common usage patterns for the desired methods from the latent correlations between the methods and related words. This manner is somehow analogous to the idea that to extract templates from the overlapping of different method descriptions. But the diversity in the method mentions brings the traditional template generation algorithms with low support and low confidence problems. Furthermore, when there are multiple methods in one document, the traditional approach would fail to figure out the latent correlations. In contrast, the generative model deals naturally with the missing data and provides a more feasible and theoretical framework.

The paper is organized as follows: in the Methods section, we present detailed descriptions about the proposed model and discuss the inference and parameter estimation procedures for the model; in the Results section, we perform extensive experiments to validate the proposed model; and in the Conclusions section, we would conclude the work and demonstrate our contributions in this paper.


Correlated Method-Word model

We present the Correlated Method-Word model (CMW model) to extract the detection methods from the biological literature. The CMW model is depicted in Figure 1 with graphical representation. In the standard graphical model formalism [19]: nodes represent the random variables and edges indicate the possible dependence. The joint probability can be obtained from the graph by taking the product of the conditional distribution of nodes given their parents, see Eq(1).

Figure 1
figure 1

Graphical model representation of the CMW model. Following the standard graphical model formalism [19]: nodes represent the random variables and edges indicate the possible dependence. The joint probability can be obtained from the graph by taking the product of the conditional distribution of nodes given their parents.

The model can be viewed in the terms of generative process that, the author should first select a set of topics for his/her manuscripts (e.g. physical protein-protein interactions); under different kind of topics, there are different choices of detection methods to confirm the findings (e.g. pull down to confirm protein interactions); the selected methods are represented by the particular word occurrences (e.g. descriptions of the experiment conditions, properties and materials), which are also governed by the selected topics. Therefore, the correlations between the detection methods and related words are characterized by the latent topic factors; and from the observed words, we are able to infer the potential methods in the given document according to such correlations.

Formally, we define a corpus consists of D documents, E methods and V words, and a given document consists of N methods and M words. To simplify the model, we have assumed the topic size k is known and fixed on the whole corpus. In the given document d, we denote θ as the document-specific topic distribution; z = {z1, z2, z3,..., z N } as the particular discrete topic assignments for each method; y = {y1, y2, y3,..., y M } as the indexing variables to indicate which topic factor generates the corresponding word and ϵ as the method distribution under the topics. These are the latent variables. e = {e1, e2, e3,..., e N } and w = {w1, w2, w3,..., w M } are the observed methods and words in document d. Besides, α and η are the parameters of k-dimensional and E-dimensional Dirichlet distributions that postulate the topic and method prior distributions on the corpus and β is a k × V matrix, which represents the word distribution under topics. These are the model parameters.

Conditioned on the model parameters (α, β, η), the CMW model assumes the following generative process of the methods and related words in one document:

  1. 1.

    Sample topic proportion θ from the Dirichlet distribution: θ ~Dir(α)

  2. 2.

    For each method e n , n {1, 2, 3,..., N}:

  3. a.

    Sample topic factor z n from the multinomial distribution : z n ~Mul(θ)

  4. b.

    Sample method e n from the multinomial distribution conditioned on z n : e n ~p(e n |ϵ, z n )

  5. 3.

    For each related word w m , m {1, 2, 3,..., M}:

  6. a.

    Sample indexing variable y m from the Uniform distribution conditioned on N: y m ~Unif (1, 2, 3,..., N)

  7. b.

    Sample word w m from the multinomial distribution conditioned on z y m MathType@MTEF@5@5@+=feaagaart1ev2aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xH8viVGI8Gi=hEeeu0xXdbba9frFj0xb9qqpG0dXdb9aspeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGaemOEaO3aaSbaaSqaaiabdMha5naaBaaameaacqWGTbqBaeqaaaWcbeaaaaa@3094@ : w m ~p(w m |β, z y m MathType@MTEF@5@5@+=feaagaart1ev2aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xH8viVGI8Gi=hEeeu0xXdbba9frFj0xb9qqpG0dXdb9aspeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGaemOEaO3aaSbaaSqaaiabdMha5naaBaaameaacqWGTbqBaeqaaaWcbeaaaaa@3094@ )

Our basic notion about each component of this model is that, the discrete occurrences of detection methods and related words in the given document are governed by the topic-specific distributions (e.g. matrix ϵ and β) respectively. We use such conditional distribution to bridge the correlation between the methods and word occurrences: under different topics, there are different choices of detection methods and the corresponding word descriptions. To formulate this notion in a probabilistic framework, we follow the general settings in the LDA model that we assume the document-specific topic proportion θ is drawn from the k-dimensional Dirichlet distribution Dir(α), which determines the topic mixture proportion. Especially, we treat the parameter of method's multinomial distribution ϵ as a k × E matrix (one row represents for each mixture component), and, to avoid over-fitting caused by the unbalanced and sparse method occurrences, we assume that each row of ϵ is independently drawn from the E-dimensional Dirichlet distribution: ϵ i ~Dir(η), which smooths the method distribution under each topic. Each row of matrix β represents the particular word distribution under the topics. Besides, since the correlation between the methods and word occurrences is underlying (a document usually associates with multiple detection methods), we use the indexing variable y to indicate such latent structure between them.

Thus, the joint probability on the observed methods, words and latent variables in one document is given as follows:

p ( e , w , θ , ϵ , y , z | α , β , η ) = p ( θ | α ) i = 1 k p ( ϵ i | η ) ( n = 1 N p ( z n | θ ) p ( e n | ϵ , z n ) ) ( m = 1 M p ( y m | N ) p ( w m | β , z y m ) ) MathType@MTEF@5@5@+=feaagaart1ev2aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xI8qiVKYPFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGaemiCaaNaeiikaGIaeCyzauMaeiilaWIaeC4DaCNaeiilaWIaeqiUdeNaeiilaWYeuvgwd1utHrhAjrxySL2yaeHbJ1wBPfdmaGabaiab=v=aYkabcYcaSiabhMha5jabcYcaSiabhQha6jabcYha8jabeg7aHjabcYcaSiabek7aIjabcYcaSiabeE7aOjabcMcaPiabg2da9iabdchaWjabcIcaOiabeI7aXjabcYha8jabeg7aHjabcMcaPmaarahabaGaemiCaaNaeiikaGIae8x9di7aaSbaaSqaaiabdMgaPbqabaGccqGG8baFcqaH3oaAcqGGPaqkaSqaaiabdMgaPjabg2da9iabigdaXaqaaiabdUgaRbqdcqGHpis1aOWaaeWaaeaadaqeWbqaaiabdchaWjabcIcaOiabdQha6naaBaaaleaacqWGUbGBaeqaaOGaeiiFaWNaeqiUdeNaeiykaKIaemiCaaNaeiikaGIaemyzau2aaSbaaSqaaiabd6gaUbqabaGccqGG8baFcqWF1pGScqGGSaalcqWG6bGEdaWgaaWcbaGaemOBa4gabeaakiabcMcaPaWcbaGaemOBa4Maeyypa0JaeGymaedabaGaemOta4eaniabg+GivdaakiaawIcacaGLPaaadaqadaqaamaarahabaGaemiCaaNaeiikaGIaemyEaK3aaSbaaSqaaiabd2gaTbqabaGccqGG8baFcqWGobGtcqGGPaqkcqWGWbaCcqGGOaakcqWG3bWDdaWgaaWcbaGaemyBa0gabeaakiabcYha8jabek7aIjabcYcaSiabdQha6naaBaaaleaacqWG5bqEcqWGTbqBaeqaaOGaeiykaKcaleaacqWGTbqBcqGH9aqpcqaIXaqmaeaacqWGnbqta0Gaey4dIunaaOGaayjkaiaawMcaaaaa@A655@

An intuitive comparison between the traditional approach (e.g. discriminative classification and template matching method) and the proposed CMW model is illustrated in Figure 2.

Figure 2
figure 2

Comparison between the traditional approach and the CMW model. In the this representation, e denotes the detection methods associating with the document, w denotes the observed words and t in the right panel denotes the latent topic factors in the CMW model.

The traditional approach (the left panel of Figure 2) simply assumes the relation between the detection methods and related words is determined by the direct mapping. On the contrary, the CMW model (the right panel of Figure 2) formulates the relationship within a more throughout consideration: via the latent topic factors, word occurrences are formulated as a finite mixture under particular methods, so that they are not restricted to any methods and multiple words could contribute to the same method. This framework is more suitable and robust to deal with the diversity in the method descriptions. Furthermore, the discriminative classification algorithms assume the methods are independent in prior and the words are also independent when observing the given methods. Thus they would neglect the latent patterns within both methods and words. But in the CMW model, different topics govern dissimilar methods and words occurrences, embedding the correlation not only between different methods but also within the related words (see the Correlation between methods and words section and the Methods correlation analysis section for the detailed experiment results).

Efficient dimensional decomposition is explicitly implemented: V-dimensional word space and E-dimensional method space are mapped into the k-dimensional topic space, in which it will be easier for us to reveal the latent correlations between the detection methods and the variant word occurrences.

Inference and parameter estimation

Variational inference

In order to utilize the CMW model, we need to compute the posterior distribution of the methods in a given document, that is:

p ( e , θ , ϵ , y , z | w , α , β , η ) = p ( e , w , θ , ϵ , y , z | α , β , η ) p ( w | α , β , η ) MathType@MTEF@5@5@+=feaagaart1ev2aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xI8qiVKYPFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGaemiCaaNaeiikaGIaeCyzauMaeiilaWIaeqiUdeNaeiilaWYeuvgwd1utHrhAjrxySL2yaeHbJ1wBPfdmaGabaiab=v=aYkabcYcaSiabhMha5jabcYcaSiabhQha6jabcYha8jabhEha3jabcYcaSiabeg7aHjabcYcaSiabek7aIjabcYcaSiabeE7aOjabcMcaPiabg2da9KqbaoaalaaabaGaemiCaaNaeiikaGIaeCyzauMaeiilaWIaeC4DaCNaeiilaWIaeqiUdeNaeiilaWIae8x9diRaeiilaWIaeCyEaKNaeiilaWIaeCOEaONaeiiFaWNaeqySdeMaeiilaWIaeqOSdiMaeiilaWIaeq4TdGMaeiykaKcabaGaemiCaaNaeiikaGIaeC4DaCNaeiiFaWNaeqySdeMaeiilaWIaeqOSdiMaeiilaWIaeq4TdGMaeiykaKcaaaaa@7503@

Unfortunately, this posterior distribution is intractable: the couples between the continuous variable θ and discrete variable β, ϵ induce a combinatorial number of terms, making it impossible to efficiently get the exact inference result.

Although the exact inference is intractable, there are a wide variety of approximate inference algorithms can serve for the propose, including: expectation propagation [20], variational inference [21] and Markov chain Monte Carlo (MCMC) [22] etc. For computational efficiency, we develop a variational inference procedure to approximate the lower bound of the desired posterior distribution of methods in a given document.

In particular, we define the following fully factorized distribution on the latent variables:

p ( θ , ϵ , y , z | γ , ϕ , λ , σ ) = p ( θ | γ ) i = 1 k q ( ϵ i | σ i ) n = 1 N q ( z n | ϕ n ) M = 1 M q ( y m | λ m ) MathType@MTEF@5@5@+=feaagaart1ev2aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xI8qiVKYPFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGaemiCaaNaeiikaGIaeqiUdeNaeiilaWYeuvgwd1utHrhAjrxySL2yaeHbJ1wBPfdmaGabaiab=v=aYkabcYcaSiabhMha5jabcYcaSiabhQha6jabcYha8jabeo7aNjabcYcaSiabew9aMjabcYcaSiabeU7aSjabcYcaSiabeo8aZjabcMcaPiabg2da9iabdchaWjabcIcaOiabeI7aXjabcYha8jabeo7aNjabcMcaPmaarahabaGaemyCaeNaeiikaGIae8x9di7aaSbaaSqaaiabdMgaPbqabaGccqGG8baFcqaHdpWCdaWgaaWcbaGaemyAaKgabeaakiabcMcaPaWcbaGaemyAaKMaeyypa0JaeGymaedabaGaem4AaSganiabg+GivdGcdaqeWbqaaiabdghaXjabcIcaOiabdQha6naaBaaaleaacqWGUbGBaeqaaOGaeiiFaWNaeqy1dy2aaSbaaSqaaiabd6gaUbqabaGccqGGPaqkaSqaaiabd6gaUjabg2da9iabigdaXaqaaiabd6eaobqdcqGHpis1aOWaaebCaeaacqWGXbqCcqGGOaakcqWG5bqEdaWgaaWcbaGaemyBa0gabeaakiabcYha8jabeU7aSnaaBaaaleaacqWGTbqBaeqaaOGaeiykaKcaleaacqWGnbqtcqGH9aqpcqaIXaqmaeaacqWGnbqta0Gaey4dIunaaaa@8A5F@

where the Dirichlet parameters γ, σ and the Multinomial parameters ϕ, λ are free variational parameters.

The meaning of the above variational distribution is that: we discard the dependence among the latent variables by assuming they are independently drawn from the respective distributions. In that case, the aim of the variational inference is to find the optimal variational parameters which could minimize the Kullback-Leibler (KL) divergence between the variational distribution and the true posterior distribution.

Following the general recipe for the variational approximation, we take derivatives with respect to the variational parameters and obtain the following coordinate ascent algorithm:

  1. 1.

    Dirichlet parameter γ:

    γ i = α i + n = 1 N ϕ n i MathType@MTEF@5@5@+=feaagaart1ev2aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xI8qiVKYPFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGaeq4SdC2aaSbaaSqaaiabdMgaPbqabaGccqGH9aqpcqaHXoqydaWgaaWcbaGaemyAaKgabeaakiabgUcaRmaaqahabaGaeqy1dy2aaSbaaSqaaiabd6gaUjabdMgaPbqabaaabaGaemOBa4Maeyypa0JaeGymaedabaGaemOta4eaniabggHiLdaaaa@3FDE@
  2. 2.

    Multinomial parameter ϕ:

    log ϕ n i m = 1 M λ m n w m s β i s + [ ψ ( γ i ) ψ ( n = 1 k γ t ) ] + e n j [ ψ ( σ i j ) ψ ( t = 1 E σ i t ) ] MathType@MTEF@5@5@+=feaagaart1ev2aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xI8qiVKYPFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGagiiBaWMaei4Ba8Maei4zaCMaeqy1dy2aaSbaaSqaaiabd6gaUjabdMgaPbqabaGccqGHDisTdaaeWbqaaiabeU7aSnaaBaaaleaacqWGTbqBcqWGUbGBaeqaaOGaem4DaC3aa0baaSqaaiabd2gaTbqaaiabdohaZbaakiabek7aInaaBaaaleaacqWGPbqAcqWGZbWCaeqaaaqaaiabd2gaTjabg2da9iabigdaXaqaaiabd2eanbqdcqGHris5aOGaey4kaSIaei4waSLaeqiYdKNaeiikaGIaeq4SdC2aaSbaaSqaaiabdMgaPbqabaGccqGGPaqkcqGHsislcqaHipqEcqGGOaakdaaeWbqaaiabeo7aNnaaBaaaleaacqWG0baDaeqaaaqaaiabd6gaUjabg2da9iabigdaXaqaaiabdUgaRbqdcqGHris5aOGaeiykaKIaeiyxa0Laey4kaSIaemyzau2aa0baaSqaaiabd6gaUbqaaiabdQgaQbaakiabcUfaBjabeI8a5jabcIcaOiabeo8aZnaaBaaaleaacqWGPbqAcqWGQbGAaeqaaOGaeiykaKIaeyOeI0IaeqiYdKNaeiikaGYaaabCaeaacqaHdpWCdaWgaaWcbaGaemyAaKMaemiDaqhabeaaaeaacqWG0baDcqGH9aqpcqaIXaqmaeaacqWGfbqra0GaeyyeIuoakiabcMcaPiabc2faDbaa@838D@

where w m s MathType@MTEF@5@5@+=feaagaart1ev2aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xH8viVGI8Gi=hEeeu0xXdbba9frFj0xb9qqpG0dXdb9aspeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGaem4DaC3aa0baaSqaaiabd2gaTbqaaiabdohaZbaaaaa@304B@ means the m th word in the document is the s th one in the vocabulary, and e n j MathType@MTEF@5@5@+=feaagaart1ev2aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xH8viVGI8Gi=hEeeu0xXdbba9frFj0xb9qqpG0dXdb9aspeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGaemyzau2aa0baaSqaaiabd6gaUbqaaiabdQgaQbaaaaa@3017@ means the n th method in the document is the j th method in the list.

  1. 3.

    Multinomial parameter λ:

    log λ m n i = 1 k ϕ n i w m s β i s MathType@MTEF@5@5@+=feaagaart1ev2aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xI8qiVKYPFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGagiiBaWMaei4Ba8Maei4zaCMaeq4UdW2aaSbaaSqaaiabd2gaTjabd6gaUbqabaGccqGHDisTdaaeWbqaaiabew9aMnaaBaaaleaacqWGUbGBcqWGPbqAaeqaaOGaem4DaC3aa0baaSqaaiabd2gaTbqaaiabdohaZbaakiabek7aInaaBaaaleaacqWGPbqAcqWGZbWCaeqaaaqaaiabdMgaPjabg2da9iabigdaXaqaaiabdUgaRbqdcqGHris5aaaa@4B2F@
  2. 4.

    Dirichlet parameter σ:

    σ i j = n j + d = 1 D n = 1 N d ϕ d n i e d n j MathType@MTEF@5@5@+=feaagaart1ev2aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xI8qiVKYPFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGaeq4Wdm3aaSbaaSqaaiabdMgaPjabdQgaQbqabaGccqGH9aqpcqWGUbGBdaWgaaWcbaGaemOAaOgabeaakiabgUcaRmaaqahabaWaaabCaeaacqaHvpGzdaWgaaWcbaGaemizaqMaemOBa4MaemyAaKgabeaakiabdwgaLnaaDaaaleaacqWGKbazcqWGUbGBaeaacqWGQbGAaaaabaGaemOBa4Maeyypa0JaeGymaedabaGaemOta40aaSbaaWqaaiabdsgaKbqabaaaniabggHiLdaaleaacqWGKbazcqGH9aqpcqaIXaqmaeaacqWGebara0GaeyyeIuoaaaa@5025@

These update equations are invoked repeatedly until the relative change in KL is small (< 0.0001%).

When we have achieved the approximate posterior probability, we could handle the conditional distribution of interest – p(e|w α, β, η) to infer the potential methods in a given document as follows:

p ( e | w , α , β , η ) n = 1 N z n q ( z n | ϕ n ) p ( e | ϵ , z n ) p ( ϵ | η ) MathType@MTEF@5@5@+=feaagaart1ev2aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xI8qiVKYPFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGaemiCaaNaeiikaGIaeCyzauMaeiiFaWNaeC4DaCNaeiilaWIaeqySdeMaeiilaWIaeqOSdiMaeiilaWIaeq4TdGMaeiykaKIaeyisIS7aaabCaeaadaaeqbqaaiabdghaXjabcIcaOiabdQha6naaBaaaleaacqWGUbGBaeqaaOGaeiiFaWNaeqy1dy2aaSbaaSqaaiabd6gaUbqabaGccqGGPaqkcqWGWbaCcqGGOaakcqWHLbqzcqGG8baFtqvzynutnfgDOLeDHXwAJbqegmwBTLwmWaaceaGae8x9diRaeiilaWIaemOEaO3aaSbaaSqaaiabd6gaUbqabaGccqGGPaqkcqWGWbaCcqGGOaakcqWF1pGScqGG8baFcqaH3oaAcqGGPaqkaSqaaiabdQha6naaBaaameaacqWGUbGBaeqaaaWcbeqdcqGHris5aaWcbaGaemOBa4Maeyypa0JaeGymaedabaGaemOta4eaniabggHiLdaaaa@7031@

Parameter estimation

Following the similar procedure in the variational inference, in this section, we utilize an empirical Bayesian method to estimate the parameters of the CMW model. This time, we are looking for the optimal model parameters to tighten the lower bound of likelihood and obtain the following update equations:

  1. 1.

    Update the Dirichlet parameter α by the Newton-Raphson algorithm:

    L ( α ) α i = d = 1 D { ψ ( t = 1 k α t ) ψ ( α i ) + ψ ( γ d i ) ψ ( t = 1 k γ d t ) } MathType@MTEF@5@5@+=feaagaart1ev2aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xI8qiVKYPFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaqcfa4aaSaaaeaacqGHciITcqWGmbatcqGGOaakcqaHXoqycqGGPaqkaeaacqGHciITcqaHXoqydaWgaaqaaiabdMgaPbqabaaaaOGaeyypa0ZaaabCaeaadaGadaqaaiabeI8a5jabcIcaOmaaqahabaGaeqySde2aaSbaaSqaaiabdsha0bqabaaabaGaemiDaqNaeyypa0JaeGymaedabaGaem4AaSganiabggHiLdGccqGGPaqkcqGHsislcqaHipqEcqGGOaakcqaHXoqydaWgaaWcbaGaemyAaKgabeaakiabcMcaPiabgUcaRiabeI8a5jabcIcaOiabeo7aNnaaBaaaleaacqWGKbazcqWGPbqAaeqaaOGaeiykaKIaeyOeI0IaeqiYdKNaeiikaGYaaabCaeaacqaHZoWzdaWgaaWcbaGaemizaqMaemiDaqhabeaaaeaacqWG0baDcqGH9aqpcqaIXaqmaeaacqWGRbWAa0GaeyyeIuoakiabcMcaPaGaay5Eaiaaw2haaaWcbaGaemizaqMaeyypa0JaeGymaedabaGaemiraqeaniabggHiLdaaaa@6F49@
    2 L ( α ) α i α j = D { ψ ( t = 1 k α t ) δ ( i , j ) ψ ( α i ) } MathType@MTEF@5@5@+=feaagaart1ev2aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xI8qiVKYPFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaqcfa4aaSaaaeaacqGHciITdaahaaqabeaacqaIYaGmaaGaemitaWKaeiikaGIaeqySdeMaeiykaKcabaGaeyOaIyRaeqySde2aaSbaaeaacqWGPbqAaeqaaiabgkGi2kabeg7aHnaaBaaabaGaemOAaOgabeaaaaGccqGH9aqpcqWGebardaGadaqaaiqbeI8a5zaafaGaeiikaGYaaabCaeaacqaHXoqydaWgaaWcbaGaemiDaqhabeaaaeaacqWG0baDcqGH9aqpcqaIXaqmaeaacqWGRbWAa0GaeyyeIuoakiabcMcaPiabgkHiTiabes7aKjabcIcaOiabdMgaPjabcYcaSiabdQgaQjabcMcaPiqbeI8a5zaafaGaeiikaGIaeqySde2aaSbaaSqaaiabdMgaPbqabaGccqGGPaqkaiaawUhacaGL9baaaaa@5D6A@

where δ (i, j) = 1 when j = k, otherwise 0.

  1. 2.

    Update the Dirichlet parameter η by the Newton-Raphson algorithm:

    L ( η ) η j = i = 1 k { ψ ( t = 1 E η t ) ψ ( η j ) + ψ ( σ i j ) ψ ( t = 1 E σ i t ) } MathType@MTEF@5@5@+=feaagaart1ev2aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xI8qiVKYPFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaqcfa4aaSaaaeaacqGHciITcqWGmbatcqGGOaakcqaH3oaAcqGGPaqkaeaacqGHciITcqaH3oaAdaWgaaqaaiabdQgaQbqabaaaaOGaeyypa0ZaaabCaeaadaGadaqaaiabeI8a5jabcIcaOmaaqahabaGaeq4TdG2aaSbaaSqaaiabdsha0bqabaaabaGaemiDaqNaeyypa0JaeGymaedabaGaemyraueaniabggHiLdGccqGGPaqkcqGHsislcqaHipqEcqGGOaakcqaH3oaAdaWgaaWcbaGaemOAaOgabeaakiabcMcaPiabgUcaRiabeI8a5jabcIcaOiabeo8aZnaaBaaaleaacqWGPbqAcqWGQbGAaeqaaOGaeiykaKIaeyOeI0IaeqiYdKNaeiikaGYaaabCaeaacqaHdpWCdaWgaaWcbaGaemyAaKMaemiDaqhabeaaaeaacqWG0baDcqGH9aqpcqaIXaqmaeaacqWGfbqra0GaeyyeIuoakiabcMcaPaGaay5Eaiaaw2haaaWcbaGaemyAaKMaeyypa0JaeGymaedabaGaem4AaSganiabggHiLdaaaa@6F8F@
    2 L ( η ) η i η j = k { ψ ( t = 1 E η t ) δ ( i , j ) ψ ( η i ) } MathType@MTEF@5@5@+=feaagaart1ev2aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xI8qiVKYPFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaqcfa4aaSaaaeaacqGHciITdaahaaqabeaacqaIYaGmaaGaemitaWKaeiikaGIaeq4TdGMaeiykaKcabaGaeyOaIyRaeq4TdG2aaSbaaeaacqWGPbqAaeqaaiabgkGi2kabeE7aOnaaBaaabaGaemOAaOgabeaaaaGccqGH9aqpcqWGRbWAdaGadaqaaiqbeI8a5zaafaGaeiikaGYaaabCaeaacqaH3oaAdaWgaaWcbaGaemiDaqhabeaaaeaacqWG0baDcqGH9aqpcqaIXaqmaeaacqWGfbqra0GaeyyeIuoakiabcMcaPiabgkHiTiabes7aKjabcIcaOiabdMgaPjabcYcaSiabdQgaQjabcMcaPiqbeI8a5zaafaGaeiikaGIaeq4TdG2aaSbaaSqaaiabdMgaPbqabaGccqGGPaqkaiaawUhacaGL9baaaaa@5DAD@
  2. 3.

    Update the Multinomial parameter β:

    β j s d = 1 D n = 1 N d m = 1 M d λ d m n w d m s ϕ d n j MathType@MTEF@5@5@+=feaagaart1ev2aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xI8qiVKYPFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGaeqOSdi2aaSbaaSqaaiabdQgaQjabdohaZbqabaGccqGHDisTdaaeWbqaamaaqahabaWaaabCaeaacqaH7oaBdaWgaaWcbaGaemizaqMaemyBa0MaemOBa4gabeaakiabdEha3naaDaaaleaacqWGKbazcqWGTbqBaeaacqWGZbWCaaGccqaHvpGzdaWgaaWcbaGaemizaqMaemOBa4MaemOAaOgabeaaaeaacqWGTbqBcqGH9aqpcqaIXaqmaeaacqWGnbqtdaWgaaadbaGaemizaqgabeaaa0GaeyyeIuoaaSqaaiabd6gaUjabg2da9iabigdaXaqaaiabd6eaonaaBaaameaacqWGKbazaeqaaaqdcqGHris5aaWcbaGaemizaqMaeyypa0JaeGymaedabaGaemiraqeaniabggHiLdaaaa@5B2C@

These update equations correspond to find the maximum likelihood estimation with the expected sufficient statistics for each document taken under the variational posterior.

We develop an alternating EM procedure to find the optimal parameters as follows:

  1. 1.

    (E-Step) For each document in the training corpus, optimizing the variational parameters (γ, ϕ, λ, σ) according to equations (3) – (6);

  2. 2.

    (M-Step) Maximizing the resulting lower bound on the variational likelihood on the whole corpus with respect to the model parameters (α, β, η) according to equations (8) – (12).

The E-Step and M-Step are repeated until the bound on the likelihood converges (relative change in likelihood is less than 0.001%). The convergency rate of the process depends on the size of parameters in the model, (e.g. number of words, methods and topics). In our experiments (3000 words, 115 methods and up to 500 topics), the algorithm terminates in less than 30 iterations in all the cases.

Results and discussion

We collect 5319 full-text documents from PubMed [23] with method annotations from another two public curated interaction databases: MINT and IntAct. We perform the following pre-processions on the data set: 1) parsing the HTML file; 2) converting the words into lower cases; 3) removing a standard list of 400 stop words, punctuations, and the terms occur less than 50 times; 4) stemming the words to its root by Porter Stemming [24]. We utilize the macro-precision, macro-recall and macro-Fscore [25] to evaluate the performance in average.

Test corpora

The whole corpus consists of 115 unique method annotations, and each document associates with 1.99 different methods in average. Unfortunately, the standard deviation of the method frequency is so large that the corpus is heavily unbalanced: the most popular method "pull down" (MI:0096) occurs 2040 times while there are 57 methods (49.6% of all) occurs less than 10 times. Figure 3 demonstrates the unbalanced method distribution on the whole corpus.

Figure 3
figure 3

Statistics of the corpus. In the whole corpus, 5 dominate detection methods take up nearly 59.3% occurrences and 86.1% (99 out of 115) methods occur in less than 10% documents.

We can discover from Figure 3: 1) the 5 dominate detection methods, i.e. pull down (MI:0096), 2 hybrid (MI:0018), coip (MI:0019), anti tag coip (MI:0007) and anti bait coip (MI:0006), take up nearly 59.3% occurrences in the whole corpus; 2) 86.1% (99 out of 115) methods occur in less than 10% documents. In this case, smoothing the estimated parameters is essential to achieve better performance.

Feature selection

The CMW model is proposed to capture the correlation between methods and the "related" words. However, no curations explicitly annotate which words or sentences are related to the curated methods. So we employ χ2 statistic [26] to select the most relevant feature words from the whole text.

Word t's χ2 value associating with the method e is calculated according to the following equation:

χ 2 ( t , e ) = N × ( A D B C ) 2 ( A + C ) × ( B + D ) × ( A + B ) × ( C + D ) MathType@MTEF@5@5@+=feaagaart1ev2aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xI8qiVKYPFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGaeq4Xdm2aaWbaaSqabeaacqaIYaGmaaGccqGGOaakcqWG0baDcqGGSaalcqWGLbqzcqGGPaqkcqGH9aqpjuaGdaWcaaqaaiabd6eaojabgEna0kabcIcaOiabdgeabjabdseaejabgkHiTiabdkeacjabdoeadjabcMcaPmaaCaaabeqaaiabikdaYaaaaeaacqGGOaakcqWGbbqqcqGHRaWkcqWGdbWqcqGGPaqkcqGHxdaTcqGGOaakcqWGcbGqcqGHRaWkcqWGebarcqGGPaqkcqGHxdaTcqGGOaakcqWGbbqqcqGHRaWkcqWGcbGqcqGGPaqkcqGHxdaTcqGGOaakcqWGdbWqcqGHRaWkcqWGebarcqGGPaqkaaaaaa@5A29@

where A is the number of times t co-occurs with e, B is the number of times t occurs without e, C is the number of times e occurs without t, D is the number of times neither e or t occurs, and N is the total number of documents.

By χ2 statistic, we approximate the dependence between word t and method e, so that we can preserve the words, which are the most relevant to the method descriptions, by the following formulation:

χ a v g 2 ( t ) = i = 1 E p ( e i ) χ 2 ( t , e i ) MathType@MTEF@5@5@+=feaagaart1ev2aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xI8qiVKYPFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGaeq4Xdm2aa0baaSqaaiabdggaHjabdAha2jabdEgaNbqaaiabikdaYaaakiabcIcaOiabdsha0jabcMcaPiabg2da9maaqahabaGaemiCaaNaeiikaGIaemyzau2aaSbaaSqaaiabdMgaPbqabaGccqGGPaqkcqaHhpWydaahaaWcbeqaaiabikdaYaaakiabcIcaOiabdsha0jabcYcaSiabdwgaLnaaBaaaleaacqWGPbqAaeqaaOGaeiykaKcaleaacqWGPbqAcqGH9aqpcqaIXaqmaeaacqWGfbqra0GaeyyeIuoaaaa@4DAF@

where p(e i ) is the prior probability of method e i .

In the following experiments, we select the top 3000 terms to build up the feature set according to Eq(13).

Effect of topic factors

We first use the perplexity as the criterion to evaluate the effect of the number of topic factors, which is the only arbitrary parameter in the CMW model. The perplexity on a set of testing documents is calculated as follows:

p e r p l e x i t y = exp { d = 1 D n = 1 N d log p ( e n | w d ) d = 1 D N d } MathType@MTEF@5@5@+=feaagaart1ev2aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xI8qiVKYPFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGaemiCaaNaemyzauMaemOCaiNaemiCaaNaemiBaWMaemyzauMaemiEaGNaemyAaKMaemiDaqNaemyEaKNaeyypa0JagiyzauMaeiiEaGNaeiiCaaNaei4EaSxcfa4aaSaaaeaacqGHsisldaaeWaqaamaaqadabaGagiiBaWMaei4Ba8Maei4zaCMaemiCaaNaeiikaGIaemyzau2aaSbaaeaacqWGUbGBaeqaaiabcYha8jabhEha3naaBaaabaGaeCizaqgabeaacqGGPaqkaeaacqWGUbGBcqGH9aqpcqaIXaqmaeaacqWGobGtdaWgaaqaaiabdsgaKbqabaaacqGHris5aaqaaiabdsgaKjabg2da9iabigdaXaqaaiabdseaebGaeyyeIuoaaeaadaaeWaqaaiabd6eaonaaBaaabaGaemizaqgabeaaaeaacqWGKbazcqGH9aqpcqaIXaqmaeaacqWGebaraiabggHiLdaaaOGaeiyFa0haaa@6972@

where D is the set of testing documents and N d is the number of methods in the document d.

Better generalization capability is indicated by a lower perplexity over the held-out testing samples. We held out 20% of collection for the testing purpose and used the remaining 80% to train the model, in accordance with 5-fold cross-validation.

Figure 4 demonstrates that the generalization power of the CMW model gets improved with more topic factors. Since with more topic factors the documents could be partitioned into finer segments, more precise correlations between the methods and words could be captured. But as the number of topics exceeds a limit, the model becomes too specific (higher perplexity). Therefore we could conclude that the topic factors could be treated as the discriminate granularity of the model, that is it operates as a tradeoff between the generality and specificity. Besides, as the number of topic factors increase, there will be more parameters to be estimated (linearly increase with the number of topics), so that more training data is needed to obtain the reliable parameters. In this sense, when the number of topic factors exceeds a limit, the quality of the estimated parameters decreases and hampers the prediction power.

Figure 4
figure 4

Methods perplexity. Lower perplexity on the testing data indicates a better generalization capability. Here we held out 20% of collection for the testing purpose and used the remaining 80% to train the model, in accordance with 5-fold cross-validation.

Besides understanding the impact of the number of topic factors on the generalization capability, we would be more interested in their explicit effect on the extraction performance. Here, we evaluate the precision and recall performance of the model under different number of topic factors. We use the same data set partition as in Figure 4.

We could discover from Figure 5 that the extraction performance peaks close to the place where the perplexity reaches the minimum. This is consistent with the foregoing perplexity result. These results give us insight about determining the proper size of topics for the CMW model.

Figure 5
figure 5

Performance on the number of topics. We use the same data set partition as in Figure 4 and evaluate the precision and recall performance of the CMW model.

Extraction performance

Since there is few work to compare with, we employ the well studied Naïve Bayes, KNN and SVM as the baseline methods to evaluate the capability of the proposed CMW model. We choose Naïve Bayes because it is the simplest generative model with complete independence assumptions, and KNN model could exploit the heterogeneity among the similar documents. These are the two basic notions in the CMW model. Besides, SVM model is the most powerful discriminative model for the classification task with decent performance [11]. All the baseline models are operating on the same feature set as the CMW model employs.

In the Naïve Bayes model, we estimate the posterior probability of the methods in a given document by Eq(14). We use a pre-estimated threshold to retrieval the most probable methods.

p ( e | w ) n p ( w n | e ) p ( e ) MathType@MTEF@5@5@+=feaagaart1ev2aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xI8qiVKYPFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGaemiCaaNaeiikaGIaemyzauMaeiiFaWNaeC4DaCNaeiykaKIaeyyhIu7aaebuaeaacqWGWbaCcqGGOaakcqWG3bWDdaWgaaWcbaGaemOBa4gabeaakiabcYha8jabdwgaLjabcMcaPiabdchaWjabcIcaOiabdwgaLjabcMcaPaWcbaGaemOBa4gabeqdcqGHpis1aaaa@45F1@

In the KNN model, we make the prediction by ranking the candidate methods in the union of the unlabeled sample's k-nearest labeled neighbors, and weight the candidate methods by the similarity between the desired unlabeled sample and its neighbors.

In the SVM model, we follow Boutell's strategy [12] to train a set of binary classifiers for each method and predict the unknown methods by the classifiers' output. We use SV Mlight[27] toolkit to implement a linear kernel SVM model with the default parameters.

We perform comparisons on different proportions of the data used for training. In this comparison, we set the size of topics in the CMW model to be 250 and k in the KNN model to be 37.

We could discover from Figure 6 that, as the training set increases, the performance of the CMW model improves rapidly. The reason for this phenomenon is that in the CMW model, there are E + k(V + 1) parameters to be estimated, when the training set is not large enough, most of the parameters cannot be fully estimated, which would directly hinder the performance of the model.

Figure 6
figure 6

Comparison with the baseline models. We compare the F-score performance of the four models on different proportions of the data used for training. In this comparison, we set the size of topics in the CMW model to be 250 and k in the KNN model to be 37.

One thing we should note is that, since the data set is unbalanced, we should attend the retrieval performance on the minor methods as well. In the method-level evaluation, the baseline models only retrieve most of the major methods (e.g. the top 5 methods) but ignoring the other minor ones, while the CMW model exhibits superior retrieve power. We demonstrate the coverage performance of each model on the testing set to compare their retrieval capability.

Figure 7 demonstrates that the CMW model possesses better retrieval capability than all the baseline methods when the training set is large enough. We contribute the nice coverage performance to the smoothing factor introduced to the method distribution. Because the whole corpus is sparse and unbalanced, the minor methods possess little proportion in the training set. However, the baseline models do not take the sparseness into account, so that they fail to retrieve the minor ones from the testing cases. In contrast, the CMW model attends the smoothing issue and overcomes the sparseness.

Figure 7
figure 7

Coverage comparison with the baseline models. We compare the coverage performance of the four models on the same data set partition as in Figure 4 and we use the same model parameter settings.

Rinaldi utilized the expert revised patterns to perform the extraction and achieved the best performance in the BioCreative II challenge evaluation [13]. To compare with their approach, we operate the CMW model on the same testing corpus (300 full text documents) and set the topic size to be 300 according to the result in the previous section. The CMW model achieved competitive results (F-Score improved 12.4%), illustrated in Table 2.

Table 2 Comparison with BioCreative II best result.

Here, we briefly conclude the performance of the CMW model. The extraction performance outperforms the discriminative baseline methods confirms that the dependence assumptions in the proposed CMW model are reasonable. Besides, the traditional discriminative classifiers fail to model the correlation within either the methods or the related words, while in the biological domain such correlations convey important domain dependent information. In this sense, the major advantage of the CMW model is that it properly exploits such informative correlations to reinforce the extraction performance. The improvements against the manually revised templates approach validate that the CMW model does exploit more precise and general patterns for the desired methods from the large-scale statistics, confirming the reasonable underlying semantic structure from another perspective.

Correlation between methods and words

To demonstrate the correlation between the different methods and words exploited by the CMW model, we utilize the method-specific distribution over words by the conditional distribution p(w|e) to retrieval the most relevant terms under each desired method:

s ( w | e ) = d D log p ( w | e ) d D M d MathType@MTEF@5@5@+=feaagaart1ev2aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xI8qiVKYPFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGaem4CamNaeiikaGIaem4DaCNaeiiFaWNaemyzauMaeiykaKIaeyypa0tcfa4aaSaaaeaadaaeqaqaaiGbcYgaSjabc+gaVjabcEgaNjabdchaWjabcIcaOiabdEha3jabcYha8jabdwgaLjabcMcaPaqaaiabdsgaKjabgIGiolabdseaebqabiabggHiLdaabaWaaabeaeaacqWGnbqtdaWgaaqaaiabdsgaKbqabaaabaGaemizaqMaeyicI4SaemiraqeabeGaeyyeIuoaaaaaaa@4EAC@

where D is the set of documents associating with the desired method e and M d is the number of words in the document d.

In Table 3, we collect top 20 terms for 6 different methods according to Eq(15) from the corpus. We can see from the table, most of the terms are appropriately gathered to the given methods. For example, "structure", "crystal", "helix" are gathered to x-ray, and "yeast", "two-hybrid", "site" are gathered to two hybrid. These are the informative terms in the MI ontology definition of these methods. From this result, we could discover that the CMW model properly selects suitable "indicators" for the given methods. From another perspective, since these "indicators" are organized in a probability framework and accordingly contribute to the desired methods, the CMW model could better overcome the issue caused by the diversity in the method mentions. The reasonable word distribution under methods confirms that the CMW model captures the in-depth correlation between the methods and related words from the literature.

Table 3 Top 20 relevant terms for methods.

Methods correlation analysis

By the CMW model, we map different methods into the latent topic space, where we are able to analyze the relationship between the different methods. Meanwhile, there are intrinsic inherit relationships between the methods, defined in the MI ontology and organized as a concept hierarchy.

To represent a given method in the latent topic space, we re-normalize the topic-specific method distribution matrix ϵ by column as follows:

r ( e i ) = ϵ i s = 1 k ϵ i s MathType@MTEF@5@5@+=feaagaart1ev2aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xI8qiVKYPFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGaemOCaiNaeiikaGIaemyzau2aaSbaaSqaaiabdMgaPbqabaGccqGGPaqkcqGH9aqpjuaGdaWcaaqaambvLH1qn1uy0Hws0fgBPngaryWyT1wAXadaiqaacqWF1pGSdaWgaaqaaiabdMgaPbqabaaabaWaaabmaeaacqWF1pGSdaWgaaqaaiabdMgaPjabdohaZbqabaaabaGaem4CamNaeyypa0JaeGymaedabaGaem4AaSgacqGHris5aaaaaaa@4AFB@

where ϵ i is the i th column of ϵ matrix.

Recall that, each row of the multinomial parameter ϵ is the method distribution under a particular topic, so that each column of ϵ represents a method in the topic space. By normalizing ϵ by column, we can represent the different methods over the latent topic factors.

Based on this representation, we employ an accumulative clustering algorithm to perform the hierarchical clustering and utilize a visualization tool gCluto [28] to demonstrate the captured "pedigree" tree. (We only illustrate part of the clustering result because of the page limit.)

From the clustering result in Figure 8, we can discover that most of the sibling nodes within the MI ontology are successfully clustered with the correct hierarchy (red circles mean the correct clusters). The promising result confirms that the CMW model captures the proper correlations not only between the methods and words but also among the different methods. The traditional discriminative classifiers are not able to figure out such relationships.

Figure 8
figure 8

Methods clustering tree. We utilize an accumulative clustering algorithm to perform the hierarchical clustering and build up the "pedigree" tree of the detection methods. Red circles in the figure mean the correct clusters according to the MI ontology definition.

Classify irrelevant documents

Although the CMW model is proposed to address the extraction problem in documents with at least one detection method, in most situation, the curators don't know whether the document is PPI related or experimentally confirmed beforehand. So it is necessary to evaluate the model's capability to classify the irrelevant documents.

We randomly select 1000 documents from PubMed, none of which are annotated by MINT nor IntAct. These documents are taken as the irrelevant documents. Meanwhile, we randomly select another 1000 documents from the evaluation corpus as relevant documents. In Eq (17), we define the relevance score of each document by the posterior probability of the most potential method in that document as follows:

r e l e v a n c e ( d ) = max e p ( e | w d ) MathType@MTEF@5@5@+=feaagaart1ev2aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xI8qiVKYPFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGaemOCaiNaemyzauMaemiBaWMaemyzauMaemODayNaemyyaeMaemOBa4Maem4yamMaemyzauMaeiikaGIaemizaqMaeiykaKIaeyypa0ZaaCbeaeaacyGGTbqBcqGGHbqycqGG4baEaSqaaiabdwgaLbqabaGccqWGWbaCcqGGOaakcqWGLbqzcqGG8baFcqWH3bWDdaWgaaWcbaGaeCizaqgabeaakiabcMcaPaaa@4B17@

This measurement indicates the maximum probability of a document containing at least one interaction detection method.

We arrange the relevance scores in a descending order in Figure 9, so that it is easy to discover that the relevance scores in the relevant document set are significantly greater than those in the irrelevant document set. If we select the threshold as the green line indicated, we would achieve a promising classification performance: in terms of precision 0.745, recall 0.676 and AUC 0.819. The result indicates that the proposed CMW model possesses the capability to reject the irrelevant documents before extracting.

Figure 9
figure 9

Relevance distribution in relevant and irrelevant documents. In the diagram, red line indicates the relevance scores in the relevant document set and the blue dots indicate the relevance scores in the irrelevant document set. If we select the classification threshold as the green line indicates, we would achieve a promising classification performance: in terms of precision 0.745, recall 0.676 and AUC 0.819.


In this paper, we propose a generative probabilistic model, the Correlated Method-Word model, to automatically extract the interaction detection methods from the biological literature. This problem is not well studied by the previous researches. By introducing the latent topic factors, the proposed model formulates the correlation between the detection methods and related words in a probabilistic framework in order to infer the potential methods from the observed words.

In our experiments, the proposed CMW model achieved competitive performance against the other well-studied discriminative classifiers on a corpus of 5319 full text documents. And it outperforms the best result reported in the BioCreative II challenge evaluation (F-Score improved 12.4%). From the promising results, we could see that the proposed CMW model overcomes the diversity in the method descriptions and appropriately solve the detection method extraction issue. Furthermore, the model captures the in-depth relationship not only between the methods and related words (see the Correlation between methods and words section), but also among the different methods (see the Methods correlation analysis section). Most of the discriminative classifiers fail to exploit such relations. The competitive performance confirms that the dependence assumptions in the model are reasonable and it is necessary to model the correlation between the different methods and words in the detection method extraction issue.

Our contributions in this paper lie in: 1) propose a generative probabilistic model with proper underlying semantics for the detection method extraction issue, and the model achieves promising performance; 2) properly model the correlation between the detection methods and related words in the biological literature, which captures the in-depth relationship not only between the methods and related words but also among the different methods.

The CMW model is now integrating to our ONBIRES system [5] to provide on-line service. And in the future work, we are planning to associate the extracted methods with the annotated interaction pairs and retrieve the evidence sentences in the documents, which would provide a more throughout annotation of the protein interactions in the biological literature.


  1. Daraselia N, Yuryev A, Egorov S, Novichkova S, Nikitin A, Mazo I: Extracting human protein interactions from MEDLINE using a full-sentence parser. Bioinformatics 2004, 20(5):604–611.

    Article  CAS  PubMed  Google Scholar 

  2. Ono T, Hishigaki H, Tanigami A, Takagi T: Automated extraction of information on protein protein interactions from the biological literature. IEEE Intelligent Systems 2001, 17(2):155–161.

    CAS  Google Scholar 

  3. Huang M, Zhu X, Hao Y, Payan DG, Qu K, Li M: Discovering patterns to extract protein protein interactions from full texts. Bioinformatics 2004, 20(18):3604–3612.

    Article  CAS  PubMed  Google Scholar 

  4. Blaschke C, Valencia A: The Frame-Based Module of the SUISEKI Information Extraction System. Bioinformatics 2002, 17(2):14–20.

    Google Scholar 

  5. Huang M, Zhu X, Ding S, Yu H, Li M: ONBIRES: ONtology-based BIological Relation Extraction System. Proceedings of the Fourth Asia Pacific Bioinformatics Conference 2006, 327–336.

    Google Scholar 

  6. Molecular INTeraction database Home[]

  7. IntAct[]

  8. Chatr-aryamontri A, Ceol A, Licata L, Cesareni G: Annotating molecular interactions in the MINT database. Proceedings of the second biocreative challenge evaluation workshop 2007, 55–59.

    Google Scholar 

  9. Krallinger M, Leitner F, Valencia A: Assessment of the Second BioCreative PPI task: Automatic Extraction of Protein-Protein Interactions. Proceedings of the second biocreative challenge evaluation workshop 2007, 41–54.

    Google Scholar 

  10. S O, L MP, H H, R A: The use of common ontologies and controlled vocabularies to enable data exchange and deposition for complex proteomic experiments. Pacific Symposium on Biocomputing 2005, 186–196.

    Google Scholar 

  11. Joachims T: Text categorization with support vector machines: learning with many relevant features. European Conference on Machine Learning 1998, 137–142.

    Google Scholar 

  12. Matthew R, Boutell XS, Luo Jiebo, M Brown C: Learning multi-label scene classification. Pattern Recognition 2004, 37(9):1757–1771.

    Article  Google Scholar 

  13. Rinaldi F, Kappeler T, Kaijurand K: OntoGene In Biocreative II. Proceedings of the second biocreative challenge evaluation workshop 2007, 193–198.

    Google Scholar 

  14. Blei DM, Ng AY, Jordan MI: Latent Dirichlet Allocation. The Journal of Machine Learning Research 2003, 3(2–3):993–1022.

    Google Scholar 

  15. Steyvers M, Griffiths T: Probabilistic Topic Models. In Handbook of Latent Semantic Analysis Edited by: Landauer TK, McNamara DS, Dennis S, Kintsch W, Routledge. 2007, 424–440.

    Google Scholar 

  16. Wei X, Croft W: LDA-based document models for ad-hoc retrieval. Proceedings of the 29th annual international ACM SIGIR 2006, 178–185.

    Google Scholar 

  17. Mimno D, McCallum A: Expertise modeling for matching papers with reviewers. Proceedings of the 13th ACM SIGKDD 2007, 500–509.

    Google Scholar 

  18. Li FF, Perona P: A Bayesian Hierarchical Model for Learning Natural Scene Categories. 2005 IEEE Computer Society Conference on Computer Vision and Pattern Recognition 2005, 524–531.

    Google Scholar 

  19. Buntine WL: Operations for Learning with Graphical Models. Journal of Artificial Intelligence Research 1994, 2: 159–225.

    Google Scholar 

  20. Minka TP: Expectation Propagation for approximate Bayesian inference. Proceedings of the 17th Conference in Uncertainty in Artificial Intelligence 2001, 362–369.

    Google Scholar 

  21. Attias H: A variational Bayesian framework for graphical models. Advances in Neural Information Processing Systems 2000, 209–215.

    Google Scholar 

  22. Andrieu C, de Freitas N, Doucet A, Jordan MI: An Introduction to MCMC for Machine Learning. Machine Learning 2003, 50(1–2):5–43.

    Article  Google Scholar 

  23. PubMed[]

  24. Porter M: An algorithm for suffix stripping. Program 1980, 14(3):130–137.

    Article  Google Scholar 

  25. Chai KMA, Chieu HL, Ng HT: Bayesian online classifiers for text classification and filtering. SIGIR '02: Proceedings of the 26th annual international ACM SIGIR 2002, 97–104.

    Chapter  Google Scholar 

  26. Yang Y, OPedersen J: A Comparative Study on Feature Selection in Text Categorization. Proceedings of the Fourteenth International Conference on Machine Learning 1997, 412–420.

    Google Scholar 

  27. Thorsten J: Learning to Classify Text Using Support Vector Machines. Heidelberg, Germany: Springer; 2002.

    Google Scholar 

  28. Matt Rasmussen, gCluto Home[]

Download references


This work was supported by the Chinese Natural Science Foundation under grant No. 60572084, National High Technology Research and Development Program of China (863 Program) under No. 2006AA02Z321, as well as Tsinghua Basic Research Foundation under grant No. 052220205 and No. 053220002.

This article has been published as part of BMC Bioinformatics Volume 10 Supplement 1, 2009: Proceedings of The Seventh Asia Pacific Bioinformatics Conference (APBC) 2009. The full contents of the supplement are available online at

Author information

Authors and Affiliations


Corresponding author

Correspondence to Xiaoyan Zhu.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors' contributions

Wang carried out the major work of the paper, proposed the model, implemented the experiments and drafted the manuscript. Huang gave directions in the process and revised the draft. Zhu supervised the whole work, gave great amount of valuable suggestions and helped to revise the manuscript. All authors have read and approved the final manuscript.

Rights and permissions

Open Access This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution License ( ), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Wang, H., Huang, M. & Zhu, X. Extract interaction detection methods from the biological literature. BMC Bioinformatics 10 (Suppl 1), S55 (2009).

Download citation

  • Published:

  • DOI: