- Methodology article
- Open Access
Improved multi-level protein–protein interaction prediction with semantic-based regularization
© Saccàet al.; licensee BioMed Central Ltd. 2014
Received: 8 August 2013
Accepted: 3 March 2014
Published: 12 April 2014
Protein–protein interactions can be seen as a hierarchical process occurring at three related levels: proteins bind by means of specific domains, which in turn form interfaces through patches of residues. Detailed knowledge about which domains and residues are involved in a given interaction has extensive applications to biology, including better understanding of the binding process and more efficient drug/enzyme design. Alas, most current interaction prediction methods do not identify which parts of a protein actually instantiate an interaction. Furthermore, they also fail to leverage the hierarchical nature of the problem, ignoring otherwise useful information available at the lower levels; when they do, they do not generate predictions that are guaranteed to be consistent between levels.
Inspired by earlier ideas of Yip et al. (BMC Bioinformatics 10:241, 2009), in the present paper we view the problem as a multi-level learning task, with one task per level (proteins, domains and residues), and propose a machine learning method that collectively infers the binding state of all object pairs. Our method is based on Semantic Based Regularization (SBR), a flexible and theoretically sound machine learning framework that uses First Order Logic constraints to tie the learning tasks together. We introduce a set of biologically motivated rules that enforce consistent predictions between the hierarchy levels.
We study the empirical performance of our method using a standard validation procedure, and compare its performance against the only other existing multi-level prediction technique. We present results showing that our method substantially outperforms the competitor in several experimental settings, indicating that exploiting the hierarchical nature of the problem can lead to better predictions. In addition, our method is also guaranteed to produce interactions that are consistent with respect to the protein–domain–residue hierarchy.
Physical interactions between proteins are the workhorse of cell life and development , and play an extremely important role both in the mechanisms of disease  and in the design of new drugs . In recent years, there has been enormous interest in reverse engineering the protein–protein interaction (PPI) networks of several species, particularly due to the availability of high-throughput experimental techniques, leading to an abundance of large databases on all aspects of PPIs .
Notwithstanding the increased availability of interaction data, the natural question of whether two arbitrary proteins interact, and why, is still open. The growing literature on protein interaction prediction [4–6] is symptomatic of the gap separating the amount of available data and the effective size of the interaction network . The present paper is a contribution towards filling this gap.
Despite the significance of low-level details in elucidating the mechanics of protein–protein interactions, most of the current experimental data comes from high-throughput screening techniques, such as yeast two-hybrid (Y2H) assays . These techniques do not provide information on domain- or residue-level interactions, which require solving the three-dimensional structure of each protein-protein complex, an expensive and time consuming task addressed by X-Ray crystallography, NMR, or electron microscopy techniques . As a consequence, protein–protein interaction data is under-characterized at the domain and residue levels: the current databases are relatively lacking when compared to the magnitude of the existing body of data about protein-level interactions . At the time of writing, the PDB hosts 84,418 structures, but merely 4,210 resolved complexes (according to http://www.rcsb.org/pdb/statistics/holdings.do, retrieved on 2013/06/20). The latter cover only a tiny fraction of the interactions stored in databases such as BioGRID and MIPS.
From a purely biological perspective, predictions at different levels have several important applications. The network topology and individual features of protein interactions are an essential component of a wide range of biological tasks: inferring protein function  and localization , reconstructing signal and metabolic pathways , discovering candidate targets for drug development . Finer granularity predictions at the domain level allow to discover affinities between domain types that can be carried over to other proteins [17, 18]; domain–domain networks have also been assessed as being typically more reliable than their protein counterparts . Finally, residue-level predictions, i.e., interface recognition, enable the detailed study of the principles of protein interactions, and are crucial for tasks such as rational drug design , metabolic reconstruction and engineering , and identification of hot-spots  in the absence of structure information.
Given the usefulness of knowing the details of protein–protein interactions at diverse levels of detail, and based on earlier ideas of Yip et al., in this paper we address the problem of collectively predicting the binding state of all proteins, domains, and residues in a network. We call this task the multi-level protein interaction prediction problem (MLPIP for short).
From a computational point of view, the most important feature of the multi-level prediction problem is its inherently relational nature. Proteins, domains and residues are organized in a hierarchy, which dictates constraints on the binding state of pairs of objects at the different levels, as follows. On the one hand, whenever two proteins are bound, at least two of their domains must also be bound, and, similarly, there must be residues in the two domains that form an interface. On the other hand, if no residues of the two proteins interact, neither do their domains, nor the proteins themselves. In other words, predictions at different levels must be consistent.
In this paper we cast the multi-level prediction problem as a statistical-relational learning task, leveraging the latest developments in the field. Our prediction method is based on Semantic Based Regularization , an elegant semi-supervised prediction framework that caters both the effectiveness of kernel machines and the expressivity of First Order Logic (FOL). The constraints described above are encoded as FOL rules, which are used to enforce consistent predictions at all levels of the interaction hierarchy. By computing multi-level predictions, our method can not only infer which protein pairs are likely to interact, but also provide details about how the interactions take place. Our empirical evaluation shows the effectiveness of this constraint-based approach in boosting predictive performance, achieving substantial improvements over both an unconstrained baseline and the only existing alternative MLPIP method .
PPI networks are most naturally formalized as graphs, where nodes represent proteins and edges represent interactions. Given a set of features describing the properties of the proteins in the network (e.g. primary structure, localization, tertiary structure —when available—, etc.), inferring the PPI network topology amounts to determining those pairs of proteins that are likely to interact. This task is often cast as a pairwise classification problem, where a binary classifier takes as input a pair of proteins (or rather their feature-based representations) and predicts whether they interact or not. Standard binary classification methods, such as Support Vector Machines , can be used to implement the pairwise classifier. In this setting, the interaction depends only on the features of the two incident nodes, and is independent of all other nodes. Interactions between domains or residues can be predicted similarly.
The most straightforward way to address the MLPIP problem is to cast the three interaction prediction problems, for proteins, domains and residues respectively, as independent pairwise classification tasks. However, as previously discussed, these problems are clearly strongly related: two proteins interact via one or more domains, which in turn contain patches of residues that constitute the interaction surface. Ignoring these relationships can lead to heavily suboptimal, inconsistent predictions, where, e.g. two proteins are predicted to interact but none of their domains are predicted to be involved in this interaction. Making these relationships explicit and forcing predictors to satisfy consistency constraints is the key contribution of this work. In the machine learning community, this kind of scenario characterized by multiple related prediction tasks is usually cast as a statistical-relational learning problem [24, 25], where the goal is to collectively classify the state of all objects of interest, taking into account the relations existing between them. The solution we adopt is grounded in this learning framework.
Overview of the proposed method
true iff the protein pair (p,p’) is bound
true iff the domain pair (d,d’) is bound
true iff the residue pair (r,r’) is bound
true iff protein p is parent of domain d
true iff domain d is parent of residue r
true iff protein p is parent of residue r
true iff protein p has at least one domain
true iff domain d has at least one residue
∀(p,p’) hasdom(p) ∧hasdom(p’) ⇒
boundp(p,p’) ⇒ ∃(d,d’) boundd(d,d’) ∧parentpd(p,d) ∧parentpd(p’,d’)
boundd(d,d’) ∧parentpd(p,d) ∧parentpd(p’,d’) ⇒boundp(p,p’)
∀(d,d’) hasres(d) ∧hasres(d’) ⇒
boundd(d,d’) ⇒ ∃ n (r,r’) boundr(r,r’) ∧parentdr(d,r) ∧parentdr(d’,r’)
boundr(r,r’) ∧parentdr(d,r) ∧parentdr(d’,r’) ⇒boundd(d,d’)
Same as D →R, with proteins in place of domains
Same as R →D, with proteins in place of domains
A kernel machine could for instance predict whether a protein is an enzyme or not (binary classification), in terms of weighted similarity to other proteins. Being similar to an enzyme x i will drive the prediction towards the positive (enzyme) class (positive weight w i ), while being similar to a non-enzyme x j will drive the prediction towards the opposite class (negative weight w j ).
Here the kernel function measures the similarity between two pairs of objects, so that, e.g. two proteins will be predicted as interacting if they are similar to protein pairs which are known to interact, and dissimilar from pairs known to not interact.
This transformation guarantees that, if the input function K is a valid kernel, so is the resulting pairwise function.
As already explained, in SBR each target predicate is implemented as a kernel machine, and the state of a predicate for an uncharacterized pair of proteins can be inferred by querying the machine. Positive predictions correspond to true predicates, i.e. bound protein pairs, and negative predictions to false ones. The confidence of the kernel machine, also called margin, embodies the confidence in the state of the predicate, that is, how strongly two proteins are believed to interact (or not).
Given the output of the kernel machines for all target predicates, SBR uses the First Order Logic rules to condition the state of the correlated predicates. It does so by first translating the FOL rules into continuous constraints, which we discuss more thoroughly in Methods. The variables coming into play into the continuous constraints are the confidences of all target predicates (and the state of all given predicates) appearing in the equivalent FOL constraint. The amount of violation is reflected by the value of the continuous constraints: if the predicted predicates satisfy a FOL rule, the corresponding constraint will have a value equal to 1; on the other hand, the closer the constraint value to zero, the more the FOL rule is violated.
where the first term accounts for consistency between inferred truth values and confidence of the individual predictions, and the second incorporates information on the degree of satisfaction of the constraints build from the FOL knowledge. Contrarily to standard kernel methods, this optimization problem is non-convex. This is commonly the case for complex statistical-relational learning tasks , and implies that we are restricted to finding local optima. SBR employs a two-stage learning process to make training effective even in presence of local optima. In particular, the first stage of SBR learning takes into account the fitting of the individual predictions to the supervised data. This learning task is convex and can be efficiently solved. The solution found in the first stage is used as starting point for a second stage, where the FOL knowledge is also considered. This optimization strategy has been experimentally proved to find high-quality solutions without adding the computational burden of other non-convex optimization techniques .
SBR is a semi-supervised method , meaning that the set of target proteins is given beforehand and can be exploited during the learning stage to fine-tune the model. Semi-supervised learning is known to enhance the prediction ability when appropriately used , and can be applied very naturally to PPI prediction, as the full set of proteins is always known.
To summarize, at each level the state of an uncharacterized pair of objects, e.g. proteins p and p′, is mainly inferred by the similarity of the pair (p,p′) to other pairs that are known to interact or not, through the pairwise kernel function K and the learned weights w. Thus the kernel allows to propagate information horizontally within the same level. At the same time, the FOL constraints allow to propagate information vertically between the levels, by keeping the interaction pattern along the protein-domain-residue hierarchy consistent.
Modeling multi-level interactions
The hasdom predicate can be computed directly by SBR using the above definition; we instead pre-compute its value for all protein pairs for run-time efficiency.
The boundp(p,p’) target predicate models the binding state of two distinct proteins. Its state is known for certain protein pairs, i.e. those in the training set, and our goal is to predict its state on the remaining ones. The boundd(d,d’) predicate plays the same role for domains. For a complete list of predicates, see Table 1. For a visualization of the predicates instantiated over a protein pair, see Figure 2-b.
In what follows we describe how to design inter-level FOL constraints to properly enforce consistency between predictions at different levels. We focus on modeling the constraints tying proteins and domains; it is easy to see that the ones between domains and residues can be modelled similarly (with one peculiar exception that will be pointed out later). Table 2 reports the complete list of rules.
Intuitively, the rule means that whenever two proteins are bound (and therefore the left-hand side (LHS) of the implication is true) then there must be at least one pair of child domains that are bound (the right-hand side (RHS) is true). In classical First Order Logic the rule would require that, whenever none of the child domains is bound (the RHS is false), then the parent proteins must not be bound (the LHS is false).
This rule is applied to all protein pairs, demanding that if there is a pair of bound children domains then the proteins must be bound too, and vice versa that if the parent proteins are unbound so are the domains. The P →D and D →P rules could be merged into a single equivalent rule using the double implication (⇔). However, the rules have been considered separately to keep their effects on the results separated and easier to analyze.
To simulate the unidirectional information propagation between levels, as done by Yip et al. (see Related work), we modified how SBR converts logic implications by using the t-norm residuum, which states that a logic implication is true if the RHS is at least as true as the LHS. This modification also removes a bias in the translation of the implication that was affecting the original formulation of SBR, whose effect is to often move the LHS toward the false value. See Methods for details.
The constraints for domains and residues can be similarly defined with one important exception. The P →D rule described above (correctly) requires at least one domain couple to be bound for each interacting protein pair. However, when two domains are bound, the interaction interface involves more than one residue pair: for instance, binding sites collected in the protein–protein docking benchmark version 3.0  consist of 25 residues on average . We integrate this observation in the D →R rule using the n-existential operator ∃ n in place of the regular existential (see Table 2 for the complete formulation), so that whenever two domains are bound, at least n pairs of their residues must be bound. Since interfaces in the employed dataset are typically 5 residues long, n=5 has been used in the experiments. Our results demonstrate that this seemingly small modification has a rather extensive impact on the prediction of domain and residue level interactions.
In this section we briefly summarize previous PPI interaction prediction approaches using methods that are most closely related to the present paper: kernel methods, semi-supervised methods, and logic-based methods. For a broader exposition of interaction prediction methods, please refer to one of the several surveys on the subject [4, 6, 9, 32].
The earliest attempt to employ kernel methods  for PPI prediction is the work of Bock et al., which casts interaction prediction as pairwise classification, using amino-acid composition and physico-chemical properties alone. Ben-Hur et al. extended the previous work by applying pairwise kernels and combining multiple data sources (primary sequence, Pfam domains, Gene Ontology annotations and interactions between orthologues). Successive publications focused primarily on aggregating more diverse sources, including phylogenetic profiles, genetic interactions, and subcellular localization and function . Kernel machines have also been applied to the prediction of binding sites from sequence, as summarized in . The appeal of supervised kernel methods is that they provide a proved and theoretically grounded set of techniques that can easily integrate various information sources, and can naturally handle noise in the data. However, they have two inherent limitations: (i) the binding state of two proteins is inferred independently from the state of all other proteins, and (ii) due to their supervised nature, they do not take advantage of unsupervised data, which is very abundant in the biological network setting.
Semi-supervised learning (SSL) techniques [28, 29] attempt to solve these issues. In the SSL setting the set of target proteins is known in advance, meaning that the learning algorithm has access to their distribution in feature space. This way the inference task can be simplified by introducing unsupervised constraints that assign the same label to proteins that are, e.g., close enough in feature space, or linked in the interaction network, instantiating a form of information propagation. There are several works in the PPI literature that embed the known network topology using SSL constraints. Qi et al. employ SSL methods to the special case of viral-host protein interactions, where supervised examples are extremely scarce. Using similar methods, You et al. attempt to detect spurious interactions in a known network by projecting it on a low-dimensional manifold. Other studies [36, 37] applied SSL techniques to the closely related problems of gene–protein and drug–protein interaction prediction. Despite the ability of SSL to integrate topology information, no study so far has applied it to highly relational problems such as the MLPIP.
An alternative strategy for interaction prediction is Inductive Logic Programming (ILP) , a group of logic-based formalisms that extract rules explaining the likely underlying causes of interactions. ILP methods were studied in the work of Tran et al. using a large number of features: SWISS-PROT keywords and enzyme properties, Gene Ontology functional annotations, gene expression, cell cycle and subcellular localization. Further advances in this direction, with a special focus on using domain information, can be found in [17, 18]. The advantage of ILP methods over purely statistical methods is that they are inherently able to deal with relational information, making them ideal candidates for solving the MLPIP problem. Alas, contrary to kernel methods, they tend to be very susceptible to noise, which is a very prominent feature of interaction dataset, and are less effective in exploiting complex feature representations, e.g. involving highly non-linear interactions between continuous features.
Recently, some works highlighted the importance of the multi-level nature of protein–protein interactions. Gonzalez et al. propose a method to infer the residue contact matrix from a known set of protein interactions using SVMs; on the contrary, our goal is to predict the interactions concurrently at all levels of the hierarchy. Another study  highlights the relevance of domain-level interactions, and the unfortunate lack of details thereof, and formulates a method to reinterpret a known PPI network in terms of its constituent domain interactions; the present work has a different focus and a more general scope.
Most relevant to this paper is the work of Yip et al., where the authors propose a procedure to solve the MLPIP problem based on a mixture of different techniques. The idea is to decompose the problem as a sequence of three prediction tasks, which are solved iteratively. Given an arbitrary order of the three levels (e.g. proteins first, then domains, then residues), their procedure involves computing putative interactions in the first level (in this case proteins), then using the most confident predictions as novel training examples at the following level (i.e., domains). The procedure is repeated until a termination criterion is met.
Intra-level predictions are obtained with Support Vector Regression (SVR) . In particular, each object has an associated SVR machine that models its propensity to bind any other object in the same level. The extrapolated values act as confidences for the predictions themselves. The mechanism for translating the most confident predictions at one level into training examples for the next level depends on the relative position of the two levels in the hierarchy. Downward propagation (e.g. from proteins to domains) simply associates to each novel example the same confidence as the parent prediction: in other words, if two proteins are predicted as bound with high confidence, all their domains will be considered bound with the same confidence. Upward propagation (e.g. from domains to proteins) is a bit more involved: the confidence assigned to the novel example (protein) is a noisy-OR combination of confidences for all the involved child objects (domains).
While this method has been shown to work reasonably well, it is afflicted by several flaws. First of all, while the iterative procedure is grounded in co-training , the specific choice of components is not as theoretically sound. For instance, the authors apply regression techniques on a classification task, which may lead to sub-optimal results. The inter-level example propagation mechanisms are ad hoc, do not exploit all the information at each level (only the most confident predictions are propagated), and are designed to merely propagate information between levels, not to enforce consistency on the predictions. In particular, the downward propagation rule is rather arbitrary: it is not clear why all domains of bound proteins should be themselves bound with the same confidence. Finally, these rules, which are intimately tied to the specific implementation, are not defined using a formal language, and are therefore difficult to extend. For instance, it would be difficult to implement in said framework something similar to an n-existential propagation rule, which is extremely useful for dealing with residue interactions.
Semantic Based Regularization seems to have many obvious advantages in this context. A first advantage is that it decouples the implementation of the functions from how consistency among levels is defined. Indeed, consistency is implemented via a set of constraints, which are applied over the output of the predictors. However, there is no limitation to which kind of predictors are used. For example, we used kernel machines as basic machinery for implementing the predictor, where different state-of-the-art kernels can be used at the single levels, while still be able to define a single optimization problem.
Furthermore, SBR allows to natively propagate the predictions of one level to the other levels. Since the predictions and not the supervisions are propagated, SBR accuracy can get advantage of the abundant unsupervised data. The availability of an efficient implementation of the n-existential quantifier is also a crucial advantage: if two proteins or domains are interacting, a small set of residues must be interacting as well. SBR does not simply propagate a generic prior to all the residues for a protein or domain, which could decrease accuracy of the reductions for the negative supervisions. SBR instead performs a search process in order to select a subset of residue candidates, where to enforce the interaction. As shown in the experimental results, this greatly improves residue prediction accuracy. Finally, the circular dependencies that make learning difficult are dealt in the context of a general and well defined framework, which implements various heuristics to make training effective.
Results and discussion
In this work we use the dataset of Yip et al., described here for completeness. The dataset represents proteins, domains and residues using features gathered from a variety of different sources:
Protein features include phylogenetic profiles derived from COG, subcellular localization, cell cycle and environmental response gene expression; protein-pair features were extracted from Y2H and TAP-MS data. The gold standard of positive interactions was constructed by aggregating experimentally verified or structurally determined interactions taken from MIPS, DIP, and iPfam.
At the domain level, the dataset includes both features for domain families and for domain instances based on frequencies of domains within one or more species and phylogenetic correlations of Pfam alignments. The gold standard of positive interactions was built from 3D structures of complexed proteins taken from PDB.
Residue features consist of sequence-based properties, namely charge complementarity, Psi-Blast  profiles, predicted secondary structure, and predicted solvent accessibility.
Kernels computed from the individual features were combined additively into a single kernel function for each level, and then transformed into pairwise kernels using Equation (1); the resulting functions were used as inputs to SBR. A visualization of the process can be found in Figure 2.
This procedure yields a dataset of 1681 proteins, 2389 domains, and 3035 residues, with a gold standard of 3201 positive (interacting) protein pairs, 422 domain pairs, and 2000 residue pairs. Since interaction experiments can not determine which pairs do not interact, the gold standard of negative pairs is built by randomly sampling, at each level, a number of pairs that are not known to interact (i.e. not positive). This is a common approach to negative labeling in the PPI prediction literature . To keep the dataset balanced, the number of sampled negative pairs is identical to the number of objects in the gold standard of positives. For more details on the dataset preparation, please refer to . We further refined the dataset by running CD-HIT  with a 20% sequence similarity threshold, identifying 23 redundant proteins. These proteins were not used when comparing the method performances.
Turning our attention to the resulting dataset, we note that most of the supervision is located at the protein level: out of all possible interactions between pairs of proteins, which are (1681×1680)/2, 0.226% are known (either positive or negative). On the contrary, the other levels hold much less information: only 0.042% of all possible residue pairs, and 0.014% of all possible domain pairs, are in the dataset. The low number of residue pairs is due to i) different requirements for experimentally determining the interactions at the three levels, i.e. whether the structure is available; and ii) sampling choices operated by Yip et al..
In this work we compare our method to that of Yip et al., where the authors evaluated their method using a 10-fold cross-validation procedure. To keep the comparison completely fair, we repeated said procedure with SBR, reusing the very same train/test splits. Since correlated objects, e.g. a protein and its domains/residues, share information, the folds were structured as to avoid such information to leak between train and test folds: this was achieved by keeping correlated objects in the same fold. In order not to bias the performance estimates, all redundant proteins were ignored, along with their domains and residues, when computing the results of both SBR and the method of Yip et al. The full experimental setup and instructions to replicate the experiments can be downloaded at http://sites.google.com/site/semanticbasedregularization/home/software/protein_interaction.
SBR has two scalar hyper-parameters that control the contribution of various parts of the objective function: λ c is the weight associated to the constraints (how much the current solution is consistent with respect to the rules) and λ r , which controls the model complexity (see the Methods section for more details). The λ r parameter was optimized on the first fold by training the model without the logic rules and it was then kept fixed for all the folds of the k-fold cross-validation. The resulting value is λ r =0.1. The λ c parameter has not been optimized and kept fixed at λ c =1. Please note that further significant gains for the proposed method could be achieved by fine-tuning this meta-parameter. However, since the dataset from Yip et al. does not include a validation split, no sound way to optimize this parameter was possible without looking at the test set, or redefining the splits (making it difficult to compare against the results of Yip et al.). Therefore, we decided to not perform any tuning for this meta-parameter.
We computed three performance metrics: the Receiver Operating Characteristic (ROC) curve, the area under the ROC (AUCROC, or AUC for short), and the F1 score. The ROC curve represents the relation between the false positive rate (FPR) and the true positive rate (TPR), and can be seen as the proportion of true positives gained by “paying” a given proportion of false positives. By definition, the ROC curve is monotonically non-decreasing; the steeper the curve, the better the predictions. The AUC measures the ability to correctly discriminate between positives and negatives, or alternatively, the ability to rank positives above negatives. It is independent of any classification threshold, and thus particularly fit to evaluate models over the whole spectrum of possible decision thresholds. The F1 score is the harmonic mean of precision and recall. Contrary to the AUC, the F1 takes into account the predicted label, but not its confidence (margin).
P → D
D → R
P → R
P ⇔ D
D ⇔ R
P ⇔ R
Results for Yip et al.
Results for SBR
Results for SBR- ∃ n
Results ( F 1 )
P → D
D → R
P → R
P ⇔ D
D ⇔ R
P ⇔ R
Results for Yip et al.
Results for SBR
Results for SBR- ∃ n
To evaluate the effects of the constraints on the performances of SBR, we performed three independent experiments using rules of increasing complexity. This setup follows closely that of Yip et al..
As a baseline, we estimate the performance of our method when constraints are ignored. This is equivalent to the method of Yip et al. when no information flow between levels is allowed. The results can be found in the “Independent” column of Tables 3 and 4.
In absence of constraints SBR reduces to standard ℓ2-regularized SVM classification: learning and inference become convex problems, and the method computes the globally optimal solution. Thus, the only differences between our method and the competitor are: (i) using classification versus regression, and (ii) using pairwise classification, instead of training a single model for each entity (protein, domain, residues). These differences alone produce a substantial increase in performance: the F1 changes by about +0.05 in all three cases. The AUC of proteins and domains is improved by about +0.09 and +0.07, respectively, while residues are less affected, with a +0.03 difference.
In the second experiment, we evaluate the effect of introducing unidirectional constraints between pairs of levels. In the P →D case only the P →D rule is active, meaning that bound protein pairs enforce positive domain pairs and negative domain pairs enforce negative protein pairs. The D →R and P →R cases are defined similarly. In all three cases, the level not appearing in the rule (e.g. the residue level in the P →D case) is predicted independently. This setup makes it easy to study the effects of propagating information from one level to the other without interferences. The results can be found in the “Unidirectional” column of Tables 3 and 4. In the same column we also show the results for Yip et al. for the unidirectional flow setting, where examples are propagated from one level to the next but not vice versa. However, since the competitor’s algorithm is iterative, information about lower levels can indeed affect the upper levels in successive iterations.
The results show that introducing unidirectional constraints in SBR improves the predictions in all cases. In particular, using (predicted and known) protein interactions helps inferring correct domain interactions, which improve by about +0.13 F1 and +0.2 AUC (P →D case). Residues improve independently of whether protein or domain-level information is used, with a +0.03 F1 in both cases, and a +0.08/ +0.07 AUC difference, respectively. Interestingly, proteins tend to help residue predictions slightly more than domains, despite the indirection between the two levels; this is likely an effect of the larger percentage of supervised pairs available.
Compared to SBR, the method of Yip et al. does not benefit as much from unidirectional information flow. Protein-level information allows to improve domain predictions only (+0.1 F1, +0.06 AUC for P →D), while residue predictions are worse than in the independent case (-0.05 and -0.04 AUC, and -0.01 F1, in the D →R and P →R cases, respectively).
In the third experiment we study the impact of using bidirectional constraints between pairs of levels; the level not appearing in the rules is predicted independently, as above. In the P ⇔D case, both the P →D and D →P rules are active, meaning that the protein and domain levels are enforced to be fully consistent; the P ⇔R and D ⇔R cases are defined analogously. This experiment is comparable to the bidirectional flow setting of Yip et al.. The results can be found in the “Bidirectional” column of Tables 3 and 4.
We observe that the new constraints have a positive effect on predictions at all three levels: proteins change from 0.808 AUC to 0.820, domains from 0.814 to 0.896 and residues from 0.671 to 0.673. In terms of F1, the changes are from 0.718 to (up to) 0.722 for proteins, from 0.693 to 0.731 for domains, and no change for residues. The change is not as marked as between the independent and unidirectional experiments. In particular, domains see the largest increase in performance (+0.08 AUC, +0.04 F1), in particular thanks to the contribution of residue-level information, which is more abundant. Proteins and residues are less affected. The result is unsurprising for proteins, which hold most of the supervision and are thus (i) more likely to be predicted correctly in the independent setting, and (ii) less likely to be assisted from hints coming from the other, less supervised levels.
As for the method of Yip et al., the bidirectional flow mostly affects the domain and residue levels, whose improvement is +0.07 AUC/ +0.04 F1 and +0.11 AUC/ +0.09 F1, respectively; the change for protein interactions is negligible. Regardless of the relative performance increase, SBR is able to largely outperform the competitor in all configurations except one (F1 of the P ⇔R case for residues).
We note that the fact that all three cases (P ⇔D, P ⇔R and D ⇔R) improve over both the independent and the unidirectional experiments shows that not only the bidirectional constraints are in fact sound, but also that, despite the increased computational complexity, SBR is still able to exploit them appropriately.
In the final experiment we activate the P →D, D →P, D →R and R →D rules, as defined in Table 2, making all levels interact. This is the most complex setting, and produces fully consistent predictions through the hierarchy. It is comparable to the “PDR” bidirectional setting of Yip et al.. The AUC scores can be found in column “Full” of Tables 3 and 4.
In this experiment the P →R and R →P constraints are not used. Direct information flow between proteins and residues is not needed, because it would be redundant: from a formal logic point of view, this corresponds to the observation that the logic rule expressing protein to residue consistency is implied by the other consistency rules. Indeed, we have experimentally verified that adding this propagation flow does not significantly affect the results.
In this experiment, protein predictions are stable with respect to the previous experiments, confirming the intuition that the abundance of supervision at this level makes it less likely to benefit from predictions at the other ones. On the contrary, domains see a large performance upgrade, from 0.896 to 0.937 AUC and from 0.731 to 0.750 F1, when made to interact with both proteins and residues. The change for residues is instead only marginal.
The results for Yip et al. are mixed, with proteins faring almost identically to the previous experiment, domains showing a slight drop in AUC but an equally slight increase in F1, and residues improving in AUC (+0.08) but not in F1 (unchanged) over the bidirectional P ⇔R case. The improvement in residue prediction (in terms of AUC) stands in contrast with the results of SBR, and is the only case in which the method of Yip and colleagues works better than SBR. The issue lies within our formulation of the D →R rule: whenever two domains are bound, the rule is satisfied when at least one residue pair is bound. As already mentioned above, this is not realistic: protein interfaces span more than two residues, typically five or more. We therefore extended SBR to support the n-existential quantifier, which allows to reformulate the D →R rule to take this observation into account (see the Methods section for more details on the n-existential quantifier). The new D →R rule, shown in Table 2, requires for each pair of bound domains at least n=5 residues to be bound. We chose the constant n=5 to be both realistic and, since the computational cost increases with n, small enough to be easily computable. We applied the same modification to the P →R rule.
Summing up, these results highlight the ability of SBR to enforce constraints even with highly complex combinations of rules, allowing the modeler to fully exploit the flexibility and performance improvement offered by non-standard FOL extensions like the n-existential operator.
The results presented in the previous section offer a clear perspective on the advantages of the proposed method. By employing appropriate classification techniques and training a single global pairwise model per level, rather than relying on the less than optimal local (per-object) regression models of Yip et al., a considerable improvement was achieved even in the unconstrained experiment. Furthermore, when enforcing consistency among the protein, domain and residue levels and using the n-existential quantifier, the experimental results are significantly better than both the unconstrained baseline and the corresponding results of Yip and colleagues, at all levels and in all experimental settings.
It is worth noting that SBR performance improves monotonically with the increase of constraint complexity in the reported experiments. This result is far from obvious, and confirms both that the biologically-motivated knowledge base is useful, and that SBR is able to effectively apply it. In contrast, the competitor’s method does not always improve in a similar manner.
In general, the performance gain brought forth by inter-level propagation is not homogeneously distributed between the three levels. We register a large improvement for domains and residues, especially when SBR is used in conjunction with the n-existential quantifier. Proteins are less affected by consistency enforcement, most likely due to the availability of more supervised examples.
We note that the FOL rules have a twofold effect. Firstly, they propagate information between the levels, enabling predicted interactions at one level to help inferring correct interactions at the other two levels. This is especially clear in the “Full” experiment with the n-existential quantifier: in this case, better residue level predictions increase the overall quality of domain predictions as well. Secondly, the rules also guarantee that the predictions are consistent along the object hierarchy.
Summarizing, SBR is able to largely outperform that of Yip and colleagues, and moreover enforces the predictions to be consistent among levels. As previously mentioned, the data taken from Yip et al. has some peculiarities worth discussing. First, it contains a low number of residue–residue interactions, partially due to design choices taken by . Second, it is artificially balanced by including an appropriate number of non-interactions, while in a real-world case all possible pairs would qualify as candidates. We decided to keep the dataset as-is in order to facilitate a fair comparison with Yip et al.. We postpone further analysis with other datasets to future work.
In this work we tackle the multi-level protein interaction prediction (MLPIP) problem, first introduced by Yip et al., which requires to establish the binding state of all uncharacterized pairs of proteins, domains and residues. Contrary to standard protein–protein interaction prediction, the MLPIP problem offers many advantages and opens up new challenges. The primary contribution of this paper is the extension and application to the MLPIP task of a state-of-the-art statistical relational learning technique called Semantic Based Regularization.
SBR is a flexible framework to inject domain knowledge into kernel machines. In this paper SBR has been used to tie together protein, domain and residue interaction predictions tasks. In particular, the domain knowledge expresses that two proteins interact if and only if there is an interaction between at least one pair of domains of the proteins. Similarly two domains can interact if and only if there are at least some residues interacting. While these tasks could be learned separately, tying them together has multiple advantages. First the predictions will be consistent and more accurate, as the predictions at one level will help the predictions at the other levels. Secondly, the domain knowledge can be enforced also on the unsupervised data (proteins, domains and residues for which interactions are unknown). Unsupervised data is typically abundant in protein interaction prediction tasks but often neglected. This methodology allows to powerfully leverage it, significantly improving the prediction accuracy. Note also that, while the resolved complexes are required during the training stage, no structural information is required for performing inference on novel proteins.
While other work in the literature has exploited the possibility of tying the predictions at multiple levels, the presented methodology employs a more principled inference process among the levels, where the domain knowledge can be exactly represented and precisely enforced. The experimental results confirm the theoretical advantages by showing significant improvements in domain and residue interaction prediction accuracy both with respect to approaches performing independent predictions and the only previous approach attempting at linking the prediction tasks.
Given the flexibility offered by SBR, the proposed method can be extended in several ways. The simplest extension involves engineering a more refined rule set, for instance by introducing (soft) constraints between the binding state of consecutive residues, which are likely to share the same state. More ambitious goals, requiring a redesign of the experimental dataset, include encoding selected information sources, such as domain types, subcellular co-localization and Gene Ontology annotations, as First Order Logic constraints rather than with kernels, to better leverage their relational nature.
Machine learning and statistical methods are very well defined for the linear case, and statistical learning theory can provide optimal solutions in terms of generalization performance. Unfortunately, non-linearity is often required in order to solve most applications, where exploiting complex dependencies is essential to predict some higher level property of the objects. Kernel methods try to combine the potential classification power of non-linear methods and the optimality and computational efficiency of linear methods by mapping the input patterns into a high dimensional feature space, where parameter optimization remains linear.
where x is the representation of the pattern, K(x,x i )=<Φ(x),Φ(x i )> is a kernel function, where Φ(·) is some mapping from the input space to the feature space. Intuitively, the kernel function measures the similarity between pairs of instances, and the prediction f(x) for a novel instance is computed as a weighted similarity to training instances x i . There is a large body of literature on kernel machines, see e.g.  for an introduction.
reduces to the formulation of hard margin ℓ2SVMs  if L(·) is the hinge loss and λ→∞. A very similar cost function has been employed to solve the protein, domain, residue interaction presented in this paper.
Propositional logic is based on the basic concept of propositions, which can assume either a true or false value. It is possible to perform operations on the propositions by connecting them via the and (∧), or (∨) and not (¬) operators. In particular, given two propositions A,B, it holds that: A∧B is true iff A=t r u e,B=t r u e, A∨B is false iff A=f a l s e,B=f a l s e and ¬A flips the current truth value of A. The operator ⇒ can be used to express a conditional statement: A⇒B expresses the fact that B must hold true if A is true. The sentence A⇒B is false iff A=t r u e and B=f a l s e, and it can be expressed in terms of other operators through the equivalence A⇒B≡¬A∨B.
First-order logic (FOL) extends propositional logic to compactly express generic properties for a class of objects, thanks to the use of predicates, variables and quantifiers. A variable can assume as value any object in some considered domain. A variable is said to be grounded once it is assigned a specific object. A predicate is a function that, taking as input some objects (or grounded variables), returns either true or false. Predicates can be connected with other predicates using the same operators defined for propositional logic. The universal quantifier (∀) expresses the fact that some proposition is true for any object, while the existential quantifier (∃) expresses the fact that some proposition is true for at least one object.
For example, let x be a variable and let Protein(x),Enzyme(x),NonEnzyme(x) indicate three predicates expressing whether, given a grounding x=PDB1a3a, PDB1a3a is a protein, an enzyme, a non-enzyme, respectively. The following FOL clause can be used to express that any protein is either an enzyme or it is not: ∀x Protein(x)⇒ Enzyme(x)∨NonEnzyme(x).
Variables and quantifiers can be combined. For example, given the predicates Protein(x) holding true if x is a protein and ResidueOf(x, y) holding true if y is a residue of x, the following clause expresses the fact that any protein has a at least one residue: ∀x Protein(x)⇒∃y ResidueOf(x,y).
Semantic Based Regularization (SBR)  is a general framework for injecting prior knowledge expressed in FOL into kernel machines for semi-supervised learning tasks. The prior knowledge is converted into a set of continuous constraints, which are enforced during training. The SBR framework is very general and allows to employ the full expressiveness of FOL in the definition of the prior knowledge. The SBR framework also allows to perform collective classification on the test set, in order to enforce the output to respect the logic knowledge.
where L(·) is a loss function measuring the distance of the function output from the desired one and is the considered sample of data points over which the functions are evaluated. In the experimental setting, L(·) has been set to be the hinge function. It is possible to extend the Representer Theorem to show that the best solution for Equation 3 can be expressed as a kernel expansion as showed in Equation 2 .
where G k , w k , f k =G k w k and y k are the gram matrix, the weights, the function values over the data sample and the desired output column vectors for the patterns in the domain of the k-th task. Evidence tasks do not need to be approximated as their are fully known.
Optimization of the w k parameters for the cost function in Equation 4 can be done using gradient descent. Constraint ϕ h are non-linear in most interesting cases like the one presented in this paper. Therefore, the cost function can present multiple local minima, making optimization difficult. SBR uses a two-step heuristic to solve this problem: first it computes the theoretically global optimum for all predicates independently (setting λ c =0), which are convex kernel machines. Then, it introduces the constraints and proceeds to find a good solution using a gradient descent.
In the following we show how to express first order logic clauses in terms of constraints ϕ h .
Translation of first-order logic clauses into real-valued constraints
Please note that the quantifier-free part of the expression is equivalent to an assertion in propositional logic for any given grounding of the quantified variables. As studied in the context of fuzzy logic and symbolic AI, different methods can be used for the conversion of a propositional expression into a continuous function with [ 0,1] input variables.
Please note that the t-norm behaves as classical logic when the variable approaches the value 0 (false) or 1 (true).
The residuum allows to relax the condition of satisfaction for the implication: an implication is satisfied if the right end side of the implication is more verified than the pre-condition on the left side. This makes the fuzzy or probabilistic inference process easier and better defined. While the original SBR formulation represents implications using modus ponens, the minimum t-norm residuum has been used in the experimental section of this paper to convert implications.
where the predicates have been substituted by the unknown corresponding functions and x1,x2 are the representations of the objects identified by the grounded variables v1,v2, respectively. The representation of the object must be compatible with what is accepted as input by the kernels used by the predicate approximations f A ,f B . For example, they can be a vector of real valued features when using a linear or Gaussian kernel, or graph or tree representations when using kernel for structures. It is also possible to use the methodology when no explicit representations are known, but only the kernel values for each pair of input objects.
where f A (x)f B (x) is the t-norm generalization of the propositional expression A(v)∧B(v) for a given grounding of v.
where indicates the n assignments of x that maximize the value of 1-t E (·) over the set . The conversion of the ∃ n operator consistently reduces to the ∀ conversion when , and to the conversion of the ∃ operator when n=1.
This research was supported by grant PRIN 2009LNP494 (Statistical Relational Learning: Algorithms and Applications) from Italian Ministry of University and Research and by a Google Faculty Research Award (Integrated Prediction of Protein Function, Interactions and Pathways with Statistical Relational Learning).
- Keskin O, Gursoy A, Ma B, Nussinov R: Principles of protein-protein interactions: what are the preferred ways for proteins to interact?. Chem Rev. 2008, 108 (4): 1225-1244.View ArticlePubMedGoogle Scholar
- Hopkins AL: Network pharmacology: the next paradigm in drug discovery. Nat Chem Biol. 2008, 4 (11): 682-690.View ArticlePubMedGoogle Scholar
- Csermely P, Korcsmáros T, Kiss HJ, London G, Nussinov R: Structure and dynamics of molecular networks: A novel paradigm of drug discovery. Pharmacol Ther. 2013, 138 (3): 333-408.View ArticlePubMed CentralPubMedGoogle Scholar
- Tuncbag N, Kar G, Keskin O, Gursoy A, Nussinov R: A survey of available tools and web servers for analysis of protein–protein interactions and interfaces. Brief Bioinform. 2009, 10 (3): 217-232.View ArticlePubMed CentralPubMedGoogle Scholar
- Lewis AC, Saeed R, Deane CM: Predicting protein–protein interactions in the context of protein evolution. Mol Biosyst. 2010, 6: 55-64.View ArticlePubMedGoogle Scholar
- Skrabanek L, Saini HK, Bader GD, Enright AJ: Computational prediction of protein–protein interactions. Mol Biotechnol. 2008, 38: 1-17.View ArticlePubMedGoogle Scholar
- Wass MN, David A, Sternberg MJ: Challenges for the prediction of macromolecular interactions. Curr Opin Struct Biol. 2011, 21 (3): 382-390.View ArticlePubMedGoogle Scholar
- Cary MP, Bader GD, Sander C: Pathway information for systems biology. FEBS Lett. 2005, 579 (8): 1815-1820.View ArticlePubMedGoogle Scholar
- Shoemaker BA, Panchenko AR: Deciphering protein–protein interactions. Part II. Computational methods to predict protein and domain interaction partners. PLoS Comput Biol. 2007, 3 (4): e43-View ArticlePubMed CentralPubMedGoogle Scholar
- Ezkurdia I, Bartoli L, Fariselli P, Casadio R, Valencia A, Tress ML: Progress and challenges in predicting protein–protein interaction sites. Brief Bioinform. 2009, 10 (3): 233-246.View ArticlePubMedGoogle Scholar
- Fields S, Song O: A novel genetic system to detect protein protein interactions. Nature. 1989, 340 (6230): 245-246.View ArticlePubMedGoogle Scholar
- Shoemaker BA, Panchenko AR: Deciphering protein–protein interactions. Part I. experimental techniques and databases. PLoS Comput Biol. 2007, 3 (3): e42-View ArticlePubMed CentralPubMedGoogle Scholar
- Memiševič V, Wallqvist A, Reifman J: Reconstituting protein interaction networks using parameter-dependent domain-domain interactions. BMC Bioinformatics. 2013, 14: 154-View ArticlePubMed CentralPubMedGoogle Scholar
- Radivojac P, Clark WT, Oron TR, Schnoes AM, Wittkop T, Sokolov A, Graim K, Funk C, Verspoor K, Ben-Hur A, Pandey G, Yunes JM, Talwalkar AS, Repo S, Souza ML, Piovesan D, Casadio R, Wang Z, Cheng J, Fang H, Gough J, Koskinen P, Törönen P, Nokso-Koivisto J, Holm L, Cozzetto D, Buchan DW, Bryson K, Jones DT, Limaye B, et al: A large-scale evaluation of computational protein function prediction. Nat Methods. 2013, 10 (3): 221-227.View ArticlePubMed CentralPubMedGoogle Scholar
- Jiang JQ, Wu M: Predicting multiplex subcellular localization of proteins using protein-protein interaction network: a comparative study. BMC Bioinformatics. 2012, 13 (Suppl 10): S20-View ArticlePubMed CentralPubMedGoogle Scholar
- Zhao XM, Wang RS, Chen L, Aihara K: Uncovering signal transduction networks from high-throughput data by integer linear programming. Nucleic Acids Res. 2008, 36 (9): e48-e48.View ArticlePubMed CentralPubMedGoogle Scholar
- Nguyen TP, Ho TB: Discovering signal transduction networks using signaling domain-domain interactions. Genome Inform. 2006, 17 (2): 35-45.PubMedGoogle Scholar
- Nguyen TP, Ho TB: An integrative domain-based approach to predicting protein–protein interactions. J Bioinform Comput Biol. 2008, 6 (06): 1115-1132.View ArticlePubMedGoogle Scholar
- Pitkänen E, Rousu J, Ukkonen E: Computational methods for metabolic reconstruction. Curr Opin Biotechnol. 2010, 21: 70-77.View ArticlePubMedGoogle Scholar
- Tuncbag N, Gursoy A, Keskin O: Identification of computational hot spots in protein interfaces: combining solvent accessibility and inter-residue potentials improves the accuracy. Bioinformatics. 2009, 25 (12): 1513-1520.View ArticlePubMedGoogle Scholar
- Yip KY, Kim PM, McDermott D, Gerstein M: Multi-level learning: improving the prediction of protein, domain and residue interactions by allowing information flow between levels. BMC Bioinformatics. 2009, 10: 241-View ArticlePubMed CentralPubMedGoogle Scholar
- Diligenti M, Gori M, Maggini M, Rigutini L: Bridging logic and kernel machines. Mach Learn. 2012, 86: 57-88.View ArticleGoogle Scholar
- Cortes C, Vapnik V: Support-vector networks. Mach Learn. 1995, 20 (3): 273-297.Google Scholar
- Getoor L, Taskar B: Introduction to Statistical Relational Learning. 2007Google Scholar
- Probabilistic Inductive Logic Programming - Theory and Applications, Volume 4911 of Lecture Notes in Computer Science. Edited by: Raedt LD, Frasconi P, Kersting K, Muggleton S. 2008Google Scholar
- Hofmann T, Schölkopf B, Smola AJ: Kernel methods in machine learning. Ann Stat. 2008, 36 (3): 1031-1508.View ArticleGoogle Scholar
- Ben-Hur A, Noble WS: Kernel methods for predicting protein–protein interactions. Bioinformatics. 2005, 21 (suppl 1): i38-i46.View ArticlePubMedGoogle Scholar
- Chapelle O, Schölkopf B, Zien A: Semi-Supervised Learning, Volume 2. 2006, MIT pressView ArticleGoogle Scholar
- Zhu X: Semi-supervised learning literature survey. Computer Science, University of Wisconsin-Madison. 2006, 2: 3-[http://pages.cs.wisc.edu/jerryzhu/research/ssl/semireview.html],Google Scholar
- Hwang H, Pierce B, Mintseris J, Janin J, Weng Z: Protein–protein docking benchmark version 3.0. Proteins: Struct, Funct, Bioinf. 2008, 73 (3): 705-709.View ArticleGoogle Scholar
- Li B, Kihara D: Protein docking prediction using predicted protein-protein interface. BMC Bioinformatics. 2012, 13: 7-View ArticlePubMed CentralPubMedGoogle Scholar
- Qi Y, Noble WS: Protein interaction networks: protein domain interaction and protein function prediction. Handbook of Statistical Bioinformatics. 2011, 427-459.View ArticleGoogle Scholar
- Bock JR, Gough DA: Predicting protein–protein interactions from primary structure. Bioinformatics. 2001, 17 (5): 455-460.View ArticlePubMedGoogle Scholar
- Qi Y, Tastan O, Carbonell JG, Klein-Seetharaman J, Weston J: Semi-supervised multi-task learning for predicting interactions between HIV-1 and human proteins. Bioinformatics. 2010, 26 (18): i645-i652.View ArticlePubMed CentralPubMedGoogle Scholar
- You ZH, Lei YK, Gui J, Huang DS, Zhou X: Using manifold embedding for assessing and predicting protein interactions from high-throughput experimental data. Bioinformatics. 2010, 26 (21): 2744-2751.View ArticlePubMed CentralPubMedGoogle Scholar
- Xia Z, Wu LY, Zhou X, Wong ST: Semi-supervised drug-protein interaction prediction from heterogeneous biological spaces. BMC Syst Biol. 2010, 4 (Suppl 2): S6-View ArticlePubMed CentralPubMedGoogle Scholar
- Nguyen TP, Ho TB: Detecting disease genes based on semi-supervised learning and protein–protein interaction networks. Artif Intell Med. 2012, 54: 63-71.View ArticlePubMedGoogle Scholar
- De Raedt L: Inductive Logic Programming. 2010View ArticleGoogle Scholar
- Tran TN, Satou K, Ho TB: Using inductive logic programming for predicting protein-protein interactions from multiple genomic data. Knowledge Discovery in Databases: PKDD 2005. 2005, 321-330.View ArticleGoogle Scholar
- González AJ, Liao L, Wu CH: Prediction of contact matrix for protein–protein interaction. Bioinformatics. 2013, 29 (8): 1018-1025.View ArticlePubMed CentralPubMedGoogle Scholar
- Smola AJ, Schölkopf B: A tutorial on support vector regression. Stat Comput. 2004, 14 (3): 199-222.View ArticleGoogle Scholar
- Blum A, Mitchell T: Combining labeled and unlabeled data with co-training. Proceedings of the Eleventh Annual Conference on Computational Learning Theory. 1998, ACM, 92-100.View ArticleGoogle Scholar
- Altschul SF, Madden TL, Schäffer AA, Zhang J, Zhang Z, Miller W, Lipman DJ: Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res. 1997, 25 (17): 3389-3402.View ArticlePubMed CentralPubMedGoogle Scholar
- Gomez SM, Noble WS, Rzhetsky A: Learning to predict protein–protein interactions from protein sequences. Bioinformatics. 2003, 19 (15): 1875-1881.View ArticlePubMedGoogle Scholar
- Li W, Godzik A: Cd-hit: a fast program for clustering and comparing large sets of protein or nucleotide sequences. Bioinformatics. 2006, 22 (13): 1658-1659.View ArticlePubMedGoogle Scholar
- Schölkopf B, Smola AJ: Learning with Kernels. 2002Google Scholar
- Klement E, Pap E, Mesiar R: Triangular Norms. 2000View ArticleGoogle Scholar
- Klir G, Yuan B: Fuzzy Sets and Fuzzy Logic: Theory and Applications. 1995Google Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.