- Research
- Open Access
Classification of protein sequences by means of irredundant patterns
- Matteo Comin^{1}Email author and
- Davide Verzotto^{1}
https://doi.org/10.1186/1471-2105-11-S1-S16
© Comin and Verzotto; licensee BioMed Central Ltd. 2010
- Published: 18 January 2010
Abstract
Background
The classification of protein sequences using string algorithms provides valuable insights for protein function prediction. Several methods, based on a variety of different patterns, have been previously proposed. Almost all string-based approaches discover patterns that are not "independent, " and therefore the associated scores overcount, a multiple number of times, the contribution of patterns that cover the same region of a sequence.
Results
In this paper we use a class of patterns, called irredundant, that is specifically designed to address this issue. Loosely speaking the set of irredundant patterns is the smallest class of "independent" patterns that can describe all common patterns in two sequences, thus they avoid overcounting. We present a novel discriminative method, called Irredundant Class, based on the statistics of irredundant patterns combined with the power of support vector machines.
Conclusion
Tests on benchmark data show that Irredundant Class outperforms most of the string algorithms previously proposed, and it achieves results as good as current state-of-the-art methods. Moreover the footprints of the most discriminative irredundant patterns can be used to guide the identification of functional regions in protein sequences.
Keywords
- Common Pattern
- Biological Sequence
- String Kernel
- Maximal Pattern
- Consensus Pattern
Background
The increasing availability of biological sequences, from proteins to entire genomes, poses the need for the automatic analysis and classification of such a huge collection of data. Alignment methods and pattern discovery techniques have been used, for quite some time, to attach various problems emerging in the field of computational biology. Unfortunately most of these methods do not scale well with the lengths of the biological sequences under examination and therefore they are unpractical for genome-wide applications. To overcome this recent obstacle a number of techniques, which do not relay on alignment, have been conceived; these methods are also called alignment-free methods (see [1] for a comprehensive review). Although several alignment methods have been proposed over the years, the development of tools to classify and digest entire genome sequences is still in its infancy.
In this paper we address the classification and analysis of biological sequences without resorting to alignment. We present a method for the classification of protein sequences by means of irredundant patterns. Our method achieves results as good as the current state-of-the-art methods. We show also that the most discriminative irredundant patterns, which are used to distinguish different protein families, provide valuable information on the consensus pattern shared by a particular protein family.
Related works
The protein classification problem can also be treated as a string classification problem. Historically this problem has been studied in the field of text documents classification. Unfortunately most of these approaches, developed for a different kind of strings, fail when applied to biological sequences. The main reasons of this failure are the different nature of biological sequences, particularly rich of regularities known as patterns, and because of their lengths difficult to digest for a general purpose application.
Thus a number of methods have been studied for protein classification based on primary structure. The main distinction is between generative methods against discriminative methods. The former class includes methods such as protein family profiles [2], hidden Markov models (HMMs) [3–5], and PSI-BLAST [6]. These methods try to derive a model for a set of proteins and then check whether a candidate protein fits the model or not. Unlike generative methods, discriminative methods such as [7–14] are supervised approaches focused on finding which sequences (including negative examples) can describe a set of proteins despite of another set. This class makes extensive use of support vector machines (SVMs, [15]) based on features of proteins.
Recent results [8, 10] suggest that the best-performing methods are discriminative string algorithms, sometimes called string kernels. The string-based learning algorithm extracts information from the sequences and computes either a feature vector for each sequence or directly a matrix with scores between pairs of sequences. Then sequences are seen as a set of labeled examples (positive if they are in the family and negative otherwise) and an SVM attempts to learn a decision boundary between the different classes. The first string kernel, called Fisher kernel [7], uses an HMM to provide the necessary means of converting proteins into fixed-length vectors. The vector summarizes how different the given sequence is from a typical member of the given protein family. In contrast in the Pairwise method [8] the feature vector, corresponding to a protein, is formed by all the E-values of the Smith-Waterman score (or BLAST score) between the sequence analyzed and each of the training sequences. Other methods, like Spectrum and Mismatch [9, 10], use as protein features the set of all possible substrings of amino acids of fixed length k (k-mer). If two sequences contain many of the same k-mers, their inner product under the k-Spectrum kernel will be large. Equivalently, the Mismatch kernel computes a large inner product between two sequences if these sequences contain many k-mers that differ by at most m mismatches. Then, the SVM-I-sites method [12] encodes into feature vectors the protein's three-dimensional structure, instead of using sequence similarity; conversely the eMOTIF database method [16] defines the kernel in terms of the sequence motifs that appear in a pair of sequences. These motifs were previously extracted from the eBLOCKs database using the eMOTIF-maker algorithm that derives patterns from sequence alignments [17, 18].
More recently, the Local Alignment method [11] mimics the behavior of the Smith-Waterman score to build a family of valid kernels. Following the work of [19] they defined a kernel to detect all local alignment between strings by convolving simpler kernels, and hence comparing, in a simple way, strings of different lengths which share common parts. The Profile-based Mismatch method [14] uses probabilistic profiles, such as those produced by the PSI-BLAST algorithm, to define a kernel established on the position-dependent mutation neighborhoods for inexact matching of k-mers in the data (like the Mismatch method). Whereas in the Profile-based Direct method [13] kernel functions are constructed by combining sequence profiles with different approaches for determining the similarity between pairs of protein sequences. However, this kernel makes an extensive use of hyperparameters that increases the risk of overfitting when no dedicated validation data set is used. Finally, in the Word Correlation Matrices method [20] the kernel is defined by average pairwise k-mer (or word) similarity between two sequences (similarly to the Spectrum kernel). The consequent analysis of discriminative words allows also the identification of characteristic regions in biological sequences.
In general, all pattern-based methods operate, directly or indirectly, two distinct steps: first extract patterns from sequences, then use this set of patterns as features for an automatic classification tool, e.g. SVMs.
Remarkably almost all string algorithms consider patterns that are not independent, and therefore the associated scores are obtained using a set of redundant features. In this paper we want to stress the idea that if the learning process has to deal with a set of redundant features, the resulting score is overcounting, a multiple number of times, the contribution of patterns that cover the same region of a sequence. Our conjecture is that the overcounting, due to redundant patterns, might mislead the classification. In the case of pattern-based method we can define if two patterns are independent through the resort of the notion of irredundancy. The class of irredundant patterns is specifically designed to address this issue and it was previously studied in a number of papers [21–24]. In [25] we applied a similar technique for the detection of transcription factor binding sites and we showed that the use of irredundancy is useful since redundant patterns can be filtered out without leading to overcounts.
In the next sections we will show that most approaches are using patterns, of various forms, that are redundant. We suppose that a set of irredundant patterns, and consequently a set of independent features, can improve the automatic learning and classification of sequences. Here we present a novel discriminative pairwise method for the classification of protein sequences, called Irredundant Class. Irredundant Class is based on the statistics of irredundant patterns between two sequences, combined with the discriminative power of SVMs. We selected for comparison some of the string-based learning approaches presented above, including the best-performing methods in the protein classification and remote homology detection [7-11, 20]. Our method outperforms most of the previous approaches, and it achieves results as good as current state-of-the-art methods. Moreover, given a protein family, we study the set of most discriminative irredundant patterns and their distribution. We found that the location of the functional consensus pattern for a family, as reported by PROSITE [26], can be identified by the irredundant patterns footprint.
Methods
Irredundant Class
Our method is based on the extraction of irredundant patterns that are common to two sequences. We use the idea of irredundant patterns, to avoid overcounting of patterns that cover the same region multiple times on a sequence. In this section we present the class of irredundant patterns, along with some properties.
Recently the notion of irredundancy was studied for the case of a single sequence [21, 22], in this paper we extend the notion of irredundancy to the case of two sequences. In particular, our approach is substantially inspired by [22], and we refer the reader for a more comprehensive treatment on the original notion. To keep the paper self-contained, we report here the basic concepts, although most of them are only slightly adapted.
Let s = σ_{1}σ_{2} ... σ_{ n }be a sequence of length |s| = n over an alphabet Σ. A character from Σ, say σ, is called a solid character, while a "don't care" character '.' represents any character. A pattern is a string over Σ·(Σ ∪{.})*·Σ, thus containing at least two solid characters. For instance, p = d.dg.g.i...e is a pattern that occurs at locations 1 and 15 in the sequence s = dadgggdistketvdedgsgtidfee.
To give an idea of the notion of irredundancy applied to two sequences s_{1} and s_{2}, let us consider s_{1} = abababab and s_{2} = babababa, and the two patterns p_{1} = abababa and p_{2} = ababa that are contained in both the sequences. We can note that the existence of p_{1} in both s_{1} and s_{2} affects the existence of p_{2} in all locations in which p_{2} appears. By simply deleting the last ba from p_{1} or right shifting twice p_{2} along p_{1} we can obtain p_{2}. Loosely speaking the two patterns p_{1} and p_{2} are not independent, or equivalently they are not irredundant. Intuitively, a pattern is irredundant if it cannot be deduced, along with its location list, by some other patterns. Consequently any redundant pattern can be derived from the set of all irredundant patterns without knowing the original sequences, thus it is not informative. We want to discard all redundant common patterns as non-informative for the learning process.
Definition 1. (Common pattern p, Location list ℒ_{ p }) Let s _{1}and s_{2} be two sequences on Σ. A string p on Σ ∪ {.} is a common pattern with location list ℒ_{ p }= (l_{1}, l_{2}, ..., l_{ q }) if all of the following hold. (i) |p| ≥ 2, (ii) p[1], p[|p|] ∈ Σ, (iii) p occurs both on s_{1} and s_{2}, and (iv) there exists no location l ∉ ℒ_{ p }such that p occurs at l either on s_{1} or s_{2} (the location list is of maximal size).
Clearly, in the two sequences context, a common pattern occurs at least twice, one per sequence. Extending the notation of location list we can then denote by (ℒ_{ p }+ i), 0 ≤ i ≤ |p| - 1, all the locations in ℒ_{ p }shifted to the right by i positions.
Definition 2. (p_{1} ≼ p_{2}) For characters σ_{1} and σ_{2} we write that σ_{1} ≼ σ_{2} if and only if σ_{1} is a don't care or σ_{1} = σ_{2}. Given two patterns p_{1} and p_{2}, with |p_{1}| ≤ |p_{2}|, p_{1} ≼ p_{2} holds if p_{1}|j| ≼ p_{2}[j + d] for all 1 ≤ j ≤ |p_{1}|, with 0 ≤ d ≤ |p_{2}| - |p_{1}|.
We also say in this case that p_{1} is a subpattern of p_{2} (or p_{1} occurs in p_{2}), and that p_{2} implies or extends p_{1}.
Definition 3. (Coverage) Given two patterns p_{1} and p_{2} we say that the occurrence at j of p_{2} is covered by p_{1} (or by a subpattern of p_{1} ) if p_{2} ≼ p_{1} and j ∈ ( + i) for an integer 0 ≤ i ≤ |p_{1}| - |p_{2}|.
For instance, the pattern p_{2} = ababa with location list over s_{1} = abababab and s_{2} = babababa, where denotes the occurrence at location j in the sequence s_{ h }, is covered at position by p_{1} = abababa with and i = 2. Note that because p_{2} is a subpattern of p_{1} obtained by deleting the last ba from p_{1} (i.e., the shift integer i is 0), and that because p_{2} occurs at location 2 in p_{1}. Another example with don't cares is the following, p_{3} = a.a.a with over s_{3} = aabababab and s_{4} = babacacac is covered at all the positions by p_{4} = aba.a.a with .
Remark 1. Let s_{1} and s_{2} be two sequences, and p_{1} and p_{2} be two common patterns with p_{1} ≼ p_{2}. Then, by definition p_{2} must cover at least an occurrence of p_{1} per sequence.
Definition 4. (Maximal common pattern) Let be the set of common patterns of the sequences s_{1} and s_{2}. A pattern is maximal in composition if and only if there exists no , l ≠ i, with p_{ i } ≼ p_{ l } and . A pattern p_{ i } maximal in composition is also maximal in length if and only if there exists no pattern , j ≠ i, such that p_{ i }≼ p_{ j }and . A maximal common pattern is a pattern that is maximal both in composition and length.
Requiring maximality in composition and length limits the number of common patterns that may be usefully extracted and accounted for in two sequences. However, the notion of maximality alone does not suffice to bound the number of such patterns. It can be shown that there are sequences that have an unusually large number of maximal common patterns without conveying extra information about the input. A maximal common pattern p is irredundant if p and the list ℒ_{ p }of its occurrences cannot be deduced by the union of a number of lists of subpatterns of other maximal common patterns. Conversely, we call a common pattern p redundant if p (and its location list ℒ_{ p }) can be deduced from the other common patterns without knowing the input sequences. More formally:
Definition 5. (Redundant/Irredundant common pattern) A maximal common pattern p, with location list ℒ_{ p }, is redundant if there exist maximal common subpatterns p_{ i }, 1 ≤ i ≤k, such that (i.e., every occurrence of p on s_{1} and s_{2} is already covered by some maximal common patterns). A maximal common pattern that is not redundant is called irredundant common pattern.
Since the redundant common patterns bring no extra information about the two sequences, the set of independent patterns is precisely the class of irredundant common patterns. For instance, we consider the sequences s_{1} = abababab and s_{2} = babababa of length 8. Then the list of all irredundant common patterns is the following. p_{1} = abababa with , p_{2} = bababab with . The other redundant maximal common patterns are. p_{3} = ababab with , p_{4} = bababa with , p_{5} = ababa with , p_{6} = babab with , p_{7} = abab with , p_{8} = baba with , p_{9} = aba with , p_{10} = bab with , p_{11} = ab with , p_{12} = ba with .
It is easy to check that p_{1} and p_{2} cannot be deduced by other common patterns, whereas p_{5} along with its location list can be simply deduced by p_{1}, and p_{7} can be derived from the union of the occurrence lists of some subpatterns of p_{1} and p_{2}. We want to emphasize that if the two sequences are identical there is only one irredundant common pattern, the sequence itself. This difference with the single sequence approach is due to a peculiarity of the original notion, in which the sequence itself is not considered as a valid pattern.
Definition 6. (Consensus, Meet) The consensus of two sequences s_{1} and s_{2} is obtained by matching the characters corresponding to the same positions of the sequences, and inserting a don't care in case of mismatch. Deleting all leading and trailing don't cares from the consensus yields the meet of s_{1} and s_{2}.
For instance, the consensus of the sequences aaaaab and baaaaa is .aaaa., while their meet is aaaa. Note that a meet is a common pattern between two sequences.
Let now s = σ_{1}σ_{2} ... σ_{ n }be a sequence of n characters over an alphabet Σ. We use to denote the suffix σ_{ j }σ_{ j }_{+1}... σ_{ n }of s; s^{(j)}the sequence where the location j appears; and the location list of p on s. Clearly ; and j ∈ ℒ_{ p }if and only if .
In the following we will briefly present the most important properties of the irredundant common patterns. Those properties are specular with respect to the single sequence approach, as in [22].
Lemma 1. Let p be a common pattern on s_{1} and s_{2}. Then, p is irredundant if and only if there exists such that the meet of s^{(min{j, k})}and is p for all , where s_{ h }is the other sequence with respect to s^{(j)}.
Proof. We can prove this Lemma using the Lemmas 4 and 5 described in [22]. They prove that a pattern p' is irredundant, relatively to the approach with a single sequence s, if and only if there exists an occurrence j' of p' in s such that the meet of s and is p' for all k' in . Following this work, we have that a pattern p' is irredundant if and only if at least one of its occurrences in s is not covered by other patterns. Similarly a common pattern p is irredundant, in the two sequences approach, if and only if it has at least an occurrence j in s_{1}, or in s_{2}, that is not covered by other maximal common patterns. Now assume, w.l.o.g., that j ∈ s_{1}, and therefore s_{ h }is s_{2}. To make the statement of the Lemma valid we have to examine the case in which a common pattern p'' covers together the occurrence j of p in s_{1} and another occurrence l of p in s_{1} (i.e., the same sequence in which j appears). Note that p'' covers some occurrences of p if p ≼ p"; hence, by definition, p" must cover also an occurrence of p in s_{2} (as seen in Remark 1). This contradicts the assumptions in which we have to obtain p whenever we intersect the sequences in correspondence of two particular occurrences of p: j in s_{1} and a in s_{2}.
By the previous facts and the proofs in [22], we can conclude that an irredundant common pattern must have at least an occurrence j in one of the two sequences such that the second part of the Lemma holds, and vice versa. □
To clarify the meaning of this lemma, we refer the reader to the general example that follows the algorithm.
Theorem 1. Every irredundant common pattern over two sequences s_{1} and s_{2} is the meet of a sequence and a suffix of another one.
Proof. From Lemma 1, an irredundant common pattern p (that occurs certainly in both the sequences) must appear at least in the meet of a sequence and a suffix of the other one. □
An immediate consequence of Theorem 1 is a linear bound for the cardinality of the set of irredundant common patterns. Thus
Theorem 2. The number of irredundant common patterns over two sequences s_{1} and s_{2} of length, respectively, m and n is O(m + n).
Proof. From Theorem 1 we have at most m + n - 1 irredundant common patterns, and so O(m + n). □
With its underlying convolutory structure, Lemma 1 suggests an immediate way for the extraction of irredundant common patterns from sequences and arrays, using available pattern matching with or without FFT [22].
The algorithm
The discovery of all irredundant common patterns over two sequences s_{1} and s_{2} is derived from Lemma 1. In this context we are interested in a proof of concept, because the aim is to verify the effectiveness of our method in the classification of protein sequences. The algorithm complexity is dominated by the most computationally intensive operation, step (b), which is the global searching of all occurrences of patterns in the sequences. We implemented this step by means of a naive string searching algorithm with don't cares that accounts for O((m + n)^{3}) time, where m and n are the sizes of s_{1} and s_{2}. The complete description of the algorithm follows.
Input: two sequences s_{1} and s_{2}, where |s_{1}| = m and |s_{2}| = n.
- 1.
Compute the m + n -1 meets between s_{1} or s_{2} and a suffix of the other sequence; then discard patterns of length < 2.
- 2.
For each meet p:
(a) for each occurrence j found in the previous step, called exposed occurrence, increment a counter, I_{1}[j] or I_{2}[j], depending on the sequence in which j appears;
(b) perform a string search over s_{1} and s_{2} to find the number of occurrences of p in s_{1} and s_{2}, called respectively q_{1} and q_{2};
(c) check if the meet p is irredundant (see Lemma 1) by finding at least an exposed occurrence j of p in s^{(j)}that has a counter value equal to the number of occurrences of p in the other sequence (with respect to s^{(j)}). Equivalently, find if there exists an occurrence j in s_{1} such that I_{1}[j] = q_{2} or an occurrence j in s_{2} such that I_{2}[j] = q_{1}.
Example of meet between a sequence and a suffix of the other sequence.
positionj | 1 | 2 | 3 | 4 | 5 | 6 | 7 | 8 | 9 | ||
---|---|---|---|---|---|---|---|---|---|---|---|
s _{2} | b | a | b | a | c | a | c | a | c | ||
s _{1} | a | a | b | a | b | a | b | a | b | ||
position j | 1 | 2 | 3 | 4 | 5 | 6 | 7 | 8 | 9 | ||
p | a | . | a | . | a |
Example of counters I_{1} and I_{2} of a meet.
position j | 1 | 2 | 3 | 4 | 5 | 6 | 7 | 8 | 9 |
---|---|---|---|---|---|---|---|---|---|
I_{1}[j] | 0 | 1 | 0 | 0 | 0 | 0 | 0 | 0 | 0 |
position j | 1 | 2 | 3 | 4 | 5 | 6 | 7 | 8 | 9 |
I_{2}[j] | 0 | 0 | 0 | 1 | 0 | 0 | 0 | 0 | 0 |
We note that and that . Then step (b) performs a string search of p over s_{1} and s_{2}. We obtain that is with cardinality q_{1} = 2, and that is with cardinality q_{2} = 2. Since z_{1} <q_{2} and z_{2} <q_{1} we can conclude by Lemma 1 that p is redundant.
Scoring irredundant patterns
Irredundant vs. maximal patterns.
No. | s _{1} | s _{2} | |s_{1}| = m | |s_{2}| = m | m+ n | Maximals | Irredundants | % Of irredundants |
---|---|---|---|---|---|---|---|---|
1. | 1alo | 1bjt | 597 | 760 | 1357 | ≫16697 | 1256 | ≪7.5 |
2. | 1qax | 1cxp | 316 | 466 | 782 | 8397 | 682 | 8.1 |
3. | 1gai | 1nmt | 472 | 227 | 699 | 7037 | 612 | 8.7 |
4. | 1cvu | 1lgr | 511 | 368 | 879 | 9014 | 787 | 8.7 |
5. | 1gpe | 1yrg | 392 | 343 | 735 | 6853 | 653 | 9.5 |
6. | 1qqj | 3pcc | 415 | 236 | 651 | 5090 | 566 | 11.1 |
7. | 1bxk | 1ofg | 352 | 220 | 572 | 3549 | 489 | 13.8 |
8. | 1ebf | 2nac | 169 | 188 | 357 | 1126 | 277 | 24.6 |
9. | 1a03 | 1mho | 90 | 88 | 178 | 257 | 108 | 42.0 |
10. | 1gpt | 1ayj | 47 | 50 | 97 | 64 | 45 | 70.3 |
where m, n, and |p| are, respectively, the lengths of the two sequences and the pattern p. Given a set of N sequences the input for the SVM learning process is the matrix of scores, i.e. Score(s_{ i }, s_{ j }) with 1 ≤ i, j ≤ N. The result of this process is a learning distance function that can be treated as an indefinite kernel. When applying SVMs with this kernel, we therefore have to resort to one of the workarounds discussed in [28]. In particular, in case of weak non-positivity of the learning function, we force the optimizer to stop after a maximum number of iterations. Despite these manageable problems, we successfully applied the matrix of scores as a kernel matrix in SVMs, and we retain for future work the task of bridging the gap in the non-positivity of the learning function.
Results and discussion
Experiments of Liao and Noble.
Training | Test | Training | Test | ||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
No. | Target family | Pos. | Neg. | Pos. | Neg. | No. | Target family | Pos. | Neg. | Pos. | Neg. |
0 | 7.3.5.2 | 12 | 2330 | 9 | 1746 | 27 | 7.3.10.1 | 11 | 423 | 95 | 3653 |
1 | 2.56.1.2 | 11 | 2509 | 8 | 1824 | 28 | 3.32.1.11 | 46 | 3880 | 5 | 421 |
2 | 3.1.8.1 | 19 | 3002 | 8 | 1263 | 29 | 3.32.1.13 | 43 | 3627 | 8 | 674 |
3 | 3.1.8.3 | 17 | 2686 | 10 | 1579 | 30 | 7.3.6.1 | 33 | 3203 | 9 | 873 |
4 | 1.27.1.1 | 12 | 2890 | 6 | 1444 | 31 | 7.3.6.2 | 16 | 1553 | 26 | 2523 |
5 | 1.27.1.2 | 10 | 2408 | 8 | 1926 | 32 | 7.3.6.4 | 37 | 3591 | 5 | 485 |
6 | 3.42.1.1 | 29 | 3208 | 10 | 1105 | 33 | 2.38.4.1 | 30 | 3682 | 5 | 613 |
7 | 1.45.1.2 | 33 | 3650 | 6 | 663 | 34 | 2.1.1.1 | 90 | 3102 | 31 | 1068 |
8 | 1.4.1.1 | 26 | 2256 | 23 | 1994 | 35 | 2.1.1.2 | 99 | 3412 | 22 | 758 |
9 | 2.9.1.2 | 17 | 2370 | 14 | 1951 | 36 | 3.32.1.1 | 42 | 3542 | 9 | 759 |
10 | 1.4.1.2 | 41 | 3557 | 8 | 693 | 37 | 2.38.4.3 | 24 | 2946 | 11 | 1349 |
11 | 2.9.1.3 | 26 | 3625 | 5 | 696 | 38 | 2.1.1.3 | 113 | 3895 | 8 | 275 |
12 | 1.4.1.3 | 40 | 3470 | 9 | 780 | 39 | 2.1.1.4 | 88 | 3033 | 33 | 1137 |
13 | 2.44.1.2 | 11 | 307 | 140 | 3894 | 40 | 2.38.4.5 | 26 | 3191 | 9 | 1104 |
14 | 2.9.1.4 | 21 | 2928 | 10 | 1393 | 41 | 2.1.1.5 | 94 | 3240 | 27 | 930 |
15 | 3.42.1.5 | 26 | 2876 | 13 | 1437 | 42 | 7.39.1.2 | 20 | 3204 | 7 | 1121 |
16 | 3.2.1.2 | 37 | 3002 | 16 | 1297 | 43 | 2.52.1.2 | 12 | 3060 | 5 | 1275 |
17 | 3.42.1.8 | 34 | 3761 | 5 | 552 | 44 | 7.39.1.3 | 13 | 2083 | 14 | 2242 |
18 | 3.2.1.3 | 44 | 3569 | 9 | 730 | 45 | 1.36.1.2 | 29 | 3477 | 7 | 839 |
19 | 3.2.1.4 | 46 | 3732 | 7 | 567 | 46 | 3.32.1.8 | 40 | 3374 | 11 | 927 |
20 | 3.2.1.5 | 46 | 3732 | 7 | 567 | 47 | 1.36.1.5 | 10 | 1199 | 26 | 3117 |
21 | 3.2.1.6 | 48 | 3894 | 5 | 405 | 48 | 7.41.5.1 | 10 | 2241 | 9 | 2016 |
22 | 2.28.1.1 | 18 | 1246 | 44 | 3044 | 49 | 7.41.5.2 | 10 | 2241 | 9 | 2016 |
23 | 3.3.1.2 | 22 | 3280 | 7 | 1043 | 50 | 1.41.1.2 | 36 | 3692 | 6 | 615 |
24 | 3.2.1.7 | 48 | 3894 | 5 | 405 | 51 | 2.5.1.1 | 13 | 2345 | 11 | 1983 |
25 | 2.28.1.3 | 56 | 3875 | 6 | 415 | 52 | 2.5.1.3 | 14 | 2525 | 10 | 1803 |
26 | 3.3.1.5 | 13 | 1938 | 16 | 2385 | 53 | 1.41.1.5 | 17 | 1744 | 25 | 2563 |
For the comparison against state-of-the-art methods we use as metric the receiver operating characteristic (ROC) score. Given a ranking of the test sequences scores in output from the SVM, ROC score is the normalized area under the curve that plots the number of positive examples correctly predicted (true positives) as a function of the number of negative examples found to be positive (false positives) for each different possible classification threshold of a specific experiment. Our results were obtained using the Gist SVM version 2.3; realized by Noble and Pavlidis [30], it implements the soft margin optimization algorithm described in [31].
Comparison of results against state-of-the-art methods
Learning algorithm | Mean ROC | Mean ROC50 | Mean mRFP |
---|---|---|---|
Irredundant Class | 0.929 | 0.524 | 0.0554 |
Local Alignment ("ekm," β = 0.5) | 0.929 | 0.600 | 0.0515 |
Local Alignment ("eig," β = 0.5) | 0.925 | 0.649 | 0.0541 |
Word Correlation Matrices (k = 6) | 0.904 | 0.447 | 0.0778 |
Pairwise | 0.896 | 0.464 | 0.0837 |
Mismatch (k = 5, m = 1) | 0.872 | 0.400 | 0.0837 |
Spectrum (k = 3) | 0.824 | 0.294 | 0.1535 |
Fisher | 0.773 | 0.250 | 0.2040 |
Analysis of irredundant patterns footprint
Although the classification of protein sequences, through an SVM, does not provide an alignment per se, one can use the footprint of irredundant patterns to detect candidate functional sites in protein sequences. Here we are not interested in aligning a set of sequences, but we want to analyze the distribution of the most discriminative irredundant common patterns.
We recall that the result of the SVM learning process, for a target protein family, is a set of weights α = (α_{1}, ..., α_{ N }) associated to the N training sequences of its superfamily. We want to study the distribution of irredundant common patterns in the test sequences using for each pattern p a weight that is proportional to its score and to the weight α_{ i }of the corresponding training sequence that generated p, with 1 ≤ i ≤ N. Consider a test sequence s_{ test }and the set of training sequences s_{1}, ..., s_{ N }; each pair (s_{ test }, s_{ i }) generates a set of irredundant patterns ℐ_{ test, i }. For each pattern p in ℐ_{ test, i }we compute its score as the product ln ( ) × α_{ i }and we associate this score to the positions of s test covered by solid characters of p. We repeat the same process for all patterns in ℐ_{ test, i }, for each 1 ≤ i ≤ N; and for each location we sum the contributions of all patterns that cover that location. We obtain a histogram of the footprint of the irredundant patterns for the test sequence s_{ test }. The interpretation of this histogram is that picks should correspond to conserved regions, thus to candidate functional sites.
Note that these results are obtained comparing sequences from a protein family and its superfamily, thus the chances to find the actual signal are weak as opposed to standard alignment methods, which consider only the protein family. Nevertheless the profile of the family functional site can be computed as a post process of our analysis by a multiple alignment of candidate regions.
This analysis does not yield to an alignment of sequences, but it is a way to interpret the distribution of irredundant patterns. In summary the most discriminative patterns contain information about the functional site of a protein family, but this information is not explicitly available by simple inspection.
Conclusion
In this paper we studied how the notion of irredundant patterns can solve an issue that is rising in the field of string-based learning algorithms. Almost all string-based approaches consider patterns that are not independent, and therefore the associated scores overcount, a multiple number of times, the contribution of patterns that cover the same region of a sequence, called redundant patterns. Thus we use the class of irredundant common patterns over two sequences, which is specifically designed to address this issue. We design a novel method, called Irredundant Class, that is a discriminative pairwise learning algorithm based on support vector machines. Results on benchmark data show that the Irredundant Class outperforms most of the approaches previously proposed, and it achieves results as good as current state-of-the-art methods. Moreover the footprint of the most discriminative irredundant patterns can be interpreted, in a biological fashion, as a guide for the identification of characteristic regions in protein sequences. Finally, the Irredundant Class approach is not limited to protein sequences, but it can also be applied to DNA or RNA sequences. The investigation of these areas will be part of future work.
Declarations
Acknowledgements
We want to thank Raffaele Giancarlo for having brought to our attention this problem and the related work. M.C. is partially supported by the Ateneo Project of University of Padova and by the CaRiPaRo Project.
This article has been published as part of BMC Bioinformatics Volume 11 Supplement 1, 2010: Selected articles from the Eighth Asia-Pacific Bioinformatics Conference (APBC 2010). The full contents of the supplement are available online at http://www.biomedcentral.com/1471-2105/11?issue=S1.
Authors’ Affiliations
References
- Vinga S, Almeida J: Alignment-free sequence comparison - a review. Bioinformatics 2003, 19(4):513–523. 10.1093/bioinformatics/btg005View ArticlePubMedGoogle Scholar
- Gribskov M, McLachlan AD, Eisenberg D: Profile analysis: detection of distantly related proteins. Proc Natl Acad Sci U S A 1987, 84(13):4355–4358. 10.1073/pnas.84.13.4355PubMed CentralView ArticlePubMedGoogle Scholar
- Krogh A, Brown M, Mian IS, Sjöolander K, Haussler D: Hidden Markov models in computational biology: applications to protein modeling. Journal of Molecular Biology 1994, 235: 1501–1531. 10.1006/jmbi.1994.1104View ArticlePubMedGoogle Scholar
- Baldi P, Chauvin Y, Hunkapiller T, McClure M: Hidden Markov models of biological primary sequence information. Proc Natl Acad Sci U S A 1994, 91: 1059–1063. 10.1073/pnas.91.3.1059PubMed CentralView ArticlePubMedGoogle Scholar
- Karplus K, Barrett C, Hughey R: Hidden Markov models for detecting remote protein homologies. Bioinformatics 1998, 14: 846–856. 10.1093/bioinformatics/14.10.846View ArticlePubMedGoogle Scholar
- Altschul SF, Madden TL, Schffer AA, Schaffer RA, Zhang J, Zhang Z, Miller W, Lipman DJ: Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Research 1997, 25: 3389–3402. 10.1093/nar/25.17.3389PubMed CentralView ArticlePubMedGoogle Scholar
- Jaakkola T, Diekhans M, Haussler D: Using the Fisher kernel method to detect remote protein homologies. In Proceedings of the 7th International Conference on Intelligent Systems for Molecular Biology (ISMB 1999). Edited by: Thomas Lengauer, Richard Schneider, Peer Bork, Douglas L Brutlag, Janice I Glasgow, Hans-Werner Mewes and Ralf Zimmer. AAAI; 1999:149–158.Google Scholar
- Liao L, Noble WS: Combining pairwise sequence similarity and support vector machines for detecting remote protein evolutionary and structural relationships. Journal of Computational Biology 2003, 10: 857–868. 10.1089/106652703322756113View ArticlePubMedGoogle Scholar
- Leslie CS, Eskin E, Noble WS: The spectrum kernel: a string kernel for SVM protein classification. In Proceedings of the 7th Pacific Symposium on Biocomputing (PSB 2002) Edited by: Russ B Altman, A Keith Dunker, Lawrence Hunter and Teri E Klein. 2002, 564–575.Google Scholar
- Leslie CS, Eskin E, Cohen A, Weston J, Noble WS: Mismatch string kernels for discriminative protein classification. Bioinformatics 2004, 20(4):467–476. 10.1093/bioinformatics/btg431View ArticlePubMedGoogle Scholar
- Saigo H, Vert JP, Ueda N, Akutsu T: Protein homology detection using string alignment kernels. Bioinformatics 2004, 20(11):1682–1689. 10.1093/bioinformatics/bth141View ArticlePubMedGoogle Scholar
- Hou Y, Hsu W, Lee ML, Bystroff C: Efficient remote homology detection using local structure. Bioinformatics 2003, 19: 2294–2301. 10.1093/bioinformatics/btg317View ArticlePubMedGoogle Scholar
- Rangwala H, Karypis G: Profile-based direct kernels for remote homology detection and fold recognition. Bioinformatics 2005, 21(23):4239–4247. 10.1093/bioinformatics/bti687View ArticlePubMedGoogle Scholar
- Kuang R, Ie E, Wang K, Wang K, Siddiqi M, Freund Y, Leslie C: Profile-based string kernels for remote homology detection and motif extraction. Journal of Bioinformatics and Computational Biology 2005, 3(3):527–550. 10.1142/S021972000500120XView ArticlePubMedGoogle Scholar
- Vapnik VN: Statistical Learning Theory. Wiley-Interscience; 1998.Google Scholar
- Ben-Hur A, Brutlag DL: Remote homology detection: a motif based approach. Bioinformatics 2003, 19(Suppl 1):i26-i33. 10.1093/bioinformatics/btg1002View ArticlePubMedGoogle Scholar
- Nevill-Manning CG, Wu TD, Brutlag DL: Highly specific protein sequence motifs for genome analysis. Proc Natl Acad Sci U S A 1998, 95: 5865–5871. 10.1073/pnas.95.11.5865PubMed CentralView ArticlePubMedGoogle Scholar
- Huang JY, Brutlag DL: The eMOTIF database. Nucleic Acids Research 2001, 29: 202–204. 10.1093/nar/29.1.202PubMed CentralView ArticlePubMedGoogle Scholar
- Haussler D: Convolution kernels on discrete structures. In Technical Report. University of California at Santa Cruz, Santa Cruz, CA, USA; 1999.Google Scholar
- Lingner T, Meinicke P: Word correlation matrices for protein sequence analysis and remote homology detection. BMC Bioinformatics 2008, 9: 259. 10.1186/1471-2105-9-259PubMed CentralView ArticlePubMedGoogle Scholar
- Apostolico A, Tagliacollo C: Incremental discovery of the irredundant motif bases for all suffixes of a string in O ( n^{2}log n ) time. Theoretical Computer Science 2008, 408(2–3):106–115. 10.1016/j.tcs.2008.08.002View ArticleGoogle Scholar
- Apostolico A, Tagliacollo C: Optimal offline extraction of irredundant motif bases. Lecture Notes in Computer Science. In Proceedings of the 13th Annual International Computing and Cobinatorics Conference (COCOON 2007). Volume 4598. Edited by: Lin G. Springer; 2007:360–371.Google Scholar
- Apostolico A, Parida L: Incremental paradigms of motif discovery. Journal of Computational Biology 2004, 11: 15–25. 10.1089/106652704773416867View ArticlePubMedGoogle Scholar
- Pisanti N, Crochemore M, Grossi R, Sagot MF: Bases of motifs for generating repeated patterns with wild cards. IEEE/ACM Transactions on Computational Biology and Bioinformatics 2005, 2: 40–50. 10.1109/TCBB.2005.5View ArticlePubMedGoogle Scholar
- Comin M, Parida L: Detection of subtle variations as consensus motifs. Theoretical Computer Science 2008, 395(2–3):158–170. 10.1016/j.tcs.2008.01.017View ArticleGoogle Scholar
- Hulo N, Bairoch A, Bulliard V, Cerutti L, Cuche B, de Castro E, Lachaize C, Langendijk-Genevaux P, Sigrist C: The 20 years of PROSITE. Nucleic Acids Research 2008, (36 Database):D245–9.Google Scholar
- Henikoff S, Henikoff JG: Amino acid substitution matrices from protein blocks. Proc Natl Acad Sci U S A 1992, 89(22):10915–10919. 10.1073/pnas.89.22.10915PubMed CentralView ArticlePubMedGoogle Scholar
- Eichhorn J: Applications of kernel machines to structured data. PhD thesis. Technische Universität Berlin; 2007.Google Scholar
- Murzin AG, Brenner SE, Hubbard T, Chothia C: SCOP: a structural classification of proteins database for the investigation of sequences and structures. Journal of Molecular Biology 1995, 247: 536–540.PubMedGoogle Scholar
- Noble WS, Pavlidis P: Gist support vector machine and kernel principal components analysis software toolkit. Columbia University; 2002.Google Scholar
- Jaakkola T, Diekhans M, Haussler D: A discriminative framework for detecting remote protein homologies. Journal of Computational Biology 2000, 7: 95–114. 10.1089/10665270050081405View ArticlePubMedGoogle Scholar
Copyright
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.