- Research article
- Open access
- Published:

# An intuitionistic approach to scoring DNA sequences against transcription factor binding site motifs

*BMC Bioinformatics*
**volumeÂ 11**, ArticleÂ number:Â 551 (2010)

## Abstract

### Background

Transcription factors (TFs) control transcription by binding to specific regions of DNA called transcription factor binding sites (TFBSs). The identification of TFBSs is a crucial problem in computational biology and includes the subtask of predicting the location of known TFBS motifs in a given DNA sequence. It has previously been shown that, when scoring matches to known TFBS motifs, interdependencies between positions within a motif should be taken into account. However, this remains a challenging task owing to the fact that sequences similar to those of known TFBSs can occur by chance with a relatively high frequency. Here we present a new method for matching sequences to TFBS motifs based on intuitionistic fuzzy sets (IFS) theory, an approach that has been shown to be particularly appropriate for tackling problems that embody a high degree of uncertainty.

### Results

We propose *SC*_{
intuit
}, a new scoring method for measuring sequence-motif affinity based on IFS theory. Unlike existing methods that consider dependencies between positions, *SC*_{
intuit
} is designed to prevent overestimation of less conserved positions of TFBSs. For a given pair of bases, *SC*_{
intuit
} is computed not only as a function of their combined probability of occurrence, but also taking into account the individual importance of each single base at its corresponding position. We used *SC*_{
intuit
} to identify known TFBSs in DNA sequences. Our method provides excellent results when dealing with both synthetic and real data, outperforming the sensitivity and the specificity of two existing methods in all the experiments we performed.

### Conclusions

The results show that *SC*_{
intuit
} improves the prediction quality for TFs of the existing approaches without compromising sensitivity. In addition, we show how *SC*_{
intuit
} can be successfully applied to real research problems. In this study the reliability of the IFS theory for motif discovery tasks is proven.

## Background

Cells control the abundance of proteins by means of diverse mechanisms. One such mechanism is the regulation of transcription, which is a continuous process whereby many factors combine to ensure appropriate rates of protein synthesis. Understanding such complex processes is one of the main objectives in computational biology. In its early stages, transcription is controlled, among other mechanisms, by the binding of proteins called transcription factors (TFs) to specific regions of a given chromosome called transcription factor binding sites (TFBSs). These interactions between proteins and DNA usually take place upstream from the gene, close to the transcription start site (TSS), in the so-called promoter region of the gene.

One of the biggest issues in identifying TFBSs is that a single binding protein can bind to different DNA sequences. Related DNA sequences to which the same TF can bind are grouped together into a TFBS motif. The identification of TFBSs within a given set of DNA sequences is an active area of research. In this context there exist two main approaches: *i)* the *de novo* discovery of motifs, and *ii)* the detection of TFBSs using motifs that are already known.

*De novo* methods aim to find significant sub-sequence patterns within a set of TFBS sequences. Some of the most popular approaches are MEME [1], Gibbs sampling [2], AlignACE [3], and more recently PRIORITY and Trawler [4, 5]. For a review see [6].

Detection methods, on the other hand, focus on inferring new TFBSs from known binding motifs. Early detection methods assumed independence between positions within a putative TFBS sequence, e.g. in Patser [7] and ConSite [8]. However, it is now well established that this assumption is wrong [9â€“11], and some methods that consider position dependency for modeling and finding TFBSs using advanced HMM and Bayesian models have appeared [12, 13]. Likewise, two recent detection methods have been developed that take into account interdependencies between TFBS positions. Tomovic and Oakeley proposed a method that incorporates a measure of positional interdependence into the overall score [14]. More recently, Zare-Mirakabad et al. developed a method based on joint information content and mutual information [15]. In this method, positional dependencies are taken into account by considering all pairwise combinations of positions (see the Methods section for more information).

The fact that TFBS sequences are usually very short means that the same or very similar sequences tend to occur by chance at a relatively high frequency. Consequently one of the main goals in the prediction of TFBSs is to reduce the false positive rate without compromising sensitivity. Methods that take into account positional dependencies tend to be significantly more effective at meeting this challenge. However, there remains room for improvement. As we will show in the Results section, existing methods have some drawbacks, such as overlearning of the training data, arbitrary threshold selection for testing dependencies, etc. The purpose of the work presented here is to provide a new method for measuring sequence-motif affinity that improves on existing approaches.

Zadeh proposed fuzzy sets theory to mathematically model the imprecision inherent in certains concepts [16]. Briefly, fuzzy sets theory allows an object to partially belong to a set with a membership degree between 0 and 1. Classical set theory is a special case of its fuzzy counterpart in which membership and certainty degrees are restricted to either 0 or 1. Atanassov proposed intuitionistic fuzzy sets (IFS) theory as an extension of the fuzzy sets theory [17]. IFSs generalize the notion of a fuzzy set representing uncertainty with respect to both the degree of membership (*Î¼*) and non-membership (*Î½*) of a set by allowing that the sum *Î¼* + *Î½* â‰¤ 1.

Owing to the fact that IFSs are capable of modelling the uncertainty present in real-life situations, they have been widely applied during the past decades to a variety of problems (see the Methods section). In recent years, it has been seen that the inherent uncertainty and noise that characterize biological data cannot always be modeled sufficiently well using probabilistic approaches and that, as a consequence, alternative approaches to modelling this uncertainty may be required [18â€“21]. In addition to the usual problems of missing values and noisy data associated with biological data, there exist some additional hidden factors that affect binding affinities in the context of sequence-motif scoring, e.g. cooperative binding and chromatin structure [22]. Furthermore, the described motifs are subject to change as new experiments confirm new binding sites. In this work we make use of IFS theory to formally model the uncertainty associated with the problem of scoring DNA sequences against TFBS motifs.

## Results

### Case studies

First, we wanted to show the ability of our proposed method, *SC*_{
intuit
}, to discriminate between the relative importance of poorly-conserved positions and well-conserved positions comparing it with the most representative scoring methods: *i)* *SC*_{
indep
}, a probabilistic method that assumes positional independence; *ii) SC*_{
dep
}, a scoring method proposed by Tomovic and Oakeley that take into account statistical interdependencies between TFBS positions [14]; and *iii)* *SC*_{
mat
}, a scoring function proposed by Zare-Mirakabad et al. based on the dependency between all pairwise combinations of binding site positions [15].

In Figure 1(A) we show the binding sequences of the motif Dof3 found in the JASPAR database. It can be observed how the first position is highly conserved while the fifth position is poorly conserved. We then considered two sequences (Figure 1(B)): *i)* a sequence with a mismatch in the conserved position; and *ii)* a sequence with a mismatch in the poorly-conserved position. As has been explained above, it would be desirable that the score obtained for the case of the mismatch in the conserved position is lower than the scoring for the other sequence, as it shares the similarities in the most conserved positions of the motif. In Table 1 we show the results for the three methods. It can be observed that our proposed scoring method discriminates between the two cases, while the other three approaches provide almost the same score for both sequences, missing the difference between the conservation level of the positions being compared.

In the majority of cases, the sequences known to belong to a given TFBS motif have very similar nucleotide compositions and highly conserved positions. However, in the databases of known motifs there are a number of examples where individual sequences differ from the majority in highly-conserved positions. Such a binding sequence can be considered an outlier with respect to the motif, i.e. a binding site that is not closely related to the other binding sites in the motif. When scoring new sequences against a given TFBS motif, we should generally tolerate small, additional variations in the sequence with respect to non-outliers, but be far less tolerant of mutations to outlier sequences. Here we evaluate the extent to which each scoring method is able to discriminate between sequences belonging to these two categories. Take, as a preliminary example, the binding sequences for motif MZF1 in the JASPAR database, as shown in Figure 2(A). It can be observed how the highlighted outlier sequence GGAGGA does not contain the higly-conserved base G at the third position, while the highlighted sequence TGGGGA is clearly a non-outlier (see motif logo in Figure 2(B)). We selected the highlighted sequences and created two new sequences by mutating its sixth position giving GGAGGG (derived from an outlier) and TGGGGT (derived from a non-outlier). In order to observe the discrimination degrees of the different scoring methods, we scored each sequence against the motif by means of the different methods.

In reality, it would be desirable that the scoring for the case of the mutated outlier sequence be lower than the scoring for the mutated non-oulier sequence. Results obtained by the *SC*_{
mat
}, *SC*_{
dep
}, and *SC*_{
indep
} methods failed to capture the expected differences, giving the incorrect impression that binding is likely to occur. On the other hand, our proposed method obtained a more realistic distance between the sequences, providing a much lower score for the mutated outlier sequence (Figure 2(C)).

These insights are confirmed in the following sections where the experiments are extended to use large datasets, and the results are measured in terms of discovery rates.

### Prediction of TFBSs

#### Synthetic sequences

In order to compare the performance of the different methods in predicting TFBSs, we used the non-redundant publicly available JASPAR motifs database for our experiments [23]. We selected all motifs for which binding sequences are available (not only matrix profiles), resulting in a dataset of 124 motifs. For each of these motifs, a random number between 2 and 6 binding sites were randomly selected and inserted in random sequences of a random length between 200 bp and 500 bp from a third-order Markov model background distribution obtained from the RSAT (Regulatory Sequence Analysis Tools) [24]. For each position of each sequence we computed the score for their corresponding motifs with an assumed known TFBS length (the length of the inserted motif).

Usually, methods have a high sensitivity (i.e. can detect true positives), so that the key difference between them is the number of false positives. Although our ultimate aim is not to rely on essentially arbitrary thresholds to assess performance, we began our analysis by following the recommendations of Tomovic and Oakeley in [14], selecting thresholds of 0.7 and 0.8 indicating a correct classification for a binding site. Table 2 shows the precision (*TP/*(*TP* + *FP* )) of the different methods. In the additional file 1: "Synthetic sequences experiment" we show thresholded results for the different methods. These indicate that our proposed scoring function performed best, giving the smallest number of false positives per TF whilst simultaneously giving a high number of true positives.

In order not to rely on the selection of an arbitrary threshold for evaluating the results, we computed a precision-recall (PR) curve for each considered method. PR curves are commonly used in information retrieval for evaluating classification performance and give a more informative picture of a method's performance than ROC (Receiver Operating Characteristic) curves [25] when dealing with highly skewed datasets as is the case here [26]. Figure 3 shows the PR graphs. *SC*_{
intuit
} produces a better PR graph than the remaining methods (see Table 3 for AUC values). In addition, ROC curves can be found in the the additional file 2: "ROC curves".

#### Mutated sequences

To further evaluate our proposed method, we obtained a set of putative binding sites that are very similar to those that are already known. This is a common scenario in motif discovery, where the set of known sequences belonging to a given binding motif is incomplete. In order to simulate this situation, we proceeded in a similar way to our previous experiment; all the steps were the same except that we gave a single base mutation at a random position within the selected binding site for each motif. PR curves and AUC values were computed to compare the performance of the different methods (Figure 4 and Table 3). ROC curves for this experiment can be obtained from the additional file 2: "ROC curves". The ROC and precision-recall graphs shows how *SC*_{
intuit
} gives consistently superior values, with a higher AUC value (Table 3). It can be observed that the improvement of the performance of our method compared to *SC*_{
indep
}, *SC*_{
mat
} and *SC*_{
dep
} grew with respect to the synthetic sequences experiment discussed in the previous section.

#### Real Data

We analyzed the performance of the proposed methods when dealing with real experimental data. In order to do so, we made use of the published ChIP-seq data on binding of TFs in embryonic stem cells from mouse by Chen et al. in [27], as provided in the supplementary material of [28]. We considered the three TFs (SMAD1, c-Myc, and STAT3) that have binding sequences available in the TRANSFAC database [29]. Thus, we obtained three sets of 200 bp sequence segments centered at TF binding locations, and we randomly selected 50 sequence segments from each set for our study (see additional file 3: "FASTA sequences"). We scanned each set of sequences using the 124 TFs from JASPAR for which binding sequences are available. The results demonstrate the superior performance of our new scoring method, as it gives the smallest number of false positives per nucleotide and per TF (Figure 5), and maintains an excellent true-positive rate (Table 4). Detailed results can be found in the additional file 4: "Motif statistics". It can be seen that our method presents consistently low false-positive rates with all three sets of sequences, whereas the performance of the other methods is much more variable.

### Study of Single Nucleotide Polymorphisms in TNFR1 Gene for the Response against Aspergillus Fumigatus

Hematological patients are typically treated by chemotherapy and/or radiation. These treatments usually produce immunosuppression and severe neutropenia. This clinical situation can be exploited by opportunistic pathogens such as *Aspergillus fumigatus* to cause a deadly infection called Invasive Pulmonary Aspergillosis (IPA) [30, 31]. The importance of finding ways to combat this pathogen is evidenced by the fact that IPA occurs in roughly 10% to 40% of hematological patients, with overall mortality rates ranging from 50% to 90% [32, 33].

Tumor necrosis factor (TNF) activates T lymphocytes in response to fungal infections through TNF receptors. One of the most important TNF receptors is TNFR1, which plays a crucial role in immune regulation and host immune responses. Experimental studies with TNFR1 knockout mice indicate that TNFR1 is indispensable in host resistance against several infections [34]. Our hypothesis is that single nucleotide polymorphisms (SNPs) in the TNFR1 gene may influence the innate immune response against *Aspergillus fumigatus*.

The gene encoding TNFR1 contain numerous polymorphisms [35, 36]. By means of different experiments, we concluded that TNFR1_{-609(G/T)} polymorphism is critical in the development of the response against *Aspergillus* because it might be regulating the cell-mediated Th1 immune response. Details on these experiments are out of the scope of this work and can be consulted in [37]. In this section, we use our proposed scoring method *SC*_{
intuit
} to investigate whether the TNFR1_{-609(G/T)} promoter polymorphism is involved in the disruption of the recognition of a potential binding site for a critical transcription factor that could influence TNFR1 transcription level.

#### TNFR1_{-609(G/T)} Polymorphism Binding Affinity

For this experiment we used TFBS motifs found in TRANSFAC database [29], which has been widely used in research involving regulatory elements [38]. In order to find interesting dependencies between the TNFR1_{-609(G/T)} SNP and TFs binding affinity we scored the human TRANSFAC TFBSs against the TNFR1_{-609(G/T)} polymorphism by means of the *SC*_{
intuit
} method.

TFs bind to short parts of the TNFR1 promoter region and, therefore, for each trial, we need to define a fragment of the promoter sequence containing the TNFR1_{-609(G/T)} SNP that might be considered as the putative TFBS. To this end, we need to determine the length of the sub-sequences and the relative offset to the position of the TNFR1_{-609(G/T)} SNP. For each of the 446 human TFs in TRANSFAC, we generated a set of putative binding sequences by using a window size of a fixed length equal to the number of position of the corresponding TF. Moving the window across the sequence in 5'-3' direction gave us the sub-sequences for the TNFR1_{-609(G/T)} SNP that we considered to be putative TFBSs (see Figure 6 for an example). Next, we scored each pair of sub-sequences (one sub-sequence for the G allele, and one for the T allele) against the given TF applying the *SC*_{
intuit
} method.

We were interested in those sub-sequences that fulfil two properties: *i)* they have a high score in one allele (G or T) so they can be considered as candidates to be binding sites, and *ii)* the score is substantially lower when considering the remaining allele so the SNP may affect to the binding affinity. For our current research we chose a conservative cut-off of 0.7, and retained TFBSs with a score above this threshold for further analysis. The scores for the selected TFBSs with respect to their corresponding alleles (sequences with a G(T) instead of a T(G) at position -609) are shown in Table 5. Subsequently we will discuss these findings from a biological perspective and show that the most interesting insights arise in the context of the ICSBP TF, which represents the highest scoring of all the human TRANSFAC motifs.

#### Functional Effect of ICSBP/IRF-8 in the TNFR1_{-609(C/T)} SNP

In the previous section, we obtained predictive results using our *SC*_{
intuit
} scoring method and TRANSFAC database (Table 5). From them, we selected four candidates according to the two properties outlined in the previous section, i.e E2A, HNF4, ICSBP, and Pax-2. We did not find described relations between IPA response for any of E2A, HNF4, and Pax-2 TFs. Logos for these TFs are provided in Figure 7.

On the other hand, we found ICSBP (also known as IRF-8) to be directly related with the purpose of our study. ICSBP/IRF-8 shows a preference for binding the T allele (see Table 5). As a member of IRF family of transcription factors it is an important modulator of IFNÎ³ signalling cascade and was identified in association on the promoter region of numerous macrophage essential genes such as IL12, IL1*Î²*, IL18, iNOS or ISG15 [39].

In addition, several genes regulated by ICSBP/IRF-8, such as MAP4K4, IL-17R, and SOCS7, are involved in different stages of the nuclear factor *Îº* B (NF*Îº* B) signaling pathway [39]. Therefore, we can hypothesize that ICSBP/IRF-8 transcription factor might be also regulating the NF*Îº* B signaling pathway through the control of the first gene of this signalling cascade, the TNFR1 gene. In support of this hypothesis, Zhao et al. established that ICSBP/IRF-8 and TNFR1 are closely related genes [40]. They found ICSBP/IRF-8 to be associated with an enhanced ubiquination of TNFR associated factor 6 (TRAF6), a protein that mediate the signal transduction from members of the TNF receptor superfamily, and the activation of AP-1 and NF*Îº* B transcription factors.

On the other hand, several studies demonstrated that ICSBP/IRF-8 promotes the differentiation and activation of dendritic cells and macrophages cells [41, 42], and that, at the same time, TNFR1 mRNA level is increased during this biological process [43].

Taken into account these observations, we hypothesize that the presence of TNFR1_{-609(G/T)} promoter polymorphisms can modify the binding affinity to ICSBP/IRF-8 (see Figure 8) and, therefore, it could be used to predict susceptibility to infection and to facilitate risk stratification of hematological patients. However, the question of whether the TNFR1 polymorphisms have biological relevance regulating mRNA TNFR1 levels through ICSBP/IRF-8 transcription factor remains unanswered. Functional analysis should be performed to demonstrate the role of TNFR1_{-609(G/T)} polymorphism mediating the binding of ICSBP/IRF-8 to TNFR1 promoter.

## Discussion

We have introduced a new IFS-based approach for scoring DNA sequences against DNA motifs called *SC*_{
intuit
}. In this work we review three scoring schemes. These approaches have several drawbacks. *SC*_{
indep
} is based on an incorrect assumption that the nucleotides of a given TFBS are independent. In that context, *SC*_{
dep
} extended the score in order to account for positional dependencies. The problems associated with unnormalized scores at each position have been pointed out [15]. In addition, the results vary depending on the choice of the method and parameters for testing the dependencies. The main drawback with *SC*_{
mat
} is that it has a tendency to overlearn the training data and consequently its performance decreases when applied to real problems. There is therefore a need for a scoring method that accounts for positional dependencies without compromising either the consistency or the accuracy of the results.

As explained above, *SC*_{
intuit
} is based on the IFS theory, which has been successfully applied to problems that suffers from noisy and imprecise data. IFS theory represents uncertainty with respect to both the degree of membership and non-membership. The uncertainty associated with the tasks of scoring DNA sequences against motifs makes intuitionistic concepts particularly suitable for handling this kind of data. Taking advantage of such properties, we define the membership and non-membership degrees of a given pair bases at a given position not only as a function of their combined probability of occurrence, but also taking into account the importance of each individual base at its corresponding position.

One of the biggest issues for this kind of scoring methods is giving high scores for the known binding sequences of the motifs without overfitting. Our proposed approach adequately solves the problem of computing the score of a given sequence against a given motif by considering the binding sequences that comprise the motif not only individually but also as part of such set of sequences. Simple experiments shows how other methods fail in capturing realistic differences, while *SC*_{
intuit
} provides good results (Figure 1, 2). Our method assigned high scores for known binding sites, disfavouring mutations in the conserved positions of the binding site.

These insights are confirmed from experiments for predicting TFBSs in large datasets. We compared the performance of the proposed scoring methods on recognizing motifs in sets of random sequences from a third-order Markov model background distribution in two circumstances: *i)* when inserting known binding sequences, and *ii)* when inserting mutated binding sequences. In both situations we found that our proposed method gave the smallest number of false positives per TF whilst simultaneously giving a high number of true positives (Figures 3, 4). More importantly, our method outperforms the other approaches when dealing with real experimental data derived from Chip-seq assays. In this case, again, the number of false positive is significantly reduced (Figure 5). Finally, we validated our method studying the effect of known SNPs of TNFR1 gene in the binding affinity of TRANSFAC TFs for the response against Aspergillus fumigatus. We found the highest scoring for ICSBP TF among all the human TRANSFAC motifs. Although functional analysis should be performed, according to several previous studies, we hypothesize that the presence of TNFR1_{-609(G/T)}polymorphisms could be used to predict susceptibility to infection of hematological patients.

In general, the obtained results on the different experiments demonstrated that the proposed intuitionistic approach provide a better and more accurate model for the detection of motifs and for the relationships between positions of the TFBSs.

## Conclusions

In the present study, we have introduced *SC*_{
intuit
}, a new scoring method for measuring sequence-motif affinity, based on IFS theory. Our main objective was to improve the prediction quality for TFs of the existing approaches, reducing the false positive rate without compromising sensitivity. We show that *SC*_{
intuit
} outperforms other approaches in motif recognition tasks, and prove how it can be successfully applied to real research problems. We have used our approach as a scanning method for the prediction of TFBSs, but it also can be incorporated with methods for *de novo* discovery of motifs. As intuitionistic theory is specially suitable for problems that deal with imprecise concepts, we are currently working on a fuzzy approach that applies the proposed scoring in an *ab initio* method to find motifs in large sets of related DNA sequences.

## Methods

### Alternative approaches

In recent years, several scoring methods for the prediction of TFBSs have been proposed. In this section we give a brief overview of those methods that take account of positional dependencies, as they have been shown to outperform methods that assume independence. Let us first introduce the notation. Let *B* = {*A*, *C*, *G*, *T*} be the set of the four DNA nucleotides. Let *D* be a set of ordered DNA sequences on *B* of length *n*. Let us suppose that we have a motif *M* = *S*_{1},..., *S*_{
t
}, where *S*_{
i
} is a DNA sequence on *D* consisting of *t* aligned binding sites of length *n*. The problem is then reduced to assigning a score to the pair formed by a given putative TFBS, *S* âˆˆ *D*, and a given motif, *M*.

In what follows we will follow the notation proposed by Wasserman and Sandelin in [44], where *F* (*b*, *i*), for *b* âˆˆ *B* and 1 â‰¤ *i* â‰¤ *n* shows the occurrences of nucleotide *b* in position *i*, and P(b,i)=\frac{F(b,i)}{t}+a(b), for *b*âˆˆ *B* and 1 â‰¤ *i* â‰¤ *n* is the corrected probability of base *b* at position *i*, where *a*(*b*) is a smoothing parameter (*a*(*b*) = 0.001). *a*(*b*) = 0.01 is usually reported but our experiments show that smaller values provide more accurate results.

#### Statistical dependencies

Tomovic and Oakeley extended the previous method that assumed positional independence [14]. The authors also followed the notation of [44] and defined *W*_{b, i}as a position weighted matrix (PWM) of base *b* in position *i* computed as:

where *P*(*b*) is the background probability of base *b*. In the case where independence is assumed, the score for a given DNA sequence *S* can be computed by summing all the values of *W*_{b, i}for every base in *S*:

The first step for extending this score involves testing the dependencies between each pair of positions *i* and *j*. The authors introduced three different methods: *i*) Ï‡^{2} test; *ii*) *G* statistics; and *iii*) Bayesian hypothesis testing. The authors used these three methods to calculate the dependencies between pairs of positions in the motifs available in the public database JASPAR [23]. The reader should note that the accurate computation of positional dependencies is still an open problem since different results are obtained depending on the method and parameters used in their computation (see Supplementary Material 2-4 in [14]). Further details about obtaining the position dependencies and multiple test corrections can be found in [14].

In order to compute the new score, the corrected probability for the bases *b*_{1}*b*_{2} ... *b*_{
m
}in the dependent positions *i*_{1}*i*_{2} ... *i*_{
m
}is defined by:

where *a*(*b*_{1},..., *b*_{
m
}) = *a*(*b*_{1}) ... *a*(*b*_{
m
}) is a smoothing parameter.

It is straightforward then to obtain values that correspond to the PWM values:

Finally, their proposed scoring function, which incorporates positional dependencies, can be computed as:

where, *k*_{1} is the number of independent positions, *k*_{2} is the number of dependent positions of order 2 (nucleotides at positions *j*_{
i
}and *j*_{i+1}) and *k*_{
m
}the number of dependent positions of order m (nucleotides at positions *j*_{
i
}, *j*_{i+1},..., *j*_{i+m-1}).

For both the *SC*_{
indep
} and *SC*_{
dep
} it is advisable to perform the following normalization:

#### Matrix based

Zare-Mirakabad et al. proposed a new scoring function based on the dependency between all pairwise combinations of binding site positions [15]. Their method is based on the mutual information matrix, defined as: (see equation (7))

and on the joint information content (JIC), defined as:

In order to compute their score, the authors defined a PWM, *W*^{PW} , containing 16 rows and (*n* Â· (*n* - 1)/2) columns for all the pairwise combinations of the positions:

where *b*_{1,} *b*_{2} âˆˆ *2* *B* and 1 â‰¤ i, *j* â‰¤ *n* and *i* â‰ *j*. For more on this method see [15].

Finally, for a given DNA sequence *S* âˆˆ *D* of length *n* the score *SC*_{
mat
} is computed as:

In order to obtain a normalized value for the score, equation (6) should be applied.

### Intuitionistic fuzzy sets

Intuitionistic fuzzy sets (IFS) theory was proposed by Atanassov [17]. It has been applied in such diverse fields as decision making [45], logic programming [46] medical diagnosis [47, 48], pattern recognition [49], etc. IFS theory is an extension of the fuzzy sets theory previously proposed by Zadeh [16] that allows the degrees of membership and non-membership to be independently uncertain, which makes the representation more flexible at capturing the current state of our understanding given inconclusive data [50, 51]. Next, we introduce some basic IFS concepts.

Let *X* be the universe of discourse. An intuitionistic fuzzy set *A* in *X* is an object having the form:

where *Î¼*_{
A
}, *Î½*_{
A
}: *X* â†’ 0[1] denote membership function and non-membership function of *A*, satisfying 0 â‰¤ *Î¼*_{
A
} + *Î½*_{
A
}â‰¤ 1 for every *x* âˆˆ *X*. Therefore, the degree of uncertainty of *x* to *A* is *Ï€*_{
A
}(*x*) = 1 - *Î¼*_{
A
} - *Î½*_{
A
}. For more on this topic please refer to [17, 50, 51].

### Intuitionistic representation of motifs

For our approach, a given motif *M* is represented as the set of IFSs of all the pairwise combinations of its positions: {I}^{M}=\left\{{I}_{i,j}^{M}\right\}, where 1 â‰¤ *i*, *j* â‰¤ *n* and *i* â‰ *j*. Each of the *i*, *j* combinations for the motif positions is then an IFS of 16 elements defined as:

where *B* Ã— *B* is the universe of discourse, i.e. the set of all 16 possible combinations of bases for two given positions *i* and *j* (AA, AC,..., TT).

#### Membership degree computation

represents the degree of membership of the pairs for the basis *b*_{1,} *b*_{2} âˆˆ *B* in a given pair of positions *i*, *j* in a motif *M*. It can be automatically computed as:

where the above notation holds. As can be seen, the membership degree is a function of the probability of the pair of bases being compared and their individual conservation. Obviously, 0\xe2\u2030\xa4{\mathrm{\xce\xbc}}_{{I}_{i,j}^{M}}({b}_{1},\phantom{\rule{0.1em}{0ex}}{b}_{2})\xe2\u2030\xa41 and the degree increases as do the corrected probabilities of bases *b*_{1} and *b*_{2} in positions *i* and *j*, as well as the individual corrected probabilities *P* (*b*_{1,} *i*) and *P* (*b*_{2,} *j*).

#### Non-membership degree computation

represents the non-membership degree of the pairs for the basis *b*_{1,} *b*_{2} âˆˆ *B* in a given pair of positions *i* and *j* in a motif *M*. It can be automatically computed as:

where I{C}_{p}^{b}=\frac{2+P(b,p){\mathrm{log}}_{2}(P(b,p))}{2} is the normalized information content of base *b* in position *p* and {\mathrm{\xce\xbd}}_{{I}_{i,j}^{M}}({b}_{1},\phantom{\rule{0.1em}{0ex}}{b}_{2}) is in the range 0\xe2\u2030\xa4{\mathrm{\xce\xbd}}_{{I}_{i,j}^{M}}({b}_{1},\phantom{\rule{0.1em}{0ex}}{b}_{2})\xe2\u2030\xa41. The reader should note that we make use of the IC of the two positions *i*, *j* to assign the allowed degree of uncertainty for such pair of positions. When the IC is high, the degree of uncertainty will be low and viceversa. Likewise, as the information content of the two basis increases, the sum of the membership degrees, {\mathrm{\xce\xbc}}_{{I}_{i,j}^{M}}({b}_{1},\phantom{\rule{0.1em}{0ex}}{b}_{2})+{\mathrm{\xce\xbd}}_{{I}_{i,j}^{M}}({b}_{1},\phantom{\rule{0.1em}{0ex}}{b}_{2}), gets closer to 1. Therefore, the non-membership degree in positions *i*, *j* is a function of the corresponding membership degree and the uncertainty level computed for *b*_{1} and *b*_{2}. It is easy to prove that {\mathrm{\xce\xbc}}_{{I}_{i,j}^{M}}({b}_{1},\phantom{\rule{0.1em}{0ex}}{b}_{2})+{\mathrm{\xce\xbd}}_{{I}_{i,j}^{M}}({b}_{1},\phantom{\rule{0.1em}{0ex}}{b}_{2})\xe2\u2030\xa41

### Scoring

In order to define our proposed score, we first introduce the simplest case of scoring a length-2 DNA subsequence *D* = *b*_{1}, *b*_{2} in the positions *i* and *j* of a motif *M*:

where \mathrm{max}({\mathrm{\xce\xbd}}_{{I}_{i,j}^{M}}) is the maximum degree of non-membership in M found in the pair of positions *i* and *j* considering all the possible combination of basis *b*_{1}, *b*_{2} âˆˆ *B*^{2}, and {\mathrm{\xce\xbc}}_{{I}_{i,j}^{M}}({b}_{1},\phantom{\rule{0.1em}{0ex}}{b}_{2}) and {\mathrm{\xce\xbd}}_{{I}_{i,j}^{M}}({b}_{1},\phantom{\rule{0.1em}{0ex}}{b}_{2}) are the membership degree and non-membership degree of the pairs for the basis (*b*_{1}, *b*_{2}) âˆˆ *B* in the pair of positions *i*, *j* of *M*, computed as stated in sections and respectively.

As with the previously defined scores, a normalization step needs to be performed in order to obtain comparable results. The source code can obtained from http://genome.ugr.es/intuit.

where are the \mathrm{min}(S{C}_{intuit}^{i,j}) and \text{max}(S{C}_{intuit}^{i,j}) are the min/max possible scores in the positions (*i*, *j*) of the motif.

Finally, for a given DNA sequence *S* âˆˆ *D* of length *n* the score *SC*_{
intuit
} is computed as:

## References

Bailey T, Elkan C: Fitting a mixture model by expectation maximization to discover motifs in biopolymers.

*Proc Int Conf Intell Syst Mol Biol*1994, 2: 28â€“36. Citeseer CiteseerLawrence C, Altschul S, Boguski M, Liu J, Neuwald A, Wootton J: Detecting subtle sequence signals: a Gibbs sampling strategy for multiple alignment.

*Chem Rev*1993, 93: 741. 10.1021/cr00018a006Hughes J, Estep P, Tavazoie S, Church G: Computational identification of cis-regulatory elements associated with groups of functionally related genes in Saccharomyces cerevisiae.

*Journal of molecular biology*2000, 296(5):1205â€“1214. 10.1006/jmbi.2000.3519Gordan R, Narlikar L, Hartemink A: Finding regulatory DNA motifs using alignment-free evolutionary conservation information.

*Nucleic Acids Research*2010, 38(6):e90. 10.1093/nar/gkp1166Haudry Y, Ramialison M, Paten B, Wittbrodt J, Ettwiller L: Using Trawler standalone to discover overrepresented motifs in DNA and RNA sequences derived from various experiments including chromatin immunoprecipitation.

*Nature Protocols*2010, 5(2):323â€“334. 10.1038/nprot.2009.158Das M, Dai H: A survey of DNA motif finding algorithms.

*BMC Bioinformatics*2007, 8(Suppl 7):S21. 10.1186/1471-2105-8-S7-S21Hertz G, Hartzell G III, Stormo G: Identification of consensus patterns in unaligned DNA sequences known to be functionally related.

*Bioinformatics*1990, 6(2):81. 10.1093/bioinformatics/6.2.81Sandelin A, Wasserman W, Lenhard B: ConSite: web-based prediction of regulatory elements using cross-species comparison.

*Nucleic acids research*2004, (32 Web Server):W249. 10.1093/nar/gkh372Benos P, Lapedes A, Stormo G: Probabilistic code for DNA recognition by proteins of the EGR family.

*Journal of molecular biology*2002, 323(4):701â€“727. 10.1016/S0022-2836(02)00917-8Bulyk M, Johnson P, Church G: Nucleotides of transcription factor binding sites exert interdependent effects on the binding affinities of transcription factors.

*Nucleic acids research*2002, 30(5):1255. 10.1093/nar/30.5.1255Eisen M: All motifs are not created equal: structural properties of transcription factor-DNA interactions and the inference of sequences specificity.

*Genome Biology*2005, 6(5):P7. 10.1186/gb-2005-6-5-p7Won K, Sandelin A, Marstrand T, Krogh A: Modeling promoter grammars with evolving hidden Markov models.

*Bioinformatics*2008, 24(15):1669. 10.1093/bioinformatics/btn254Kim N, Tharakaraman K, MariÃ±o-RamÃrez L, Spouge J: Finding sequence motifs with Bayesian models incorporating positional information: an application to transcription factor binding sites.

*BMC bioinformatics*2008, 9: 262. 10.1186/1471-2105-9-262Tomovic A, Oakeley E: Position dependencies in transcription factor binding sites.

*Bioinformatics*2007, 23(8):933. 10.1093/bioinformatics/btm055Zare-Mirakabad F, Ahrabian H, Sadeghi M, Nowzari-Dalini A, Goliaei B: New scoring schema for finding motifs in DNA Sequences.

*BMC Bioinformatics*2009, 10: 93. 10.1186/1471-2105-10-93Zadeh L: Fuzzy sets*.

*Information and control*1965, 8(3):338â€“353. 10.1016/S0019-9958(65)90241-XAtanassov K: Intuitionistic fuzzy sets.

*Fuzzy Sets and Systems*1986, 20: 87â€“96. 10.1016/S0165-0114(86)80034-3Garcia F, Lopez F, Cano C, Blanco A: FISim: A new similarity measure between transcription factor binding sites based on the fuzzy integral.

*BMC Bioinformatics*2009, 10: 224. 10.1186/1471-2105-10-224Lopez F, Blanco A, Garcia F, Cano C, Marin A: Fuzzy association rules for biological data analysis: a case study on yeast.

*BMC Bioinformatics*2008, 9: 107. 10.1186/1471-2105-9-107Liang L, Mandal V, Lu Y, Kumar D: MCM-test: a fuzzy-set-theory-based approach to differential analysis of gene pathways.

*BMC Bioinformatics*2008, 9(6):S16. 10.1186/1471-2105-9-S6-S16Chen J, Bardes E, Aronow B, Jegga A: ToppGene Suite for gene list enrichment analysis and candidate gene prioritization.

*Nucleic Acids Research*2009, (37 Web Server):W305. 10.1093/nar/gkp427Lam F, Steger D, O'Shea E: Chromatin decouples promoter threshold from dynamic range.

*Nature*2008, 453(7192):246â€“250. 10.1038/nature06867Sandelin A, Alkema W, Engstrom P, Wasserman W, Lenhard B: JASPAR: an open-access database for eukaryotic transcription factor binding profiles.

*Nucleic Acids Research*2004, (32 Database):D91. 10.1093/nar/gkh012Van Helden J: Regulatory sequence analysis tools.

*Nucleic acids research*2003, 31(13):3593. 10.1093/nar/gkg567Hanley J, McNeil B: The meaning and use of the area under a receiver operating characteristic (ROC) curve.

*Radiology*1982, 143: 29â€“36.Fawcett T: An introduction to ROC analysis.

*Pattern recognition letters*2006, 27(8):861â€“874. 10.1016/j.patrec.2005.10.010Chen X, Xu H, Yuan P, Fang F, Huss M, Vega V, Wong E, Orlov Y, Zhang W, Jiang J,

*et al*.: Integration of external signaling pathways with the core transcriptional network in embryonic stem cells.*Cell*2008, 133(6):1106â€“1117. 10.1016/j.cell.2008.04.043Sharov A, Ko M: Exhaustive Search for Over-represented DNA Sequence Motifs with CisFinder.

*DNA Research*2009, 16:(5):261â€“73. 10.1093/dnares/dsp014Matys V, Kel-Margoulis O, Fricke E, Liebich I, Land S, Barre-Dirrie A, Reuter I, Chekmenev D, Krull M, Hornischer K,

*et al*.: TRANSFAC (R) and its module TRANSCompel (R): transcriptional gene regulation in eukaryotes.*Nucleic acids research*2006, (34 Database):D108. 10.1093/nar/gkj143Denning D: Invasive aspergillosis.

*Clinical infectious diseases*1998, 781â€“803. 10.1086/513943Offner F, Cordonnier C, Ljungman P, Prentice H, Engelhard D, Bacquer D, Pauw B: Impact of previous aspergillosis on the outcome of bone marrow transplantation.

*Clinical infectious diseases*1998, 26(5):1098â€“1103. 10.1086/520274Chamilos G, Luna M, Lewis R, Bodey G, Chemaly R, Tarrand J, Safdar A, Raad I, Kontoyiannis D: Invasive fungal infections in patients with hematologic malignancies in a tertiary care cancer center: an autopsy study over a 15-year period (1989â€“2003).

*Haematologica*2006, 91(7):986.Diop G, Spadoni J, Do H, Hirtzig T, Coulonges C, Labib T, Issing W, Rappaport J, Therwath A, Lathrop M,

*et al*.: Genomic approach of AIDS pathogenesis: exhaustive genotyping of the TNFR1 gene in a French AIDS cohort.*Biomedicine & Pharmacotherapy*2005, 59(8):474â€“480.Hehlgans T, Pfeffer K: The intriguing biology of the tumour necrosis factor/tumour necrosis factor receptor superfamily: players, rules and the games.

*Immunology*2005, 115: 1. 10.1111/j.1365-2567.2005.02143.xBochud P, Chien J, Marr K, Leisenring W, Upton A, Janer M, Rodrigues S, Li S, Hansen J, Zhao L,

*et al*.: Toll-like receptor 4 polymorphisms and aspergillosis in stem-cell transplantation.*New England Journal of Medicine*2008, 359(17):1766. 10.1056/NEJMoa0802629Baker E, Chen L, Smith C, Callen D, Goodwin R, Sutherland G: Chromosomal location of the human tumor necrosis factor receptor genes.

*Cytogenetics and cell genetics*1991, 57:(2â€“3):117.Sainz J, Salas I, Lopez E, Olmedo C, Comino A, Garcia F, Blanco A, Oyonarte S, Bueno P, Jurado M: TNFR1 mRNA Expression Level and TNFR1 Gene Polymorphisms Are Predictive Markers for Susceptibility to Develop Invasive Pulmonary Aspergillosis.

*International Journal of Immunopathology and Pharmacology*2009, 22(3):557â€“565.Wingender E: TheTRANSFAC project as an example of framework technology that supports the analysis of genomic regulation.

*Brief Bioinform*2008, 9: 326â€“332. 10.1093/bib/bbn016Dror N, Alter-Koltunoff M, Azriel A, Amariglio N, Jacob-Hirsch J, Zeligson S, Morgenstern A, Tamura T, Hauser H, Rechavi G,

*et al*.: Identification of IRF-8 and IRF-1 target genes in activated macrophages.*Molecular immunology*2007, 44(4):338â€“346. 10.1016/j.molimm.2006.02.026Zhao J, Kong H, Li H, Huang B, Yang M, Zhu C, Bogunovic M, Zheng F, Mayer L, Ozato K,

*et al*.: IRF-8/interferon (IFN) consensus sequence-binding protein is involved in Toll-like receptor (TLR) signaling and contributes to the cross-talk between TLR and IFN-Î³ signaling pathways.*Journal of Biological Chemistry*2006, 281(15):10073. 10.1074/jbc.M507788200Tamura T, Ozato K: Review: ICSBP/IRF-8: its regulatory roles in the development of myeloid cells.

*Journal of Interferon & Cytokine Research*2002, 22: 145â€“152.Tamura T, Nagamura-Inoue T, Shmeltzer Z, Kuwata T, Ozato K: ICSBP directs bipotential myeloid progenitor cells to differentiate into mature macrophages.

*Immunity*2000, 13(2):155â€“165. 10.1016/S1074-7613(00)00016-9Schling P, Rudolph C, Heimerl S, Fruth S, Schmitz G: Expression of tumor necrosis factor alpha and its receptors during cellular differentiation.

*Cytokine*2006, 33(5):239â€“245. 10.1016/j.cyto.2006.02.007Wasserman W, Sandelin A: Applied bioinformatics for the identification of regulatory elements.

*Nature Reviews Genetics*2004, 5(4):276â€“287. 10.1038/nrg1315Szmidt E, Kacprzyk J: Intuitionistic fuzzy sets in group decision making.

*Notes on IFS*1996, 2: 11â€“14.Atanassov K, Gargov G: Intuitionistic fuzzy logic.

*Compt Rend Acad Bulg Sci*1990, 43: 9â€“12.De S, Biswas R, Roy A: An application of intuitionistic fuzzy sets in medical diagnosis.

*Fuzzy Sets and Systems*2001, 117(2):209â€“213. 10.1016/S0165-0114(98)00235-8Khatibi V, Montazer G: Intuitionistic fuzzy set vs. fuzzy set application in medical pattern recognition.

*Artificial Intelligence In Medicine*2009, 47: 43â€“52. 10.1016/j.artmed.2009.03.002Hung W, Yang M: Similarity measures of intuitionistic fuzzy sets based on Hausdorff distance.

*Pattern Recognition Letters*2004, 25(14):1603â€“1611. 10.1016/j.patrec.2004.06.006Atanassov K: New operations defined over the intuitionistic fuzzy sets.

*Fuzzy Sets and Systems*1994, 61(2):142. 10.1016/0165-0114(94)90229-1Atanassov K: Intuitionistic Fuzzy Sets: Theory and Applications.

*Physica-Verlag, Heidelberg, New York*1999.

## Acknowledgements

We would like to thank the anonymous reviewers for their helpful comments which significantly improved the quality of the paper. This work has been carried out as part of projects P08-TIC-4299 of J. A., Sevilla and TIN2009-13489 of DGICT, Madrid.

## Author information

### Authors and Affiliations

### Corresponding author

## Additional information

### Authors' contributions

FG designed the study, designed and implemented the intuitionistic scoring, performed the experiments, helped with the analysis of the results and drafted the paper. AB assisted with the design of the study and helped to draft the paper. AS provided help with the analysis of the results and assisted in drafting the paper. All authors read and approved the final manuscript.

## Electronic supplementary material

### 12859_2010_4134_MOESM1_ESM.XLS

Additional file 1: **Synthetic sequences experiment**. This files contains thresholded results for the different methods for the synthetic sequences experiment. (XLS 34 KB)

### 12859_2010_4134_MOESM2_ESM.PDF

Additional file 2: **ROC curves**. This file contains the ROC curves associated to the synthetic and mutated sequences experiments. (PDF 60 KB)

### 12859_2010_4134_MOESM3_ESM.ZIP

Additional file 3: **FASTA sequences**. This file contains the 50 sequence segments in FASTA format for each one of the motifs SMAD1, Myc, and STAT3. (ZIP 11 KB)

### 12859_2010_4134_MOESM4_ESM.XLS

Additional file 4: **Motif statistics**. This file contains some statistics for the experiment discussed in the section *Real data*. (XLS 9 KB)

## Authorsâ€™ original submitted files for images

Below are the links to the authorsâ€™ original submitted files for images.

## Rights and permissions

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

## About this article

### Cite this article

Garcia-Alcalde, F., Blanco, A. & Shepherd, A.J. An intuitionistic approach to scoring DNA sequences against transcription factor binding site motifs.
*BMC Bioinformatics* **11**, 551 (2010). https://doi.org/10.1186/1471-2105-11-551

Received:

Accepted:

Published:

DOI: https://doi.org/10.1186/1471-2105-11-551